Fragmented and scrambled mitochondrial ribosomal RNA coding regions among green algae: a model for their origin and evolution

Mitochondrial ribosomal RNA coding regions in the only three green algal taxa investigated to date are fundamentally different in that they are continuous in Prototheca wickerhamii, but highly fragmented and scrambled in Chlamydomonas reinhardtii and Chlamydomonas eugametos. To gain more insight into the mode of evolution of fragmented and scrambled mitochondrial ribosomal RNA (rRNA) genes within the green algal group, this work (1) provides additional information on fragmentation patterns of mitochondrial small- and large-subunit (SSU and LSU) rRNAs that strongly supports the concept of a gradual increase in the extent of discontinuity of mitochondrial rRNAs among chlorophycean green algae and (2) reports the first example of fragmented and scrambled mitochondrial LSU rRNA coding regions in a green algal taxon outside the Chlamydomonas group. The present study (1) suggests that the scrambling of the mitochondrial rRNA coding regions may have occurred early in the evolution of fragmented and scrambled mitochondrial rRNA genes within the chlorophycean green algal group, most likely in parallel with the fragmentation events, (2) proposes recombination as a possible mechanism involved in the evolution of these mitochondrial rRNA genes, and (3) presents a hypothetical pathway for converting continuous mitochondrial rRNA genes into the highly fragmented and scrambled rRNA coding regions of Chlamydomonas through a series of recombinatorial events between short repeated sequences.      

DISCONTINUOUS MITOCHONDRIAL AND CHLOROPLAST LARGE SUBUNIT RIBOSOMAL RNAs AMONG GREEN ALGAE: PHYLOGENETIC IMPLICATIONS1

To provide insights into the occurrence, evolution, and phylogenetic distribution of discontinuous mitochondrial and chloroplast large subunit ribosomal RNAs (LSU rRNAs) among green algae, we surveyed 12 taxa representing three classes of green algae: the Chlorophyceae, Pleurastrophyceae, and Micromonadophyceae (sensu Mattox and Stewart 1984). We present evidence that discontinuous mitochondrial and chloroplast LSU rRNAs are quite widespread among green algae. Mitochondrial LSU rRNAs appear discontinuous in zoosporic chlorophycean lineages displaying a clockwise or directly opposed configuration in their flagellar apparatus, as well as in chlorococcalean autosporic taxa phylogenetically related to them, but are continuous among zoosporic green algal lineages with a counterclockwise flagellar apparatus configuration, as well as among chlorococcalean autosporic taxa phylogenetically related to them. Chloroplast LSU rRNAs appear discontinuous in all of the lineages investigated. Discontinuous mitochondrial LSU rRNA represents a molecular trait that might have originated at or near the base of Chlorophyceae, whereas discontinuous chloroplast LSU rRNA might have developed very early in the evolutionary history of the green algal group itself. We suggest, therefore, that the presence of discontinuous mitochondrial but not chloroplast LSU rRNA can be used as an additional character in assessing phylogenetic affiliations among green algae.     

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

In the protist Euglena gracilis, the cytosolic small subunit (SSU) rRNA is a single, covalently continuous species typical of most eukaryotes; in contrast, the large subunit (LSU) rRNA is naturally fragmented, comprising 14 separate RNA molecules instead of the bipartite (28S + 5.8S) eukaryotic LSU rRNA typically seen. We present extensively revised secondary structure models of the E. gracilis SSU and LSU rRNAs and have mapped the positions of all of the modified nucleosides in these rRNAs (88 in SSU rRNA and 262 in LSU rRNA, with only 3 LSU rRNA modifications incompletely characterized). The relative proportions of ribose-methylated nucleosides and pseudouridine (∼ 60% and ∼ 35%, respectively) are closely similar in the two rRNAs; however, whereas the Euglena SSU rRNA has about the same absolute number of modifications as its human counterpart, the Euglena LSU rRNA has twice as many modifications as the corresponding human LSU rRNA. The increased levels of rRNA fragmentation and modification in E. gracilis LSU rRNA are correlated with a 3-fold increase in the level of mispairing in helical regions compared to the human LSU rRNA. In contrast, no comparable increase in mispairing is seen in helical regions of the SSU rRNA compared to its homologs in other eukaryotes. In view of the reported effects of both ribose-methylated nucleoside and pseudouridine residues on RNA structure, these correlations lead us to suggest that increased modification in the LSU rRNA may play a role in stabilizing a ‘looser’ structure promoted by elevated helical mispairing and a high degree of fragmentation.     

A 7.1 kb linear DNA molecule of Theileria parva has scrambled rDNA sequences and open reading frames for mitochondrially encoded proteins.

Theileria parva, an intralymphocytic protozoan parasite of cattle, contains a linear 7.1 kb DNA element with terminal inverted repeat sequences. The molecule is transcribed into low molecular weight RNA, and both DNA strands encode short stretches of unique sequences, usually < 100 nucleotides, which are similar to large (LSU) or small (SSU) ribosomal subunit RNA. Phylogenetically conserved conformational rRNA domains were assembled from the discontinuous rDNA sequences using comparative secondary structure modelling. For example, a minimum of four predicted sequences, two derived from each DNA strand, is required to assemble domain V of LSU rRNA which participates in peptidyl transferase activity. The discontinuities in the identified rRNA domains fall within regions of no known functional significance. Hence, it is likely that the element encodes fragmented rDNA genes and the mature rRNA is unconventional, consisting of several fragments of RNA, primarily held together by intermolecular and intramolecular base pairing. The element also has ORFs for components of the last two mitochondrial electron transport enzyme complexes. The structure of the parasite DNA element, its protein coding capacity and scrambled rDNA gene sequences, are reminiscent of the mitochondrial genome of Chlamydomonas reinhardtii. We propose that the 7.1 kb element is equivalent to the mitochondrial DNA of T. parva, although a number of its features are unusual for this family of extrachromosomal DNA molecules.     

DNA sequence and secondary structure of the mitochondrial small subunit ribosomal RNA coding region including a group-IC2 intron from the cultivated basidiomycete Agrocybe aegerita

Due to their structural complexity and their evolutionary dimension, rRNAs are the most investigated nucleic acids in prokaryotes, eukaryotes and organelles. However, no complete sequence of a mitochondrial small subunit (SSU) rRNA was available in the basidiomycotina subdivision. The mitochondrial gene encoding the SSU rRNA of the cultivated basidiomycete Agrocybe aegerita was cloned and its complete nucleotide sequence achieved; the 5′- and 3′-ends were localized by nuclease S1 mapping, leading to a size of 3277 nt. The secondary structure of the SSU rRNA (1906 nt in size) possessed all the helices and loops of the prokaryotic model; a unique modification was found in a conserved nucleotide predicted by the model: the nt 487 was A instead of C. The same modification, has been found in all the partial basidiomycete mitochondrial sequences available in databases. The Agrocybe aegerita SSU rRNA was characterized by large and unusual extensions leading to additional helices in the variable domains V4, V6 and V9, which were the longest of the known prokaryotic or mitochondrial SSU rRNAs. Nucleotide sequence analysis indicated a 1371-bp intron, belonging to subgroup-IC2, located in a conserved loop in the 3′-part of the SSU rRNA. This intron, which is the second example reported in a fungal mitochondrial SSU rDNA, encoded a putative protein (407 aa) sharing homologies with endonucleases involved in group-I intron mobility. This report constitutes the first complete mitochondrial SSU rRNA sequence and secondary structure of any member of the basidiomycotina subdivision.     

Identities and phylogenetic comparisons of posttranscriptional modifications in 16 S ribosomal RNA from Haloferax volcanii

Small subunit (16 S) rRNA from the archaeon Haloferax volcanii, for which sites of modification were previously reported, was examined using mass spectrometry. A census of all modified residues was taken by liquid chromatography/electrospray ionization-mass spectrometry analysis of a total nucleoside digest of the rRNA. Following rRNA hydrolysis by RNase T(1), accurate molecular mass values of oligonucleotide products were measured using liquid chromatography/electrospray ionization-mass spectrometry and compared with values predicted from the corresponding gene sequence. Three modified nucleosides, distributed over four conserved sites in the decoding region of the molecule, were characterized: 3-(3-amino-3-carboxypropyl)uridine-966, N(6)-methyladenosine-1501, and N(6),N(6)-dimethyladenosine-1518 and -1519 (all Escherichia coli numbering). Nucleoside 3-(3-amino-3-carboxypropyl)uridine, previously unknown in rRNA, occurs at a highly conserved site of modification in all three evolutionary domains but for which no structural assignment in archaea has been previously reported. Nucleoside N(6)-methyladenosine, not previously placed in archaeal rRNAs, frequently occurs at the analogous location in eukaryotic small subunit rRNA but not in bacteria. H. volcanii small subunit rRNA appears to reflect the phenotypically low modification level in the Crenarchaeota kingdom and is the only cytoplasmic small subunit rRNA shown to lack pseudouridine.     

Sequence and organization of large subunit rRNA genes from the extrachromosomal 35 kb circular DNA of the malaria parasite Plasmodium falciparum.

The malaria parasite Plasmodium falciparum carries an extrachromosomal 35 kb circular DNA molecule of unknown provenance. A striking feature of the circle is a palindromic sequence of genes for subunit rRNAs and several tRNAs, spanning ca. 10.5 kb. The palindrome has an intriguing resemblance to the inverted repeat of plastid genomes, and the sequence and putative secondary structure of the malarial large subunit (LSU) rRNA described in this report were used as the basis of a phylogenetic study. The malarial rRNA was found to be highly divergent in comparison with a selected group of chloroplast LSU rRNAs but was more closely related to them than to mitochondrial LSU rRNA genes.     

Hyaloraphidium curvatum: a linear mitochondrial genome, tRNA editing, and an evolutionary link to lower fungi

We have sequenced the mitochondrial DNA (mtDNA) of Hyaloraphidium curvatum, an organism previously classified as a colorless green alga but now recognized as a lower fungus based on molecular data. The 29.97-kbp mitochondrial chromosome is maintained as a monomeric, linear molecule with identical, inverted repeats (1.43 kbp) at both ends, a rare genome architecture in mitochondria. The genome encodes only 14 known mitochondrial proteins, 7 tRNAs, the large subunit rRNA and small subunit rRNA (SSU rRNA), and 3 ORFs. The SSU rRNA is encoded in two gene pieces that are located 8 kbp apart on the mtDNA. Scrambled and fragmented mitochondrial rRNAs are well known from green algae and alveolate protists but are unprecedented in fungi. Protein genes code for apocytochrome b; cytochrome oxidase 1, 2, and 3, NADH dehydrogenase 1, 2, 3, 4, 4L, 5, and 6, and ATP synthase 6, 8, and 9 subunits, and several of these genes are organized in operon-like clusters. The set of seven mitochondrially encoded tRNAs is insufficient to recognize all codons that occur in the mitochondrial protein genes. When taking into account the pronounced codon bias, at least 16 nuclear-encoded tRNAs are assumed to be imported into the mitochondria. Three of the seven predicted mitochondria-encoded tRNA sequences carry mispairings in the first three positions of the acceptor stem. This strongly suggests that these tRNAs are edited by a mechanism similar to the one seen in the fungus Spizellomyces punctatus and the rhizopod amoeba Acanthamoeba castellanii. Our phylogenetic analysis confirms with overwhelming support that H. curvatum is a member of the chytridiomycete fungi, specifically related to the Monoblepharidales.     

Mitochondrial genome of the colorless green alga Polytomella parva: two linear DNA molecules with homologous inverted repeat Termini

Most of the well-characterized mitochondrial genomes from diverse green algal lineages are circular mapping DNA molecules; however, Chlamydomonas reinhardtii has a linear 15.8 kb unit mitochondrial genome with 580 or 581 bp inverted repeat ends. In mitochondrial-enriched fractions prepared from Polytomella parva (=P. agilis), a colorless, naturally wall-less relative of C. reinhardtii, we have detected two linear mitochondrial DNA (mtDNA) components with sizes of 13.5 and 3.5 kb. Sequences spanning 97% and 86% of the 13.5- and 3.5-kb mtDNAs, respectively, reveal that these molecules contain long, at least 1.3 kb, homologous inverted repeat sequences at their termini. The 3.5-kb mtDNA has only one coding region (nad6), the functionality of which is supported by both the relative rate at which it has accumulated nonsynonymous nucleotide substitutions and its absence from the 13.5-kb mtDNA which encodes nine genes (i.e., large and small subunit rRNA [LSU and SSU rRNA] genes, one tRNA gene, and six protein-coding genes). On the basis of DNA sequence data, we propose that a variant start codon, GTG, is utilized by the P. parva 13.5-kb mtDNA-encoded gene, nad5. Using the relative rate test with Chlamydomonas moewusii (=C. eugametos) as the outgroup, we conclude that the nonsynonymous nucleotide substitution rate in the mitochondrial protein-coding genes of P. parva is on an average about 3.3 times that of the C. reinhardtii counterparts.     

NSUN4 Is a Dual Function Mitochondrial Protein Required for Both Methylation of 12S rRNA and Coordination of Mitoribosomal Assembly

Biogenesis of mammalian mitochondrial ribosomes requires a concerted maturation of both the small (SSU) and large subunit (LSU). We demonstrate here that the m⁵C methyltransferase NSUN4, which forms a complex with MTERF4, is essential in mitochondrial ribosomal biogenesis as mitochondrial translation is abolished in conditional Nsun4 mouse knockouts. Deep sequencing of bisulfite-treated RNA shows that NSUN4 methylates cytosine 911 in 12S rRNA (m5C911) of the SSU. Surprisingly, NSUN4 does not need MTERF4 to generate this modification. Instead, the NSUN4/MTERF4 complex is required to assemble the SSU and LSU to form a monosome. NSUN4 is thus a dual function protein, which on the one hand is needed for 12S rRNA methylation and, on the other hand interacts with MTERF4 to facilitate monosome assembly. The presented data suggest that NSUN4 has a key role in controlling a final step in ribosome biogenesis to ensure that only the mature SSU and LSU are assembled.     

Characterization of mitochondrial ribosomal RNA genes in gadiformes: sequence variations,secondary structural features,and phylogenetic implications

Secondary structure features of mitochondrial ribosomal RNAs (mt-rRNAs) of bony fishes were investigated by a DNA sequence alignment approach. The small subunit (SSU) and large subunit (LSU) mt-rRNA genes were found to contain several additional variable regions compared to their mammalian counterparts. Fish mt-LSU rRNA genes were found to be longer than the mammalians due to increased length of some of the variable regions. The 5' and 3' ends of Atlantic cod mt-rRNAs were precisely mapped. The 3' ends of mt-SSU rRNAs were found to be homogenous and mono-adenylated, whereas that of the mt-LSU rRNAs were heterogenous and oligo-adenylated. The 5' ends of mt-SSU rRNAs appeared to be heterogenous, corresponding to the presumed first and second positions of the gene. Sequences of the central domain and the D-domain of the mt-SSU and mt-LSU rRNA genes, respectively, were determined and characterized for 11 gadiform species (representing the families Gadidae, Lotidae, Ranicipitidae, Merlucciidae, Phycidae, and Macrouridae) and one Lophiidae species. Detailed secondary structure models of the RNA regions are presented for the Atlantic cod (Gadus morhua) and Roundnose grenadier (Coryphaeonides rupestris). Saturation plots revealed that DNA nucleotide positions corresponding to unpaired RNA regions become saturated with transitions at sequence divergence levels about 0.15. Phylogenetic analyses revealed some aspects of gadiform relationships. Gadidae was identified as the most derived of the gadiform families. Lotidae was found to be the family closest related to Gadidae, and Ranicipitidae was also recognized as a derived gadiform taxon.     

Combined heat shock protein 90 and ribosomal RNA sequence phylogeny supports multiple replacements of dinoflagellate plastids

Dinoflagellates harbour diverse plastids obtained from several algal groups, including haptophytes, diatoms, cryptophytes, and prasinophytes. Their major plastid type with the accessory pigment peridinin is found in the vast majority of photosynthetic species. Some species of dinoflagellates have other aberrantly pigmented plastids. We sequenced the nuclear small subunit (SSU) ribosomal RNA (rRNA) gene of the "green" dinoflagellate Gymnodinium chlorophorum and show that it is sister to Lepidodinium viride, indicating that their common ancestor obtained the prasinophyte (or other green alga) plastid in one event. As the placement of dinoflagellate species that acquired green algal or haptophyte plastids is unclear from small and large subunit (LSU) rRNA trees, we tested the usefulness of the heat shock protein (Hsp) 90 gene for dinoflagellate phylogeny by sequencing it from four species with aberrant plastids (G. chlorophorum, Karlodinium micrum, Karenia brevis, and Karenia mikimotoi) plus Alexandrium tamarense, and constructing phylogenetic trees for Hsp90 and rRNAs, separately and together. Analyses of the Hsp90 and concatenated data suggest an ancestral origin of the peridinin-containing plastid, and two independent replacements of the peridinin plastid soon after the early radiation of the dinoflagellates. Thus, the Hsp90 gene seems to be a promising phylogenetic marker for dinoflagellate phylogeny.     

metaxa2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data

The ribosomal rRNA genes are widely used as genetic markers for taxonomic identification of microbes. Particularly the small subunit (SSU; 16S/18S) rRNA gene is frequently used for species‐ or genus‐level identification, but also the large subunit (LSU; 23S/28S) rRNA gene is employed in taxonomic assignment. The metaxa software tool is a popular utility for extracting partial rRNA sequences from large sequencing data sets and assigning them to an archaeal, bacterial, nuclear eukaryote, mitochondrial or chloroplast origin. This study describes a comprehensive update to metaxa – metaxa 2 – that extends the capabilities of the tool, introducing support for the LSU rRNA gene, a greatly improved classifier allowing classification down to genus or species level, as well as enhanced support for short‐read (100 bp) and paired‐end sequences, among other changes. The performance of metaxa 2 was compared to other commonly used taxonomic classifiers, showing that metaxa 2 often outperforms previous methods in terms of making correct predictions while maintaining a low misclassification rate. metaxa 2 is freely available from http://microbiology.se/software/metaxa2/ .