Fragmentation of mitochondrial large subunit rRNA in the dinoflagellate Alexandrium catenella and the evolution of rRNA structure in alveolate mitochondria

Complete Mitochondrial Genomes of Ancyromonads Provide Clues for the Gene Content and Genome Structures of Ancestral Mitochondria

Mitochondria of eukaryotic cells are direct descendants of an endosymbiotic bacterium related to Alphaproteobacteria. These organelles retain their own genomes, which are highly reduced and divergent when compared to those of their bacterial relatives. To better understand the trajectory of mitochondrial genome evolution from the last eukaryotic common ancestor (LECA) to extant species, mitochondrial genome sequences from phylogenetically diverse lineages of eukaryotes-particularly protists-are essential. For this reason, we focused on the mitochondrial genomes of Ancyromonadida, an independent and understudied protist lineage in the eukaryote tree of life. Here we report the mitochondrial genomes from three Ancyromonadida: Ancyromonas sigmoides, Nutomonas longa, and Fabomonas tropica. Our analyses reveal that these mitochondrial genomes are circularly mapping molecules with inverted repeats that carry genes. This inverted repeat structure has been observed in other mitochondrial genomes but is patchily distributed over the tree of eukaryotes. Ancyromonad mitochondrial genomes possess several protein-coding genes, which have not been detected from any other mitochondrial genomes of eukaryotes sequenced to date, thereby extending the known mitochondrial gene repertoire of ancestral eukaryotes, including LECA. These findings significantly expand our understanding of mitochondrial genome diversity across eukaryotes, shedding light on the early phases of mitochondrial genome evolution.

Long-read DNA sequencing reveals the organization of the mitochondrial genome in the early-branching dinoflagellate Oxyrrhis marina

The mitochondrial genomes of dinoflagellate protists are remarkable for their highly fragmented and heterogeneous organization. Early attempts to determine their structure without 'next-generation' DNA sequencing failed to recover a defined genome. Still, it coincided in showing that the proteins coding genes, three in total, and parts of the ribosomal RNA genes were spread across a diffuse assortment of small linear fragments. In contrast, a recent study employed Illumina sequencing to assemble a 326 kbp long single-molecule, circular mitochondrial genome in the symbiotic dinoflagellate Breviolum minutum. Here, we used a combination of short- and long-read massively-parallel DNA sequencing to analyze further the mitochondrial DNA of the early-branching dinoflagellate Oxyrrhis marina. We found that the mitochondrial genome of O. marina consists of 3 linear chromosomes sized 15.9, 33.8 and 40.6 kbp for a total of 90.3 kbp. It contains the cox1, cox3 and cob genes, the same three proteins encoded in the mitochondrion of all myzozoans (Apicomplexa and Dinophyceae), some fragments of ribosomal RNA genes as well as many non-functional gene fragments and extensive noncoding DNA. Our analysis unveiled segments syntenic patterns and rearrangements encompassing coding and non-coding regions, suggesting that recombination is a pervasive process driving the evolution of these genomes.

Mitochondrial genomes revisited: why do different lineages retain different genes?

The mitochondria contain their own genome derived from an alphaproteobacterial endosymbiont. From thousands of protein-coding genes originally encoded by their ancestor, only between 1 and about 70 are encoded on extant mitochondrial genomes (mitogenomes). Thanks to a dramatically increasing number of sequenced and annotated mitogenomes a coherent picture of why some genes were lost, or relocated to the nucleus, is emerging. In this review, we describe the characteristics of mitochondria-to-nucleus gene transfer and the resulting varied content of mitogenomes across eukaryotes. We introduce a 'burst-upon-drift' model to best explain nuclear-mitochondrial population genetics with flares of transfer due to genetic drift.

Unique Properties of Apicomplexan Mitochondria

Apicomplexan parasites constitute more than 6,000 species infecting a wide range of hosts. These include important pathogens such as those causing malaria and toxoplasmosis. Their evolutionary emergence coincided with the dawn of animals. Mitochondrial genomes of apicomplexan parasites have undergone dramatic reduction in their coding capacity, with genes for only three proteins and ribosomal RNA genes present in scrambled fragments originating from both strands. Different branches of the apicomplexans have undergone rearrangements of these genes, with Toxoplasma having massive variations in gene arrangements spread over multiple copies. The vast evolutionary distance between the parasite and the host mitochondria has been exploited for the development of antiparasitic drugs, especially those used to treat malaria, wherein inhibition of the parasite mitochondrial respiratory chain is selectively targeted with little toxicity to the host mitochondria. We describe additional unique characteristics of the parasite mitochondria that are being investigated and provide greater insights into these deep-branching eukaryotic pathogens.

Intraspecies variation of the mitochondrial genome: An evaluation for phylogenetic approaches based on the conventional choices of genes and segments on mitogenome

Intraspecies nucleotide sequence variation is a key to understanding the evolutionary history of a species, such as the geographic distribution and population structure. To date, numerous phylogenetic and population genetics studies have been conducted based on the sequences of a gene or an intergenic region on the mitochondrial genome (mtDNA), such as cytochrome c oxidase subunits or the D-loop. To evaluate the credibility of the usage of such 'classic' markers, we compared the phylogenetic inferences based on the analyses of the partial and entire mtDNA sequences. Importantly, the phylogenetic reconstruction based on the short marker sequences did not necessarily reproduce the tree topologies based on the analyses of the entire mtDNA. In addition, analyses on the datasets of various organisms revealed that the analyses based on the classic markers yielded phylogenetic trees with poor confidence in all tested cases compared to the results based on full-length mtDNA. These results demonstrated that phylogenetic analyses based on complete mtDNA sequences yield more insightful results compared to those based on mitochondrial genes and segments. To ameliorate the shortcomings of the classic markers, we identified a segment of mtDNA that may be used as an 'approximate marker' to closely reproduce the phylogenetic inference obtained from the entire mtDNA in the case of mammalian species, which can be utilized to design amplicon-seq-based studies. Our study demonstrates the importance of the choice of mitochondrial markers for phylogenetic analyses and proposes a novel approach to choosing appropriate markers for mammalian mtDNA that reproduces the phylogenetic inferences obtained from full-length mtDNA.

The Elusive Mitochondrial Genomes of Apicomplexa: Where Are We Now?

Mitochondria are vital organelles of eukaryotic cells, participating in key metabolic pathways such as cellular respiration, thermogenesis, maintenance of cellular redox potential, calcium homeostasis, cell signaling, and cell death. The phylum Apicomplexa is entirely composed of obligate intracellular parasites, causing a plethora of severe diseases in humans, wild and domestic animals. These pathogens include the causative agents of malaria, cryptosporidiosis, neosporosis, East Coast fever and toxoplasmosis, among others. The mitochondria in Apicomplexa has been put forward as a promising source of undiscovered drug targets, and it has been validated as the target of atovaquone, a drug currently used in the clinic to counter malaria. Apicomplexans present a single tubular mitochondria that varies widely both in structure and in genomic content across the phylum. The organelle is characterized by massive gene migrations to the nucleus, sequence rearrangements and drastic functional reductions in some species. Recent third generation sequencing studies have reignited an interest for elucidating the extensive diversity displayed by the mitochondrial genomes of apicomplexans and their intriguing genomic features. The underlying mechanisms of gene transcription and translation are also ill-understood. In this review, we present the state of the art on mitochondrial genome structure, composition and organization in the apicomplexan phylum revisiting topological and biochemical information gathered through classical techniques. We contextualize this in light of the genomic insight gained by second and, more recently, third generation sequencing technologies. We discuss the mitochondrial genomic and mechanistic features found in evolutionarily related alveolates, and discuss the common and distinct origins of the apicomplexan mitochondria peculiarities.

Identification and phylogenetic analysis of RNA binding domain abundant in apicomplexans or RAP proteins

The RNA binding domain abundant in apicomplexans (RAP) is a protein domain identified in a diverse group of proteins, called RAP proteins, many of which have been shown to be involved in RNA binding. To understand the expansion and potential function of the RAP proteins, we conducted a hidden Markov model based screen among the proteomes of 54 eukaryotes, 17 bacteria and 12 archaea. We demonstrated that the domain is present in closely and distantly related organisms with particular expansions in Alveolata and Chlorophyta, and are not unique to Apicomplexa as previously believed. All RAP proteins identified can be decomposed into two parts. In the N-terminal region, the presence of variable helical repeats seems to participate in the specific targeting of diverse RNAs, while the RAP domain is mostly identified in the C-terminal region and is highly conserved across the different phylogenetic groups studied. Several conserved residues defining the signature motif could be crucial to ensure the function(s) of the RAP proteins. Modelling of RAP domains in apicomplexan parasites confirmed an ⍺/β structure of a restriction endonuclease-like fold. The phylogenetic trees generated from multiple alignment of RAP domains and full-length proteins from various distantly related eukaryotes indicated a complex evolutionary history of this family. We further discuss these results to assess the potential function of this protein family in apicomplexan parasites.

Intercompartmental Piecewise Gene Transfer

Gene relocation from the residual genomes of organelles to the nuclear genome still continues, although as a scaled down evolutionary phenomenon, limited in occurrence mostly to protists (sensu lato) and land plants. During this process, the structural integrity of transferred genes is usually preserved. However, the relocation of mitochondrial genes that code for respiratory chain and ribosomal proteins is sometimes associated with their fragmentation into two complementary genes. Herein, this review compiles cases of piecewise gene transfer from the mitochondria to the nucleus, and discusses hypothesized mechanistic links between the fission and relocation of those genes.

Group II Intron-Mediated Trans-Splicing in the Gene-Rich Mitochondrial Genome of an Enigmatic Eukaryote, Diphylleia rotans

Although mitochondria have evolved from a single endosymbiotic event, present day mitochondria of diverse eukaryotes display a great range of genome structures, content and features. Group I and group II introns are two features that are distributed broadly but patchily in mitochondrial genomes across branches of the tree of eukaryotes. While group I intron-mediated trans-splicing has been reported from some lineages distantly related to each other, findings of group II intron-mediated trans-splicing has been restricted to members of the Chloroplastida. In this study, we found the mitochondrial genome of the unicellular eukaryote Diphylleia rotans possesses currently the second largest gene repertoire. On the basis of a probable phylogenetic position of Diphylleia, which is located within Amorphea, current mosaic gene distribution in Amorphea must invoke parallel gene losses from mitochondrial genomes during evolution. Most notably, although the cytochrome c oxidase subunit (cox) 1 gene was split into four pieces which located at a distance to each other, we confirmed that a single mature mRNA that covered the entire coding region could be generated by group II intron-mediated trans-splicing. This is the first example of group II intron-mediated trans-splicing outside Chloroplastida. Similar trans-splicing mechanisms likely work for bipartitely split cox2 and nad3 genes to generate single mature mRNAs. We finally discuss origin and evolution of this type of trans-splicing in D. rotans as well as in eukaryotes.

The Large Mitochondrial Genome of Symbiodinium minutum Reveals Conserved Noncoding Sequences between Dinoflagellates and Apicomplexans

Even though mitochondrial genomes, which characterize eukaryotic cells, were first discovered more than 50 years ago, mitochondrial genomics remains an important topic in molecular biology and genome sciences. The Phylum Alveolata comprises three major groups (ciliates, apicomplexans, and dinoflagellates), the mitochondrial genomes of which have diverged widely. Even though the gene content of dinoflagellate mitochondrial genomes is reportedly comparable to that of apicomplexans, the highly fragmented and rearranged genome structures of dinoflagellates have frustrated whole genomic analysis. Consequently, noncoding sequences and gene arrangements of dinoflagellate mitochondrial genomes have not been well characterized. Here we report that the continuous assembled genome (∼326 kb) of the dinoflagellate, Symbiodinium minutum, is AT-rich (∼64.3%) and that it contains three protein-coding genes. Based upon in silico analysis, the remaining 99% of the genome comprises transcriptomic noncoding sequences. RNA edited sites and unique, possible start and stop codons clarify conserved regions among dinoflagellates. Our massive transcriptome analysis shows that almost all regions of the genome are transcribed, including 27 possible fragmented ribosomal RNA genes and 12 uncharacterized small RNAs that are similar to mitochondrial RNA genes of the malarial parasite, Plasmodium falciparum. Gene map comparisons show that gene order is only slightly conserved between S. minutum and P. falciparum. However, small RNAs and intergenic sequences share sequence similarities with P. falciparum, suggesting that the function of noncoding sequences has been preserved despite development of very different genome structures.

Divergent mitochondrial respiratory chains in phototrophic relatives of apicomplexan parasites

Four respiratory complexes and ATP-synthase represent central functional units in mitochondria. In some mitochondria and derived anaerobic organelles, a few or all of these respiratory complexes have been lost during evolution. We show that the respiratory chain of Chromera velia, a phototrophic relative of parasitic apicomplexans, lacks complexes I and III, making it a uniquely reduced aerobic mitochondrion. In Chromera, putative lactate:cytochrome c oxidoreductases are predicted to transfer electrons from lactate to cytochrome c, rendering complex III unnecessary. The mitochondrial genome of Chromera has the smallest known protein-coding capacity of all mitochondria, encoding just cox1 and cox3 on heterogeneous linear molecules. In contrast, another photosynthetic relative of apicomplexans, Vitrella brassicaformis, retains the same set of genes as apicomplexans and dinoflagellates (cox1, cox3, and cob).

The fragmented mitochondrial ribosomal RNAs of Plasmodium falciparum

Background

The mitochondrial genome in the human malaria parasite Plasmodium falciparum is most unusual. Over half the genome is composed of the genes for three classic mitochondrial proteins: cytochrome oxidase subunits I and III and apocytochrome b. The remainder encodes numerous small RNAs, ranging in size from 23 to 190 nt. Previous analysis revealed that some of these transcripts have significant sequence identity with highly conserved regions of large and small subunit rRNAs, and can form the expected secondary structures. However, these rRNA fragments are not encoded in linear order; instead, they are intermixed with one another and the protein coding genes, and are coded on both strands of the genome. This unorthodox arrangement hindered the identification of transcripts corresponding to other regions of rRNA that are highly conserved and/or are known to participate directly in protein synthesis.

Principal Findings

The identification of 14 additional small mitochondrial transcripts from P. falcipaurm and the assignment of 27 small RNAs (12 SSU RNAs totaling 804 nt, 15 LSU RNAs totaling 1233 nt) to specific regions of rRNA are supported by multiple lines of evidence. The regions now represented are highly similar to those of the small but contiguous mitochondrial rRNAs of Caenorhabditis elegans. The P. falciparum rRNA fragments cluster on the interfaces of the two ribosomal subunits in the three-dimensional structure of the ribosome.

Significance

All of the rRNA fragments are now presumed to have been identified with experimental methods, and nearly all of these have been mapped onto the SSU and LSU rRNAs. Conversely, all regions of the rRNAs that are known to be directly associated with protein synthesis have been identified in the P. falciparum mitochondrial genome and RNA transcripts. The fragmentation of the rRNA in the P. falciparum mitochondrion is the most extreme example of any rRNA fragmentation discovered.

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.

Systematically fragmented genes in a multipartite mitochondrial genome

Arguably, the most bizarre mitochondrial DNA (mtDNA) is that of the euglenozoan eukaryote Diplonema papillatum. The genome consists of numerous small circular chromosomes none of which appears to encode a complete gene. For instance, the cox1 coding sequence is spread out over nine different chromosomes in non-overlapping pieces (modules), which are transcribed separately and joined to a contiguous mRNA by trans-splicing. Here, we examine how many genes are encoded by Diplonema mtDNA and whether all are fragmented and their transcripts trans-spliced. Module identification is challenging due to the sequence divergence of Diplonema mitochondrial genes. By employing most sensitive protein profile search algorithms and comparing genomic with cDNA sequence, we recognize a total of 11 typical mitochondrial genes. The 10 protein-coding genes are systematically chopped up into three to 12 modules of 60–350 bp length. The corresponding mRNAs are all trans-spliced. Identification of ribosomal RNAs is most difficult. So far, we only detect the 3′-module of the large subunit ribosomal RNA (rRNA); it does not trans-splice with other pieces. The small subunit rRNA gene remains elusive. Our results open new intriguing questions about the biochemistry and evolution of mitochondrial trans-splicing in Diplonema.

Fragmentation of the large subunit ribosomal RNA gene in oyster mitochondrial genomes

Background

Discontinuous genes have been observed in bacteria, archaea, and eukaryotic nuclei, mitochondria and chloroplasts. Gene discontinuity occurs in multiple forms: the two most frequent forms result from introns that are spliced out of the RNA and the resulting exons are spliced together to form a single transcript, and fragmented gene transcripts that are not covalently attached post-transcriptionally. Within the past few years, fragmented ribosomal RNA (rRNA) genes have been discovered in bilateral metazoan mitochondria, all within a group of related oysters.

Results

In this study, we have characterized this fragmentation with comparative analysis and experimentation. We present secondary structures, modeled using comparative sequence analysis of the discontinuous mitochondrial large subunit rRNA genes of the cupped oysters C. virginica, C. gigas, and C. hongkongensis. Comparative structure models for the large subunit rRNA in each of the three oyster species are generally similar to those for other bilateral metazoans. We also used RT-PCR and analyzed ESTs to determine if the two fragmented LSU rRNAs are spliced together. The two segments are transcribed separately, and not spliced together although they still form functional rRNAs and ribosomes.

Conclusions

Although many examples of discontinuous ribosomal genes have been documented in bacteria and archaea, as well as the nuclei, chloroplasts, and mitochondria of eukaryotes, oysters are some of the first characterized examples of fragmented bilateral animal mitochondrial rRNA genes. The secondary structures of the oyster LSU rRNA fragments have been predicted on the basis of previous comparative metazoan mitochondrial LSU rRNA structure models.

Substitutional editing of Heterocapsa triquetra chloroplast transcripts and a folding model for its divergent chloroplast 16S rRNA

Substitutional editing increases genomic plasticity by changing or modifying bases at the RNA level. In this study we sequenced 10 mature chloroplast mRNAs, the chloroplast 16S rRNA and a partial chloroplast 23S rRNA from the dinoflagellate Heterocapsa triquetra, and found multiple types of substitutional editing, with A-to-G editing predominating. A-to-G editing of mRNAs converts two unusual AUA start codons into conventional AUG start codons, but three AUA start codons are not edited, showing that this dinoflagellate chloroplast has three possible start codons: AUG, AUA and UUG. To analyze the editing effects on rRNAs, we computationally predicted the secondary structure of the 16S rRNA based on the E. coli model. There are twenty editing sites in well-conserved regions of the secondary structure and eleven out of them restore conservation with other models. Moreover, A-to-G editing sites are frequently found in loop regions rather than double-stranded regions, suggesting that the A-to-G editing mechanism in dinoflagellate chloroplasts is different from that responsible for animal nuclear A-to-I(G) editing. The model of the edited 16S rRNA derived by the comparative method shares conserved secondary structural elements with other 16S rRNAs in spite of its very divergent primary sequence, supporting its role as a functional component of the chloroplast ribosome.

Dinoflagellate mitochondrial genomes: stretching the rules of molecular biology

Mitochondrial genomes represent relict bacterial genomes derived from a progenitor alpha-proteobacterium that gave rise to all mitochondria through an ancient endosymbiosis. Evolution has massively reduced these genomes, yet despite relative simplicity their organization and expression has developed considerable novelty throughout eukaryotic evolution. Few organisms have reengineered their mitochondrial genomes as thoroughly as the protist lineage of dinoflagellates. Recent work reveals dinoflagellate mitochondrial genomes as likely the most gene-impoverished of any free-living eukaryote, encoding only two to three proteins. The organization and expression of these genomes, however, is far from the simplicity their gene content would suggest. Gene duplication, fragmentation, and scrambling have resulted in an inflated and complex genome organization. Extensive RNA editing then recodes gene transcripts, and trans-splicing is required to assemble full-length transcripts for at least one fragmented gene. Even after these processes, messenger RNAs (mRNAs) lack canonical start codons and most transcripts have abandoned stop codons altogether.

Dinoflagellates: a mitochondrial genome all at sea

Dinoflagellate algae are notorious for their highly unusual organization of nuclear and chloroplast genomes. Early studies on the dinoflagellate mitochondrial genome indicated that it encodes the same three protein-coding genes found in Plasmodium spp., but with a complex organization and transcript editing. Recent work has extended this view, showing that the dinoflagellate mitochondrial genome contains a wide array of gene fragments and genes interspersed with noncoding inverted repeats. The genome seems to require noncanonical start and stop codons, as well as high levels of editing, trans-splicing and the addition of oligonucleotide caps at the 5' and 3' ends of transcripts. Despite its small coding content, the dinoflagellate mitochondrial genome is one of the most complex known.

Broad genomic and transcriptional analysis reveals a highly derived genome in dinoflagellate mitochondria

Background

Dinoflagellates comprise an ecologically significant and diverse eukaryotic phylum that is sister to the phylum containing apicomplexan endoparasites. The mitochondrial genome of apicomplexans is uniquely reduced in gene content and size, encoding only three proteins and two ribosomal RNAs (rRNAs) within a highly compacted 6 kb DNA. Dinoflagellate mitochondrial genomes have been comparatively poorly studied: limited available data suggest some similarities with apicomplexan mitochondrial genomes but an even more radical type of genomic organization. Here, we investigate structure, content and expression of dinoflagellate mitochondrial genomes.

Results

From two dinoflagellates, Crypthecodinium cohnii and Karlodinium micrum, we generated over 42 kb of mitochondrial genomic data that indicate a reduced gene content paralleling that of mitochondrial genomes in apicomplexans, i.e., only three protein-encoding genes and at least eight conserved components of the highly fragmented large and small subunit rRNAs. Unlike in apicomplexans, dinoflagellate mitochondrial genes occur in multiple copies, often as gene fragments, and in numerous genomic contexts. Analysis of cDNAs suggests several novel aspects of dinoflagellate mitochondrial gene expression. Polycistronic transcripts were found, standard start codons are absent, and oligoadenylation occurs upstream of stop codons, resulting in the absence of termination codons. Transcripts of at least one gene, cox3, are apparently trans-spliced to generate full-length mRNAs. RNA substitutional editing, a process previously identified for mRNAs in dinoflagellate mitochondria, is also implicated in rRNA expression.

Conclusion

The dinoflagellate mitochondrial genome shares the same gene complement and fragmentation of rRNA genes with its apicomplexan counterpart. However, it also exhibits several unique characteristics. Most notable are the expansion of gene copy numbers and their arrangements within the genome, RNA editing, loss of stop codons, and use of trans-splicing.

The dinoflagellates Durinskia baltica and Kryptoperidinium foliaceum retain functionally overlapping mitochondria from two evolutionarily distinct lineages

Background

The dinoflagellates Durinskia baltica and Kryptoperidinium foliaceum are distinguished by the presence of a tertiary plastid derived from a diatom endosymbiont. The diatom is fully integrated with the host cell cycle and is so altered in structure as to be difficult to recognize it as a diatom, and yet it retains a number of features normally lost in tertiary and secondary endosymbionts, most notably mitochondria. The dinoflagellate host is also reported to retain mitochondrion-like structures, making these cells unique in retaining two evolutionarily distinct mitochondria. This redundancy raises the question of whether the organelles share any functions in common or have distributed functions between them.

Results

We show that both host and endosymbiont mitochondrial genomes encode genes for electron transport proteins. We have characterized cytochrome c oxidase 1 (cox1), cytochrome oxidase 2 (cox2), cytochrome oxidase 3 (cox3), cytochrome b (cob), and large subunit of ribosomal RNA (LSUrRNA) of endosymbiont mitochondrial ancestry, and cox1 and cob of host mitochondrial ancestry. We show that all genes are transcribed and that those ascribed to the host mitochondrial genome are extensively edited at the RNA level, as expected for a dinoflagellate mitochondrion-encoded gene. We also found evidence for extensive recombination in the host mitochondrial genes and that recombination products are also transcribed, as expected for a dinoflagellate.

Conclusion

Durinskia baltica and K. foliaceum retain two mitochondria from evolutionarily distinct lineages, and the functions of these organelles are at least partially overlapping, since both express genes for proteins in electron transport.

The highly reduced and fragmented mitochondrial genome of the early-branching dinoflagellate Oxyrrhis marina shares characteristics with both apicomplexan and dinoflagellate mitochondrial genomes

The mitochondrial genome and the expression of the genes within it have evolved to be highly unusual in several lineages. Within alveolates, apicomplexans and dinoflagellates share the most reduced mitochondrial gene content on record, but differ from one another in organisation and function. To clarify how these characteristics originated, we examined mitochondrial genome form and expression in a key lineage that arose close to the divergence of apicomplexans and dinoflagellates, Oxyrrhis marina. We show that Oxyrrhis is a basal member of the dinoflagellate lineage whose mitochondrial genome has some unique characteristics while sharing others with apicomplexans or dinoflagellates. Specifically, Oxyrrhis has the smallest gene complement known, with several rRNA fragments and only two protein coding genes, cox1 and a cob-cox3 fusion. The genome appears to be highly fragmented, like that of dinoflagellates, but genes are frequently arranged as tandem copies, reminiscent of the repeating nature of the Plasmodium genome. In dinoflagellates and Oxyrrhis, genes are found in many arrangements, but the Oxyrrhis genome appears to be more structured, since neighbouring genes or gene fragments are invariably the same: cox1 and the cob-cox3 fusion were never found on the same genomic fragment. Analysing hundreds of cDNAs for both genes and circularized mRNAs from cob-cox3 showed that neither uses canonical start or stop codons, although a UAA terminator is created in the cob-cox3 fusion mRNA by post-transcriptional oligoadenylation. mRNAs from both genes also use a novel 5' oligo(U) cap. Extensive RNA editing is characteristic of dinoflagellates, but we find no editing in Oxyrrhis. Overall, the combination of characteristics found in the Oxyrrhis genome allows us to plot the sequence of many events that led to the extreme organisation of apicomplexan and dinoflalgellate mitochondrial genomes.