The enigmatic RNase MRP of kinetoplastids.

The ribonucleoprotein RNase MRP is responsible for the processing of ribosomal RNA precursors. It is found in virtually all eukaryotes that have been examined. In the Euglenozoa, including the genera Euglena, Diplonema and kinetoplastids, MRP RNA and protein subunits have so far escaped detection using bioinformatic methods. However, we now demonstrate that the RNA component is widespread among the Euglenozoa and that these RNAs have secondary structures that conform to the structure of all other phylogenetic groups. In Euglena, we identified the same set of P/MRP protein subunits as in many other protists. However, we failed to identify any of these proteins in the kinetoplastids. This finding poses interesting questions regarding the structure and function of RNase MRP in these species.

Minimal residual disease assessed with deep sequencing of NPM1 mutations predicts relapse after allogeneic stem cell transplant in AML.

Mutations in NPM1 can be used for minimal residual disease (MRD) analysis in acute myeloid leukemia (AML). We here applied a newly introduced method, deep sequencing, allowing for simultaneous analysis of all recurrent NPM1 insertions and thus constituting an attractive alternative to multiple PCRs for the clinical laboratory. We retrospectively used deep sequencing for measurement of MRD pre- and post-allogeneic hematopoietic stem cell transplantation (alloHCT). For 29 patients in morphological remission at the time of alloHCT, the effect of deep sequencing MRD on outcome was assessed. MRD positivity was defined as variant allele frequency ≥0.02%. Post-transplant MRD status was significantly and independently associated with clinical outcome; 3-year relapse-free survival 20% vs 85% (p < .001), HR 45 (95% CI 2-1260), and overall survival 20% vs 89% (p < .001), HR 49 (95% CI 2-1253). Thus, the new methodology deep sequencing is an applicable and predictive tool for MRD assessment in AML.

Accurate and Sensitive Analysis of Minimal Residual Disease in Acute Myeloid Leukemia Using Deep Sequencing of Single Nucleotide Variations.

Minimal residual disease (MRD) in acute myeloid leukemia (AML) is of major prognostic importance. The genetic landscape of AML is characterized by numerous somatic mutations, which constitute potential MRD markers. Leukemia-specific mutations can be identified with exome sequencing at diagnosis and assessed during follow-up at low frequencies by using targeted deep sequencing. Our aim was to further validate this patient-tailored assay for substitution mutations. By applying a statistical model, which corrects for position-specific errors, a limit of detection for single nucleotide variations of variant allele frequency (VAF) of 0.02% was achieved. The assay was linear in MRD range (0.03% to 1%) with good precision [CV, 4.1% (2.2% to 5.7%) at VAF 1% and 13.3% (8.8% to 19.4%) at VAF 0.1%], and low relative bias [7.9% (2.5% to 15.3%) at VAF 1%]. When applied to six childhood AML cases and compared with multiparameter flow cytometry for MRD analysis, deep sequencing showed concordance and superior sensitivity. Further high concordance was found with expression of fusion transcripts RUNX1-RUNX1T1 and KMT2A-MLLT10. The deep sequencing assay also detected mutations in blood when VAF in bone marrow exceeded 0.1% (n = 19). In conclusion, deep sequencing enables reliable detection of low levels of residual leukemic cells. Introduction of this method in patient care will allow for highly sensitive MRD surveillance in virtually every patient with AML.

Helicobacter suis binding to carbohydrates on human and porcine gastric mucins and glycolipids occurs via two modes.

Helicobacter suis colonizes the stomach of most pigs and is the most prevalent non-Helicobacter pylori Helicobacter species found in the human stomach. In the human host, H. suis contributes to the development of chronic gastritis, peptic ulcer disease and MALT lymphoma, whereas in pigs it is associated with gastritis, decreased growth and ulcers. Here, we demonstrate that the level of H. pylori and H. suis binding to human and pig gastric mucins varies between individuals with species dependent specificity. The binding optimum of H. pylori is at neutral pH whereas that of H. suis has an acidic pH optimum, and the mucins that H. pylori bind to are different than those that H. suis bind to. Mass spectrometric analysis of mucin O-glycans from the porcine mucin showed that individual variation in binding is reflected by a difference in glycosylation; of 109 oligosaccharide structures identified, only 14 were present in all examined samples. H. suis binding to mucins correlated with glycans containing sulfate, sialic acid and terminal galactose. Among the glycolipids present in pig stomach, binding to lactotetraosylceramide (Galß3GlcNAcß3Galß4Glcß1Cer) was identified, and adhesion to Galß3GlcNAcß3Galß4Glc at both acidic and neutral pH was confirmed using other glycoconjugates. Together with that H. suis bound to DNA (used as a proxy for acidic charge), we conclude that H. suis has two binding modes: one to glycans terminating with Galß3GlcNAc, and one to negatively charged structures. Identification of the glycan structures H. suis interacts with can contribute to development of therapeutic strategies alternative to antibiotics.

CDK9 and SPT5 proteins are specifically required for expression of herpes simplex virus 1 replication-dependent late genes.

DNA replication greatly enhances expression of the herpes simplex virus 1 (HSV-1) γ2 late genes by still unknown mechanisms. Here, we demonstrate that 5,6-dichloro-1-ß-d-ribofuranosylbenzimidazole (DRB), an inhibitor of CDK9, suppresses expression of γ2 late genes with an IC50 of 5 µm, which is at least 10 times lower than the IC50 value required for inhibition of expression of early genes. The effect of DRB could not be explained by inhibition of DNA replication per se or loading of RNA polymerase II to late promoters and subsequent reduction of transcription. Instead, DRB reduces accumulation of γ2 late mRNA in the cytoplasm. In addition, we show that siRNA-mediated knockdown of the transcription factor SPT5, but not NELF-E, also gives rise to a specific inhibition of HSV-1 late gene expression. Finally, addition of DRB reduces co-immunoprecipitation of ICP27 using an anti-SPT5 antibody. Our results suggest that efficient expression of replication-dependent γ2 late genes is, at least in part, regulated by CDK9 dependent co- and/or post-transcriptional events involving SPT5 and ICP27.

Aeromonas salmonicida Growth in Response to Atlantic Salmon Mucins Differs between Epithelial Sites, Is Governed by Sialylated and N-Acetylhexosamine-Containing O-Glycans, and Is Affected by Ca2.

Aeromonas salmonicida causes furunculosis in salmonids and is a threat to Atlantic salmon aquaculture. The epithelial surfaces that the pathogen colonizes are covered by a mucus layer predominantly comprised of secreted mucins. By using mass spectrometry to identify mucin glycan structures with and without enzymatic removal of glycan residues, coupled to measurements of bacterial growth, we show here that the complex Atlantic salmon intestinal mucin glycans enhance A. salmonicida growth, whereas the more simple skin mucin glycans do not. Of the glycan residues present terminally on the salmon mucins, only N-acetylglucosamine (GlcNAc) enhances growth. Sialic acids, which have an abundance of 75% among terminal glycans from skin and of <50% among intestinal glycans, cannot be removed or used by A. salmonicida for growth-enhancing purposes, and they shield internal GlcNAc from utilization. A Ca2+ concentration above 0.1 mM is needed for A. salmonicida to be able to utilize mucins for growth-promoting purposes, and 10 mM further enhances both A. salmonicida growth in response to mucins and binding of the bacterium to mucins. In conclusion, GlcNAc and sialic acids are important determinants of the A. salmonicida interaction with its host at the mucosal surface. Furthermore, since the mucin glycan repertoire affects pathogen growth, the glycan repertoire may be a factor to take into account during breeding and selection of strains for aquaculture.

Patient-tailored analysis of minimal residual disease in acute myeloid leukemia using next-generation sequencing.

Next-generation sequencing techniques have revealed that leukemic cells in acute myeloid leukemia often are characterized by a limited number of somatic mutations. These mutations can be the basis for the detection of leukemic cells in follow-up samples. The aim of this study was to identify leukemia-specific mutations in cells from patients with acute myeloid leukemia and to use these mutations as markers for minimal residual disease. Leukemic cells and normal lymphocytes were simultaneously isolated at diagnosis from 17 patients with acute myeloid leukemia using fluorescence-activated cell sorting. Exome sequencing of these cells identified 240 leukemia-specific single nucleotide variations and 22 small insertions and deletions. Based on estimated allele frequencies and their accuracies, 191 of these mutations qualified as candidates for minimal residual disease analysis. Targeted deep sequencing with a significance threshold of 0.027% for single nucleotide variations and 0.006% for NPM1 type A mutation was developed for quantification of minimal residual disease. When tested on follow-up samples from a patient with acute myeloid leukemia, targeted deep sequencing of single nucleotide variations as well as NPM1 was more sensitive than minimal residual disease quantification with multiparameter flow cytometry. In conclusion, we here describe how exome sequencing can be used for identification of leukemia-specific mutations in samples already at diagnosis of acute myeloid leukemia. We also show that targeted deep sequencing of such mutations, including single nucleotide variations, can be used for high-sensitivity quantification of minimal residual disease in a patient-tailored manner.

Searching the Evolutionary Origin of Epithelial Mucus Protein Components-Mucins and FCGBP.

The gel-forming mucins are large glycosylated proteins that are essential components of the mucus layers covering epithelial cells. Using novel methods of identifying mucins based on profile hidden Markov models, we have found a large number of such proteins in Metazoa, aiding in their classification and allowing evolutionary studies. Most vertebrates have 5-6 gel-forming mucin genes and the genomic arrangement of these genes is well conserved throughout vertebrates. An exception is the frog Xenopus tropicalis with an expanded repertoire of at least 26 mucins of this type. Furthermore, we found that the ovomucin protein, originally identified in chicken, is characteristic of reptiles, birds, and amphibians. Muc6 is absent in teleost fish, but we now show that it is present in animals such as ghost sharks, demonstrating an early origin in vertebrate evolution. Public RNA-Seq data were analyzed with respect to mucins in zebrafish, frog, and chicken, thus allowing comparison in regard of tissue and developmental specificity. Analyses of invertebrate proteins reveal that gel-forming-mucin type of proteins is widely distributed also in this group. Their presence in Cnidaria, Porifera, and in Ctenophora (comb jellies) shows that these proteins were present early in metazoan evolution. Finally, we examined the evolution of the FCGBP protein, abundant in mucus and related to gel-forming mucins in terms of structure and localization. We demonstrate that FCGBP, ubiquitous in vertebrates, has a conserved N-terminal domain. Interestingly, this domain is also present as an N-terminal sequence in a number of bacterial proteins.

Loss of the Mediator subunit Med20 affects transcription of tRNA and other non-coding RNA genes in fission yeast.

Mediator is a co-regulator of RNA polymerase II (Pol II), transducing signals from regulatory elements and transcription factors to the general transcription machinery at the promoter. We here demonstrate that Med20 influences ribosomal protein expression in fission yeast. In addition, loss of Med20 leads to an accumulation of aberrant, readthrough tRNA transcripts. These transcripts are polyadenylated and targeted for degradation by the exosome. Similarly, other non-coding RNA molecules, such as snRNA, snoRNA and rRNA, are also enriched in the polyadenylate preparations in the absence of Med20. We suggest that fission yeast Mediator takes part in a regulatory pathway that affects Pol III-dependent transcripts.

Evolution, diversification, and expression of KNOX proteins in plants.

The KNOX (KNOTTED1-like homeobox) transcription factors play a pivotal role in leaf and meristem development. The majority of these proteins are characterized by the KNOX1, KNOX2, ELK, and homeobox domains whereas the proteins of the KNATM family contain only the KNOX domains. We carried out an extensive inventory of these proteins and here report on a total of 394 KNOX proteins from 48 species. The land plant proteins fall into two classes (I and II) as previously shown where the class I family seems to be most closely related to the green algae homologs. The KNATM proteins are restricted to Eudicots and some species have multiple paralogs of this protein. Certain plants are characterized by a significant increase in the number of KNOX paralogs; one example is Glycine max. Through the analysis of public gene expression data we show that the class II proteins of this plant have a relatively broad expression specificity as compared to class I proteins, consistent with previous studies of other plants. In G. max, class I protein are mainly distributed in axis tissues and KNATM paralogs are overall poorly expressed; highest expression is in the early plumular axis. Overall, analysis of gene expression in G. max demonstrates clearly that the expansion in gene number is associated with functional diversification.

Noncanoncial signal recognition particle RNAs in a major eukaryotic phylum revealed by purification of SRP from the human pathogen Cryptococcus neoformans.

Despite conservation of the signal recognition particle (SRP) from bacteria to man, computational approaches have failed to identify SRP components from genomes of many lower eukaryotes, raising the possibility that they have been lost or altered in those lineages. We report purification and analysis of SRP in the human pathogen Cryptococcus neoformans, providing the first description of SRP in basidiomycetous yeast. The C. neoformans SRP RNA displays a predicted structure in which the universally conserved helix 8 contains an unprecedented stem-loop insertion. Guided by this sequence, we computationally identified 152 SRP RNAs throughout the phylum Basidiomycota. This analysis revealed additional helix 8 alterations including single and double stem-loop insertions as well as loop diminutions affecting RNA structural elements that are otherwise conserved from bacteria to man. Strikingly, these SRP RNA features in Basidiomycota are accompanied by phylum-specific alterations in the RNA-binding domain of Srp54, the SRP protein subunit that directly interacts with helix 8. Our findings reveal unexpected fungal SRP diversity and suggest coevolution of the two most conserved SRP features-SRP RNA helix 8 and Srp54-in basidiomycetes. Because members of this phylum include important human and plant pathogens, these noncanonical features provide new targets for antifungal compound development.

Regulation of DNA replication at the end of the mitochondrial D-loop involves the helicase TWINKLE and a conserved sequence element.

The majority of mitochondrial DNA replication events are terminated prematurely. The nascent DNA remains stably associated with the template, forming a triple-stranded displacement loop (D-loop) structure. However, the function of the D-loop region of the mitochondrial genome remains poorly understood. Using a comparative genomics approach we here identify two closely related 15 nt sequence motifs of the D-loop, strongly conserved among vertebrates. One motif is at the D-loop 5'-end and is part of the conserved sequence block 1 (CSB1). The other motif, here denoted coreTAS, is at the D-loop 3'-end. Both these sequences may prevent transcription across the D-loop region, since light and heavy strand transcription is terminated at CSB1 and coreTAS, respectively. Interestingly, the replication of the nascent D-loop strand, occurring in a direction opposite to that of heavy strand transcription, is also terminated at coreTAS, suggesting that coreTAS is involved in termination of both transcription and replication. Finally, we demonstrate that the loading of the helicase TWINKLE at coreTAS is reversible, implying that this site is a crucial component of a switch between D-loop formation and full-length mitochondrial DNA replication.

Simultaneous DNA and RNA Mapping of Somatic Mitochondrial Mutations across Diverse Human Cancers.

Somatic mutations in the nuclear genome are required for tumor formation, but the functional consequences of somatic mitochondrial DNA (mtDNA) mutations are less understood. Here we identify somatic mtDNA mutations across 527 tumors and 14 cancer types, using an approach that takes advantage of evidence from both genomic and transcriptomic sequencing. We find that there is selective pressure against deleterious coding mutations, supporting that functional mitochondria are required in tumor cells, and also observe a strong mutational strand bias, compatible with endogenous replication-coupled errors as the major source of mutations. Interestingly, while allelic ratios in general were consistent in RNA compared to DNA, some mutations in tRNAs displayed strong allelic imbalances caused by accumulation of unprocessed tRNA precursors. The effect was explained by altered secondary structure, demonstrating that correct tRNA folding is a major determinant for processing of polycistronic mitochondrial transcripts. Additionally, the data suggest that tRNA clusters are preferably processed in the 3' to 5' direction. Our study gives insights into mtDNA function in cancer and answers questions regarding mitochondrial tRNA biogenesis that are difficult to address in controlled experimental systems.

Zebrafish as a model to study live mucus physiology.

Dysfunctional mucus barriers can result in important pulmonary and gastrointestinal conditions, but model systems to study the underlying causes are largely missing. We identified and characterized five mucin homologues in zebrafish, and demonstrated a strategy for fluorescence labeling of one selected mucin. These tools can be used for in vivo experiments and in pharmacological and genetic screens to study the dynamics and mechanisms of mucosal physiology.

The UbL protein UBTD1 stably interacts with the UBE2D family of E2 ubiquitin conjugating enzymes.

UBTD1 is a previously uncharacterized ubiquitin-like (UbL) domain containing protein with high homology to the mitochondrial Dc-UbP/UBTD2 protein. Here we show that UBTD1 and UBTD2 belong to a family of proteins that is conserved through evolution and found in metazoa, funghi, and plants. To gain further insight into the function of UBTD1, we screened for interacting proteins. In a yeast-2-hybrid (Y2H) screen, we identified several proteins involved in the ubiquitylation pathway, including the UBE2D family of E2 ubiquitin conjugating enzymes. An affinity capture screen for UBTD1 interacting proteins in whole cell extracts also identified members of the UBE2D family. Biochemical characterization of recombinant UBTD1 and UBE2D demonstrated that the two proteins form a stable, stoichiometric complex that can be purified to near homogeneity. We discuss the implications of these findings in light of the ubiquitin proteasome system (UPS).

The landscape of viral expression and host gene fusion and adaptation in human cancer.

Viruses cause 10-15% of all human cancers. Massively parallel sequencing has recently proved effective for uncovering novel viruses and virus-tumour associations, but this approach has not yet been applied to comprehensive patient cohorts. Here we screen a diverse landscape of human cancer, encompassing 4,433 tumours and 19 cancer types, for known and novel expressed viruses based on >700 billion transcriptome sequencing reads from The Cancer Genome Atlas Research Network. The resulting map confirms and extends current knowledge. We observe recurrent fusion events, including human papillomavirus insertions in RAD51B and ERBB2. Patterns of coadaptation between host and viral gene expression give clues to papillomavirus oncogene function. Importantly, our analysis argues strongly against viral aetiology in several cancers where this has frequently been proposed. We provide a virus-tumour map of unprecedented scale that constitutes a reference for future studies of tumour-associated viruses using transcriptome sequencing data.

In vivo mutagenesis reveals that OriL is essential for mitochondrial DNA replication.

The mechanisms of mitochondrial DNA replication have been hotly debated for a decade. The strand-displacement model states that lagging-strand DNA synthesis is initiated from the origin of light-strand DNA replication (OriL), whereas the strand-coupled model implies that OriL is dispensable. Mammalian mitochondria cannot be transfected and the requirements of OriL in vivo have therefore not been addressed. We here use in vivo saturation mutagenesis to demonstrate that OriL is essential for mtDNA maintenance in the mouse. Biochemical and bioinformatic analyses show that OriL is functionally conserved in vertebrates. Our findings strongly support the strand-displacement model for mtDNA replication.

RNAcentral: A vision for an international database of RNA sequences.

During the last decade there has been a great increase in the number of noncoding RNA genes identified, including new classes such as microRNAs and piRNAs. There is also a large growth in the amount of experimental characterization of these RNA components. Despite this growth in information, it is still difficult for researchers to access RNA data, because key data resources for noncoding RNAs have not yet been created. The most pressing omission is the lack of a comprehensive RNA sequence database, much like UniProt, which provides a comprehensive set of protein knowledge. In this article we propose the creation of a new open public resource that we term RNAcentral, which will contain a comprehensive collection of RNA sequences and fill an important gap in the provision of biomedical databases. We envision RNA researchers from all over the world joining a federated RNAcentral network, contributing specialized knowledge and databases. RNAcentral would centralize key data that are currently held across a variety of databases, allowing researchers instant access to a single, unified resource. This resource would facilitate the next generation of RNA research and help drive further discoveries, including those that improve food production and human and animal health. We encourage additional RNA database resources and research groups to join this effort. We aim to obtain international network funding to further this endeavor.