Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes.

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

Systematic comparison of unilamellar vesicles reveals that archaeal core lipid membranes are more permeable than bacterial membranes.

One of the deepest branches in the tree of life separates the Archaea from the Bacteria. These prokaryotic groups have distinct cellular systems including fundamentally different phospholipid membrane bilayers. This dichotomy has been termed the lipid divide and possibly bestows different biophysical and biochemical characteristics on each cell type. Classic experiments suggest that bacterial membranes (formed from lipids extracted from Escherichia coli, for example) show permeability to key metabolites comparable to archaeal membranes (formed from lipids extracted from Halobacterium salinarum), yet systematic analyses based on direct measurements of membrane permeability are absent. Here, we develop a new approach for assessing the membrane permeability of approximately 10 µm unilamellar vesicles, consisting of an aqueous medium enclosed by a single lipid bilayer. Comparing the permeability of 18 metabolites demonstrates that diether glycerol-1-phosphate lipids with methyl branches, often the most abundant membrane lipids of sampled archaea, are permeable to a wide range of compounds useful for core metabolic networks, including amino acids, sugars, and nucleobases. Permeability is significantly lower in diester glycerol-3-phosphate lipids without methyl branches, the common building block of bacterial membranes. To identify the membrane characteristics that determine permeability, we use this experimental platform to test a variety of lipid forms bearing a diversity of intermediate characteristics. We found that increased membrane permeability is dependent on both the methyl branches on the lipid tails and the ether bond between the tails and the head group, both of which are present on the archaeal phospholipids. These permeability differences must have had profound effects on the cell physiology and proteome evolution of early prokaryotic forms. To explore this further, we compare the abundance and distribution of transmembrane transporter-encoding protein families present on genomes sampled from across the prokaryotic tree of life. These data demonstrate that archaea tend to have a reduced repertoire of transporter gene families, consistent with increased membrane permeation. These results demonstrate that the lipid divide demarcates a clear difference in permeability function with implications for understanding some of the earliest transitions in cell origins and evolution.

New Parabasalia symbionts Snyderella spp. and Daimonympha gen. nov. from South American Rugitermes termites and the parallel evolution of a cell with a rotating "head".

Most Parabasalia are symbionts in the hindgut of "lower" (non-Termitidae) termites, where they widely vary in morphology and degree of morphological complexity. Large and complex cells in the class Cristamonadea evolved by replicating a fundamental unit, the karyomastigont, in various ways. We describe here four new species of Calonymphidae (Cristamonadea) from Rugitermes hosts, assigned to the genus Snyderella based on diagnostic features (including the karyomastigont pattern) and molecular phylogeny. We also report a new genus of Calonymphidae, Daimonympha, from Rugitermes laticollis. Daimonympha's morphology does not match that of any known Parabasalia, and its SSU rRNA gene sequence corroborates this distinction. Daimonympha does however share a puzzling feature with a few previously described, but distantly related, Cristamonadea: a rapid, smooth, and continuous rotation of the anterior end of the cell, including the many karyomastigont nuclei. The function of this rotatory movement, the cellular mechanisms enabling it, and the way the cell deals with the consequent cell membrane shear, are all unknown. "Rotating wheel" structures are famously rare in biology, with prokaryotic flagella being the main exception; these mysterious spinning cells found only among Parabasalia are another, far less understood, example.

A distinct lineage of giant viruses brings a rhodopsin photosystem to unicellular marine predators.

Giant viruses are remarkable for their large genomes, often rivaling those of small bacteria, and for having genes thought exclusive to cellular life. Most isolated to date infect nonmarine protists, leaving their strategies and prevalence in marine environments largely unknown. Using eukaryotic single-cell metagenomics in the Pacific, we discovered a Mimiviridae lineage of giant viruses, which infects choanoflagellates, widespread protistan predators related to metazoans. The ChoanoVirus genomes are the largest yet from pelagic ecosystems, with 442 of 862 predicted proteins lacking known homologs. They are enriched in enzymes for modifying organic compounds, including degradation of chitin, an abundant polysaccharide in oceans, and they encode 3 divergent type-1 rhodopsins (VirR) with distinct evolutionary histories from those that capture sunlight in cellular organisms. One (VirRDTS) is similar to the only other putative rhodopsin from a virus (PgV) with a known host (a marine alga). Unlike the algal virus, ChoanoViruses encode the entire pigment biosynthesis pathway and cleavage enzyme for producing the required chromophore, retinal. We demonstrate that the rhodopsin shared by ChoanoViruses and PgV binds retinal and pumps protons. Moreover, our 1.65-Å resolved VirRDTS crystal structure and mutational analyses exposed differences from previously characterized type-1 rhodopsins, all of which come from cellular organisms. Multiple VirR types are present in metagenomes from across surface oceans, where they are correlated with and nearly as abundant as a canonical marker gene from Mimiviridae Our findings indicate that light-dependent energy transfer systems are likely common components of giant viruses of photosynthetic and phagotrophic unicellular marine eukaryotes.

Taxonomy of the Apicomplexan Symbionts of Coral, including Corallicolida ord. nov., Reassignment of the Genus Gemmocystis, and Description of New Species Corallicola aquarius gen. nov. sp. nov. and Anthozoaphila gnarlus gen. nov. sp. nov.

Corals (Metazoa; Cnidaria; Anthozoa) have recently been shown to play host to a widespread and diverse group of intracellular symbionts of the phylum Apicomplexa. These symbionts, colloquially called "corallicolids", are mostly known through molecular analyses, and no formal taxonomy has been proposed. Another apicomplexan, Gemmocystis cylindrus (described from the coral Dendrogyra cylindrus), may be related to corallicolids, but lacks molecular data. Here, we isolate and describe motile trophozoite (feeding) corallicolids cells using microscopic (light, SEM, and TEM) and molecular phylogenetic analysis to provide the basis for a formal description. Phylogenetic analyses using nuclear and plastid rRNA operons, and three mitochondrial protein sequences derived from single cell transcriptomes, all confirm that these organisms fall into a discrete deep-branching clade within the Apicomplexa not closely related to any known species or major subgroup. As a result, we assign this clade to a new order, Corallicolida ord. nov., and family, Corallicolidae fam. nov. We describe a type species, Corallicola aquarius gen. nov. sp. nov. from its Rhodactis sp. host, and also describe a second species, Anthozoaphila gnarlus gen. nov. sp. nov., from the coral host Madracis mirabilis. Finally, we propose reassigning the incertae sedis taxon G. cylindrus from the order Agamococcidiorida to the Corallicolida, based on similarities in morphology and host localization to that of the corallicolids.

Phylogenomics supports the monophyly of the Cercozoa.

The phylum Cercozoa consists of a diverse assemblage of amoeboid and flagellated protists that forms a major component of the supergroup, Rhizaria. However, despite its size and ubiquity, the phylogeny of the Cercozoa remains unclear as morphological variability between cercozoan species and ambiguity in molecular analyses, including phylogenomic approaches, have produced ambiguous results and raised doubts about the monophyly of the group. Here we sought to resolve these ambiguities using a 161-gene phylogenetic dataset with data from newly available genomes and deeply sequenced transcriptomes, including three new transcriptomes from Aurigamonas solis, Abollifer prolabens, and a novel species, Lapot gusevi n. gen. n. sp. Our phylogenomic analysis strongly supported a monophyletic Cercozoa, and approximately-unbiased tests rejected the paraphyletic topologies observed in previous studies. The transcriptome of L. gusevi represents the first transcriptomic data from the large and recently characterized Aquavolonidae-Treumulida-'Novel Clade 12' group, and phylogenomics supported its position as sister to the cercozoan subphylum, Endomyxa. These results provide insights into the phylogeny of the Cercozoa and the Rhizaria as a whole.

What Can Environmental Sequences Tell Us About the Distribution of Low-Rank Taxa? The Case of Euplotes (Ciliophora, Spirotrichea), Including a Description of Euplotes enigma sp. nov.

Environmental sequences have become a major source of information. High-throughput sequencing (HTS) surveys have been used to infer biogeographic patterns and distribution of broad taxa of protists. This approach is, however, more questionable for addressing low-rank (less inclusive) taxa such as species and genera, because of the increased chance of errors in identification due to blurry taxonomic boundaries, low sequence divergence, or sequencing errors. The specious ciliate genus Euplotes partially escapes these limitations. It is a ubiquitous, monophyletic taxon, clearly differentiated from related genera, and with a relatively well-developed internal systematics. It has also been the focus of several ecological studies. We present an update on Euplotes biogeography, taking into consideration for the first time environmental sequences, both traditional (Sanger) and HTS. We inferred a comprehensive small subunit rRNA gene phylogeny of the genus including a newly described marine species, Euplotes enigma, characterized by a unique question mark-shaped macronucleus. We then added available environmental sequences to the tree, mapping associated metadata. The resulting scenario conflicts on many accounts with previously held views, suggesting, for example, that a large diversity of anaerobic Euplotes species exist, and that marine representatives of mainly freshwater lineages (and vice-versa) might be more common than previously thought.

Molecular characterization and phylogeny of four new species of the genus Trichonympha (Parabasalia, Trichonymphea) from lower termite hindguts.

Members of the genus Trichonympha are among the most well-known, recognizable and widely distributed parabasalian symbionts of lower termites and the wood-eating cockroach species of the genus Cryptocercus. Nevertheless, the species diversity of this genus is largely unknown. Molecular data have shown that the superficial morphological similarities traditionally used to identify species are inadequate, and have challenged the view that the same species of the genus Trichonympha can occur in many different host species. Ambiguities in the literature, uncertainty in identification of both symbiont and host, and incomplete samplings are limiting our understanding of the systematics, ecology and evolution of this taxon. Here we describe four closely related novel species of the genus Trichonympha collected from South American and Australian lower termites: Trichonympha hueyi sp. nov. from Rugitermes laticollis, Trichonympha deweyi sp. nov. from Glyptotermes brevicornis, Trichonympha louiei sp. nov. from Calcaritermes temnocephalus and Trichonympha webbyae sp. nov. from Rugitermes bicolor. We provide molecular barcodes to identify both the symbionts and their hosts, and infer the phylogeny of the genus Trichonympha based on small subunit rRNA gene sequences. The analysis confirms the considerable divergence of symbionts of members of the genus Cryptocercus, and shows that the two clades of the genus Trichonympha harboured by termites reflect only in part the phylogeny of their hosts.

Molecular Phylogeny of Mobilid and Sessilid Ciliates Symbiotic in Eastern Pacific Limpets (Mollusca: Patellogastropoda).

The phylogenetic relationships of the ciliate subclass Peritrichia, composed of the orders Mobilida and Sessilida, have recently come under debate as morphological and molecular analyses have struck contrasting conclusions as to the monophyly of the group. We provide additional molecular data to assess the monophyly of the Peritrichia by sequencing the small subunit ribosomal RNA genes of two symbiotic peritrichs, Urceolaria korschelti and Scyphidia ubiquita, found inhabiting the mantle cavity of limpets. Although phylogenetic analyses indicated a nonmonophyletic Peritrichia, approximately unbiased tests revealed that the monophyletic hypothesis could not be rejected. With regard to the Mobilida, our analysis showed divergence within the family Trichodinidae related to host taxa-a molluscan clade and a fish clade. For the Sessilida, the family Scyphidiidae was sister to the Astylozoidae. In our sampling of U. korschelti and S. ubiquita, both species showed significant genetic divergence among geographically isolated, yet morphologically indistinguishable populations. We hypothesize that cryptic speciation has produced these morphologically identical species and argue that more extensive genomic analyses are required to fully assess the monophyly, biogeography, and ultimately biodiversity of the peritrichs.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes.

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

Multiple parallel origins of parasitic Marine Alveolates.

Microbial eukaryotes are important components of marine ecosystems, and the Marine Alveolates (MALVs) are consistently both abundant and diverse in global environmental sequencing surveys. MALVs are dinoflagellates that are thought to be parasites of other protists and animals, but the lack of data beyond ribosomal RNA gene sequences from all but a few described species means much of their biology and evolution remain unknown. Using single-cell transcriptomes from several MALVs and their free-living relatives, we show that MALVs evolved independently from two distinct, free-living ancestors and that their parasitism evolved in parallel. Phylogenomics shows one subgroup (MALV-II and -IV, or Syndiniales) is related to a novel lineage of free-living, eukaryovorous predators, the eleftherids, while the other (MALV-I, or Ichthyodinida) is related to the free-living predator Oxyrrhis and retains proteins targeted to a non-photosynthetic plastid. Reconstructing the evolution of photosynthesis, plastids, and parasitism in early-diverging dinoflagellates shows a number of parallels with the evolution of their apicomplexan sisters. In both groups, similar forms of parasitism evolved multiple times and photosynthesis was lost many times. By contrast, complete loss of the plastid organelle is infrequent and, when this does happen, leaves no residual genes.

Systematic evaluation of horizontal gene transfer between eukaryotes and viruses.

Gene exchange between viruses and their hosts acts as a key facilitator of horizontal gene transfer and is hypothesized to be a major driver of evolutionary change. Our understanding of this process comes primarily from bacteria and phage co-evolution, but the mode and functional importance of gene transfers between eukaryotes and their viruses remain anecdotal. Here we systematically characterized viral-eukaryotic gene exchange across eukaryotic and viral diversity, identifying thousands of transfers and revealing their frequency, taxonomic distribution and projected functions. Eukaryote-derived viral genes, abundant in the Nucleocytoviricota, highlighted common strategies for viral host-manipulation, including metabolic reprogramming, proteolytic degradation and extracellular modification. Furthermore, viral-derived eukaryotic genes implicate genetic exchange in the early evolution and diversification of eukaryotes, particularly through viral-derived glycosyltransferases, which have impacted structures as diverse as algal cell walls, trypanosome mitochondria and animal tissues. These findings illuminate the nature of viral-eukaryotic gene exchange and its impact on the evolution of viruses and their eukaryotic hosts.

All essential endosymbionts of the ciliate Euplotes are cyclically replaced.

Symbiotic systems vary in the degree to which the partners are bound to each other1. At one extreme, there are intracellular endosymbionts in mutually obligate relationships with their host, often interpreted as mutualistic. The symbiosis between the betaproteobacterium Polynucleobacter and the ciliate Euplotes (clade B) challenges this view2: although freshwater Euplotes species long ago became dependent on endosymbionts, the many extant Polynucleobacter lineages they harbour arose recently and in parallel from different free-living ancestors2. The host requires the endosymbionts for reproduction and survival3, but each newly established symbiont is ultimately driven to extinction in a cycle of establishment, degeneration, and replacement. Similar replacement events have been observed in sap-feeding insects4-6, a model for bacteria-eukaryote symbioses7, but usually only affect a small subset of the host populations. Most insects retain an ancient coevolving symbiont, suggesting that long-term mutualism and permanent integration remain the rule and symbiont turnovers are mere evolutionary side-stories. Here we show that this is not the case for Euplotes. We examined all known essential Euplotes symbionts and found that none are ancient or coevolving; rather, all are recently established and continuously replaced over relatively short evolutionary time spans, making the symbiosis ancient for the host but not for any bacterial lineage.

A diversified and segregated mRNA spliced-leader system in the parasitic Perkinsozoa.

Spliced-leader trans-splicing (SLTS) has been described in distantly related eukaryotes and acts to mark mRNAs with a short 5' exon, giving different mRNAs identical 5' sequence-signatures. The function of these systems is obscure. Perkinsozoa encompasses a diversity of parasitic protists that infect bivalves, toxic-tide dinoflagellates, fish and frog tadpoles. Here, we report considerable sequence variation in the SLTS-system across the Perkinsozoa and find that multiple variant SLTS-systems are encoded in parallel in the ecologically important Perkinsozoa parasite Parvilucifera sinerae. These results demonstrate that the transcriptome of P. sinerae is segregated based on the addition of different spliced-leader (SL) exons. This segregation marks different gene categories, suggesting that SL-segregation relates to functional differentiation of the transcriptome. By contrast, both sets of gene categories are present in the single SL-transcript type sampled from Maranthos, implying that the SL-segregation of the Parvilucifera transcriptome is a recent evolutionary innovation. Furthermore, we show that the SLTS-system marks a subsection of the transcriptome with increased mRNA abundance and includes genes that encode the spliceosome system necessary for SLTS-function. Collectively, these data provide a picture of how the SLTS-systems can vary within a major evolutionary group and identify how additional transcriptional-complexity can be achieved through SL-segregation.

Microscopic marine invertebrates are reservoirs for cryptic and diverse protists and fungi.

BACKGROUND: Microbial symbioses in marine invertebrates are commonplace. However, characterizations of invertebrate microbiomes are vastly outnumbered by those of vertebrates. Protists and fungi run the gamut of symbiosis, yet eukaryotic microbiome sequencing is rarely undertaken, with much of the focus on bacteria. To explore the importance of microscopic marine invertebrates as potential symbiont reservoirs, we used a phylogenetic-focused approach to analyze the host-associated eukaryotic microbiomes of 220 animal specimens spanning nine different animal phyla. RESULTS: Our data expanded the traditional host range of several microbial taxa and identified numerous undescribed lineages. A lack of comparable reference sequences resulted in several cryptic clades within the Apicomplexa and Ciliophora and emphasized the potential for microbial invertebrates to harbor novel protistan and fungal diversity. CONCLUSIONS: Microscopic marine invertebrates, spanning a wide range of animal phyla, host various protist and fungal sequences and may therefore serve as a useful resource in the detection and characterization of undescribed symbioses. Video Abstract.

The molecular phylogeny of Chionaster nivalis reveals a novel order of psychrophilic and globally distributed Tremellomycetes (Fungi, Basidiomycota).

Snow and ice present challenging substrates for cellular growth, yet microbial snow communities not only exist, but are diverse and ecologically impactful. These communities are dominated by green algae, but additional organisms, such as fungi, are also abundant and may be important for nutrient cycling, syntrophic interactions, and community structure in general. However, little is known about these non-algal community members, including their taxonomic affiliations. An example of this is Chionaster nivalis, a unicellular fungus that is morphologically enigmatic and frequently observed in snow communities globally. Despite being described over one hundred years ago, the phylogeny and higher-level taxonomic classifications of C. nivalis remain unknown. Here, we isolated and sequenced the internal transcribed spacer (ITS) and the D1-D2 region of the large subunit ribosomal RNA gene of C. nivalis, providing a molecular barcode for future studies. Phylogenetic analyses using the ITS and D1-D2 region revealed that C. nivalis is part of a novel lineage in the class Tremellomycetes (Basidiomycota, Agaricomycotina) for which a new order Chionasterales ord. nov. (MB838717) and family Chionasteraceae fam. nov. (MB838718) are proposed. Comparisons between C. nivalis and sequences generated from environmental surveys revealed that the Chionasterales are globally distributed and probably psychrophilic, as they appear to be limited to the high alpine and arctic regions. These results highlight the unexplored diversity that exists within these extreme habitats and emphasize the utility of single-cell approaches in characterizing these complex algal-dominated communities.

Phylogenomics Identifies a New Major Subgroup of Apicomplexans, Marosporida class nov., with Extreme Apicoplast Genome Reduction.

The phylum Apicomplexa consists largely of obligate animal parasites that include the causative agents of human diseases such as malaria. Apicomplexans have also emerged as models to study the evolution of nonphotosynthetic plastids, as they contain a relict chloroplast known as the apicoplast. The apicoplast offers important clues into how apicomplexan parasites evolved from free-living ancestors and can provide insights into reductive organelle evolution. Here, we sequenced the transcriptomes and apicoplast genomes of three deep-branching apicomplexans, Margolisiella islandica, Aggregata octopiana, and Merocystis kathae. Phylogenomic analyses show that these taxa, together with Rhytidocystis, form a new lineage of apicomplexans that is sister to the Coccidia and Hematozoa (the lineages including most medically significant taxa). Members of this clade retain plastid genomes and the canonical apicomplexan plastid metabolism. However, the apicoplast genomes of Margolisiella and Rhytidocystis are the most reduced of any apicoplast, are extremely GC-poor, and have even lost genes for the canonical plastidial RNA polymerase. This new lineage of apicomplexans, for which we propose the class Marosporida class nov., occupies a key intermediate position in the apicomplexan phylogeny, and adds a new complexity to the models of stepwise reductive evolution of genome structure and organelle function in these parasites.

The Function and Evolution of Motile DNA Replication Systems in Ciliates.

DNA replication is a ubiquitous and conserved cellular process. However, regulation of DNA replication is only understood in a small fraction of organisms that poorly represent the diversity of genetic systems in nature. Here we used computational and experimental approaches to examine the function and evolution of one such system, the replication band (RB) in spirotrich ciliates, which is a localized, motile hub that traverses the macronucleus while replicating DNA. We show that the RB can take unique forms in different species, from polar bands to a "replication envelope," where replication initiates at the nuclear periphery before advancing inward. Furthermore, we identify genes involved in cellular transport, including calcium transporters and cytoskeletal regulators, that are associated with the RB and may be involved in its function and translocation. These findings highlight the evolution and diversity of DNA replication systems and provide insights into the regulation of nuclear organization and processes.

Characterization of new cristamonad species from kalotermitid termites including a novel genus, Runanympha.

Cristamonadea is a large class of parabasalian protists that reside in the hindguts of wood-feeding insects, where they play an essential role in the digestion of lignocellulose. This group of symbionts boasts an impressive array of complex morphological characteristics, many of which have evolved multiple times independently. However, their diversity is understudied and molecular data remain scarce. Here we describe seven new species of cristamonad symbionts from Comatermes, Calcaritermes, and Rugitermes termites from Peru and Ecuador. To classify these new species, we examined cells by light and scanning electron microscopy, sequenced the symbiont small subunit ribosomal RNA (rRNA) genes, and carried out barcoding of the mitochondrial large subunit rRNA gene of the hosts to confirm host identification. Based on these data, five of the symbionts characterized here represent new species within described genera: Devescovina sapara n. sp., Devescovina aymara n. sp., Macrotrichomonas ashaninka n. sp., Macrotrichomonas secoya n. sp., and Macrotrichomonas yanesha n. sp. Additionally, two symbionts with overall morphological characteristics similar to the poorly-studied and probably polyphyletic 'joeniid' Parabasalia are classified in a new genus Runanympha n. gen.: Runanympha illapa n. sp., and Runanympha pacha n. sp.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists.

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.