Long read sequencing characterises a novel structural variant, revealing underactive AKR1C1 with overactive AKR1C2 as a possible cause of severe chronic fatigue.

BACKGROUND: Causative genetic variants cannot yet be found for many disorders with a clear heritable component, including chronic fatigue disorders like myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). These conditions may involve genes in difficult-to-align genomic regions that are refractory to short read approaches. Structural variants in these regions can be particularly hard to detect or define with short reads, yet may account for a significant number of cases. Long read sequencing can overcome these difficulties but so far little data is available regarding the specific analytical challenges inherent in such regions, which need to be taken into account to ensure that variants are correctly identified. Research into chronic fatigue disorders faces the additional challenge that the heterogeneous patient populations likely encompass multiple aetiologies with overlapping symptoms, rather than a single disease entity, such that each individual abnormality may lack statistical significance within a larger sample. Better delineation of patient subgroups is needed to target research and treatment. METHODS: We use nanopore sequencing in a case of unexplained severe fatigue to identify and fully characterise a large inversion in a highly homologous region spanning the AKR1C gene locus, which was indicated but could not be resolved by short-read sequencing. We then use GC-MS/MS serum steroid analysis to investigate the functional consequences. RESULTS: Several commonly used bioinformatics tools are confounded by the homology but a combined approach including visual inspection allows the variant to be accurately resolved. The DNA inversion appears to increase the expression of AKR1C2 while limiting AKR1C1 activity, resulting in a relative increase of inhibitory GABAergic neurosteroids and impaired progesterone metabolism which could suppress neuronal activity and interfere with cellular function in a wide range of tissues. CONCLUSIONS: This study provides an example of how long read sequencing can improve diagnostic yield in research and clinical care, and highlights some of the analytical challenges presented by regions containing tandem arrays of genes. It also proposes a novel gene associated with a novel disease aetiology that may be an underlying cause of complex chronic fatigue. It reveals biomarkers that could now be assessed in a larger cohort, potentially identifying a subset of patients who might respond to treatments suggested by the aetiology.

Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes.

DNA transfer from cytoplasmic organelles to the cell nucleus is a legacy of the endosymbiotic event-the majority of nuclear-mitochondrial segments (NUMTs) are thought to be ancient, preceding human speciation1-3. Here we analyse whole-genome sequences from 66,083 people-including 12,509 people with cancer-and demonstrate the ongoing transfer of mitochondrial DNA into the nucleus, contributing to a complex NUMT landscape. More than 99% of individuals had at least one of 1,637 different NUMTs, with 1 in 8 individuals having an ultra-rare NUMT that is present in less than 0.1% of the population. More than 90% of the extant NUMTs that we evaluated inserted into the nuclear genome after humans diverged from apes. Once embedded, the sequences were no longer under the evolutionary constraint seen within the mitochondrion, and NUMT-specific mutations had a different mutational signature to mitochondrial DNA. De novo NUMTs were observed in the germline once in every 104 births and once in every 103 cancers. NUMTs preferentially involved non-coding mitochondrial DNA, linking transcription and replication to their origin, with nuclear insertion involving multiple mechanisms including double-strand break repair associated with PR domain zinc-finger protein 9 (PRDM9) binding. The frequency of tumour-specific NUMTs differed between cancers, including a probably causal insertion in a myxoid liposarcoma. We found evidence of selection against NUMTs on the basis of size and genomic location, shaping a highly heterogenous and dynamic human NUMT landscape.

An ancestral 10-bp repeat expansion in VWA1 causes recessive hereditary motor neuropathy.

The extracellular matrix comprises a network of macromolecules such as collagens, proteoglycans and glycoproteins. VWA1 (von Willebrand factor A domain containing 1) encodes a component of the extracellular matrix that interacts with perlecan/collagen VI, appears to be involved in stabilizing extracellular matrix structures, and demonstrates high expression levels in tibial nerve. Vwa1-deficient mice manifest with abnormal peripheral nerve structure/function; however, VWA1 variants have not previously been associated with human disease. By interrogating the genome sequences of 74 180 individuals from the 100K Genomes Project in combination with international gene-matching efforts and targeted sequencing, we identified 17 individuals from 15 families with an autosomal-recessive, non-length dependent, hereditary motor neuropathy and rare biallelic variants in VWA1. A single disease-associated allele p.(G25Rfs*74), a 10-bp repeat expansion, was observed in 14/15 families and was homozygous in 10/15. Given an allele frequency in European populations approaching 1/1000, the seven unrelated homozygote individuals ascertained from the 100K Genomes Project represents a substantial enrichment above expected. Haplotype analysis identified a shared 220 kb region suggesting that this founder mutation arose >7000 years ago. A wide age-range of patients (6-83 years) helped delineate the clinical phenotype over time. The commonest disease presentation in the cohort was an early-onset (mean 2.0 ± 1.4 years) non-length-dependent axonal hereditary motor neuropathy, confirmed on electrophysiology, which will have to be differentiated from other predominantly or pure motor neuropathies and neuronopathies. Because of slow disease progression, ambulation was largely preserved. Neurophysiology, muscle histopathology, and muscle MRI findings typically revealed clear neurogenic changes with single isolated cases displaying additional myopathic process. We speculate that a few findings of myopathic changes might be secondary to chronic denervation rather than indicating an additional myopathic disease process. Duplex reverse transcription polymerase chain reaction and immunoblotting using patient fibroblasts revealed that the founder allele results in partial nonsense mediated decay and an absence of detectable protein. CRISPR and morpholino vwa1 modelling in zebrafish demonstrated reductions in motor neuron axonal growth, synaptic formation in the skeletal muscles and locomotive behaviour. In summary, we estimate that biallelic variants in VWA1 may be responsible for up to 1% of unexplained hereditary motor neuropathy cases in Europeans. The detailed clinical characterization provided here will facilitate targeted testing on suitable patient cohorts. This novel disease gene may have previously evaded detection because of high GC content, consequential low coverage and computational difficulties associated with robustly detecting repeat-expansions. Reviewing previously unsolved exomes using lower QC filters may generate further diagnoses.

Profiling of Small Ribosomal Subunits Reveals Modes and Regulation of Translation Initiation.

Translation initiation is often attributed as the rate-determining step of eukaryotic protein synthesis and key to gene expression control. Despite this centrality, the series of steps involved in this process is poorly understood. Here, we capture the transcriptome-wide occupancy of ribosomes across all stages of translation initiation, enabling us to characterize the transcriptome-wide dynamics of ribosome recruitment to mRNAs, scanning across 5' UTRs and stop codon recognition, in a higher eukaryote. We provide mechanistic evidence for ribosomes attaching to the mRNA by threading the mRNA through the small subunit. Moreover, we identify features that regulate the recruitment and processivity of scanning ribosomes and redefine optimal initiation contexts. Our approach enables deconvoluting translation initiation into separate stages and identifying regulators at each step.

REPARATION: ribosome profiling assisted (re-)annotation of bacterial genomes.

Prokaryotic genome annotation is highly dependent on automated methods, as manual curation cannot keep up with the exponential growth of sequenced genomes. Current automated methods depend heavily on sequence composition and often underestimate the complexity of the proteome. We developed RibosomeE Profiling Assisted (re-)AnnotaTION (REPARATION), a de novo machine learning algorithm that takes advantage of experimental protein synthesis evidence from ribosome profiling (Ribo-seq) to delineate translated open reading frames (ORFs) in bacteria, independent of genome annotation (https://github.com/Biobix/REPARATION). REPARATION evaluates all possible ORFs in the genome and estimates minimum thresholds based on a growth curve model to screen for spurious ORFs. We applied REPARATION to three annotated bacterial species to obtain a more comprehensive mapping of their translation landscape in support of experimental data. In all cases, we identified hundreds of novel (small) ORFs including variants of previously annotated ORFs and >70% of all (variants of) annotated protein coding ORFs were predicted by REPARATION to be translated. Our predictions are supported by matching mass spectrometry proteomics data, sequence composition and conservation analysis. REPARATION is unique in that it makes use of experimental translation evidence to intrinsically perform a de novo ORF delineation in bacterial genomes irrespective of the sequence features linked to open reading frames.

Ribosome signatures aid bacterial translation initiation site identification.

BACKGROUND: While methods for annotation of genes are increasingly reliable, the exact identification of translation initiation sites remains a challenging problem. Since the N-termini of proteins often contain regulatory and targeting information, developing a robust method for start site identification is crucial. Ribosome profiling reads show distinct patterns of read length distributions around translation initiation sites. These patterns are typically lost in standard ribosome profiling analysis pipelines, when reads from footprints are adjusted to determine the specific codon being translated. RESULTS: Utilising these signatures in combination with nucleotide sequence information, we build a model capable of predicting translation initiation sites and demonstrate its high accuracy using N-terminal proteomics. Applying this to prokaryotic translatomes, we re-annotate translation initiation sites and provide evidence of N-terminal truncations and extensions of previously annotated coding sequences. These re-annotations are supported by the presence of structural and sequence-based features next to N-terminal peptide evidence. Finally, our model identifies 61 novel genes previously undiscovered in the Salmonella enterica genome. CONCLUSIONS: Signatures within ribosome profiling read length distributions can be used in combination with nucleotide sequence information to provide accurate genome-wide identification of translation initiation sites.

Maternal colonization with Streptococcus agalactiae and associated stillbirth and neonatal disease in coastal Kenya.

Streptococcus agalactiae (group B streptococcus, GBS) causes neonatal disease and stillbirth, but its burden in sub-Saharan Africa is uncertain. We assessed maternal recto-vaginal GBS colonization (7,967 women), stillbirth and neonatal disease. Whole-genome sequencing was used to determine serotypes, sequence types and phylogeny. We found low maternal GBS colonization prevalence (934/7,967, 12%), but comparatively high incidence of GBS-associated stillbirth and early onset neonatal disease (EOD) in hospital (0.91 (0.25-2.3)/1,000 births and 0.76 (0.25-1.77)/1,000 live births, respectively). However, using a population denominator, EOD incidence was considerably reduced (0.13 (0.07-0.21)/1,000 live births). Treated cases of EOD had very high case fatality (17/36, 47%), especially within 24 h of birth, making under-ascertainment of community-born cases highly likely, both here and in similar facility-based studies. Maternal GBS colonization was less common in women with low socio-economic status, HIV infection and undernutrition, but when GBS-colonized, they were more probably colonized by the most virulent clone, CC17. CC17 accounted for 267/915 (29%) of maternal colonizing (265/267 (99%) serotype III; 2/267 (0.7%) serotype IV) and 51/73 (70%) of neonatal disease cases (all serotype III). Trivalent (Ia/II/III) and pentavalent (Ia/Ib/II/III/V) vaccines would cover 71/73 (97%) and 72/73 (99%) of disease-causing serotypes, respectively. Serotype IV should be considered for inclusion, with evidence of capsular switching in CC17 strains.

Characterisation of Invasive Streptococcus pneumoniae Isolated from Cambodian Children between 2007 - 2012.

BACKGROUND: The 13-valent pneumococcal vaccine (PCV13) was introduced in Cambodia in January 2015. There are limited data concerning the common serotypes causing invasive pneumococcal disease (IPD). Knowledge of the circulating pneumococcal serotypes is important to monitor epidemiological changes before and after vaccine implementation. METHODS: All episodes of IPD defined by the isolation of Streptococcus pneumoniae from blood, cerebrospinal fluid or other sterile site in Cambodian children admitted to the Angkor Hospital for Children in Siem Reap, Northwestern Cambodia, between 1st January 2007 and 1st July 2012 were retrospectively studied. Streptococcus pneumoniae isolates that could be retrieved underwent phenotypic typing and whole genome sequencing. RESULTS: There were 90 Cambodian children hospitalized with IPD with a median (IQR) age of 2.3 years (0.9-6.2). The case fatality was 15.6% (95% CI 8-23). Of 50 Streptococcus pneumoniae isolates available for further testing, 46% were penicillin non-susceptible and 8% were ceftriaxone non-susceptible, 78% were cotrimoxazole resistant, 30% were erythromycin resistant and 30% chloramphenicol resistant. There were no significant changes in resistance levels over the five-year period. The most common serotypes were 1 (11/50; 22%), 23F (8/50; 16%), 14 (6/50; 12%), 5 (5/50; 10%) and 19A (3/50; 6%). Coverage by PCV7, PCV10 and PCV13 was 44%, 76% and 92% respectively. We identified novel multilocus sequence types and resistotypes using whole genome sequencing. CONCLUSIONS: This study suggests IPD is an important disease in Cambodian children and can have a significant mortality. PCV13 coverage of the serotypes determined in studied strains was high and consistent with another recent study. The phenotypic resistance patterns observed were similar to other regional studies. The use of whole genome sequencing in the present study provides additional typing and resistance information together with the description of novel sequence types and resistotypes.

Nested Russian Doll-Like Genetic Mobility Drives Rapid Dissemination of the Carbapenem Resistance Gene blaKPC.

The recent widespread emergence of carbapenem resistance in Enterobacteriaceae is a major public health concern, as carbapenems are a therapy of last resort against this family of common bacterial pathogens. Resistance genes can mobilize via various mechanisms, including conjugation and transposition; however, the importance of this mobility in short-term evolution, such as within nosocomial outbreaks, is unknown. Using a combination of short- and long-read whole-genome sequencing of 281 blaKPC-positive Enterobacteriaceae isolates from a single hospital over 5 years, we demonstrate rapid dissemination of this carbapenem resistance gene to multiple species, strains, and plasmids. Mobility of blaKPC occurs at multiple nested genetic levels, with transmission of blaKPC strains between individuals, frequent transfer of blaKPC plasmids between strains/species, and frequent transposition of blaKPC transposon Tn4401 between plasmids. We also identify a common insertion site for Tn4401 within various Tn2-like elements, suggesting that homologous recombination between Tn2-like elements has enhanced the spread of Tn4401 between different plasmid vectors. Furthermore, while short-read sequencing has known limitations for plasmid assembly, various studies have attempted to overcome this by the use of reference-based methods. We also demonstrate that, as a consequence of the genetic mobility observed in this study, plasmid structures can be extremely dynamic, and therefore these reference-based methods, as well as traditional partial typing methods, can produce very misleading conclusions. Overall, our findings demonstrate that nonclonal resistance gene dissemination can be extremely rapid, presenting significant challenges for public health surveillance and achieving effective control of antibiotic resistance.

Low Dose Iron Treatments Induce a DNA Damage Response in Human Endothelial Cells within Minutes.

BACKGROUND: Spontaneous reports from patients able to report vascular sequelae in real time, and recognition that serum non transferrin bound iron may reach or exceed 10µmol/L in the blood stream after iron tablets or infusions, led us to hypothesize that conventional iron treatments may provoke acute vascular injury. This prompted us to examine whether a phenotype could be observed in normal human endothelial cells treated with low dose iron. METHODOLOGY: Confluent primary human endothelial cells (EC) were treated with filter-sterilized iron (II) citrate or fresh media for RNA sequencing and validation studies. RNA transcript profiles were evaluated using directional RNA sequencing with no pre-specification of target sequences. Alignments were counted for exons and junctions of the gene strand only, blinded to treatment types. PRINCIPAL FINDINGS: Rapid changes in RNA transcript profiles were observed in endothelial cells treated with 10µmol/L iron (II) citrate, compared to media-treated cells. Clustering for Gene Ontology (GO) performed on all differentially expressed genes revealed significant differences in biological process terms between iron and media-treated EC, whereas 10 sets of an equivalent number of randomly selected genes from the respective EC gene datasets showed no significant differences in any GO terms. After 1 hour, differentially expressed genes clustered to vesicle mediated transport, protein catabolism, and cell cycle (Benjamini p = 0.0016, 0.0024 and 0.0032 respectively), and by 6 hours, to cellular response to DNA damage stimulus most significantly through DNA repair genes FANCG, BLM, and H2AFX. Comet assays demonstrated that 10µM iron treatment elicited DNA damage within 1 hour. This was accompanied by a brisk DNA damage response pulse, as ascertained by the development of DNA damage response (DDR) foci, and p53 stabilization. SIGNIFICANCE: These data suggest that low dose iron treatments are sufficient to modify the vascular endothelium, and induce a DNA damage response.

Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae at a single institution: insights into endemicity from whole-genome sequencing.

The global emergence of Klebsiella pneumoniae carbapenemase-producing K. pneumoniae (KPC-Kp) multilocus sequence type ST258 is widely recognized. Less is known about the molecular and epidemiological details of non-ST258 K. pneumoniae in the setting of an outbreak mediated by an endemic plasmid. We describe the interplay of blaKPC plasmids and K. pneumoniae strains and their relationship to the location of acquisition in a U.S. health care institution. Whole-genome sequencing (WGS) analysis was applied to KPC-Kp clinical isolates collected from a single institution over 5 years following the introduction of blaKPC in August 2007, as well as two plasmid transformants. KPC-Kp from 37 patients yielded 16 distinct sequence types (STs). Two novel conjugative blaKPC plasmids (pKPC_UVA01 and pKPC_UVA02), carried by the hospital index case, accounted for the presence of blaKPC in 21/37 (57%) subsequent cases. Thirteen (35%) isolates represented an emergent lineage, ST941, which contained pKPC_UVA01 in 5/13 (38%) and pKPC_UVA02 in 6/13 (46%) cases. Seven (19%) isolates were the epidemic KPC-Kp strain, ST258, mostly imported from elsewhere and not carrying pKPC_UVA01 or pKPC_UVA02. Using WGS-based analysis of clinical isolates and plasmid transformants, we demonstrate the unexpected dispersal of blaKPC to many non-ST258 lineages in a hospital through spread of at least two novel blaKPC plasmids. In contrast, ST258 KPC-Kp was imported into the institution on numerous occasions, with other blaKPC plasmid vectors and without sustained transmission. Instead, a newly recognized KPC-Kp strain, ST941, became associated with both novel blaKPC plasmids and spread locally, making it a future candidate for clinical persistence and dissemination.

Assessment of Mycobacterium tuberculosis transmission in Oxfordshire, UK, 2007-12, with whole pathogen genome sequences: an observational study.

BACKGROUND: Patients born outside the UK have contributed to a 20% rise in the UK's tuberculosis incidence since 2000, but their effect on domestic transmission is not known. Here we use whole-genome sequencing to investigate the epidemiology of tuberculosis transmission in an unselected population over 6 years. METHODS: We identified all residents with Oxfordshire postcodes with a Mycobacterium tuberculosis culture or a clinical diagnosis of tuberculosis between Jan 1, 2007, and Dec 31, 2012, using local databases and checking against the national Enhanced Tuberculosis Surveillance database. We used Illumina technology to sequence all available M tuberculosis cultures from identified cases. Sequences were clustered by genetic relatedness and compared retrospectively with contact investigations. The first patient diagnosed in each cluster was defined as the index case, with links to subsequent cases assigned first by use of any epidemiological linkage, then by genetic distance, and then by timing of diagnosis. FINDINGS: Although we identified 384 patients with a diagnosis of tuberculosis, country of birth was known for 380 and we sequenced isolates from 247 of 269 cases with culture-confirmed disease. 39 cases were genomically linked within 13 clusters, implying 26 local transmission events. Only 11 of 26 possible transmissions had been previously identified through contact tracing. Of seven genomically confirmed household clusters, five contained additional genomic links to epidemiologically unidentified non-household members. 255 (67%) patients were born in a country with high tuberculosis incidence, conferring a local incidence of 109 cases per 100,000 population per year in Oxfordshire, compared with 3·5 cases per 100,000 per year for those born in low-incidence countries. However, patients born in the low-incidence countries, predominantly UK, were more likely to have pulmonary disease (adjusted odds ratio 1·8 [95% CI 1·2-2·9]; p=0·009), social risk factors (4·4 [2·0-9·4]; p<0·0001), and be part of a local transmission cluster (4·8 [1·6-14·8]; p=0·006). INTERPRETATION: Although inward migration has contributed to the overall tuberculosis incidence, our findings suggest that most patients born in high-incidence countries reactivate latent infection acquired abroad and are not involved in local onward transmission. Systematic screening of new entrants could further improve tuberculosis control, but it is important that health care remains accessible to all individuals, especially high-risk groups, if tuberculosis control is not to be jeopardised. FUNDING: UK Clinical Research Collaboration (Wellcome Trust, Medical Research Council, National Institute for Health Research [NIHR]), and NIHR Oxford Biomedical Research Centre.

Whole genome sequencing and de novo assembly identifies Sydney-like variant noroviruses and recombinants during the winter 2012/2013 outbreak in England.

BACKGROUND: Norovirus is the commonest cause of epidemic gastroenteritis among people of all ages. Outbreaks frequently occur in hospitals and the community, costing the UK an estimated £110 m per annum. An evolutionary explanation for periodic increases in norovirus cases, despite some host-specific post immunity is currently limited to the identification of obvious recombinants. Our understanding could be significantly enhanced by full length genome sequences for large numbers of intensively sampled viruses, which would also assist control and vaccine design. Our objective is to develop rapid, high-throughput, end-to-end methods yielding complete norovirus genome sequences. We apply these methods to recent English outbreaks, placing them in the wider context of the international norovirus epidemic of winter 2012. METHOD: Norovirus sequences were generated from 28 unique clinical samples by Illumina RNA sequencing (RNA-Seq) of total faecal RNA. A range of de novo sequence assemblers were attempted. The best assembler was identified by validation against three replicate samples and two norovirus qPCR negative samples, together with an additional 20 sequences determined by PCR and fractional capillary sequencing. Phylogenetic methods were used to reconstruct evolutionary relationships from the whole genome sequences. RESULTS: Full length norovirus genomes were generated from 23/28 samples. 5/28 partial norovirus genomes were associated with low viral copy numbers. The de novo assembled sequences differed from sequences determined by capillary sequencing by <0.003%. Intra-host nucleotide sequence diversity was rare, but detectable by mapping short sequence reads onto its de novo assembled consensus. Genomes similar to the Sydney 2012 strain caused 78% (18/23) of cases, consistent with its previously documented association with the winter 2012 global outbreak. Interestingly, phylogenetic analysis and recombination detection analysis of the consensus sequences identified two related viruses as recombinants, containing sequences in prior circulation to Sydney 2012 in open reading frame (ORF) 2. CONCLUSION: Our approach facilitates the rapid determination of complete norovirus genomes. This method provides high resolution of full norovirus genomes which, when coupled with detailed epidemiology, may improve the understanding of evolution and control of this important healthcare-associated pathogen.

Cdc14 phosphatase promotes segregation of telomeres through repression of RNA polymerase II transcription.

Kinases and phosphatases regulate messenger RNA synthesis through post-translational modification of the carboxy-terminal domain (CTD) of the largest subunit of RNA polymerase II (ref. 1). In yeast, the phosphatase Cdc14 is required for mitotic exit(2,3) and for segregation of repetitive regions(4). Cdc14 is also a subunit of the silencing complex RENT (refs 5,6), but no roles in transcriptional repression have been described. Here we report that inactivation of Cdc14 causes silencing defects at the intergenic spacer sequences of ribosomal genes during interphase and at Y' repeats in subtelomeric regions during mitosis. We show that the role of Cdc14 in silencing is independent of the RENT deacetylase subunit Sir2. Instead, Cdc14 acts directly on RNA polymerase II by targeting CTD phosphorylation at Ser 2 and Ser 5. We also find that the role of Cdc14 as a CTD phosphatase is conserved in humans. Finally, telomere segregation defects in cdc14 mutants(4) correlate with the presence of subtelomeric Y' elements and can be rescued by transcriptional inhibition of RNA polymerase II.