Nanopore sequencers can be used to selectively sequence certain DNA molecules in a pool by reversing the voltage across individual nanopores to reject specific sequences, enabling enrichment and depletion to address biological questions. Previously, we achieved this using dynamic time warping to map the signal to a reference genome, but the method required substantial computational resources and did not scale to gigabase-sized references. Here we overcome this limitation by using GPU base calling. We show enrichment of specific chromosomes from the human genome and of low-abundance organisms in mixed populations without a priori knowledge of sample composition. Finally, we enrich targeted panels comprising 25,600 exons from 10,000 human genes and 717 genes implicated in cancer, identifying PML-RARA fusions in the NB4 cell line in <15 hours sequencing. These methods can be used to efficiently screen any target panel of genes without specialised sample preparation using any computer and suitable GPU. Our toolkit, readfish, is available at https://www.github.com/looselab/readfish .

Killer-cell immunoglobulin-like receptors (KIRs) are expressed predominantly on natural killer cells, where they play a key role in the regulation of innate immune responses. Recent studies show that inhibitory KIRs can also impact adaptive T cell-mediated immunity. In mice and in human T cells in vitro, inhibitory KIR ligation enhanced CD8+ T cell survival. To investigate the clinical relevance of these observations, we conducted an extensive immunogenetic analysis of multiple, independent cohorts of HIV-1, hepatitis C virus (HCV) and human T cell leukemia virus (HTLV-1)-infected individuals in conjunction with in vitro assays of T cell survival, analysis of ex vivo KIR expression and mathematical modeling of host-virus dynamics. Our data suggest that functional engagement of inhibitory KIRs enhances the CD8+ T cell response against HIV-1, HCV and HTLV-1 and is a significant determinant of clinical outcome in all three viral infections.

Nanopore sequencers enable selective sequencing of single molecules in real time by individually reversing the voltage across specific nanopores. Thus DNA molecules can be rejected and replaced with new molecules enabling targeted sequencing to enrich, deplete or achieve specific coverage in a set of reads to address a biological question. We previously demonstrated this method worked using dynamic time warping mapping signal to reference, but required significant compute and did not scale to gigabase references. Using direct base calling with GPU we can now scale to gigabase references. We enrich for specific chromosomes mapping against the human genome and we develop pipelines enriching low abundance organisms from mixed populations without prior knowledge of sample composition. Finally, we enrich panels including 25,600 exon targets from 10,000 human genes and 717 genes implicated in cancer. Using this approach we identify PML-RARA fusions in the NB4 cell line in under 15 hours sequencing. These methods can be used to efficiently screen any target panel of genes without specialised sample preparation using a single computer and suitably powerful GPU.

Variation in the risk and severity of many autoimmune diseases, malignancies and infections is strongly associated with polymorphisms at the HLA class I loci. These genetic associations provide a powerful opportunity for understanding the etiology of human disease. HLA class I associations are often interpreted in the light of ‘protective’ or ‘detrimental’ CD8+ T cell responses which are restricted by the host HLA class I allotype. However, given the diverse receptors which are bound by HLA class I molecules, alternative interpretations are possible. As well as binding T cell receptors on CD8+ T cells, HLA class I molecules are important ligands for inhibitory and activating killer immunoglobulin-like receptors (KIRs) which are found on natural killer cells and some T cells; for the CD94:NKG2 family of receptors also expressed mainly by NK cells and for leukocyte immunoglobulin-like receptors (LILRs) on myeloid cells. The aim of this study is to develop an immunogenetic approach for identifying and quantifying the relative contribution of different receptor-ligand interactions to a given HLA class I disease association and then to use this approach to investigate the immune interactions underlying HLA class I disease associations in three viral infections: Human T cell Leukemia Virus type 1, Human Immunodeficiency Virus type 1 and Hepatitis C Virus as well as in the inflammatory condition Crohn’s disease.

The processes governing lymphocyte fate (division, differentiation, and death), are typically assumed to be independent of cell age. This assumption has been challenged by a series of elegant studies which clearly show that, for murine cells in vitro, lymphocyte fate is age-dependent and that younger cells (i.e., cells which have recently divided) are less likely to divide or die. Here we investigate whether the same rules determine human T cell fate in vivo. We combined data from in vivo stable isotope labeling in healthy humans with stochastic, agent-based mathematical modeling. We show firstly that the choice of model paradigm has a large impact on parameter estimates obtained using stable isotope labeling i.e., different models fitted to the same data can yield very different estimates of T cell lifespan. Secondly, we found no evidence in humans in vivo to support the model in which younger T cells are less likely to divide or die. This age-dependent model never provided the best description of isotope labeling; this was true for naïve and memory, CD4 + and CD8 + T cells. Furthermore, this age-dependent model also failed to predict an independent data set in which the link between division and death was explored using Annexin V and deuterated glucose. In contrast, the age-independent model provided the best description of both naïve and memory T cell dynamics and was also able to predict the independent dataset.

In the early phases of the SARS coronavirus type 2 (SARS-CoV-2) pandemic, testing focused on individuals fitting a strict case definition involving a limited set of symptoms together with an identified epidemiological risk, such as contact with an infected individual or travel to a high-risk area. To assess whether this impaired our ability to detect and control early introductions of the virus into the UK, we PCR-tested archival specimens collected on admission to a large UK teaching hospital who retrospectively were identified as having a clinical presentation compatible with COVID-19. In addition, we screened available archival specimens submitted for respiratory virus diagnosis, and dating back to early January 2020, for the presence of SARS-CoV-2 RNA. Our data provides evidence for widespread community circulation of SARS-CoV-2 in early February 2020 and into March that was undetected at the time due to restrictive case definitions informing testing policy. Genome sequence data showed that many of these early cases were infected with a distinct lineage of the virus. Sequences obtained from the first officially recorded case in Nottinghamshire - a traveller returning from Daegu, South Korea – also clustered with these early UK sequences suggesting acquisition of the virus occurred in the UK and not Daegu. Analysis of a larger sample of sequences obtained in the Nottinghamshire area revealed multiple viral introductions, mainly in late February and through March. These data highlight the importance of timely and extensive community testing to prevent future widespread transmission of the virus.

Solventogenic clostridia represent a diverse group of anaerobic, spore-forming bacteria capable of producing acetone, butanol and ethanol through their unique biphasic metabolism. An intrinsic problem with these organisms however is their tendency to degenerate when repeatedly subcultured or when grown continuously. This phenomenon sees cells lose their ability to produce solvents and spores, posing a significant problem for industrial applications. To investigate the mechanistic and evolutionary basis of degeneration we combined comparative genomics, ultra-deep sequencing, and concepts of sociomicrobiology using Clostridium beijerinckii NCIMB 8052 as our model organism. These approaches revealed spo0A, the master regulator gene involved in spore and solvent formation, to be key to the degeneration process in this strain. Comparative genomics of 71 degenerate variants revealed four distinct hotspot regions that contained considerably more mutations than the rest of the genome. These included spo0A as well as genes suspected to regulate its expression and activity. Ultra-deep sequencing of populations during the subculturing process showed transient increases in mutations we believe linked to the spo0A network, however, these were ultimately dominated by mutations in the master regulator itself. Through frequency-dependent fitness assays, we found that spo0A mutants gained a fitness advantage, relative to the wild type, presumably allowing for propagation throughout the culture. Combined, our data provides new insights into the phenomenon of clostridial strain degeneration and the C. beijerinckii NCIMB 8052 solvent and spore regulation network.

No abstract