Ribosomal RNA genes in Euglena gracilis mitochondrial DNA: fragmented genes in a seemingly fragmented genome

Circular extrachromosomal DNA in Euglena gracilis under normal and stress conditions

Extrachromosomal circular DNA (eccDNA) enhances genomic plasticity, augmenting its coding and regulatory potential. Advances in high-throughput sequencing have enabled the investigation of these structural variants. Although eccDNAs have been investigated in numerous taxa, they remained understudied in euglenids. Therefore, we examined eccDNAs predicted from Illumina sequencing data of Euglena gracilis Z SAG 1224-5/25, grown under optimal photoperiod and exposed to UV irradiation. We identified approximately 1000 unique eccDNA candidates, about 20% of which were shared across conditions. We also observed a significant enrichment of mitochondrially encoded eccDNA in the UV-irradiated sample. Furthermore, we found that the heterogeneity of eccDNA was reduced in UV-exposed samples compared to cells that were grown in optimal conditions. Hence, eccDNA appears to play a role in the response to oxidative stress in Euglena, as it does in other studied organisms. In addition to contributing to the understanding of Euglena genomes, our results contribute to the validation of bioinformatics pipelines on a large, non-model genome.

Origin of minicircular mitochondrial genomes in red algae

Eukaryotic organelle genomes are generally of conserved size and gene content within phylogenetic groups. However, significant variation in genome structure may occur. Here, we report that the Stylonematophyceae red algae contain multipartite circular mitochondrial genomes (i.e., minicircles) which encode one or two genes bounded by a specific cassette and a conserved constant region. These minicircles are visualized using fluorescence microscope and scanning electron microscope, proving the circularity. Mitochondrial gene sets are reduced in these highly divergent mitogenomes. Newly generated chromosome-level nuclear genome assembly of Rhodosorus marinus reveals that most mitochondrial ribosomal subunit genes are transferred to the nuclear genome. Hetero-concatemers that resulted from recombination between minicircles and unique gene inventory that is responsible for mitochondrial genome stability may explain how the transition from typical mitochondrial genome to minicircles occurs. Our results offer inspiration on minicircular organelle genome formation and highlight an extreme case of mitochondrial gene inventory reduction.

Dissecting Trypanosoma brucei RRP44 function in the maturation of segmented ribosomal RNA using a regulated genetic complementation system

Trypanosoma brucei belongs to a group of protozoans presenting fragmented large subunit rRNA. Its LSU rRNA equivalent to the 25S/28S rRNA of other eukaryotes is split into six fragments, requiring additional processing for removal of the extra spacer sequences. We have used a genetic complementation strategy to further investigate the T. brucei RRP44 nuclease in pre-rRNA maturation. TbRRP44 contains both a PIN and a RNB domain whose homologues are found in association with the exosome complex. We found that the exonucleolytic activity of the RNB domain as well as the physical presence of the PIN domain are essential for TbRRP44 function, while a catalytic site mutation in the PIN domain has no detectable effect on cell growth. A new endonucleolytic cleavage site in ITS1 was identified. In addition to the 5.8S rRNA 3'-end maturation, TbRRP44 is required for degradation of the excised 5'-ETS and for removal of part of ITS1 during maturation of the 18S rRNA 3'-end. TbRRP44 deficiency leads to accumulation of many LSU intermediate precursors, most of them not detected in control cells. TbRRP44 is also required for U3 snoRNA and spliced leader processing, indicating that TbRRP44 may have a wide role in RNA processing in T. brucei.

De Novo Transcriptome Meta-Assembly of the Mixotrophic Freshwater Microalga Euglena gracilis

Euglena gracilis is a well-known photosynthetic microeukaryote considered as the product of a secondary endosymbiosis between a green alga and a phagotrophic unicellular belonging to the same eukaryotic phylum as the parasitic trypanosomatids. As its nuclear genome has proven difficult to sequence, reliable transcriptomes are important for functional studies. In this work, we assembled a new consensus transcriptome by combining sequencing reads from five independent studies. Based on a detailed comparison with two previously released transcriptomes, our consensus transcriptome appears to be the most complete so far. Remapping the reads on it allowed us to compare the expression of the transcripts across multiple culture conditions at once and to infer a functionally annotated network of co-expressed genes. Although the emergence of meaningful gene clusters indicates that some biological signal lies in gene expression levels, our analyses confirm that gene regulation in euglenozoans is not primarily controlled at the transcriptional level. Regarding the origin of E. gracilis, we observe a heavily mixed gene ancestry, as previously reported, and rule out sequence contamination as a possible explanation for these observations. Instead, they indicate that this complex alga has evolved through a convoluted process involving much more than two partners.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans

Background

The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered.

Results

We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome.

Conclusions

Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01035-y.

Phage Origin of Mitochondrion-Localized Family A DNA Polymerases in Kinetoplastids and Diplonemids

Mitochondria retain their own genomes as other bacterial endosymbiont-derived organelles. Nevertheless, no protein for DNA replication and repair is encoded in any mitochondrial genomes (mtDNAs) assessed to date, suggesting that the nucleus primarily governs the maintenance of mtDNA. As the proteins of diverse evolutionary origins occupy a large proportion of the current mitochondrial proteomes, we anticipate finding the same evolutionary trend in the nucleus-encoded machinery for mtDNA maintenance. Indeed, none of the DNA polymerases (DNAPs) in the mitochondrial endosymbiont, a putative α-proteobacterium, seemingly had been inherited by their descendants (mitochondria), as none of the known types of mitochondrion-localized DNAP showed a specific affinity to the α-proteobacterial DNAPs. Nevertheless, we currently have no concrete idea of how and when the known types of mitochondrion-localized DNAPs emerged. We here explored the origins of mitochondrion-localized DNAPs after the improvement of the samplings of DNAPs from bacteria and phages/viruses. Past studies have revealed that a set of mitochondrion-localized DNAPs in kinetoplastids and diplonemids, namely PolIB, PolIC, PolID, PolI-Perk1/2, and PolI-dipl (henceforth designated collectively as “PolIBCD+”) have emerged from a single DNAP. In this study, we recovered an intimate connection between PolIBCD+ and the DNAPs found in a particular group of phages. Thus, the common ancestor of kinetoplastids and diplonemids most likely converted a laterally acquired phage DNAP into a mitochondrion-localized DNAP that was ancestral to PolIBCD+. The phage origin of PolIBCD+ hints at a potentially large contribution of proteins acquired via nonvertical processes to the machinery for mtDNA maintenance in kinetoplastids and diplonemids.

Inventory and Evolution of Mitochondrion-localized Family A DNA Polymerases in Euglenozoa

The order Trypanosomatida has been well studied due to its pathogenicity and the unique biology of the mitochondrion. In Trypanosoma brucei, four DNA polymerases, namely PolIA, PolIB, PolIC, and PolID, related to bacterial DNA polymerase I (PolI), were shown to be localized in mitochondria experimentally. These mitochondrion-localized DNA polymerases are phylogenetically distinct from other family A DNA polymerases, such as bacterial PolI, DNA polymerase gamma (Polγ) in human and yeasts, “plant and protist organellar DNA polymerase (POP)” in diverse eukaryotes. However, the diversity of mitochondrion-localized DNA polymerases in Euglenozoa other than Trypanosomatida is poorly understood. In this study, we discovered putative mitochondrion-localized DNA polymerases in broad members of three major classes of Euglenozoa—Kinetoplastea, Diplonemea, and Euglenida—to explore the origin and evolution of trypanosomatid PolIA-D. We unveiled distinct inventories of mitochondrion-localized DNA polymerases in the three classes: (1) PolIA is ubiquitous across the three euglenozoan classes, (2) PolIB, C, and D are restricted in kinetoplastids, (3) new types of mitochondrion-localized DNA polymerases were identified in a prokinetoplastid and diplonemids, and (4) evolutionarily distinct types of POP were found in euglenids. We finally propose scenarios to explain the inventories of mitochondrion-localized DNA polymerases in Kinetoplastea, Diplonemea, and Euglenida.

Piece by piece: Building a ribozyme

The ribosome and RNase P are cellular ribonucleoprotein complexes that perform peptide bond synthesis and phosphodiester bond cleavage, respectively. Both are ancient biological assemblies that were already present in the last universal common ancestor of all life. The large subunit rRNA in the ribosome and the RNA subunit of RNase P are the ribozyme components required for catalysis. Here, we explore the idea that these two large ribozymes may have begun their evolutionary odyssey as an assemblage of RNA "fragments" smaller than the contemporary full-length versions and that they transitioned through distinct stages along a pathway that may also be relevant for the evolution of other non-coding RNAs.

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.

Perfection of eccentricity: Mitochondrial genomes of diplonemids

Mitochondria are the sandbox of evolution as exemplified most particularly by the diplonemids, a group of marine microeukaryotes. These protists are uniquely characterized by their highly multipartite mitochondrial genome and systematically fragmented genes whose pieces are spread out over several dozens of chromosomes. The type species Diplonema papillatum was the first member of this group in which the expression of fragmented mitochondrial genes was investigated experimentally. We now know that gene expression involves separate transcription of gene pieces (modules), RNA editing of module transcripts, and module joining to mature mRNAs and rRNAs. The mechanism of cognate module recognition and ligation is distinct from known intron splicing and remains to be uncovered. Here, we review the current status of research on mitochondrial genome architecture, as well as gene complement, structure, and expression modes in diplonemids. Further, we discuss the potential molecular mechanisms of posttranscriptional processing, and finally reflect on the evolutionary trajectories and trends of mtDNA evolution as seen in this protist group. © 2018 IUBMB Life, 70(12):1197-1206, 2018.

The Mitochondrion of Euglena gracilis

In the presence of oxygen, Euglena gracilis mitochondria function much like mammalian mitochondria. Under anaerobiosis, E. gracilis mitochondria perform a malonyl-CoA independent synthesis of fatty acids leading to accumulation of wax esters, which serve as the sink for electrons stemming from glycolytic ATP synthesis and pyruvate oxidation. Some components (enzymes and cofactors) of Euglena's anaerobic energy metabolism are found among the anaerobic mitochondria of invertebrates, others are found among hydrogenosomes, the H₂-producing anaerobic mitochondria of protists.

Mitochondrial Genome Evolution and a Novel RNA Editing System in Deep-Branching Heteroloboseids

Discoba (Excavata) is an evolutionarily important group of eukaryotes that includes Jakobida, with the most bacterial-like mitochondrial genomes known, and Euglenozoa, many of which have extensively fragmented mitochondrial genomes. However, little is known about the mitochondrial genomes of Heterolobosea, the third main group of Discoba. Here, we studied two heteroloboseids—an undescribed amoeba “BB2” and Pharyngomonas kirbyi. Phylogenomic analysis revealed that they form a clade that is a sister group to all other Heterolobosea. We characterized the mitochondrial genomes of BB2 and P. kirbyi, which encoded 44 and 48 putative protein-coding genes respectively. Their gene contents were similar to that of Naegleria. In BB2, mitochondrially encoded RNAs were heavily edited, with ∼500 mononucleotide insertion events, mostly guanosines. These insertions always have the same identity as an adjacent nucleotide. Editing occurs in all ribosomal RNAs and protein-coding transcripts except one, and half of the transfer RNAs. Analysis of Illumina deep-sequencing data suggested that this RNA editing is very accurate and efficient, and most likely co-transcriptional. The dissimilarity of this editing process to other RNA editing phenomena in discobids, as well as its apparent absence in P. kirbyi, suggest that this remarkably extensive system of insertional editing evolved independently in the BB2 lineage, after its divergence from the P. kirbyi lineage.

Euglena gracilis Genome and Transcriptome: Organelles, Nuclear Genome Assembly Strategies and Initial Features

Euglena gracilis is a major component of the aquatic ecosystem and together with closely related species, is ubiquitous worldwide. Euglenoids are an important group of protists, possessing a secondarily acquired plastid and are relatives to the Kinetoplastidae, which themselves have global impact as disease agents. To understand the biology of E. gracilis, as well as to provide further insight into the evolution and origins of the Kinetoplastidae, we embarked on sequencing the nuclear genome; the plastid and mitochondrial genomes are already in the public domain. Earlier studies suggested an extensive nuclear DNA content, with likely a high degree of repetitive sequence, together with significant extrachromosomal elements. To produce a list of coding sequences we have combined transcriptome data from both published and new sources, as well as embarked on de novo sequencing using a combination of 454, Illumina paired end libraries and long PacBio reads. Preliminary analysis suggests a surprisingly large genome approaching 2 Gbp, with a highly fragmented architecture and extensive repeat composition. Over 80% of the RNAseq reads from E. gracilis maps to the assembled genome sequence, which is comparable with the well assembled genomes of T. brucei and T. cruzi. In order to achieve this level of assembly we employed multiple informatics pipelines, which are discussed here. Finally, as a preliminary view of the genome architecture, we discuss the tubulin and calmodulin genes, which highlight potential novel splicing mechanisms.

Genes in Hiding

Unrecognizable genes are an unsettling problem in genomics. Here, we survey the various types of cryptic genes and the corresponding deciphering strategies employed by cells. Encryption that renders genes substantially different from homologs in other species includes sequence substitution, insertion, deletion, fragmentation plus scrambling, and invasion by mobile genetic elements. Cells decode cryptic genes at the DNA, RNA or protein level. We will focus on a recently discovered case of unparalleled encryption involving massive gene fragmentation and nucleotide deletions and substitutions, occurring in the mitochondrial genome of a poorly understood protist group, the diplonemids. This example illustrates that comprehensive gene detection requires not only auxiliary sequence information - transcriptome and proteome data - but also knowledge about a cell's deciphering arsenal.

Protists and the Wild, Wild West of Gene Expression: New Frontiers, Lawlessness, and Misfits

The DNA double helix has been called one of life's most elegant structures, largely because of its universality, simplicity, and symmetry. The expression of information encoded within DNA, however, can be far from simple or symmetric and is sometimes surprisingly variable, convoluted, and wantonly inefficient. Although exceptions to the rules exist in certain model systems, the true extent to which life has stretched the limits of gene expression is made clear by nonmodel systems, particularly protists (microbial eukaryotes). The nuclear and organelle genomes of protists are subject to the most tangled forms of gene expression yet identified. The complicated and extravagant picture of the underlying genetics of eukaryotic microbial life changes how we think about the flow of genetic information and the evolutionary processes shaping it. Here, we discuss the origins, diversity, and growing interest in noncanonical protist gene expression and its relationship to genomic architecture.

From simple to supercomplex: mitochondrial genomes of euglenozoan protists

Mitochondria are double membrane organelles of endosymbiotic origin, best known for constituting the centre of energetics of a eukaryotic cell. They contain their own mitochondrial genome, which as a consequence of gradual reduction during evolution typically contains less than two dozens of genes. In this review, we highlight the extremely diverse architecture of mitochondrial genomes and mechanisms of gene expression between the three sister groups constituting the phylum Euglenozoa - Euglenida, Diplonemea and Kinetoplastea. The earliest diverging euglenids possess a simplified mitochondrial genome and a conventional gene expression, whereas both are highly complex in the two other groups. The expression of their mitochondrial-encoded proteins requires extensive post-transcriptional modifications guided by complex protein machineries and multiple small RNA molecules. Moreover, the least studied diplonemids, which have been recently discovered as a highly abundant component of the world ocean plankton, possess one of the most complicated mitochondrial genome organisations known to date.

From simple to supercomplex: mitochondrial genomes of euglenozoan protists

Mitochondria are double membrane organelles of endosymbiotic origin, best known for constituting the centre of energetics of a eukaryotic cell. They contain their own mitochondrial genome, which as a consequence of gradual reduction during evolution typically contains less than two dozens of genes. In this review, we highlight the extremely diverse architecture of mitochondrial genomes and mechanisms of gene expression between the three sister groups constituting the phylum Euglenozoa - Euglenida, Diplonemea and Kinetoplastea. The earliest diverging euglenids possess a simplified mitochondrial genome and a conventional gene expression, whereas both are highly complex in the two other groups. The expression of their mitochondrial-encoded proteins requires extensive post-transcriptional modifications guided by complex protein machineries and multiple small RNA molecules. Moreover, the least studied diplonemids, which have been recently discovered as a highly abundant component of the world ocean plankton, possess one of the most complicated mitochondrial genome organisations known to date.

Novel modes of RNA editing in mitochondria

Gene structure and expression in diplonemid mitochondria are unparalleled. Genes are fragmented in pieces (modules) that are separately transcribed, followed by the joining of module transcripts to contiguous RNAs. Some instances of unique uridine insertion RNA editing at module boundaries were noted, but the extent and potential occurrence of other editing types remained unknown. Comparative analysis of deep transcriptome and genome data from Diplonema papillatum mitochondria reveals ∼220 post-transcriptional insertions of uridines, but no insertions of other nucleotides nor deletions. In addition, we detect in total 114 substitutions of cytosine by uridine and adenosine by inosine, amassed into unusually compact clusters. Inosines in transcripts were confirmed experimentally. This is the first report of adenosine-to-inosine editing of mRNAs and ribosomal RNAs in mitochondria. In mRNAs, editing causes mostly amino-acid additions and non-synonymous substitutions; in ribosomal RNAs, it permits formation of canonical secondary structures. Two extensively edited transcripts were compared across four diplonemids. The pattern of uridine-insertion editing is strictly conserved, whereas substitution editing has diverged dramatically, but still rendering diplonemid proteins more similar to other eukaryotic orthologs. We posit that RNA editing not only compensates but also sustains, or even accelerates, ultra-rapid evolution of genome structure and sequence in diplonemid mitochondria.

Aerobic mitochondria of parasitic protists: Diverse genomes and complex functions

In this review the main features of the mitochondria of aerobic parasitic protists are discussed. While the best characterized organelles are by far those of kinetoplastid flagellates and Plasmodium, we also consider amoebae Naegleria and Acanthamoeba, a ciliate Ichthyophthirius and related lineages. The simplistic view of the mitochondrion as just a power house of the cell has already been abandoned in multicellular organisms and available data indicate that this also does not apply for protists. We discuss in more details the following mitochondrial features: genomes, post-transcriptional processing, translation, biogenesis of iron-sulfur complexes, heme metabolism and the electron transport chain. Substantial differences in all these core mitochondrial features between lineages are compatible with the view that aerobic protists harbor organelles that are more complex and flexible than previously appreciated.

Unexpectedly Streamlined Mitochondrial Genome of the Euglenozoan Euglena gracilis

In this study, we describe the mitochondrial genome of the excavate flagellate Euglena gracilis. Its gene complement is reduced as compared with the well-studied sister groups Diplonemea and Kinetoplastea. We have identified seven protein-coding genes: Three subunits of respiratory complex I (nad1, nad4, and nad5), one subunit of complex III (cob), and three subunits of complex IV (cox1, cox2, and a highly divergent cox3). Moreover, fragments of ribosomal RNA genes have also been identified. Genes encoding subunits of complex V, ribosomal proteins and tRNAs were missing, and are likely located in the nuclear genome. Although mitochondrial genomes of diplonemids and kinetoplastids possess the most complex RNA processing machineries known, including trans-splicing and editing of the uridine insertion/deletion type, respectively, our transcriptomic data suggest their total absence in E. gracilis. This finding supports a scenario in which the complex mitochondrial processing machineries of both sister groups evolved relatively late in evolution from a streamlined genome and transcriptome of their common predecessor.

Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes

Mitochondrial and plastid genomes show a wide array of architectures, varying immensely in size, structure, and content. Some organelle DNAs have even developed elaborate eccentricities, such as scrambled coding regions, nonstandard genetic codes, and convoluted modes of posttranscriptional modification and editing. Here, we compare and contrast the breadth of genomic complexity between mitochondrial and plastid chromosomes. Both organelle genomes have independently evolved many of the same features and taken on similar genomic embellishments, often within the same species or lineage. This trend is most likely because the nuclear-encoded proteins mediating these processes eventually leak from one organelle into the other, leading to a high likelihood of processes appearing in both compartments in parallel. However, the complexity and intensity of genomic embellishments are consistently more pronounced for mitochondria than for plastids, even when they are found in both compartments. We explore the evolutionary forces responsible for these patterns and argue that organelle DNA repair processes, mutation rates, and population genetic landscapes are all important factors leading to the observed convergence and divergence in organelle genome architecture.

How amino acids and peptides shaped the RNA world

The “RNA world” hypothesis is seen as one of the main contenders for a viable theory on the origin of life. Relatively small RNAs have catalytic power, RNA is everywhere in present-day life, the ribosome is seen as a ribozyme, and rRNA and tRNA are crucial for modern protein synthesis. However, this view is incomplete at best. The modern protein-RNA ribosome most probably is not a distorted form of a “pure RNA ribosome” evolution started out with. Though the oldest center of the ribosome seems “RNA only”, we cannot conclude from this that it ever functioned in an environment without amino acids and/or peptides. Very small RNAs (versatile and stable due to basepairing) and amino acids, as well as dipeptides, coevolved. Remember, it is the amino group of aminoacylated tRNA that attacks peptidyl-tRNA, destroying the bond between peptide and tRNA. This activity of the amino acid part of aminoacyl-tRNA illustrates the centrality of amino acids in life. With the rise of the “RNA world” view of early life, the pendulum seems to have swung too much towards the ribozymatic part of early biochemistry. The necessary presence and activity of amino acids and peptides is in need of highlighting. In this article, we try to bring the role of the peptide component of early life back into focus. We argue that an RNA world completely independent of amino acids never existed.

Gene content evolution in Discobid mitochondria deduced from the phylogenetic position and complete mitochondrial genome of Tsukubamonas globosa

The unicellular eukaryotic assemblage Discoba (Excavata) comprises four lineages: the Heterolobosea, Euglenozoa, Jakobida, and Tsukubamonadida. Discoba has been considered as a key assemblage for understanding the early evolution of mitochondrial (mt) genomes, as jakobids retain the most gene-rich (i.e., primitive) genomes compared with any other eukaryotes determined to date. However, to date, mt genome sequences have been completed for only a few groups within Discoba, including jakobids, two closely related heteroloboseans, and kinetoplastid euglenozoans. The Tsukubamonadida is the least studied lineage, as the order was only recently established with the description of a sole representative species, Tsukubamonas globosa. The evolutionary relationship between T. globosa and other discobids has yet to be resolved, and no mt genome data are available for this particular organism. Here, we use a “phylogenomic” approach to resolve the relationship between T. globosa, heteroloboseans, euglenozoans, and jakobids. In addition, we have characterized the mt genome of T. globosa (48,463 bp in length), which encodes 52 putative protein-coding and 29 RNA genes. By mapping the gene repertoires of discobid mt genomes onto the well-resolved Discoba tree, we model gene loss events during the evolution of discobid mt genomes.

Organelle evolution, fragmented rRNAs, and Carl

I am honored to have been asked to contribute to this memorial issue, although I cannot claim to have known Carl Woese well. Carl's insights and the discoveries that his research group made over the years certainly stimulated my own research program, and at several points early on, interactions with him were pivotal in my career. Here I comment on these personal dealings with Carl and emphasize his influence in two areas of long-standing interest in my lab: organelle evolution and rRNA evolution.

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.

RNA-dependent genome processing during nuclear differentiation: the model systems of stichotrichous ciliates

We introduce ciliated protozoa, and more specifically the stichotrichous ciliates Oxytricha and Stylonychia, as biological model systems for the analysis of programmed DNA-reorganization processes during nuclear differentiation. These include DNA excision, DNA elimination, reordering of gene segments and specific gene amplification. We show that small nuclear RNAs specify DNA sequences to be excised or retained, but also discuss the need for a RNA template molecule derived from the parental nucleus for these processes. This RNA template guides reordering of gene segments to become functional genes and determines gene copy number in the differentiated nucleus. Since the template is derived from the parental macronucleus, gene reordering and DNA amplification are inherited in a non-Mendelian epigenetic manner.

Diversity of mitochondrial genome organization

In this review, we discuss types of mitochondrial genome structural organization (architecture), which includes the following characteristic features: size and the shape of DNA molecule, number of encoded genes, presence of cryptogenes, and editing of primary transcripts.

The mitochondrial genome of the arbuscular mycorrhizal fungus Gigaspora margarita reveals two unsuspected trans-splicing events of group I introns

• Arbuscular mycorrhizal fungi (AMF) are ubiquitous organisms that benefit ecosystems through the establishment of an association with the roots of most plants: the mycorrhizal symbiosis. Despite their ecological importance, however, these fungi have been poorly studied at the genome level. • In this study, total DNA from the AMF Gigaspora margarita was subjected to a combination of 454 and Illumina sequencing, and the resulting reads were used to assemble its mitochondrial genome de novo. This genome was annotated and compared with those of other relatives to better comprehend the evolution of the AMF lineage. • The mitochondrial genome of G. margarita is unique in many ways, exhibiting a large size (97 kbp) and elevated GC content (45%). This genome also harbors molecular events that were previously unknown to occur in fungal mitochondrial genomes, including trans-splicing of group I introns from two different genes coding for the first subunit of the cytochrome oxidase and for the small subunit of the rRNA. • This study reports the second published genome from an AMF organelle, resulting in relevant DNA sequence information from this poorly studied fungal group, and providing new insights into the frequency, origin and evolution of trans-spliced group I introns found across the mitochondrial genomes of distantly related organisms.

Unscrambling genetic information at the RNA level

Genomics aims at unraveling the blueprint of life; however, DNA sequence alone does not always reveal the proteins and structural RNAs encoded by the genome. The reason is that genetic information is often encrypted. Recognizing the logic of encryption, and understanding how living cells decode hidden information--at the level of DNA, RNA or protein--is challenging. RNA-level decryption includes topical RNA editing and more 'macroscopic' transcript rearrangements. The latter events involve the four types of introns recognized to date, notably spliceosomal, group I, group II, and archaeal/tRNA splicing. Intricate variants, such as alternative splicing and trans-splicing, have been reported for each intron type, but the biological significance has not always been confirmed. Novel RNA-level unscrambling processes were recently discovered in mitochondria of dinoflagellates and diplonemids, and potentially euglenids. These processes seem not to rely on known introns, and the corresponding molecular mechanisms remain to be elucidated.

How a neutral evolutionary ratchet can build cellular complexity

Complex cellular machines and processes are commonly believed to be products of selection, and it is typically understood to be the job of evolutionary biologists to show how selective advantage can account for each step in their origin and subsequent growth in complexity. Here, we describe how complex machines might instead evolve in the absence of positive selection through a process of "presuppression," first termed constructive neutral evolution (CNE) more than a decade ago. If an autonomously functioning cellular component acquires mutations that make it dependent for function on another, pre-existing component or process, and if there are multiple ways in which such dependence may arise, then dependence inevitably will arise and reversal to independence is unlikely. Thus, CNE is a unidirectional evolutionary ratchet leading to complexity, if complexity is equated with the number of components or steps necessary to carry out a cellular process. CNE can explain "functions" that seem to make little sense in terms of cellular economy, like RNA editing or splicing, but it may also contribute to the complexity of machines with clear benefit to the cell, like the ribosome, and to organismal complexity overall. We suggest that CNE-based evolutionary scenarios are in these and other cases less forced than the selectionist or adaptationist narratives that are generally told.

The mitochondrial genome and transcriptome of the basal dinoflagellate Hematodinium sp.: character evolution within the highly derived mitochondrial genomes of dinoflagellates

The sister phyla dinoflagellates and apicomplexans inherited a drastically reduced mitochondrial genome (mitochondrial DNA, mtDNA) containing only three protein-coding (cob, cox1, and cox3) genes and two ribosomal RNA (rRNA) genes. In apicomplexans, single copies of these genes are encoded on the smallest known mtDNA chromosome (6 kb). In dinoflagellates, however, the genome has undergone further substantial modifications, including massive genome amplification and recombination resulting in multiple copies of each gene and gene fragments linked in numerous combinations. Furthermore, protein-encoding genes have lost standard stop codons, trans-splicing of messenger RNAs (mRNAs) is required to generate complete cox3 transcripts, and extensive RNA editing recodes most genes. From taxa investigated to date, it is unclear when many of these unusual dinoflagellate mtDNA characters evolved. To address this question, we investigated the mitochondrial genome and transcriptome character states of the deep branching dinoflagellate Hematodinium sp. Genomic data show that like later-branching dinoflagellates Hematodinium sp. also contains an inflated, heavily recombined genome of multicopy genes and gene fragments. Although stop codons are also lacking for cox1 and cob, cox3 still encodes a conventional stop codon. Extensive editing of mRNAs also occurs in Hematodinium sp. The mtDNA of basal dinoflagellate Hematodinium sp. indicates that much of the mtDNA modification in dinoflagellates occurred early in this lineage, including genome amplification and recombination, and decreased use of standard stop codons. Trans-splicing, on the other hand, occurred after Hematodinium sp. diverged. Only RNA editing presents a nonlinear pattern of evolution in dinoflagellates as this process occurs in Hematodinium sp. but is absent in some later-branching taxa indicating that this process was either lost in some lineages or developed more than once during the evolution of the highly unusual dinoflagellate mtDNA.

RNA-Mediated Epigenetic Programming of Genome Rearrangements

RNA, normally thought of as a conduit in gene expression, has a novel mode of action in ciliated protozoa. Maternal RNA templates provide both an organizing guide for DNA rearrangements and a template that can transport somatic mutations to the next generation. This opportunity for RNA-mediated genome rearrangement and DNA repair is profound in the ciliate Oxytricha, which deletes 95% of its germline genome during development in a process that severely fragments its chromosomes and then sorts and reorders the hundreds of thousands of pieces remaining. Oxytricha's somatic nuclear genome is therefore an epigenome formed through RNA templates and signals arising from the previous generation. Furthermore, this mechanism of RNA-mediated epigenetic inheritance can function across multiple generations, and the discovery of maternal template RNA molecules has revealed new biological roles for RNA and has hinted at the power of RNA molecules to sculpt genomic information in cells.

Gene fragmentation: a key to mitochondrial genome evolution in Euglenozoa?

Phylum Euglenozoa comprises three groups of eukaryotic microbes (kinetoplastids, diplonemids, and euglenids), the mitochondrial (mt) genomes of which exhibit radically different modes of organization and expression. Gene fragmentation is a striking feature of both euglenid and diplonemid mtDNAs. To rationalize the emergence of these highly divergent mtDNA types and the existence of insertion/deletion RNA editing (in kinetoplastids) and trans-splicing (in diplonemids), we propose that in the mitochondrion of the common evolutionary ancestor of Euglenozoa, small expressed gene fragments promoted a rampant neutral evolutionary pathway. Interactions between small antisense transcripts of these gene fragments and full-length transcripts, assisted by RNA-processing enzymes, permitted the emergence of RNA editing and/or trans-splicing activities, allowing the system to tolerate indel mutations and further gene fragmentation, respectively, and leading to accumulation of additional mutations. In this way, dramatically different mitochondrial genome structures and RNA-processing machineries were able to evolve. The paradigm of constructive neutral evolution acting on the widely different mitochondrial genetic systems in Euglenozoa posits the accretion of initially neutral molecular interactions by genetic drift, leading inevitably to the observed 'irremediable complexity'.

Evolutionarily conserved cox1 trans-splicing without cis-motifs

In the protist Diplonema papillatum (Diplonemea, Euglenozoa), mitochondrial genes are systematically fragmented with each nonoverlapping piece (module) encoded individually on a distinct circular chromosome. Gene modules are transcribed separately, and precursor transcripts are assembled to mature mRNA by a trans-splicing process of yet unknown mechanism. Expression of the cox1 gene that consists of nine modules, also involves RNA editing by which six uridines are added between Modules 4 and 5. Here, we investigate whether the unusual features of cox1 are shared by all Diplonemea and what the mechanism of trans-splicing might be. We examine three additional species representing both Diplonemea genera, namely D. papillatum described before, and D. ambulator, Diplonema sp.2, and Rhynchopus euleeides and discover that in all Diplonemea, the cox1 gene is discontinuous and split up into nine modules that each reside on a distinct chromosome. Positions of gene breakpoints vary by up to two nucleotides. Further, all taxa have six nonencoded uridines inserted in cox1 mRNA at exactly the same position as D. papillatum. In silico searches do not detect signatures of introns known to engage in trans-splicing, in particular Group I, Group II, spliceosomal, and transfer RNA introns. Nor did we find statistically significant reverse-complementary motifs between adjacent modules and their flanking regions, or residues conserved within or across species. This provides compelling evidence that trans-splicing in Diplonemea mitochondria does not rely on sequence elements in cis but rather proceeds by a mechanism employing matchmaking trans factors, such as RNAs or proteins.