Three-Dimensional Model of Human Nicotinamide Nucleotide Transhydrogenase (NNT) and Sequence-Structure Analysis of its Disease-Causing Variations

Medical Research Council. Grant Numbers: MR/K020455/1, MR/K021613/

PolyTB: a genomic variation map for Mycobacterium tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) is the second major cause of death from an infectious disease worldwide. Recent advances in DNA sequencing are leading to the ability to generate whole genome information in clinical isolates of M. tuberculosis complex (MTBC). The identification of informative genetic variants such as phylogenetic markers and those associated with drug resistance or virulence will help barcode Mtb in the context of epidemiological, diagnostic and clinical studies. Mtb genomic datasets are increasingly available as raw sequences, which are potentially difficult and computer intensive to process, and compare across studies. Here we have processed the raw sequence data (>1500 isolates, eight studies) to compile a catalogue of SNPs (n = 74,039, 63% non-synonymous, 51.1% in more than one isolate, i.e. non-private), small indels (n = 4810) and larger structural variants (n = 800). We have developed the PolyTB web-based tool (http://pathogenseq.lshtm.ac.uk/polytb) to visualise the resulting variation and important meta-data (e.g. in silico inferred strain-types, location) within geographical map and phylogenetic views. This resource will allow researchers to identify polymorphisms within candidate genes of interest, as well as examine the genomic diversity and distribution of strains. PolyTB source code is freely available to researchers wishing to develop similar tools for their pathogen of interest

BCCTBbp: the Breast Cancer Campaign Tissue Bank bioinformatics portal

This article has been accepted for publication in Nucleic Acids Research 2015 Jan;43(Database issue):D831-6. doi: 10.1093/nar/gku984. Epub 2014 Oct 20. Published by Oxford University Press.Breast Cancer Campaign Tissue Bank [09TBBAR]; Cancer Research UK [15310 to A.Z.D.U.]. Funding for open access charge: Breast Cancer Campaign [09TBBAR

Identification of Recurrent Mutations in the microRNA-Binding Sites of B-Cell Lymphoma-Associated Genes in Follicular Lymphoma

Follicular lymphoma (FL) is a common indolent B-cell lymphoma that can transform into the more aggressive transformed FL (tFL). However, the molecular process driving this transformation is uncertain. In this work, we aimed to identify microRNA (miRNA)-binding sites recurrently mutated in follicular lymphoma patients, as well as in transformed FL patients. Using whole-genome sequencing data from FL tumors, we discovered 544 mutations located in bioinformatically predicted microRNA-binding sites. We then studied these specific regions using targeted sequencing in a cohort of 55 FL patients, found 16 recurrent mutations, and identified a further 69 variants. After filtering for QC, we identified 21 genes with mutated miRNA-binding sites that were also enriched for B-cell-associated genes by Gene Ontology. Over 40% of mutations identified in these genes were present exclusively in tFL patients. We validated the predicted miRNA-binding sites of five of the genes by luciferase assay and demonstrated that the identified mutations in BCL2 and EZH2 genes impaired the binding efficiency of miR-5008 and miR-144 and regulated the endogenous levels of messenger RNA (mRNA)

Using Whole Genome Sequences to Investigate Adenovirus Outbreaks in a Hematopoietic Stem Cell Transplant Unit

A recent surge in human mastadenovirus (HAdV) cases, including five deaths, amongst a haematopoietic stem cell transplant population led us to use whole genome sequencing (WGS) to investigate. We compared sequences from 37 patients collected over a 20-month period with sequences from GenBank and our own database of HAdVs. Maximum likelihood trees and pairwise differences were used to evaluate genotypic relationships, paired with the epidemiological data from routine infection prevention and control (IPC) records and hospital activity data. During this time period, two formal outbreaks had been declared by IPC, while WGS detected nine monophyletic clusters, seven were corroborated by epidemiological evidence and by comparison of single-nucleotide polymorphisms. One of the formal outbreaks was confirmed, and the other was not. Of the five HAdV-associated deaths, three were unlinked and the remaining two considered the source of transmission. Mixed infection was frequent (10%), providing a sentinel source of recombination and superinfection. Immunosuppressed patients harboring a high rate of HAdV positivity require comprehensive surveillance. As a consequence of these findings, HAdV WGS is being incorporated routinely into clinical practice to influence IPC policy contemporaneously

Evaluation of live attenuated Streptococcus pneumoniae vaccine strains on the epithelial response to colonisation using a human challenge model

Background: Streptococcus pneumoniae frequently colonises the human nasopharynx but causes over 500,000 deaths each year from Pneumonia, Sepsis and Meningitis. Nasopharyngeal carriage is required for transmission and is a pre-requisite for disease. The conjugate polysaccharide vaccines have proven effective in decreasing disease. However, replacement of vaccine serotypes with non-vaccine types in carriage threatens the future of the vaccines’ efficiency. Using an Experimental Human Pneumococcal Challenge model (EHPC) and epithelial cell culture models, we have previously shown that pneumococcal colonisation involves both direct epithelial association and micro-invasion, inducing innate immunity and clearance without overt disease. Repeated challenge in the EHPC with the same strain decreases subsequent carriage efficiency and diminishes transmission potential and/or progression to disease, suggesting active mucosal immunity in the nasopharynx. Methods and materials: We have generated live attenuated strains of 6B S. pneumoniae (AS1 and AS2) that have double virulence deletions and cannot revert to cause disease. Here, we have explored the hypothesis that despite their attenuation in a mouse model of disease, these attenuated strains retain their ability to invade the epithelium and induce epithelial-derived innate immunity in humans. Colonisation were measured by confocal microscopy and microbiology density by CFU counts. Epithelial activation was measured by flow cytometry, ELISA and RNAseq. Results: We found that both mutants colonised the human nasopharynx and formed epithelial associations with microinvasion in the EHPC model. In vitro, both mutants adhered, invaded and transmigrated across the epithelium 4-fold less than wild type. However, both mutants still resulted in secretion of IL-8, IL-6 and ICAM-1 secretion and barrier integrity was maintained. PCA analyses revealed that epithelial transcriptomic responses between wild type and the mutants generally overlapped, indicating overall similar stimulation of signaling pathways following exposure. Conclusion: The results reveal that attenuation of these pneumococcal strains has not led to loss of their ability to elicit amucosal immune/inflammatory response. This approach provides an exciting new pipeline for the development and testing of novel vaccines. The application of these attenuated strains in the EHPC also has the potential to provide important new knowledge on the mechanisms behind bacterial clearance, transmission and disease progression during colonisation

Rapid, comprehensive, and affordable mycobacterial diagnosis with whole-genome sequencing: a prospective study

Background Slow and cumbersome laboratory diagnostics for Mycobacterium tuberculosis complex (MTBC) risk delayed treatment and poor patient outcomes. Whole-genome sequencing (WGS) could potentially provide a rapid and comprehensive diagnostic solution. In this prospective study, we compare real-time WGS with routine MTBC diagnostic workflows. Methods We compared sequencing mycobacteria from all newly positive liquid cultures with routine laboratory diagnostic workflows across eight laboratories in Europe and North America for diagnostic accuracy, processing times, and cost between Sept 6, 2013, and April 14, 2014. We sequenced specimens once using local Illumina MiSeq platforms and processed data centrally using a semi-automated bioinformatics pipeline. We identified species or complex using gene presence or absence, predicted drug susceptibilities from resistance-conferring mutations identified from reference-mapped MTBC genomes, and calculated genetic distance to previously sequenced UK MTBC isolates to detect outbreaks. WGS data processing and analysis was done by staff masked to routine reference laboratory and clinical results. We also did a microcosting analysis to assess the financial viability of WGS-based diagnostics. Findings Compared with routine results, WGS predicted species with 93% (95% CI 90–96; 322 of 345 specimens; 356 mycobacteria specimens submitted) accuracy and drug susceptibility also with 93% (91–95; 628 of 672 specimens; 168 MTBC specimens identified) accuracy, with one sequencing attempt. WGS linked 15 (16% [95% CI 10–26]) of 91 UK patients to an outbreak. WGS diagnosed a case of multidrug-resistant tuberculosis before routine diagnosis was completed and discovered a new multidrug-resistant tuberculosis cluster. Full WGS diagnostics could be generated in a median of 9 days (IQR 6–10), a median of 21 days (IQR 14–32) faster than final reference laboratory reports were produced (median of 31 days [IQR 21–44]), at a cost of £481 per culture-positive specimen, whereas routine diagnosis costs £518, equating to a WGS-based diagnosis cost that is 7% cheaper annually than are present diagnostic workflows. Interpretation We have shown that WGS has a scalable, rapid turnaround, and is a financially feasible method for full MTBC diagnostics. Continued improvements to mycobacterial processing, bioinformatics, and analysis will improve the accuracy, speed, and scope of WGS-based diagnosis. Funding National Institute for Health Research, Department of Health, Wellcome Trust, British Colombia Centre for Disease Control Foundation for Population and Public Health, Department of Clinical Microbiology, Trinity College Dublin

Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle

Background Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals. Results Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle. Conclusions This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought

High through-put sequencing of the Parhyale hawaiensis mRNAs and microRNAs to aid comparative developmental studies

Understanding the genetic and evolutionary basis of animal morphological diversity will require comparative developmental studies that use new model organisms. This necessitates development of tools for the study of genetics and also the generation of sequence information of the organism to be studied. The development of next generation sequencing technology has enabled quick and cost effective generation of sequence information. Parhyale hawaiensis has emerged as a model organism of choice due to the development of advanced molecular tools, thus P. hawaiensis genetic information will help drive functional studies in this organism. Here we present a transcriptome and miRNA collection generated using next generation sequencing platforms. We generated approximately 1.7 million reads from a P. hawaiensis cDNA library constructed from embryos up to the germ band stage. These reads were assembled into a dataset comprising 163,501 transcripts. Using the combined annotation of Annot8r and pfam2go, Gene Ontology classifications was assigned to 20,597 transcripts. Annot8r was used to provide KEGG orthology to our transcript dataset. A total of 25,292 KEGG pathway assignments were defined and further confirmed with reciprocal blast against the NCBI nr protein database. This has identified many P. hawaiensis gene orthologs of key conserved signalling pathways involved in development. We also generated small RNA sequences from P. hawaiensis, identifying 55 conserved miRNAs. Sequenced small RNAs that were not annotated by stringent comparison to mirBase were used to search the Daphnia pulex for possible novel miRNAs. Using a conservative approach, we have identified 51 possible miRNA candidates conserved in the Daphnia pulex genome, which could be potential crustacean/arthropod specific miRNAs. Our study presents gene and miRNA discovery in a new model organism that does not have a sequenced genome. The data provided by our work will be valuable for the P. hawaiensis community as well as the wider evolutionary developmental biology community

Randomized controlled trial of molnupiravir SARS-CoV-2 viral and antibody response in at-risk adult outpatients

Viral clearance, antibody response and the mutagenic effect of molnupiravir has not 77 been elucidated in at-risk populations. Non-hospitalised participants within 5 days of 78 SARS-CoV-2 symptoms randomised to receive molnupiravir (n=253) or Usual Care 79 (n=324) were recruited to study viral and antibody dynamics and the effect of molnupiravir on viral whole genome sequence from 1437 viral genomes. Molnupiravir accelerates viral load decline, but virus is detectable by Day 5 in most cases. At Day 14 (9 days post-treatment), molnupiravir is associated with significantly higher viral persistence and significantly lower anti-SARS-CoV-2 spike antibody titres compared to Usual Care. Serial sequencing reveals increased mutagenesis with molnupiravir treatment. Persistence of detectable viral RNA at Day 14 in the molnupiravir group is associated with higher transition mutations following treatment cessation. Viral viability at Day 14 is similar in both groups with post-molnupiravir treated samples cultured up to 9 days post cessation of treatment. The current 5-day molnupiravir course is too short. Longer courses should be tested to reduce the risk of potentially transmissible molnupiravir-mutated variants being generated