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.

A functional bacteria-derived restriction modification system in the mitochondrion of a heterotrophic protist.

The overarching trend in mitochondrial genome evolution is functional streamlining coupled with gene loss. Therefore, gene acquisition by mitochondria is considered to be exceedingly rare. Selfish elements in the form of self-splicing introns occur in many organellar genomes, but the wider diversity of selfish elements, and how they persist in the DNA of organelles, has not been explored. In the mitochondrial genome of a marine heterotrophic katablepharid protist, we identify a functional type II restriction modification (RM) system originating from a horizontal gene transfer (HGT) event involving bacteria related to flavobacteria. This RM system consists of an HpaII-like endonuclease and a cognate cytosine methyltransferase (CM). We demonstrate that these proteins are functional by heterologous expression in both bacterial and eukaryotic cells. These results suggest that a mitochondrion-encoded RM system can function as a toxin-antitoxin selfish element, and that such elements could be co-opted by eukaryotic genomes to drive biased organellar inheritance.

Depletion of a Toxoplasma porin leads to defects in mitochondrial morphology and contacts with the endoplasmic reticulum.

The voltage-dependent anion channel (VDAC) is a ubiquitous channel in the outer membrane of the mitochondrion with multiple roles in protein, metabolite and small molecule transport. In mammalian cells, VDAC protein, as part of a larger complex including the inositol triphosphate receptor, has been shown to have a role in mediating contacts between the mitochondria and endoplasmic reticulum (ER). We identify VDAC of the pathogenic apicomplexan Toxoplasma gondii and demonstrate its importance for parasite growth. We show that VDAC is involved in protein import and metabolite transfer to mitochondria. Further, depletion of VDAC resulted in significant morphological changes in the mitochondrion and ER, suggesting a role in mediating contacts between these organelles in T. gondii. This article has an associated First Person interview with the first author of the paper.

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.

Constructive Neutral Evolution 20 Years Later.

Evolution has led to a great diversity that ranges from elegant simplicity to ornate complexity. Many complex features are often assumed to be more functional or adaptive than their simpler alternatives. However, in 1999, Arlin Stolzfus published a paper in the Journal of Molecular Evolution that outlined a framework in which complexity can arise through a series of non-adaptive steps. He called this framework Constructive Neutral Evolution (CNE). Despite its two-decade-old roots, many evolutionary biologists still appear to be unaware of this explanatory framework for the origins of complexity. In this perspective piece, we explain the theory of CNE and how it changes the order of events in narratives that describe the evolution of complexity. We also provide an extensive list of cellular features that may have become more complex through CNE. We end by discussing strategies to determine whether complexity arose through neutral or adaptive processes.

The Origin of Mitochondrial Cristae from Alphaproteobacteria.

Mitochondria are the respiratory organelles of eukaryotes and their evolutionary history is deeply intertwined with that of eukaryotes. The compartmentalization of respiration in mitochondria occurs within cristae, whose evolutionary origin has remained unclear. Recent discoveries, however, have revived the old notion that mitochondrial cristae could have had a pre-endosymbiotic origin. Mitochondrial cristae are likely homologous to the intracytoplasmic membranes (ICMs) used by diverse alphaproteobacteria for harnessing energy. Because the Mitochondrial Contact site and Cristae Organizing System (MICOS) that controls the development of cristae evolved from a simplified version that is phylogenetically restricted to Alphaproteobacteria (alphaMICOS), ICMs most probably transformed into cristae during the endosymbiotic origin of mitochondria. This inference is supported by the sequence and structural similarities between MICOS and alphaMICOS, and the expression pattern and cellular localization of alphaMICOS. Given that cristae and ICMs develop similarly, alphaMICOS likely functions analogously to mitochondrial MICOS by culminating ICM development with the creation of tubular connections and membrane contact sites at the alphaproteobacterial envelope. Mitochondria thus inherited a pre-existing ultrastructure adapted to efficient energy transduction from their alphaproteobacterial ancestors. The widespread nature of purple bacteria among alphaproteobacteria raises the possibility that cristae evolved from photosynthetic ICMs.

Evolutionary conservation of a core fungal phosphate homeostasis pathway coupled to development in Blastocladiella emersonii.

The model yeast Saccharomyces cerevisiae elicits a transcriptional response to phosphate (Pi) depletion. To determine the origins of the phosphate response (PHO) system, we bioinformatically identified putative PHO components in the predicted proteomes of diverse fungi. Our results suggest that the PHO system is ancient; however, components have been expanded or lost in different fungal lineages. To show that a similar physiological response is present in deeply-diverging fungi we examined the transcriptional and physiological response of PHO genes to Pi depletion in the blastocladiomycete Blastocladiella emersonii. Our physiological experiments indicate that B. emersonii relies solely on high-affinity Na+-independent Pho84-like transporters. In response to Pi depletion, BePho84 paralogues were 4-8-fold transcriptionally upregulated, whereas several other PHO homologues like phosphatases and vacuolar transporter chaperone (VTC) complex components show 2-3-fold transcriptional upregulation. Since Pi has been shown to be important during the development of B. emersonii, we sought to determine if PHO genes are differentially regulated at different lifecycle stages. We demonstrate that a similar set of PHO transporters and phosphatases are upregulated at key points during B. emersonii development. Surprisingly, some genes upregulated during Pi depletion, including VTC components, are repressed at these key stages of development indicating that PHO genes are regulated by different pathways in different developmental and environmental situations. Overall, our findings indicate that a complex PHO network existed in the ancient branches of the fungi, persists in diverse extant fungi, and that this ancient network is likely to be involved in development and cell cycle regulation.

The evolution of ERMIONE in mitochondrial biogenesis and lipid homeostasis: An evolutionary view from comparative cell biology.

The ER-mitochondria organizing network (ERMIONE) in Saccharomyces cerevisiae is involved in maintaining mitochondrial morphology and lipid homeostasis. ERMES and MICOS are two scaffolding complexes of ERMIONE that contribute to these processes. ERMES is ancient but has been lost in several lineages including animals, plants, and SAR (stramenopiles, alveolates and rhizaria). On the other hand, MICOS is ancient and has remained present in all organisms bearing mitochondrial cristae. The ERMIONE precursor evolved in the α-proteobacterial ancestor of mitochondria which had the central subunit of MICOS, Mic60. The subsequent evolution of ERMIONE and its interactors in eukaryotes reflects the integrative co-evolution of mitochondria and their hosts and the adaptive paths that some lineages have followed in their specialization to certain environments. By approaching the ERMIONE from a perspective of comparative evolutionary cell biology, we hope to shed light on not only its evolutionary history, but also how ERMIONE components may function in organisms other than S. cerevisiae. This article is part of a Special Issue entitled: The cellular lipid landscape edited by Tim P. Levine and Anant K. Menon.

The ancient and widespread nature of the ER-mitochondria encounter structure.

Mitochondria are the result of a billion years of integrative evolution, converting a once free-living bacterium to an organelle deeply linked to diverse cellular processes. One way in which mitochondria are integrated with nonendosymbiotically derived organelles is via endoplasmic reticulum (ER)-mitochondria contact sites. The ER membrane is physically tethered to the mitochondrial outer membrane by the ER-mitochondria encounter structure (ERMES). However, to date, ERMES has only ever been found in the fungal lineage. Here, we bioinformatically demonstrate that ERMES is present in lineages outside Fungi and validate this inference by mass spectrometric identification of ERMES components in Acanthamoeba castellanii mitochondria. We further demonstrate that ERMES is retained in hydrogenosome-bearing but not mitosome-bearing organisms, yielding insight into the process of reductive mitochondrial evolution. Finally, we find that the taxonomic distribution of ERMES is most consistent with rooting the eukaryotic tree between Amorphea (Animals + Fungi + Amoebozoa) + Excavata and all other eukaryotes (Diaphoratickes).

A dynamin superfamily-like pseudoenzyme coordinates with MICOS to promote cristae architecture.

Mitochondrial cristae architecture is crucial for optimal respiratory function of the organelle. Cristae shape is maintained in part by the mitochondrial contact site and cristae organizing system (MICOS) complex. While MICOS is required for normal cristae morphology, the precise mechanistic role of each of the seven human MICOS subunits, and how the complex coordinates with other cristae-shaping factors, has not been fully determined. Here, we examine the MICOS complex in Schizosaccharomyces pombe, a minimal model whose genome only encodes for four core subunits. Using an unbiased proteomics approach, we identify a poorly characterized inner mitochondrial membrane protein that interacts with MICOS and is required to maintain cristae morphology, which we name Mmc1. We demonstrate that Mmc1 works in concert with MICOS to promote normal mitochondrial morphology and respiratory function. Mmc1 is a distant relative of the dynamin superfamily of proteins (DSPs), GTPases, which are well established to shape and remodel membranes. Similar to DSPs, Mmc1 self-associates and forms high-molecular-weight assemblies. Interestingly, however, Mmc1 is a pseudoenzyme that lacks key residues required for GTP binding and hydrolysis, suggesting that it does not dynamically remodel membranes. These data are consistent with the model that Mmc1 stabilizes cristae architecture by acting as a scaffold to support cristae ultrastructure on the matrix side of the inner membrane. Our study reveals a new class of proteins that evolved early in fungal phylogeny and is required for the maintenance of cristae architecture. This highlights the possibility that functionally analogous proteins work with MICOS to establish cristae morphology in metazoans.

The persistent homology of mitochondrial ATP synthases.

Relatively little is known about ATP synthase structure in protists, and the investigated ones exhibit divergent structures distinct from yeast or animals. To clarify the subunit composition of ATP synthases across all eukaryotic lineages, we used homology detection techniques and molecular modeling tools to identify an ancestral set of 17 ATP synthase subunits. Most eukaryotes possess an ATP synthase comparable to those of animals and fungi, while some have undergone drastic divergence (e.g., ciliates, myzozoans, euglenozoans). Additionally, a ∼1 billion-year-old gene fusion between ATP synthase stator subunits was identified as a synapomorphy of the SAR (Stramenopila, Alveolata, Rhizaria) supergroup (stramenopile, alveolate, rhizaria). Our comparative approach highlights the persistence of ancestral subunits even amidst major structural changes. We conclude by urging that more ATP synthase structures (e.g., from jakobids, heteroloboseans, stramenopiles, rhizarians) are needed to provide a complete picture of the evolution of the structural diversity of this ancient and essential complex.

Evolutionary History of Oxysterol-Binding Proteins Reveals Complex History of Duplication and Loss in Animals and Fungi.

Cells maintain the specific lipid composition of distinct organelles by vesicular transport as well as non-vesicular lipid trafficking via lipid transport proteins. Oxysterol-binding proteins (OSBPs) are a family of lipid transport proteins that transfer lipids at various membrane contact sites (MCSs). OSBPs have been extensively investigated in human and yeast cells where 12 have been identified in Homo sapiens and 7 in Saccharomyces cerevisiae. The evolutionary relationship between these well-characterized OSBPs is still unclear. By reconstructing phylogenies of eukaryote OSBPs, we show that the ancestral Saccharomycotina had four OSBPs, the ancestral fungus had five OSBPs, and the ancestral animal had six OSBPs, whereas the shared ancestor of animals and fungi as well as the ancestral eukaryote had only three OSBPs. Our analyses identified three undescribed ancient OSBP orthologues, one fungal OSBP (Osh8) lost in the lineage leading to yeast, one animal OSBP (ORP12) lost in the lineage leading to vertebrates, and one eukaryotic OSBP (OshEu) lost in both the animal and fungal lineages.

A DRP-like pseudoenzyme coordinates with MICOS to promote cristae architecture.

Mitochondrial cristae architecture is crucial for optimal respiratory function of the organelle. Cristae shape is maintained in part by the mitochondrial inner membrane-localized MICOS complex. While MICOS is required for normal cristae morphology, the precise mechanistic role of each of the seven human MICOS subunits, and how the complex coordinates with other cristae shaping factors, has not been fully determined. Here, we examine the MICOS complex in Schizosaccharomyces pombe, a minimal model whose genome only encodes for four core subunits. Using an unbiased proteomics approach, we identify a poorly characterized inner mitochondrial membrane protein that interacts with MICOS and is required to maintain cristae morphology, which we name Mmc1. We demonstrate that Mmc1 works in concert with MICOS complexes to promote normal mitochondrial morphology and respiratory function. Bioinformatic analyses reveal that Mmc1 is a distant relative of the Dynamin-Related Protein (DRP) family of GTPases, which are well established to shape and remodel membranes. We find that, like DRPs, Mmc1 self-associates and forms high molecular weight assemblies. Interestingly, however, Mmc1 is a pseudoenzyme that lacks key residues required for GTP binding and hydrolysis, suggesting it does not dynamically remodel membranes. These data are consistent with a model in which Mmc1 stabilizes cristae architecture by acting as a scaffold to support cristae ultrastructure on the matrix side of the inner membrane. Our study reveals a new class of proteins that evolved early in fungal phylogeny and is required for the maintenance of cristae architecture. This highlights the possibility that functionally analogous proteins work with MICOS to establish cristae morphology in metazoans.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria).

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

Evolution: Divergent trajectories predate the origins of animals and fungi.

Animals, fungi, and their closest protist relatives comprise the clade Opisthokonta. Although they are comparatively closely related, animals and fungi have diverged greatly from one another. A new study demonstrates that the genomic features that are characteristic of animals and fungi arose even before the origin of these two kingdoms.

A Eukaryote-Wide Perspective on the Diversity and Evolution of the ARF GTPase Protein Family.

The evolution of eukaryotic cellular complexity is interwoven with the extensive diversification of many protein families. One key family is the ARF GTPases that act in eukaryote-specific processes, including membrane traffic, tubulin assembly, actin dynamics, and cilia-related functions. Unfortunately, our understanding of the evolution of this family is limited. Sampling an extensive set of available genome and transcriptome sequences, we have assembled a data set of over 2,000 manually curated ARF family genes from 114 eukaryotic species, including many deeply diverged protist lineages, and carried out comprehensive molecular phylogenetic analyses. These reconstructed as many as 16 ARF family members present in the last eukaryotic common ancestor, nearly doubling the previously inferred ancient system complexity. Evidence for the wide occurrence and ancestral origin of Arf6, Arl13, and Arl16 is presented for the first time. Moreover, Arl17, Arl18, and SarB, newly described here, are absent from well-studied model organisms and as a result their function(s) remain unknown. Analyses of our data set revealed a previously unsuspected diversity of membrane association modes and domain architectures within the ARF family. We detail the step-wise expansion of the ARF family in the metazoan lineage, including discovery of several new animal-specific family members. Delving back to its earliest evolution in eukaryotes, the resolved relationship observed between the ARF family paralogs sets boundaries for scenarios of vesicle coat origins during eukaryogenesis. Altogether, our work fundamentally broadens the understanding of the diversity and evolution of a protein family underpinning the structural and functional complexity of the eukaryote cells.

First report of mitochondrial COI in foraminifera and implications for DNA barcoding.

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in many marine environments and play a major role in global carbon cycling. Species recognition in Foraminifera is mainly based on morphological characters and nuclear 18S ribosomal RNA barcoding. The 18S rRNA contains variable sequence regions that allow for the identification of most foraminiferal species. Still, some species show limited variability, while others contain high levels of intragenomic polymorphisms, thereby complicating species identification. The use of additional, easily obtainable molecular markers other than 18S rRNA will enable more detailed investigation of evolutionary history, population genetics and speciation in Foraminifera. Here we present the first mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences ("barcodes") of Foraminifera. We applied shotgun sequencing to single foraminiferal specimens, assembled COI, and developed primers that allow amplification of COI in a wide range of foraminiferal species. We obtained COI sequences of 49 specimens from 17 species from the orders Rotaliida and Miliolida. Phylogenetic analysis showed that the COI tree is largely congruent with previously published 18S rRNA phylogenies. Furthermore, species delimitation with ASAP and ABGD algorithms showed that foraminiferal species can be identified based on COI barcodes.

Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae.

The endosymbiotic origin of plastids from cyanobacteria gave eukaryotes photosynthetic capabilities and launched the diversification of countless forms of algae. These primary plastids are found in members of the eukaryotic supergroup Archaeplastida. All known archaeplastids still retain some form of primary plastids, which are widely assumed to have a single origin. Here, we use single-cell genomics from natural samples combined with phylogenomics to infer the evolutionary origin of the phylum Picozoa, a globally distributed but seemingly rare group of marine microbial heterotrophic eukaryotes. Strikingly, the analysis of 43 single-cell genomes shows that Picozoa belong to Archaeplastida, specifically related to red algae and the phagotrophic rhodelphids. These picozoan genomes support the hypothesis that Picozoa lack a plastid, and further reveal no evidence of an early cryptic endosymbiosis with cyanobacteria. These findings change our understanding of plastid evolution as they either represent the first complete plastid loss in a free-living taxon, or indicate that red algae and rhodelphids obtained their plastids independently of other archaeplastids.

Independent accretion of TIM22 complex subunits in the animal and fungal lineages.

Background: The mitochondrial protein import complexes arose early in eukaryogenesis. Most of the components of the protein import pathways predate the last eukaryotic common ancestor. For example, the carrier-insertase TIM22 complex comprises the widely conserved Tim22 channel core. However, the auxiliary components of fungal and animal TIM22 complexes are exceptions to this ancient conservation. Methods: Using comparative genomics and phylogenetic approaches, we identified precisely when each TIM22 accretion occurred. Results: In animals, we demonstrate that Tim29 and Tim10b arose early in the holozoan lineage. Tim29 predates the metazoan lineage being present in the animal sister lineages, choanoflagellate and filastereans, whereas the erroneously named Tim10b arose from a duplication of Tim9 at the base of metazoans. In fungi, we show that Tim54 has representatives present in every holomycotan lineage including microsporidians and fonticulids, whereas Tim18 and Tim12 appeared much later in fungal evolution. Specifically, Tim18 and Tim12 arose from duplications of Sdh3 and Tim10, respectively, early in the Saccharomycotina. Surprisingly, we show that Tim54 is distantly related to AGK suggesting that AGK and Tim54 are extremely divergent orthologues and the origin of AGK/Tim54 interaction with Tim22 predates the divergence of animals and fungi. Conclusions: We argue that the evolutionary history of the TIM22 complex is best understood as the neutral structural divergence of an otherwise strongly functionally conserved protein complex. This view suggests that many of the differences in structure/subunit composition of multi-protein complexes are non-adaptive. Instead, most of the phylogenetic variation of functionally conserved molecular machines, which have been under stable selective pressures for vast phylogenetic spans, such as the TIM22 complex, is most likely the outcome of the interplay of random genetic drift and mutation pressure.

Homologue replacement in the import motor of the mitochondrial inner membrane of trypanosomes.

Many mitochondrial proteins contain N-terminal presequences that direct them to the organelle. The main driving force for their translocation across the inner membrane is provided by the presequence translocase-associated motor (PAM) which contains the J-protein Pam18. Here, we show that in the PAM of Trypanosoma brucei the function of Pam18 has been replaced by the non-orthologous euglenozoan-specific J-protein TbPam27. TbPam27 is specifically required for the import of mitochondrial presequence-containing but not for carrier proteins. Similar to yeast Pam18, TbPam27 requires an intact J-domain to function. Surprisingly, T. brucei still contains a bona fide Pam18 orthologue that, while essential for normal growth, is not involved in protein import. Thus, during evolution of kinetoplastids, Pam18 has been replaced by TbPam27. We propose that this replacement is linked to the transition from two ancestral and functionally distinct TIM complexes, found in most eukaryotes, to the single bifunctional TIM complex present in trypanosomes.