Complete modification maps for the cytosolic small and large subunit rRNAs of Euglena gracilis: functional and evolutionary implications of contrasting patterns between the two rRNA components

Typical structure of rRNA coding genes in diplonemids points to two independent origins of the bizarre rDNA structures of euglenozoans

Background

Members of Euglenozoa (Discoba) are known for unorthodox rDNA organization. In Euglenida rDNA is located on extrachromosomal circular DNA. In Kinetoplastea and Euglenida the core of the large ribosomal subunit, typically formed by the 28S rRNA, consists of several smaller rRNAs. They are the result of the presence of additional internal transcribed spacers (ITSs) in the rDNA. Diplonemea is the third of the main groups of Euglenozoa and its members are known to be among the most abundant and diverse protists in the oceans. Despite that, the rRNA of only one diplonemid species, Diplonema papillatum, has been examined so far and found to exhibit continuous 28S rRNA. Currently, the rDNA organization has not been researched for any diplonemid. Herein we investigate the structure of rRNA genes in classical (Diplonemidae) and deep-sea diplonemids (Eupelagonemidae), representing the majority of known diplonemid diversity. The results fill the gap in knowledge about diplonemid rDNA and allow better understanding of the evolution of the fragmented structure of the rDNA in Euglenozoa.

Results

We used available genomic (culture and single-cell) sequencing data to assemble complete or almost complete rRNA operons for three classical and six deep-sea diplonemids. The rDNA sequences acquired for several euglenids and kinetoplastids were used to provide the background for the analysis. In all nine diplonemids, 28S rRNA seems to be contiguous, with no additional ITSs detected. Similarly, no additional ITSs were detected in basal prokinetoplastids. However, we identified five additional ITSs in the 28S rRNA of all analysed metakinetoplastids, and up to twelve in euglenids. Only three of these share positions, and they cannot be traced back to their common ancestor.

Conclusions

Presented results indicate that independent origin of additional ITSs in euglenids and kinetoplastids seems to be the most likely. The reason for such unmatched fragmentation remains unknown, but for some reason euglenozoan ribosomes appear to be prone to 28S rRNA fragmentation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02014-9.

The Peptidyl Transferase Center: a Window to the Past

In his 2001 article, "Translation: in retrospect and prospect," the late Carl Woese made a prescient observation that there was a need for the then-current view of translation to be "reformulated to become an all-embracing perspective about which 21st century Biology can develop" (RNA 7:1055-1067, 2001, https://doi.org/10.1017/s1355838201010615). The quest to decipher the origins of life and the road to the genetic code are both inextricably linked with the history of the ribosome. After over 60 years of research, significant progress in our understanding of how ribosomes work has been made. Particularly attractive is a model in which the ribosome may facilitate an ∼180° rotation of the CCA end of the tRNA from the A-site to the P-site while the acceptor stem of the tRNA would then undergo a translation from the A-site to the P-site. However, the central question of how the ribosome originated remains unresolved. Along the path from a primitive RNA world or an RNA-peptide world to a proto-ribosome world, the advent of the peptidyl transferase activity would have been a seminal event. This functionality is now housed within a local region of the large-subunit (LSU) rRNA, namely, the peptidyl transferase center (PTC). The PTC is responsible for peptide bond formation during protein synthesis and is usually considered to be the oldest part of the modern ribosome. What is frequently overlooked is that by examining the origins of the PTC itself, one is likely going back even further in time. In this regard, it has been proposed that the modern PTC originated from the association of two smaller RNAs that were once independent and now comprise a pseudosymmetric region in the modern PTC. Could such an association have survived? Recent studies have shown that the extant PTC is largely depleted of ribosomal protein interactions. It is other elements like metallic ion coordination and nonstandard base/base interactions that would have had to stabilize the association of RNAs. Here, we present a detailed review of the literature focused on the nature of the extant PTC and its proposed ancestor, the proto-ribosome.

The ribosome epitranscriptome: inert-or a platform for functional plasticity?

A universal property of all rRNAs explored to date is the prevalence of post-transcriptional ("epitranscriptional") modifications, which expand the chemical and topological properties of the four standard nucleosides. Are these modifications an inert, constitutive part of the ribosome? Or could they, in part, also regulate the structure or function of the ribosome? In this review, we summarize emerging evidence that rRNA modifications are more heterogeneous than previously thought, and that they can also vary from one condition to another, such as in the context of a cellular response or a developmental trajectory. We discuss the implications of these results and key open questions on the path toward connecting such heterogeneity with function.

The 23S Ribosomal RNA From Pyrococcus furiosus Is Circularly Permuted

Synthesis and assembly of ribosomal components are fundamental cellular processes and generally well-conserved within the main groups of organisms. Yet, provocative variations to the general schemes exist. We have discovered an unusual processing pathway of pre-rRNA in extreme thermophilic archaea exemplified by Pyrococcus furiosus. The large subunit (LSU) rRNA is produced as a circularly permuted form through circularization followed by excision of Helix 98. As a consequence, the terminal domain VII that comprise the binding site for the signal recognition particle is appended to the 5´ end of the LSU rRNA that instead terminates in Domain VI carrying the Sarcin-Ricin Loop, the primary interaction site with the translational GTPases. To our knowledge, this is the first example of a true post-transcriptional circular permutation of a main functional molecule and the first example of rRNA fragmentation in archaea.

Cryo-EM structure of the highly atypical cytoplasmic ribosome of Euglena gracilis

Ribosomal RNA is the central component of the ribosome, mediating its functional and architectural properties. Here, we report the cryo-EM structure of a highly divergent cytoplasmic ribosome from the single-celled eukaryotic alga Euglena gracilis. The Euglena large ribosomal subunit is distinct in that it contains 14 discrete rRNA fragments that are assembled non-covalently into the canonical ribosome structure. The rRNA is substantially enriched in post-transcriptional modifications that are spread far beyond the catalytic RNA core, contributing to the stabilization of this highly fragmented ribosome species. A unique cluster of five adenosine base methylations is found in an expansion segment adjacent to the protein exit tunnel, such that it is positioned for interaction with the nascent peptide. As well as featuring distinctive rRNA expansion segments, the Euglena ribosome contains four novel ribosomal proteins, localized to the ribosome surface, three of which do not have orthologs in other eukaryotes.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

RNA ribose methylation (2'-O-methylation): Occurrence, biosynthesis and biological functions

Methylation of riboses at 2'-OH group is one of the most common RNA modifications found in number of cellular RNAs from almost any species which belong to all three life domains. This modification was extensively studied for decades in rRNAs and tRNAs, but recent data revealed the presence of 2'-O-methyl groups also in low abundant RNAs, like mRNAs. Ribose methylation is formed in RNA by two alternative enzymatic mechanisms: either by stand-alone protein enzymes or by complex assembly of proteins associated with snoRNA guides (sno(s)RNPs). In that case one catalytic subunit acts at various RNA sites, the specificity is provided by base pairing of the sno(s)RNA guide with the target RNA. In this review we compile available information on 2'-OH ribose methylation in different RNAs, enzymatic machineries involved in their biosynthesis and dynamics, as well as on the physiological functions of these modified residues.

RNA-Seq employing a novel rRNA depletion strategy reveals a rich repertoire of snoRNAs in Euglena gracilis including box C/D and Ψ-guide RNAs targeting the modification of rRNA extremities

Previous mRNA transcriptome studies of Euglena gracilis have shown that this organism possesses a large and diverse complement of protein coding genes; however, the study of non-coding RNA classes has been limited. The natural extensive fragmentation of the E. gracilis large subunit ribosomal RNA presents additional barriers to the identification of non-coding RNAs as size-selected small RNA libraries will be dominated by rRNA sequences. In this study we have developed a strategy to significantly reduce rRNA amplification prior to RNA-Seq analysis thereby producing a ncRNA library allowing for the identification of many new E. gracilis small RNAs. Library analysis reveals 113 unique new small nucleolar (sno) RNAs and a large collection of snoRNA isoforms, as well as the first significant collection of nuclear tRNAs in this organism. A 3' end AGAUGN consensus motif and conserved structural features can now be defined for E. gracilis pseudouridine guide RNAs. snoRNAs of both classes were identified that target modification of the 3' extremities of rRNAs utilizing predicted base-pairing interactions with internally transcribed spacers (ITS), providing insight into the timing of steps in rRNA maturation. Cumulatively, this represents the most comprehensive analysis of small ncRNAs in Euglena gracilis to date.

A tRNA methyltransferase paralog is important for ribosome stability and cell division in Trypanosoma brucei

Most eukaryotic ribosomes contain 26/28S, 5S, and 5.8S large subunit ribosomal RNAs (LSU rRNAs) in addition to the 18S rRNA of the small subunit (SSU rRNA). However, in kinetoplastids, a group of organisms that include medically important members of the genus Trypanosoma and Leishmania, the 26/28S large subunit ribosomal RNA is uniquely composed of 6 rRNA fragments. In addition, recent studies have shown the presence of expansion segments in the large ribosomal subunit (60S) of Trypanosoma brucei. Given these differences in structure, processing and assembly, T. brucei ribosomes may require biogenesis factors not found in other organisms. Here, we show that one of two putative 3-methylcytidine methyltransferases, TbMTase37 (a homolog of human methyltransferase-like 6, METTL6), is important for ribosome stability in T. brucei. TbMTase37 localizes to the nucleolus and depletion of the protein results in accumulation of ribosomal particles lacking srRNA 4 and reduced levels of polysome associated ribosomes. We also find that TbMTase37 plays a role in cytokinesis, as loss of the protein leads to multi-flagellated and multi-nucleated cells.

Identification of the methyltransferase targeting C2499 in Deinococcus radiodurans 23S ribosomal RNA

The bacterium Deinococcus radiodurans-like all other organisms-introduces nucleotide modifications into its ribosomal RNA. We have previously found that the bacterium contains a Carbon-5 methylation on cytidine 2499 of its 23S ribosomal RNA, which is so far the only modified version of cytidine 2499 reported. Using homology search, we identified the open reading frame DR_0049 as the primary candidate gene for the methyltransferase that modifies cytidine 2499. Mass spectrometric analysis demonstrated that recombinantly expressed DR0049 protein methylates E. coli cytidine 2499 both in vitro and in vivo. We also inactivated the DR_0049 gene in D. radiodurans through insertion of a chloramphenicol resistance cassette. This resulted in complete absence of the cytidine 2499 methylation, which all together demonstrates that DR_0049 encodes the methyltransferase producing m(5)C2499 in D. radiodurans 23S rRNA. Growth experiments disclosed that inactivation of DR_0049 is associated with a severe growth defect, but available ribosome structures show that cytidine 2499 is positioned very similar in D. radiodurans harbouring the modification and E. coli without the modification. Hence there is no obvious structure-based explanation for the requirement for the C2499 posttranscriptional modification in D. radiodurans.

Genome-wide analysis of small nucleolar RNAs of Leishmania major reveals a rich repertoire of RNAs involved in modification and processing of rRNA

Trypanosomatids are protozoan parasites and the causative agent of infamous infectious diseases. These organisms regulate their gene expression mainly at the post-transcriptional level and possess characteristic RNA processing mechanisms. In this study, we analyzed the complete repertoire of Leishmania major small nucleolar (snoRNA) RNAs by performing RNA-seq analysis on RNAs that were affinity-purified using the C/D snoRNA core protein, SNU13, and the H/ACA core protein, NHP2. This study revealed a large collection of C/D and H/ACA snoRNAs, organized in gene clusters generally containing both snoRNA types. Abundant snoRNAs were identified and predicted to guide trypanosome-specific rRNA cleavages. The repertoire of snoRNAs was compared to that of the closely related Trypanosoma brucei, and 80% of both C/D and H/ACA molecules were found to have functional homologues. The comparative analyses elucidated how snoRNAs evolved to generate molecules with analogous functions in both species. Interestingly, H/ACA RNAs have great flexibility in their ability to guide modifications, and several of the RNA species can guide more than one modification, compensating for the presence of single hairpin H/ACA snoRNA in these organisms. Placing the predicted modifications on the rRNA secondary structure revealed hypermodification regions mostly in domains which are modified in other eukaryotes, in addition to trypanosome-specific modifications.

Molecular Phylogeny and Cryptic Diversity of the Genus Phacus (Phacaceae, Euglenophyceae) and the Descriptions of Seven New Species

The photosynthetic euglenoid genus Phacus is commonly found in freshwater; it is characterized by a rigid to semi-rigid cell, usually flat with numerous small discoid chloroplasts without pyrenoids. To understand the phylogenetic relationships among Phacus species, we used combined cytoplasmic SSU and LSU rDNA and plastid-encoded SSU and LSU rDNA sequence data from 82 strains, including seven Lepocinclis, three Discoplastis, one Eutreptia, and two Eutreptiella strains, as well as morphological data. The combined molecular dataset was analyzed using Bayesian and maximum likelihood methods. The resulting tree revealed that the genus Phacus was not monophyletic and fully resolved the phylogenetic relationships among eight lineages that were congruent with unique morphological characters in each clade. Molecular phylogeny and detailed morphological data led to the descriptions of seven new species: P. brevisulca, P. claviformis, P. hordei-formis, P. longisulca, P. minimus, P. paraorbicularis, and P. viridioryza. The new species were well supported as independent species and formed close relationships with small Phacus species and P. orbicularis in the tree. In addition, the new species had unique molecular signatures and showed high genetic diversity. Although the strains of P. orbicularis sensu Hübner were morphologically very similar, the phylogenetic analyses and genetic diversity suggested that P. orbicularis sensu Hübner should be divided into two subclades.

Functional implications of ribosomal RNA methylation in response to environmental stress

The study of post-transcriptional RNA modifications has long been focused on the roles these chemical modifications play in maintaining ribosomal function. The field of ribosomal RNA modification has reached a milestone in recent years with the confirmation of the final unknown ribosomal RNA methyltransferase in Escherichia coli in 2012. Furthermore, the last 10 years have brought numerous discoveries in non-coding RNAs and the roles that post-transcriptional modification play in their functions. These observations indicate the need for a revitalization of this field of research to understand the role modifications play in maintaining cellular health in a dynamic environment. With the advent of high-throughput sequencing technologies, the time is ripe for leaps and bounds forward. This review discusses ribosomal RNA methyltransferases and their role in responding to external stress in Escherichia coli, with a specific focus on knockout studies and on analysis of transcriptome data with respect to rRNA methyltransferases.

Trans-splicing and RNA editing of LSU rRNA in Diplonema mitochondria

Mitochondrial ribosomal RNAs (rRNAs) often display reduced size and deviant secondary structure, and sometimes are fragmented, as are their corresponding genes. Here we report a mitochondrial large subunit rRNA (mt-LSU rRNA) with unprecedented features. In the protist Diplonema, the rnl gene is split into two pieces (modules 1 and 2, 534- and 352-nt long) that are encoded by distinct mitochondrial chromosomes, yet the rRNA is continuous. To reconstruct the post-transcriptional maturation pathway of this rRNA, we have catalogued transcript intermediates by deep RNA sequencing and RT-PCR. Gene modules are transcribed separately. Subsequently, transcripts are end-processed, the module-1 transcript is polyuridylated and the module-2 transcript is polyadenylated. The two modules are joined via trans-splicing that retains at the junction ∼26 uridines, resulting in an extent of insertion RNA editing not observed before in any system. The A-tail of trans-spliced molecules is shorter than that of mono-module 2, and completely absent from mitoribosome-associated mt-LSU rRNA. We also characterize putative antisense transcripts. Antisense-mono-modules corroborate bi-directional transcription of chromosomes. Antisense-mt-LSU rRNA, if functional, has the potential of guiding concomitantly trans-splicing and editing of this rRNA. Together, these findings open a window on the investigation of complex regulatory networks that orchestrate multiple and biochemically diverse post-transcriptional events.

Clustered organization, polycistronic transcription, and evolution of modification-guide snoRNA genes in Euglena gracilis

Previous studies have shown that the eukaryotic microbe Euglena gracilis contains an unusually large assortment of small nucleolar RNAs (snoRNAs) and ribosomal RNA (rRNA) modification sites. However, little is known about the evolutionary mechanisms contributing to this situation. In this study, we have examined the organization and evolution of snoRNA genes in Euglena with the additional objective of determining how these properties relate to the rRNA modification pattern in this protist. We have identified and extensively characterized a clustered pattern of genes encoding previously biochemically isolated snoRNA sequences in E. gracilis. We show that polycistronic transcription is a prevalent snoRNA gene expression strategy in this organism. Further, we have identified 121 new snoRNA coding regions through sequence analysis of these clusters. We have identified an E. gracilis U14 snoRNA homolog clustered with modification-guide snoRNA genes. The U14 snoRNAs in other eukaryotic organisms examined to date typically contain both a modification and a processing domain. E. gracilis U14 lacks the modification domain but retains the processing domain. Our analysis of U14 structure and evolution in Euglena and other eukaryotes allows us to propose a model for its evolution and suggest its processing role may be its more important function, explaining its conservation in many eukaryotes. The preponderance of apparent small and larger-scale duplication events in the genomic regions we have characterized in Euglena provides a mechanism for the generation of the unusually diverse collection and abundance of snoRNAs and modified rRNA sites. Our findings provide the framework for more extensive whole genome analysis to elucidate whether these snoRNA gene clusters are spread across multiple chromosomes and/or form dense "arrays" at a limited number of chromosomal loci.