RNA-level unscrambling of fragmented genes in Diplonema mitochondria

Recent expansion of metabolic versatility in Diplonema papillatum, the model species of a highly speciose group of marine eukaryotes

BACKGROUND

Diplonemid flagellates are among the most abundant and species-rich of known marine microeukaryotes, colonizing all habitats, depths, and geographic regions of the world ocean. However, little is known about their genomes, biology, and ecological role.

RESULTS

We present the first nuclear genome sequence from a diplonemid, the type species Diplonema papillatum. The ~ 280-Mb genome assembly contains about 32,000 protein-coding genes, likely co-transcribed in groups of up to 100. Gene clusters are separated by long repetitive regions that include numerous transposable elements, which also reside within introns. Analysis of gene-family evolution reveals that the last common diplonemid ancestor underwent considerable metabolic expansion. D. papillatum-specific gains of carbohydrate-degradation capability were apparently acquired via horizontal gene transfer. The predicted breakdown of polysaccharides including pectin and xylan is at odds with reports of peptides being the predominant carbon source of this organism. Secretome analysis together with feeding experiments suggest that D. papillatum is predatory, able to degrade cell walls of live microeukaryotes, macroalgae, and water plants, not only for protoplast feeding but also for metabolizing cell-wall carbohydrates as an energy source. The analysis of environmental barcode samples shows that D. papillatum is confined to temperate coastal waters, presumably acting in bioremediation of eutrophication.

CONCLUSIONS

Nuclear genome information will allow systematic functional and cell-biology studies in D. papillatum. It will also serve as a reference for the highly diverse diplonemids and provide a point of comparison for studying gene complement evolution in the sister group of Kinetoplastida, including human-pathogenic taxa.

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.

First finding of free-living representatives of Prokinetoplastina and their nuclear and mitochondrial genomes

Kinetoplastids are heterotrophic flagellated protists, including important parasites of humans and animals (trypanosomatids), and ecologically important free-living bacterial consumers (bodonids). Phylogenies have shown that the earliest-branching kinetoplastids are all parasites or obligate endosymbionts, whose highly-derived state makes reconstructing the ancestral state of the group challenging. We have isolated new strains of unusual free-living flagellates that molecular phylogeny shows to be most closely related to endosymbiotic and parasitic Perkinsela and Ichthyobodo species that, together with unidentified environmental sequences, form the clade at the base of kinetoplastids. These strains are therefore the first described free-living prokinetoplastids, and potentially very informative in understanding the evolution and ancestral states of morphological and molecular characteristics described in other kinetoplastids. Overall, we find that these organisms morphologically and ultrastructurally resemble some free-living bodonids and diplonemids, and possess nuclear genomes with few introns, polycistronic mRNA expression, high coding density, and derived traits shared with other kinetoplastids. Their genetic repertoires are more diverse than the best-studied free-living kinetoplastids, which is likely a reflection of their higher metabolic potential. Mitochondrial RNAs of these new species undergo the most extensive U insertion/deletion editing reported so far, and limited deaminative C-to-U and A-to-I editing, but we find no evidence for mitochondrial trans-splicing.

Targeted integration by homologous recombination enables in situ tagging and replacement of genes in the marine microeukaryote Diplonema papillatum

Diplonemids are a group of highly diverse and abundant marine microeukaryotes that belong to the phylum Euglenozoa and form a sister clade to the well-studied, mostly parasitic kinetoplastids. Very little is known about the biology of diplonemids, as few species have been formally described and just one, Diplonema papillatum, has been studied to a decent extent at the molecular level. Following up on our previous results showing stable but random integration of delivered extraneous DNA, we demonstrate here homologous recombination in D. papillatum. Targeting various constructs to the intended position in the nuclear genome was successful when 5' and 3' homologous regions longer than 1 kbp were used, achieving N-terminal tagging with mCherry and gene replacement of α- and β-tubulins. For more convenient genetic manipulation, we designed a modular plasmid, pDP002, which bears a protein-A tag and used it to generate and express a C-terminally tagged mitoribosomal protein. Lastly, we developed an improved transformation protocol for broader applicability across laboratories. Our robust methodology allows the replacement, integration as well as endogenous tagging of D. papillatum genes, thus opening the door to functional studies in this species and establishing a basic toolkit for reverse genetics of diplonemids in general.

Gene fragmentation and RNA editing without borders: eccentric mitochondrial genomes of diplonemids

Diplonemids are highly abundant heterotrophic marine protists. Previous studies showed that their strikingly bloated mitochondrial genome is unique because of systematic gene fragmentation and manifold RNA editing. Here we report a comparative study of mitochondrial genome architecture, gene structure and RNA editing of six recently isolated, phylogenetically diverse diplonemid species. Mitochondrial gene fragmentation and modes of RNA editing, which include cytidine-to-uridine (C-to-U) and adenosine-to-inosine (A-to-I) substitutions and 3′ uridine additions (U-appendage), are conserved across diplonemids. Yet as we show here, all these features have been pushed to their extremes in the Hemistasiidae lineage. For example, Namystynia karyoxenos has its genes fragmented into more than twice as many modules than other diplonemids, with modules as short as four nucleotides. Furthermore, we detected in this group multiple A-appendage and guanosine-to-adenosine (G-to-A) substitution editing events not observed before in diplonemids and found very rarely elsewhere. With >1,000 sites, C-to-U and A-to-I editing in Namystynia is nearly 10 times more frequent than in other diplonemids. The editing density of 12% in coding regions makes Namystynia’s the most extensively edited transcriptome described so far. Diplonemid mitochondrial genome architecture, gene structure and post-transcriptional processes display such high complexity that they challenge all other currently known systems.

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 sequence features that define efficient and specific hAGO2-dependent miRNA silencing guides

MicroRNAs (miRNAs) are ribonucleic acids (RNAs) of ∼21 nucleotides that interfere with the translation of messenger RNAs (mRNAs) and play significant roles in development and diseases. In bilaterian animals, the specificity of miRNA targeting is determined by sequence complementarity involving the seed. However, the role of the remaining nucleotides (non-seed) is only vaguely defined, impacting negatively on our ability to efficiently use miRNAs exogenously to control gene expression. Here, using reporter assays, we deciphered the role of the base pairs formed between the non-seed region and target mRNA. We used molecular modeling to reveal that this mechanism corresponds to the formation of base pairs mediated by ordered motions of the miRNA-induced silencing complex. Subsequently, we developed an algorithm based on this distinctive recognition to predict from sequence the levels of mRNA downregulation with high accuracy (r2 > 0.5, P-value < 10-12). Overall, our discovery improves the design of miRNA-guide sequences used to simultaneously downregulate the expression of multiple predetermined target genes.

Respiratory chain Complex I of unparalleled divergence in diplonemids

Mitochondrial genes of Euglenozoa (Kinetoplastida, Diplonemea, and Euglenida) are notorious for being barely recognizable, raising the question of whether such divergent genes actually code for functional proteins. Here we demonstrate the translation and identify the function of five previously unassigned y genes encoded by mitochondrial DNA (mtDNA) of diplonemids. As is the rule in diplonemid mitochondria, y genes are fragmented, with gene pieces transcribed separately and then trans-spliced to form contiguous mRNAs. Further, y transcripts undergo massive RNA editing, including uridine insertions that generate up to 16-residue-long phenylalanine tracts, a feature otherwise absent from conserved mitochondrial proteins. By protein sequence analyses, MS, and enzymatic assays in Diplonema papillatum, we show that these y genes encode the subunits Nad2, -3, -4L, -6, and -9 of the respiratory chain Complex I (CI; NADH:ubiquinone oxidoreductase). The few conserved residues of these proteins are essentially those involved in proton pumping across the inner mitochondrial membrane and in coupling ubiquinone reduction to proton pumping (Nad2, -3, -4L, and -6) and in interactions with subunits containing electron-transporting Fe-S clusters (Nad9). Thus, in diplonemids, 10 CI subunits are mtDNA-encoded. Further, MS of D. papillatum CI allowed identification of 26 conventional and 15 putative diplonemid-specific nucleus-encoded components. Most conventional accessory subunits are well-conserved but unusually long, possibly compensating for the streamlined mtDNA-encoded components and for missing, otherwise widely distributed, conventional subunits. Finally, D. papillatum CI predominantly exists as a supercomplex I:III:IV that is exceptionally stable, making this protist an organism of choice for structural studies.

A Revised Taxonomy of Diplonemids Including the Eupelagonemidae n. fam. and a Type Species, Eupelagonema oceanica n. gen. & sp

Recent surveys of marine microbial diversity have identified a previously unrecognized lineage of diplonemid protists as being among the most diverse heterotrophic eukaryotes in global oceans. Despite their monophyly (and assumed importance), they lack a formal taxonomic description, and are informally known as deep-sea pelagic diplonemids (DSPDs) or marine diplonemids. Recently, we documented morphology and molecular sequences from several DSPDs, one of which is particularly widespread and abundant in environmental sequence data. To simplify the communication of future work on this important group, here we formally propose to erect the family Eupelagonemidae to encompass this clade, as well as a formal genus and species description for the apparently most abundant phylotype, Eupelagonema oceanica, for which morphological information and single-cell amplified genome data are currently available.

Phylogeny and Morphology of New Diplonemids from Japan

Diplonemids were recently found to be the most species-rich group of marine planktonic protists. Based on phylogenetic analysis of 18S rRNA gene sequences and morphological observations, we report the description of new members of the genus Rhynchopus - R. humris sp. n. and R. serpens sp. n., and the establishment of two new genera - Lacrimia gen. n. and Sulcionema gen. n., represented by L. lanifica sp. n. and S. specki sp. n., respectively. In addition, we describe the organism formerly designated as Diplonema sp. 2 (ATCC 50224) as Flectonema neradi gen. n., sp. n. The newly described diplonemids share a common set of traits. Cells are sac-like but variable in shape and size, highly metabolic, and surrounded by a naked cell membrane, which is supported by a tightly packed corset of microtubules. They carry a single highly reticulated peripheral mitochondrion containing a large amount of mitochondrial DNA, with lamellar cristae. The cytopharyngeal complex and flagellar pocket are contiguous and have separate openings. Two parallel flagella are inserted sub-apically into a pronounced flagellar pocket. Rhynchopus species have their flagella concealed in trophic stages and fully developed in swimming stages, while they permanently protrude in all other known diplonemid species.

Transformation of Diplonema papillatum, the type species of the highly diverse and abundant marine microeukaryotes Diplonemida (Euglenozoa)

Diplonema papillatum is the type species of diplonemids, which are among the most abundant and diverse heterotrophic microeukaryotes in the world's oceans. Diplonemids are also known for a unique form of post-transcriptional processing in mitochondria. However, the lack of reverse genetics methodologies in these protists has hampered elucidation of their cellular and molecular biology. Here we report a protocol for D. papillatum transformation. We have identified several antibiotics to which D. papillatum is sensitive and thus are suitable selectable markers, and focus in particular on puromycin. Constructs were designed encoding antibiotic resistance markers, fluorescent tags, and additional genomic sequences from D. papillatum to facilitate vector integration into chromosomes. We established conditions for effective electroporation, and demonstrate that electroporated constructs can be stably integrated in the D. papillatum nuclear genome. In D. papillatum transformants, the heterologous puromycin resistance gene is transcribed into mRNA and translated into protein, as determined by Southern hybridization, reverse transcription, and Western blot analyses. This is the first documented case of transformation in a euglenozoan protist outside the well-studied kinetoplastids, making D. papillatum a genetically tractable organism and potentially a model system for marine microeukaryotes.

Keeping it complicated: Mitochondrial genome plasticity across diplonemids

Chromosome rearrangements are important drivers in genome and gene evolution, with implications ranging from speciation to development to disease. In the flagellate Diplonema papillatum (Euglenozoa), mitochondrial genome rearrangements have resulted in nearly hundred chromosomes and a systematic dispersal of gene fragments across the multipartite genome. Maturation into functional RNAs involves separate transcription of gene pieces, joining of precursor RNAs via trans-splicing, and RNA editing by substitution and uridine additions both reconstituting crucial coding sequence. How widespread these unusual features are across diplonemids is unclear. We have analyzed the mitochondrial genomes and transcriptomes of four species from the Diplonema/Rhynchopus clade, revealing a considerable genomic plasticity. Although gene breakpoints, and thus the total number of gene pieces (~80), are essentially conserved across this group, the number of distinct chromosomes varies by a factor of two, with certain chromosomes combining up to eight unrelated gene fragments. Several internal protein-coding gene pieces overlap substantially, resulting, for example, in a stretch of 22 identical amino acids in cytochrome c oxidase subunit 1 and NADH dehydrogenase subunit 5. Finally, the variation of post-transcriptional editing patterns across diplonemids indicates compensation of two adverse trends: rapid sequence evolution and loss of genetic information through unequal chromosome segregation.

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.

Post-transcriptional mending of gene sequences: Looking under the hood of mitochondrial gene expression in diplonemids

The instructions to make proteins and structural RNAs are laid down in gene sequences. Yet, in certain instances, these primary instructions need to be modified considerably during gene expression, most often at the transcript level. Here we review a case of massive post-transcriptional revisions via trans-splicing and RNA editing, a phenomenon occurring in mitochondria of a recently recognized protist group, the diplonemids. As of now, the various post-transcriptional steps have been cataloged in detail, but how these processes function is still unknown. Since genetic manipulation techniques such as gene replacement and RNA interference have not yet been established for these organisms, alternative strategies have to be deployed. Here, we discuss the experimental and bioinformatics approaches that promise to unravel the molecular machineries of trans-splicing and RNA editing in Diplonema mitochondria.

Hyper-eccentric structural genes in the mitochondrial genome of the algal parasite Hemistasia phaeocysticola

Diplonemid mitochondria are considered to have very eccentric structural genes. Coding regions of individual diplonemid mitochondrial genes are fragmented into small pieces and found on different circular DNAs. Short RNAs transcribed from each DNA molecule mature through a unique RNA maturation process involving assembly and three types of RNA editing (i.e., U insertion and A-to-I and C-to-U substitutions), although the molecular mechanism(s) of RNA maturation and the evolutionary history of these eccentric structural genes still remain to be understood. Since the gene fragmentation pattern is generally conserved among the diplonemid species studied to date, it was considered that their structural complexity has plateaued and further gene fragmentation could not occur. Here, we show the mitochondrial gene structure of Hemistasia phaeocysticola, which was recently identified as a member of a novel lineage in diplonemids, by comparison of the mitochondrial DNA sequences with cDNA sequences synthesized from mature mRNA. The genes of H. phaeocysticola are fragmented much more finely than those of other diplonemids studied to date. Furthermore, in addition to all known types of RNA editing, it is suggested that a novel processing step (i.e., secondary RNA insertion) is involved in the RNA maturation in the mitochondria of H. phaeocysticola. Our findings demonstrate the tremendous plasticity of mitochondrial gene structures.

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.

Non-functional genes repaired at the RNA level

Genomes and genes continuously evolve. Gene sequences undergo substitutions, deletions or nucleotide insertions; mobile genetic elements invade genomes and interleave in genes; chromosomes break, even within genes, and pieces reseal in reshuffled order. To maintain functional gene products and assure an organism's survival, two principal strategies are used - either repair of the gene itself or of its product. I will introduce common types of gene aberrations and how gene function is restored secondarily, and then focus on systematically fragmented genes found in a poorly studied protist group, the diplonemids. Expression of their broken genes involves restitching of pieces at the RNA-level, and substantial RNA editing, to compensate for point mutations. I will conclude with thoughts on how such a grotesquely unorthodox system may have evolved, and why this group of organisms persists and thrives since tens of millions of years.

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.

Gene Loss and Error-Prone RNA Editing in the Mitochondrion of Perkinsela, an Endosymbiotic Kinetoplastid

Perkinsela is an enigmatic early-branching kinetoplastid protist that lives as an obligate endosymbiont inside Paramoeba (Amoebozoa). We have sequenced the highly reduced mitochondrial genome of Perkinsela, which possesses only six protein-coding genes (cox1, cox2, cox3, cob, atp6, and rps12), despite the fact that the organelle itself contains more DNA than is present in either the host or endosymbiont nuclear genomes. An in silico analysis of two Perkinsela strains showed that mitochondrial RNA editing and processing machineries typical of kinetoplastid flagellates are generally conserved, and all mitochondrial transcripts undergo U-insertion/deletion editing. Canonical kinetoplastid mitochondrial ribosomes are also present. We have developed software tools for accurate and exhaustive mapping of transcriptome sequencing (RNA-seq) reads with extensive U-insertions/deletions, which allows detailed investigation of RNA editing via deep sequencing. With these methods, we show that up to 50% of reads for a given edited region contain errors of the editing system or, less likely, correspond to alternatively edited transcripts.

Uridine insertion/deletion-type RNA editing, which occurs in the mitochondrion of kinetoplastid protists, has been well-studied in the model parasite genera Trypanosoma, Leishmania, and Crithidia. Perkinsela provides a unique opportunity to broaden our knowledge of RNA editing machinery from an evolutionary perspective, as it represents the earliest kinetoplastid branch and is an obligatory endosymbiont with extensive reductive trends. Interestingly, up to 50% of mitochondrial transcripts in Perkinsela contain errors. Our study was complemented by use of newly developed software designed for accurate mapping of extensively edited RNA-seq reads obtained by deep sequencing.

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.

RiboFSM: frequent subgraph mining for the discovery of RNA structures and interactions

Frequent subgraph mining is a useful method for extracting meaningful patterns from a set of graphs or a single large graph. Here, the graph represents all possible RNA structures and interactions. Patterns that are significantly more frequent in this graph over a random graph are extracted. We hypothesize that these patterns are most likely to represent biological mechanisms. The graph representation used is a directed dual graph, extended to handle intermolecular interactions. The graph is sampled for subgraphs, which are labeled using a canonical labeling method and counted. The resulting patterns are compared to those created from a randomized dataset and scored. The algorithm was applied to the mitochondrial genome of the kinetoplastid species Trypanosoma brucei, which has a unique RNA editing mechanism. The most significant patterns contain two stem-loops, indicative of gRNA, and represent interactions of these structures with target mRNA.

A large number of nuclear genes in the human parasite blastocystis require mRNA polyadenylation to create functional termination codons

Termination codons in mRNA molecules are typically specified directly by the sequence of the corresponding gene. However, in mitochondria of a few eukaryotic groups, some mRNAs contain the termination codon UAA deriving one or both adenosines from transcript polyadenylation. Here, we show that a similar phenomenon occurs for a substantial number of nuclear genes in Blastocystis spp., divergent unicellular eukaryote gut parasites. Our analyses of published genomic data from Blastocystis sp. subtype 7 revealed that polyadenylation-mediated creation of termination codons occurs in approximately 15% of all nuclear genes. As this phenomenon has not been noticed before, the procedure previously employed to annotate the Blastocystis nuclear genome sequence failed to correctly define the structure of the 3'-ends of hundreds of genes. From sequence data we have obtained from the distantly related Blastocystis sp. subtype 1 strain, we show that this phenomenon is widespread within the Blastocystis genus. Polyadenylation in Blastocystis appears to be directed by a conserved GU-rich element located four nucleotides downstream of the polyadenylation site. Thus, the highly precise positioning of the polyadenylation in Blastocystis has allowed reduction of the 3'-untranslated regions to the point that, in many genes, only one or two nucleotides of the termination codon are left.

The pre-endosymbiont hypothesis: a new perspective on the origin and evolution of mitochondria

Mitochondrial DNA (mtDNA) is unquestionably the remnant of an α-proteobacterial genome, yet only ~10%-20% of mitochondrial proteins are demonstrably α-proteobacterial in origin (the "α-proteobacterial component," or APC). The evolutionary ancestry of the non-α-proteobacterial component (NPC) is obscure and not adequately accounted for in current models of mitochondrial origin. I propose that in the host cell that accommodated an α-proteobacterial endosymbiont, much of the NPC was already present, in the form of a membrane-bound metabolic organelle (the premitochondrion) that compartmentalized many of the non-energy-generating functions of the contemporary mitochondrion. I suggest that this organelle also possessed a protein import system and various ion and small-molecule transporters. In such a scenario, an α-proteobacterial endosymbiont could have been converted relatively directly and rapidly into an energy-generating organelle that incorporated the extant metabolic functions of the premitochondrion. This model (the "pre-endosymbiont hypothesis") effectively represents a synthesis of previous, contending mitochondrial origin hypotheses, with the bulk of the mitochondrial proteome (much of the NPC) having an endogenous origin and the minority component (the APC) having a xenogenous origin.

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.