Variable region V1 of Saccharomyces cerevisiae 18S rRNA participates in biogenesis and function of the small ribosomal subunit

Ribosome Structural Changes Dynamically Affect Ribosome Function

Ribosomes were known to be multicomponent complexes as early as the 1960s. Nonetheless, the prevailing view for decades considered active ribosomes to be a monolithic population, in which all ribosomes are identical in composition and function. This implied that ribosomes themselves did not actively contribute to the regulation of protein synthesis. In this perspective, I review evidence for a different model, based on results showing that ribosomes can harbor different types of ribosomal RNA (rRNA) and ribosomal proteins (r-proteins) and, furthermore, need not contain a complete set of r-proteins. I also summarize recent results favoring the notion that such distinct types of ribosomes have different affinities for specific messenger RNAs and may execute the translation process differently. Thus, ribosomes should be considered active contributors to the regulation of protein synthesis.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

With hundreds of copies of rDNA, it is unknown whether they possess sequence variations that form different types of ribosomes. Here, we developed an algorithm for long-read variant calling, termed RGA, which revealed that variations in human rDNA loci are predominantly insertion-deletion (indel) variants. We developed full-length rRNA sequencing (RIBO-RT) and in situ sequencing (SWITCH-seq), which showed that translating ribosomes possess variation in rRNA. Over 1,000 variants are lowly expressed. However, tens of variants are abundant and form distinct rRNA subtypes with different structures near indels as revealed by long-read rRNA structure probing coupled to dimethyl sulfate sequencing. rRNA subtypes show differential expression in endoderm/ectoderm-derived tissues, and in cancer, low-abundance rRNA variants can become highly expressed. Together, this study identifies the diversity of ribosomes at the level of rRNA variants, their chromosomal location, and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

Ribosomal DNA and RNA (rDNA and rRNA) sequences are usually discarded from sequencing analyses. But with hundreds of copies of rDNA genes it is unknown whether they possess sequence variations that form different types of ribosomes that affect human physiology and disease. Here, we developed an algorithm for variant-calling between paralog genes (termed RGA) and compared rDNA variations found in short- and long-read sequencing data from the 1,000 Genomes Project (1KGP) and Genome In A Bottle (GIAB). We additionally developed a novel protocol for long-read sequencing full-length rRNA (RIBO-RT) from actively translating ribosomes. Our analyses identified hundreds of rDNA variants, most of which, surprisingly, are short insertion-deletions (indels) and dozens of highly abundant rRNA variants that are incorporated into translationally active ribosomes. To visualize variant ribosomes at the single cell level, we developed an in-situ rRNA sequencing method (SWITCH-seq) which revealed that variants are co-expressed within individual cells. Strikingly, by analyzing rDNA, we found that variants assemble into distinct ribosome subtypes. We discovered that these subtypes acquire different rRNA structures by successfully employing dimethyl sulfate (DMS) probing of full length rRNA. With this atlas we investigated rRNA variation changes across human tissues and cancer types. This revealed tissue-specific rRNA subtype expression in endoderm/ectoderm-derived tissues. In cancer, low abundant rRNA variants can become highly expressed, which suggests the presence of cancer-specific ribosomes. Together, this study identifies and comprehensively characterizes the diversity of ribosomes at the level of rRNA variants which is dominated by indel variants, their chromosomal location and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

rRNA expansion segment 27Lb modulates the factor recruitment capacity of the yeast ribosome and shapes the proteome

Fine-tuned regulation of protein biosynthesis is crucial for cellular fitness and became even more vital when cellular and organismal complexity increased during the course of evolution. In order to cope with this augmented demand for translation control, eukaryal ribosomes have gained extensions both at the ribosomal protein and rRNA levels. Here we analyze the functional role of ES27L, an rRNA expansion segment in the large ribosomal subunit of Saccharomyces cerevisiae. Deletion of the b-arm of this expansion segment, called ES27Lb, did not hamper growth during optimal conditions, thus demonstrating that this 25S rRNA segment is not inherently crucial for ribosome functioning. However, reductive stress results in retarded growth and rendered unique protein sets prone to aggregation. Lack of ES27Lb negatively affects ribosome-association of known co-translational N-terminal processing enzymes which in turn contributes to the observed protein aggregation. Likely as a compensatory response to these challenges, the truncated ribosomes showed re-adjusted translation of specific sets of mRNAs and thus fine-tune the translatome in order to re-establish proteostasis. Our study gives comprehensive insight into how a highly conserved eukaryal rRNA expansion segment defines ribosomal integrity, co-translational protein maturation events and consequently cellular fitness.

Separated Siamese Twins: Intronic Small Nucleolar RNAs and Matched Host Genes May be Altered in Conjunction or Separately in Multiple Cancer Types

Small nucleolar RNAs (snoRNAs) are non-coding RNAs involved in RNA modification and processing. Approximately half of the so far identified snoRNA genes map within the intronic regions of host genes, and their expression, as well as the expression of their host genes, is dependent on transcript splicing and maturation. Growing evidence indicates that mutations and/or deregulations that affect snoRNAs, as well as host genes, play a significant role in oncogenesis. Among the possible factors underlying snoRNA/host gene expression deregulation is copy number alteration (CNA). We analyzed the data available in The Cancer Genome Atlas database, relative to CNA and expression of 295 snoRNA/host gene couples in 10 cancer types, to understand whether the genetic or expression alteration of snoRNAs and their matched host genes would have overlapping trends. Our results show that, counterintuitively, copy number and expression alterations of snoRNAs and matched host genes are not necessarily coupled. In addition, some snoRNA/host genes are mutated and overexpressed recurrently in multiple cancer types. Our findings suggest that the differential contribution to cancer development of both snoRNAs and host genes should always be considered, and that snoRNAs and their host genes may contribute to cancer development in conjunction or independently.

The Expansion Segments of 28S Ribosomal RNA Extensively Match Human Messenger RNAs

Eukaryote ribosomal RNAs (rRNAs) have expanded in the course of phylogeny by addition of nucleotides in specific insertion areas, the expansion segments. These number about 40 in the larger (25–28S) rRNA (up to 2,400 nucleotides), and about 12 in the smaller (18S) rRNA (<700 nucleotides). Expansion of the larger rRNA shows a clear phylogenetic increase, with a dramatic rise in mammals and especially in hominids. Substantial portions of expansion segments in this RNA are not bound to ribosomal proteins, and may engage extraneous interactants, including messenger RNAs (mRNAs). Studies on the ribosome-mRNA interaction have focused on proteins of the smaller ribosomal subunit, with some examination of 18S rRNA. However, the expansion segments of human 28S rRNA show much higher density and numbers of mRNA matches than those of 18S rRNA, and also a higher density and match numbers than its own core parts. We have studied that with frequent and potentially stable matches containing 7–15 nucleotides. The expansion segments of 28S rRNA average more than 50 matches per mRNA even assuming only 5% of their sequence as available for such interaction. Large expansion segments 7, 15, and 27 of 28S rRNA also have copious long (≥10-nucleotide) matches to most human mRNAs, with frequencies much higher than in other 28S rRNA parts. Expansion segments 7 and 27 and especially segment 15 of 28S rRNA show large size increase in mammals compared to other metazoans, which could reflect a gain of function related to interaction with non-ribosomal partners. The 28S rRNA expansion segment 15 shows very high increments in size, guanosine, and cytidine nucleotide content and mRNA matching in mammals, and especially in hominids. With these segments (but not with other 28S rRNA or any 18S rRNA expansion segments) the density and number of matches are much higher in 5′-terminal than in 3′-terminal untranslated mRNA regions, which may relate to mRNA mobilization via 5′ termini. Matches in the expansion segments 7, 15, and 27 of human 28S rRNA appear as candidates for general interaction with mRNAs, especially those associated with intracellular matrices such as the endoplasmic reticulum.

Molecular species identification of Central European ground beetles (Coleoptera: Carabidae) using nuclear rDNA expansion segments and DNA barcodes

BACKGROUND

The identification of vast numbers of unknown organisms using DNA sequences becomes more and more important in ecological and biodiversity studies. In this context, a fragment of the mitochondrial cytochrome c oxidase I (COI) gene has been proposed as standard DNA barcoding marker for the identification of organisms. Limitations of the COI barcoding approach can arise from its single-locus identification system, the effect of introgression events, incomplete lineage sorting, numts, heteroplasmy and maternal inheritance of intracellular endosymbionts. Consequently, the analysis of a supplementary nuclear marker system could be advantageous.

RESULTS

We tested the effectiveness of the COI barcoding region and of three nuclear ribosomal expansion segments in discriminating ground beetles of Central Europe, a diverse and well-studied invertebrate taxon. As nuclear markers we determined the 18S rDNA: V4, 18S rDNA: V7 and 28S rDNA: D3 expansion segments for 344 specimens of 75 species. Seventy-three species (97%) of the analysed species could be accurately identified using COI, while the combined approach of all three nuclear markers provided resolution among 71 (95%) of the studied Carabidae.

CONCLUSION

Our results confirm that the analysed nuclear ribosomal expansion segments in combination constitute a valuable and efficient supplement for classical DNA barcoding to avoid potential pitfalls when only mitochondrial data are being used. We also demonstrate the high potential of COI barcodes for the identification of even closely related carabid species.

The secondary structure of the unusually long 18S ribosomal RNA of the myxozoan Sphaerospora truttae and structural evolutionary trends in the Myxozoa

The nearly complete 18S rRNA sequence of the myxozoan parasite Sphaerospora truttae shows an extraordinary length (2,552bp) in comparison with other myxozoans and with metazoans in general (average 1,800-1,900bp). The sequence shows nucleotide insertions in most variable regions of the 18S rRNA (V2, V4, V5 and V7), with especially large expansion segments in V4 and V7. In the myxozoans, nucleotide insertions and specific secondary structures in these regions of the gene were found to be strongly related to large scale phylogenetic clustering and thus with the invertebrate host type. Whereas expansion segments were generally found to be absent in the malacasporeans and the clade of primary marine myxozoan species, they occur in all taxa of the clade containing freshwater species, where they showed a consistent secondary structure throughout. The longest expansion segments occur in S. truttae, Sphaerospora elegans and Leptotheca ranae, which represent a clade that has emerged after the malacosporeans and before the radiation of all other myxozoan genera. These three species demonstrate structural links to the malacosporeans as well as other unique features. A smaller number of nucleotide insertions in different subhelices and specific secondary structures appear to have evolved independently in two marine genera, i.e. Ceratomyxa and Parvicapsula. The secondary structural elements of V4 and V7 of the myxozoan 18S rRNAs were found to be highly informative and revealed evolutionary trends of various regions of the gene hitherto unknown, since previous analyses have been based on primary sequence data excluding these regions. Furthermore, the unique features of the V4 region in S. truttae allowed for the design of a highly specific PCR assay for this species.

Identification of conserved secondary structures and expansion segments in enod40 RNAs reveals new enod40 homologues in plants

enod40 is a plant gene that participates in the regulation of symbiotic interaction between leguminous plants and bacteria or fungi. Furthermore, it has been suggested to play a general role in non-symbiotic plant development. Although enod40 seems to have multiple functions, being present in many land plants, the molecular mechanisms of its activity are unclear; they may be determined though, by short peptides and/or RNA structures encoded in the enod40 genes. We utilized conserved RNA structures in enod40 sequences to search nucleotide sequence databases and identified a number of new enod40 homologues in plant species that belong to known, but also, to yet unknown enod40-containing plant families. RNA secondary structure predictions and comparative sequence analysis of enod40 RNAs allowed us to determine the most conserved structural features, present in all known enod40 genes. Remarkably, the topology and evolution of one of the conserved structural domains are similar to those of the expansion segments found in structural RNAs such as rRNAs, RNase P and SRP RNAs. Surprisingly, the enod40 RNA structural elements are much more stronger conserved than the encoded peptides. This finding suggests that some general functions of enod40 gene could be determined by the encoded RNA structure, whereas short peptides may be responsible for more diverse functions found only in certain plant families.

Structural implications of novel diversity in eucaryal RNase P RNA

Previous eucaryotic RNase P RNA secondary structural models have been based on limited diversity, representing only two of the approximately 30 phylogenetic kingdoms of the domain Eucarya. To elucidate a more generally applicable structure, we used biochemical, bioinformatic, and molecular approaches to obtain RNase P RNA sequences from diverse organisms including representatives of six additional kingdoms of eucaryotes. Novel sequences were from acanthamoeba (Acathamoeba castellanii, Balamuthia mandrillaris, Filamoeba nolandi), animals (Caenorhabditis elegans, Drosophila melanogaster), alveolates (Theileria annulata, Babesia bovis), conosids (Dictyostelium discoideum, Physarum polycephalum), trichomonads (Trichomonas vaginalis), microsporidia (Encephalitozoon cuniculi), and diplomonads (Giardia intestinalis). An improved alignment of eucaryal RNase P RNA sequences was assembled and used for statistical and comparative structural analysis. The analysis identifies a conserved core structure of eucaryal RNase P RNA that has been maintained throughout evolution and indicates that covariation in size occurs between some structural elements of the RNA. Eucaryal RNase P RNA contains regions of highly variable length and structure reminiscent of expansion segments found in rRNA. The eucaryal RNA has been remodeled through evolution as a simplified version of the structure found in bacterial and archaeal RNase P RNAs.

Selection on the structural stability of a ribosomal RNA expansion segment in Daphnia obtusa

The high rate of sequence divergence in nuclear ribosomal RNA (rRNA) expansion segments offers a unique opportunity to study the importance of natural selection in their evolution. To this end, we polymerase chain reaction amplified and cloned a 589-nt fragment of the 18S rRNA gene containing expansion segments 43/e1 and 43/e4 from six individual Daphnia obtusa from four populations. We screened 2,588 clones using single-stranded conformation polymorphism analysis and identified 103 unique haplotype sequences. We detected two pairs of indel sites in segment 43/e4 that complement each other when the secondary structure of the linear sequence is formed. Seven of the 12 observed combinations of length variants at these four sites (haplotypes) are shared between individuals from different populations, which may suggest that some of the length variation was present in their common ancestor. Haplotypes with uncompensated indels were only observed at low frequencies, while compensated indel haplotypes were found at a wide range of frequencies, supporting the hypothesis that the energetic stability of expansion segments is a trait under natural selection. In addition, there was strong linkage disequilibrium between the four complementary indel sites, particularly those that pair with one another in the secondary structure. Despite selection against unpaired bulges at these four indel sites, some nucleotides that form unpaired bulges are highly conserved in segment 43/e4, indicating that they are under a different selective constraint, possibly due to their role in higher level structural interactions.

Secondary structure of two regions in expansion segments ES3 and ES6 with the potential of forming a tertiary interaction in eukaryotic 40S ribosomal subunits

The 18S rRNA of the small eukaryotic ribosomal subunit contains several expansion segments. Electron microscopy data indicate that two of the largest expansion segments are juxtaposed in intact 40S subunits, and data from phylogenetic sequence comparisons indicate that these two expansion segments contain complementary sequences that could form a direct tertiary interaction on the ribosome. We have investigated the secondary structure of the two expansion segments in the region around the putative tertiary interaction. Ribosomes from yeast, wheat, and mouse-three organisms representing separate eukaryotic kingdoms-were isolated, and the structure of ES3 and part of the ES6 region were analyzed using the single-strand-specific chemical reagents CMCT and DMS and the double-strand-specific ribonuclease V1. The modification patterns were analyzed by primer extension and gel electrophoresis on an ABI 377 automated DNA sequencer. The investigated sequences were relatively exposed to chemical and enzymatic modification. This is in line with their indicated location on the surface at the solvent side of the subunit. The complementary ES3 and ES6 sequences were clearly inaccessible to single-strand modification, but available for cleavage by double-strand-specific RNase V1. The results are compatible with a direct helical interaction between bases in ES3 and ES6. Almost identical results were obtained with ribosomes from the three organisms investigated.

Antisense in abundance: the ribosome as a vehicle for antisense RNA

Insertions at some sites within rRNA variable regions can be tolerated without affecting rRNA function. Antisense RNAs inserted at such sites in the T. thermophila rRNA can eliminate phenotypically or immunologically detectable gene expression of three genes tested. This unusually effective antisense activity is probably due to the abundance, stability and favourable intracellular localization of these antisense rRNAs with respect to mRNAs. Since antisense RNAs function very well as a part of the rRNA, rRNA might also be useful as a vehicle for other RNAs that might affect cell function such as protein binding sites or trans-acting ribozymes. The robust function of the antisense ribosome system in T. thermophila should allow the use of this system to specifically suppress gene expression and to clone genes by their null or hypomorphic phenotypes. The use of the antisense ribosome in other eukaryotes has yet to be explored, but the realization of this goal is well within the realm of possibility.