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.

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.

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.

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.

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.

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.

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.

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.

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).

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.

eGOB: eukaryotic Gene Order Browser.

UNLABELLED: A large number of genomes have been sequenced, allowing a range of comparative studies. Here, we present the eukaryotic Gene Order Browser with information on the order of protein and non-coding RNA (ncRNA) genes of 74 different eukaryotic species. The browser is able to display a gene of interest together with its genomic context in all species where that gene is present. Thereby, questions related to the evolution of gene organization and non-random gene order may be examined. The browser also provides access to data collected on pairs of adjacent genes that are evolutionarily conserved. AVAILABILITY: eGOB as well as underlying data are freely available at http://egob.biomedicine.gu.se SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. CONTACT: tore.samuelsson@medkem.gu.se.

U12 type introns were lost at multiple occasions during evolution.

BACKGROUND: Two categories of introns are known, a common U2 type and a rare U12 type. These two types of introns are removed by distinct spliceosomes. The phylogenetic distribution of spliceosomal RNAs that are characteristic of the U12 spliceosome, i.e. the U11, U12, U4atac and U6atac RNAs, suggest that U12 spliceosomes were lost in many phylogenetic groups. We have now examined the distribution of U2 and U12 introns in many of these groups. RESULTS: U2 and U12 introns were predicted by making use of available EST and genomic sequences. The results show that in species or branches where U12 spliceosomal components are missing, also U12 type of introns are lacking. Examples are the choanoflagellate Monosiga brevicollis, Entamoeba histolytica, green algae, diatoms, and the fungal lineage Basidiomycota. Furthermore, whereas U12 splicing does not occur in Caenorhabditis elegans, U12 introns as well as U12 snRNAs are present in Trichinella spiralis, which is deeply branching in the nematode tree. A comparison of homologous genes in T. spiralis and C. elegans revealed different mechanisms whereby U12 introns were lost. CONCLUSIONS: The phylogenetic distribution of U12 introns and spliceosomal RNAs give further support to an early origin of U12 dependent splicing. In addition, this distribution identifies a large number of instances during eukaryotic evolution where such splicing was lost.

Structure of mitochondrial transcription termination factor 3 reveals a novel nucleic acid-binding domain.

In mammalian cells, a family of mitochondrial transcription termination factors (MTERFs) regulates mitochondrial gene expression. MTERF family members share a approximately 270 residues long MTERF-domain required for DNA binding and transcription regulation. However, the structure of this widely conserved domain is unknown. Here, we show that the MTERF-domain of human MTERF3 forms a half-doughnut-shaped right-handed superhelix. The superhelix is built from alpha-helical tandem repeats that display a novel triangular three-helix motif. This repeat motif, which we denote the MTERF-motif, is a conserved structural element present in proteins from metazoans, plants, and protozoans. Furthermore, a narrow, strongly positively charged nucleic acid-binding path is found in the middle of the concave side of the half-doughnut. This arrangement suggests a half clamp nucleic acid-binding mode for MTERF-domains.

LRPPRC is a mitochondrial matrix protein that is conserved in metazoans.

LRPPRC (also called LRP130) is an RNA-binding pentatricopeptide repeat protein. LRPPRC has been recognized as a mitochondrial protein, but has also been shown to regulate nuclear gene transcription and to bind specific RNA molecules in both the nucleus and the cytoplasm. We here present a bioinformatic analysis of the LRPPRC primary sequence, which reveals that orthologs to the LRPPRC gene are restricted to metazoan cells and that all of the corresponding proteins contain mitochondrial targeting signals. To address the subcellular localization further, we have carefully analyzed LRPPRC in mammalian cells and identified a single isoform that is exclusively localized to mitochondria. The LRPPRC protein is imported to the mitochondrial matrix and its mitochondrial targeting sequence is cleaved upon entry.

Early evolution of histone mRNA 3' end processing.

The replication-dependent histone mRNAs in metazoa are not polyadenylated, in contrast to the bulk of mRNA. Instead, they contain an RNA stem-loop (SL) structure close to the 3' end of the mature RNA, and this 3' end is generated by cleavage using a machinery involving the U7 snRNP and protein factors such as the stem-loop binding protein (SLBP). This machinery of 3' end processing is related to that of polyadenylation as protein components are shared between the systems. It is commonly believed that histone 3' end processing is restricted to metazoa and green algae. In contrast, polyadenylation is ubiquitous in Eukarya. However, using computational approaches, we have now identified components of histone 3' end processing in a number of protozoa. Thus, the histone mRNA stem-loop structure as well as the SLBP protein are present in many different protozoa, including Dictyostelium, alveolates, Trypanosoma, and Trichomonas. These results show that the histone 3' end processing machinery is more ancient than previously anticipated and can be traced to the root of the eukaryotic phylogenetic tree. We also identified histone mRNAs from both metazoa and protozoa that are polyadenylated but also contain the signals characteristic of histone 3' end processing. These results provide further evidence that some histone genes are regulated at the level of 3' end processing to produce either polyadenylated RNAs or RNAs with the 3' end characteristic of replication-dependent histone mRNAs.

Computational screen for spliceosomal RNA genes aids in defining the phylogenetic distribution of major and minor spliceosomal components.

The RNA molecules of the spliceosome are critical for specificity and catalysis during splicing of eukaryotic pre-mRNA. In order to examine the evolution and phylogenetic distribution of these RNAs, we analyzed 149 eukaryotic genomes representing a broad range of phylogenetic groups. RNAs were predicted using high-sensitivity local alignment methods and profile HMMs in combination with covariance models. The results provide the most comprehensive view so far of the phylogenetic distribution of spliceosomal RNAs. RNAs were predicted in many phylogenetic groups where these RNA were not previously reported. Examples are RNAs of the major (U2-type) spliceosome in all fungal lineages, in lower metazoa and many protozoa. We also identified the minor (U12-type) spliceosomal U11 and U6atac RNAs in Acanthamoeba castellanii, where U12 spliceosomal RNA as well as minor introns were reported recently. In addition, minor-spliceosome-specific RNAs were identified in a number of phylogenetic groups where previously such RNAs were not observed, including the nematode Trichinella spiralis, the slime mold Physarum polycephalum and the fungal lineages Zygomycota and Chytridiomycota. The detailed map of the distribution of the U12-type RNA genes supports an early origin of the minor spliceosome and points to a number of occasions during evolution where it was lost.

Conserved and variable domains of RNase MRP RNA.

Ribonuclease MRP is a eukaryotic ribonucleoprotein complex consisting of one RNA molecule and 7-10 protein subunits. One important function of MRP is to catalyze an endonucleolytic cleavage during processing of rRNA precursors. RNase MRP is evolutionary related to RNase P which is critical for tRNA processing. A large number of MRP RNA sequences that now are available have been used to identify conserved primary and secondary structure features of the molecule. MRP RNA has structural features in common with P RNA such as a conserved catalytic core, but it also has unique features and is characterized by a domain highly variable between species. Information regarding primary and secondary structure features is of interest not only in basic studies of the function of MRP RNA, but also because mutations in the RNA give rise to human genetic diseases such as cartilage-hair hypoplasia.

Kinship in the SRP RNA family.

The signal recognition particle (SRP) is a ribonucleoprotein complex which participates in the targeting of protein to cellular membranes. The RNA component of the SRP has been found in all domains of life, but the size of the molecule and the number of RNA secondary structure elements vary considerably between the different phylogenetic groups. We continued our efforts to identify new SRP RNAs, compare their sequences, discover new secondary structure elements, conserved motifs, and other properties. We found additional proof for the variability in the apical loop of helix 8, and we identified several bacteria which lack all of their SRP components. Based on the distribution of SRP RNA features within the taxonomy, we suggest seven alignment groups: Bacteria with a small (4.5S) SRP RNA, Bacteria with a large (6S) SRP RNA, Archaea, Fungi (Ascomycota), Metazoa group, Protozoa group, and Plants. The proposed divisions improve the prediction of more distantly related SRP RNAs and provide a more inclusive representation of the SRP RNA family. Updates of the Rfam SRP RNA sequence collection are expected to benefit from the suggested groupings.