Comparative genomics of Ascetosporea gives new insight into the evolutionary basis for animal parasitism in Rhizaria.

BACKGROUND: Ascetosporea (Endomyxa, Rhizaria) is a group of unicellular parasites infecting aquatic invertebrates. They are increasingly being recognized as widespread and important in marine environments, causing large annual losses in invertebrate aquaculture. Despite their importance, little molecular data of Ascetosporea exist, with only two genome assemblies published to date. Accordingly, the evolutionary origin of these parasites is unclear, including their phylogenetic position and the genomic adaptations that accompanied the transition from a free-living lifestyle to parasitism. Here, we sequenced and assembled three new ascetosporean genomes, as well as the genome of a closely related amphizoic species, to investigate the phylogeny, origin, and genomic adaptations to parasitism in Ascetosporea. RESULTS: Using a phylogenomic approach, we confirm the monophyly of Ascetosporea and show that Paramyxida group with Mikrocytida, with Haplosporida being sister to both groups. We report that the genomes of these parasites are relatively small (12-36 Mb) and gene-sparse (~ 2300-5200 genes), while containing surprisingly high amounts of non-coding sequence (~ 70-90% of the genomes). Performing gene-tree aware ancestral reconstruction of gene families, we demonstrate extensive gene losses at the origin of parasitism in Ascetosporea, primarily of metabolic functions, and little gene gain except on terminal branches. Finally, we highlight some functional gene classes that have undergone expansions during evolution of the group. CONCLUSIONS: We present important new genomic information from a lineage of enigmatic but important parasites of invertebrates and illuminate some of the genomic innovations accompanying the evolutionary transition to parasitism in this lineage. Our results and data provide a genetic basis for the development of control measures against these parasites.

PhyloFisher: A phylogenomic package for resolving eukaryotic relationships.

Phylogenomic analyses of hundreds of protein-coding genes aimed at resolving phylogenetic relationships is now a common practice. However, no software currently exists that includes tools for dataset construction and subsequent analysis with diverse validation strategies to assess robustness. Furthermore, there are no publicly available high-quality curated databases designed to assess deep (>100 million years) relationships in the tree of eukaryotes. To address these issues, we developed an easy-to-use software package, PhyloFisher (https://github.com/TheBrownLab/PhyloFisher), written in Python 3. PhyloFisher includes a manually curated database of 240 protein-coding genes from 304 eukaryotic taxa covering known eukaryotic diversity, a novel tool for ortholog selection, and utilities that will perform diverse analyses required by state-of-the-art phylogenomic investigations. Through phylogenetic reconstructions of the tree of eukaryotes and of the Saccharomycetaceae clade of budding yeasts, we demonstrate the utility of the PhyloFisher workflow and the provided starting database to address phylogenetic questions across a large range of evolutionary time points for diverse groups of organisms. We also demonstrate that undetected paralogy can remain in phylogenomic "single-copy orthogroup" datasets constructed using widely accepted methods such as all vs. all BLAST searches followed by Markov Cluster Algorithm (MCL) clustering and application of automated tree pruning algorithms. Finally, we show how the PhyloFisher workflow helps detect inadvertent paralog inclusions, allowing the user to make more informed decisions regarding orthology assignments, leading to a more accurate final dataset.

Identification of a new gregarine parasite associated with mass mortality events of freshwater pearl mussels (Margaritifera margaritifera) in Sweden.

Freshwater bivalves play key ecological roles in lakes and rivers, largely contributing to healthy ecosystems. The freshwater pearl mussel, Margaritifera margaritifera, is found in Europe and on the East coast of North America. Once common in oxygenated streams, M. margaritifera is rapidly declining and consequently assessed as a threatened species worldwide. Deterioration of water quality has been considered the main factor for the mass mortality events affecting this species. Yet, the role of parasitic infections has not been investigated. Here, we report the discovery of three novel protist lineages found in Swedish populations of M. margaritifera belonging to one of the terrestrial groups of gregarines (Eugregarinorida, Apicomplexa). These lineages are closely related-but clearly separated-from the tadpole parasite Nematopsis temporariae. In one lineage, which is specifically associated with mortality events of M. margaritifera, we found cysts containing single vermiform zoites in the gills and other organs of diseased individuals using microscopy and in situ hybridization. This represents the first report of a parasitic infection in M. margaritifera that may be linked to the decline of this mussel species. We propose a tentative life cycle with the distribution of different developmental stages and potential exit from the host into the environment.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians.

BACKGROUND: Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution. RESULTS: At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths. CONCLUSIONS: Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

Phylogenomic Insights into the Origin of Primary Plastids.

The origin of plastids was a major evolutionary event that paved the way for an astonishing diversification of photosynthetic eukaryotes. Plastids originated by endosymbiosis between a heterotrophic eukaryotic host and cyanobacteria, presumably in a common ancestor of the primary photosynthetic eukaryotes (Archaeplastida). A single origin of primary plastids is well supported by plastid evidence but not by nuclear phylogenomic analyses, which have consistently failed to recover the monophyly of Archaeplastida hosts. Importantly, plastid monophyly and nonmonophyletic hosts could be explained under scenarios of independent or serial eukaryote-to-eukaryote endosymbioses. Here, we assessed the strength of the signal for the monophyly of Archaeplastida hosts in four available phylogenomic data sets. The effect of phylogenetic methodology, data quality, alignment trimming strategy, gene and taxon sampling, and the presence of outlier genes were investigated. Our analyses revealed a lack of support for host monophyly in the shorter individual data sets. However, when analyzed together under rigorous data curation and complex mixture models, the combined nuclear data sets supported the monophyly of primary photosynthetic eukaryotes (Archaeplastida) and recovered a putative association with plastid-lacking Picozoa. This study represents an important step toward better understanding deep eukaryotic evolution and the origin of plastids. [Archaeplastida; Bayesian; chloroplast; maximum likelihood; mixture model; ortholog; outlier loci; paralog; protist.].

The long-time orphan protist Meringosphaera mediterranea Lohmann, 1902 [1903] is a centrohelid heliozoan.

Meringosphaera is an enigmatic marine protist without clear phylogenetic affiliation, but it has long been suggested to be a chrysophyte-related autotroph. Microscopy-based reports indicate that it has a worldwide distribution, but no sequence data exist so far. We obtained the first 18S rDNA sequence for M. mediterranea (identified using light and electron microscopy) from the west coast of Sweden. Observations of living cells revealed granulated axopodia and up to 6 globular photosynthesizing bodies about 2 µm in diameter, the nature of which requires further investigation. The ultrastructure of barbed undulating spine scales and patternless plate scales with a central thickening is in agreement with previous reports. Molecular phylogenetic analysis placed M. mediterranea inside the NC5 environmental clade of Centroplasthelida (Haptista) along with additional environmental sequences, together closely related to Choanocystidae. This placement is supported by similar scales in Meringosphaera and Choanocystidae. We searched the Tara Oceans 18S V9 metabarcoding dataset, which revealed four OTUs with 94.8%-98.2% similarity, with oceanic distribution similar to that based on morphological observations. The current taxonomic position and species composition of the genus are discussed. The planktonic lifestyle of M. mediterranea contradicts the view of some authors that centrohelids enter the plankton only temporarily.

New Phylogenomic Analysis of the Enigmatic Phylum Telonemia Further Resolves the Eukaryote Tree of Life.

The resolution of the broad-scale tree of eukaryotes is constantly improving, but the evolutionary origin of several major groups remains unknown. Resolving the phylogenetic position of these "orphan" groups is important, especially those that originated early in evolution, because they represent missing evolutionary links between established groups. Telonemia is one such orphan taxon for which little is known. The group is composed of molecularly diverse biflagellated protists, often prevalent although not abundant in aquatic environments. Telonemia has been hypothesized to represent a deeply diverging eukaryotic phylum but no consensus exists as to where it is placed in the tree. Here, we established cultures and report the phylogenomic analyses of three new transcriptome data sets for divergent telonemid lineages. All our phylogenetic reconstructions, based on 248 genes and using site-heterogeneous mixture models, robustly resolve the evolutionary origin of Telonemia as sister to the Sar supergroup. This grouping remains well supported when as few as 60% of the genes are randomly subsampled, thus is not sensitive to the sets of genes used but requires a minimal alignment length to recover enough phylogenetic signal. Telonemia occupies a crucial position in the tree to examine the origin of Sar, one of the most lineage-rich eukaryote supergroups. We propose the moniker "TSAR" to accommodate this new mega-assemblage in the phylogeny of eukaryotes.

Predatory colponemids are the sister group to all other alveolates.

Alveolates are a major supergroup of eukaryotes encompassing more than ten thousand free-living and parasitic species, including medically, ecologically, and economically important apicomplexans, dinoflagellates, and ciliates. These three groups are among the most widespread eukaryotes on Earth, and their environmental success can be linked to unique innovations that emerged early in each group. Understanding the emergence of these well-studied and diverse groups and their innovations has relied heavily on the discovery and characterization of early-branching relatives, which allow ancestral states to be inferred with much greater confidence. Here we report the phylogenomic analyses of 313 eukaryote protein-coding genes from transcriptomes of three members of one such group, the colponemids (Colponemidia), which support their monophyly and position as the sister lineage to all other known alveolates. Colponemid-related sequences from environmental surveys and our microscopical observations show that colponemids are not common in nature, but they are diverse and widespread in freshwater habitats around the world. Studied colponemids possess two types of extrusive organelles (trichocysts or toxicysts) for active hunting of other unicellular eukaryotes and potentially play an important role in microbial food webs. Colponemids have generally plesiomorphic morphology and illustrate the ancestral state of Alveolata. We further discuss their importance in understanding the evolution of alveolates and the origin of myzocytosis and plastids.

Perspectives from Ten Years of Protist Studies by High-Throughput Metabarcoding.

During the last decade, high-throughput metabarcoding became routine for analyzing protistan diversity and distributions in nature. Amid a multitude of exciting findings, scientists have also identified and addressed technical and biological limitations, although problems still exist for inference of meaningful taxonomic and ecological knowledge based on short DNA sequences. Given the extensive use of this approach, it is critical to settle our understanding on its strengths and weaknesses and to synthesize up-to-date methodological and conceptual trends. This article summarizes key scientific and technical findings, and identifies current and future directions in protist research that uses metabarcoding.

Major transitions in dinoflagellate evolution unveiled by phylotranscriptomics.

Dinoflagellates are key species in marine environments, but they remain poorly understood in part because of their large, complex genomes, unique molecular biology, and unresolved in-group relationships. We created a taxonomically representative dataset of dinoflagellate transcriptomes and used this to infer a strongly supported phylogeny to map major morphological and molecular transitions in dinoflagellate evolution. Our results show an early-branching position of Noctiluca, monophyly of thecate (plate-bearing) dinoflagellates, and paraphyly of athecate ones. This represents unambiguous phylogenetic evidence for a single origin of the group's cellulosic theca, which we show coincided with a radiation of cellulases implicated in cell division. By integrating dinoflagellate molecular, fossil, and biogeochemical evidence, we propose a revised model for the evolution of thecal tabulations and suggest that the late acquisition of dinosterol in the group is inconsistent with dinoflagellates being the source of this biomarker in pre-Mesozoic strata. Three distantly related, fundamentally nonphotosynthetic dinoflagellates, Noctiluca, Oxyrrhis, and Dinophysis, contain cryptic plastidial metabolisms and lack alternative cytosolic pathways, suggesting that all free-living dinoflagellates are metabolically dependent on plastids. This finding led us to propose general mechanisms of dependency on plastid organelles in eukaryotes that have lost photosynthesis; it also suggests that the evolutionary origin of bioluminescence in nonphotosynthetic dinoflagellates may be linked to plastidic tetrapyrrole biosynthesis. Finally, we use our phylogenetic framework to show that dinoflagellate nuclei have recruited DNA-binding proteins in three distinct evolutionary waves, which included two independent acquisitions of bacterial histone-like proteins.

PREQUAL: detecting non-homologous characters in sets of unaligned homologous sequences.

Summary: Phylogenomic datasets invariably contain undetected stretches of non-homologous characters due to poor-quality sequences or erroneous gene models. The large-scale multi-gene nature of these datasets renders impractical or impossible detailed manual curation of sequences, but few tools exist that can automate this task. To address this issue, we developed a new method that takes as input a set of unaligned homologous sequences and uses an explicit probabilistic approach to identify and mask regions with non-homologous adjacent characters. These regions are defined as sharing no statistical support for homology with any other sequence in the set, which can result from e.g. sequencing errors or gene prediction errors creating frameshifts. Our methodology is implemented in the program PREQUAL, which is a fast and accurate tool for high-throughput filtering of sequences. The program is primarily aimed at amino acid sequences, although it can handle protein coding DNA sequences as well. It is fully customizable to allow fine-tuning of the filtering sensitivity. Availability and implementation: The program PREQUAL is written in C/C++ and available through a GNU GPL v3.0 at https://github.com/simonwhelan/prequal. Supplementary information: Supplementary data are available at Bioinformatics online.

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.

Revisions to the Classification, Nomenclature, and Diversity of Eukaryotes.

This revision of the classification of eukaryotes follows that of Adl et al., 2012 [J. Euk. Microbiol. 59(5)] and retains an emphasis on protists. Changes since have improved the resolution of many nodes in phylogenetic analyses. For some clades even families are being clearly resolved. As we had predicted, environmental sampling in the intervening years has massively increased the genetic information at hand. Consequently, we have discovered novel clades, exciting new genera and uncovered a massive species level diversity beyond the morphological species descriptions. Several clades known from environmental samples only have now found their home. Sampling soils, deeper marine waters and the deep sea will continue to fill us with surprises. The main changes in this revision are the confirmation that eukaryotes form at least two domains, the loss of monophyly in the Excavata, robust support for the Haptista and Cryptista. We provide suggested primer sets for DNA sequences from environmental samples that are effective for each clade. We have provided a guide to trophic functional guilds in an appendix, to facilitate the interpretation of environmental samples, and a standardized taxonomic guide for East Asian users.

Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote-eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph.

Factors mediating plastid dependency and the origins of parasitism in apicomplexans and their close relatives.

Apicomplexans are a major lineage of parasites, including causative agents of malaria and toxoplasmosis. How such highly adapted parasites evolved from free-living ancestors is poorly understood, particularly because they contain nonphotosynthetic plastids with which they have a complex metabolic dependency. Here, we examine the origin of apicomplexan parasitism by resolving the evolutionary distribution of several key characteristics in their closest free-living relatives, photosynthetic chromerids and predatory colpodellids. Using environmental sequence data, we describe the diversity of these apicomplexan-related lineages and select five species that represent this diversity for transcriptome sequencing. Phylogenomic analysis recovered a monophyletic lineage of chromerids and colpodellids as the sister group to apicomplexans, and a complex distribution of retention versus loss for photosynthesis, plastid genomes, and plastid organelles. Reconstructing the evolution of all plastid and cytosolic metabolic pathways related to apicomplexan plastid function revealed an ancient dependency on plastid isoprenoid biosynthesis, predating the divergence of apicomplexan and dinoflagellates. Similarly, plastid genome retention is strongly linked to the retention of two genes in the plastid genome, sufB and clpC, altogether suggesting a relatively simple model for plastid retention and loss. Lastly, we examine the broader distribution of a suite of molecular characteristics previously linked to the origins of apicomplexan parasitism and find that virtually all are present in their free-living relatives. The emergence of parasitism may not be driven by acquisition of novel components, but rather by loss and modification of the existing, conserved traits.

Functional Relationship between a Dinoflagellate Host and Its Diatom Endosymbiont.

While we know much about the evolutionary patterns of endosymbiotic organelle origins, we know less about how the actual process unfolded within each system. This is partly due to the massive changes endosymbiosis appears to trigger, and partly because most organelles evolved in the distant past. The dinotoms are dinoflagellates with diatom endosymbionts, and they represent a relatively recent but nevertheless obligate endosymbiotic association. We have carried out deep sequencing of both the host and endosymbiont transcriptomes from two dinotoms, Durinskia baltica and Glenodinium foliaceum, to examine how the nucleocytosolic compartments have functionally integrated. This analysis showed little or no functional reduction in either the endosymbiont or host, and no evidence for genetic integration. Rather, host and endosymbiont seem to be bound to each other via metabolites, such as photosynthate exported from the endosymbiont to the host as indicated by the presence of plastidic phosphate translocators in the host transcriptome. The host is able to synthesize starch, using plant-specific starch synthases, as a way to store imported photosynthate.

Evolutionary Origins of Rhizarian Parasites.

The SAR group (Stramenopila, Alveolata, Rhizaria) is one of the largest clades in the tree of eukaryotes and includes a great number of parasitic lineages. Rhizarian parasites are obligate and have devastating effects on commercially important plants and animals but despite this fact, our knowledge of their biology and evolution is limited. Here, we present rhizarian transcriptomes from all major parasitic lineages in order to elucidate their evolutionary relationships using a phylogenomic approach. Our results suggest that Ascetosporea, parasites of marine invertebrates, are sister to the novel clade Apofilosa. The phytomyxean plant parasites branch sister to the vampyrellid algal ectoparasites in the novel clade Phytorhiza. They also show that Ascetosporea + Apofilosa + Retaria + Filosa + Phytorhiza form a monophyletic clade, although the branching pattern within this clade is difficult to resolve and appears to be model-dependent. Our study does not support the monophyly of the rhizarian parasitic lineages (Endomyxa), suggesting independent origins for rhizarian animal and plant parasites.

Untangling the early diversification of eukaryotes: a phylogenomic study of the evolutionary origins of Centrohelida, Haptophyta and Cryptista.

Assembling the global eukaryotic tree of life has long been a major effort of Biology. In recent years, pushed by the new availability of genome-scale data for microbial eukaryotes, it has become possible to revisit many evolutionary enigmas. However, some of the most ancient nodes, which are essential for inferring a stable tree, have remained highly controversial. Among other reasons, the lack of adequate genomic datasets for key taxa has prevented the robust reconstruction of early diversification events. In this context, the centrohelid heliozoans are particularly relevant for reconstructing the tree of eukaryotes because they represent one of the last substantial groups that was missing large and diverse genomic data. Here, we filled this gap by sequencing high-quality transcriptomes for four centrohelid lineages, each corresponding to a different family. Combining these new data with a broad eukaryotic sampling, we produced a gene-rich taxon-rich phylogenomic dataset that enabled us to refine the structure of the tree. Specifically, we show that (i) centrohelids relate to haptophytes, confirming Haptista; (ii) Haptista relates to SAR; (iii) Cryptista share strong affinity with Archaeplastida; and (iv) Haptista + SAR is sister to Cryptista + Archaeplastida. The implications of this topology are discussed in the broader context of plastid evolution.

Endosymbiotic gene transfer in tertiary plastid-containing dinoflagellates.

Plastid establishment involves the transfer of endosymbiotic genes to the host nucleus, a process known as endosymbiotic gene transfer (EGT). Large amounts of EGT have been shown in several photosynthetic lineages but also in present-day plastid-lacking organisms, supporting the notion that endosymbiotic genes leave a substantial genetic footprint in the host nucleus. Yet the extent of this genetic relocation remains debated, largely because the long period that has passed since most plastids originated has erased many of the clues to how this process unfolded. Among the dinoflagellates, however, the ancestral peridinin-containing plastid has been replaced by tertiary plastids on several more recent occasions, giving us a less ancient window to examine plastid origins. In this study, we evaluated the endosymbiotic contribution to the host genome in two dinoflagellate lineages with tertiary plastids. We generated the first nuclear transcriptome data sets for the "dinotoms," which harbor diatom-derived plastids, and analyzed these data in combination with the available transcriptomes for kareniaceans, which harbor haptophyte-derived plastids. We found low level of detectable EGT in both dinoflagellate lineages, with only 9 genes and 90 genes of possible tertiary endosymbiotic origin in dinotoms and kareniaceans, respectively, suggesting that tertiary endosymbioses did not heavily impact the host dinoflagellate genomes.

Gain and loss of multiple functionally related, horizontally transferred genes in the reduced genomes of two microsporidian parasites.

Microsporidia of the genus Encephalitozoon are widespread pathogens of animals that harbor the smallest known nuclear genomes. Complete sequences from Encephalitozoon intestinalis (2.3 Mbp) and Encephalitozoon cuniculi (2.9 Mbp) revealed massive gene losses and reduction of intergenic regions as factors leading to their drastically reduced genome size. However, microsporidian genomes also have gained genes through horizontal gene transfers (HGT), a process that could allow the parasites to exploit their hosts more fully. Here, we describe the complete sequences of two intermediate-sized genomes (2.5 Mbp), from Encephalitozoon hellem and Encephalitozoon romaleae. Overall, the E. hellem and E. romaleae genomes are strikingly similar to those of Encephalitozoon cuniculi and Encephalitozoon intestinalis in both form and content. However, in addition to the expected expansions and contractions of known gene families in subtelomeric regions, both species also were found to harbor a number of protein-coding genes that are not found in any other microsporidian. All these genes are functionally related to the metabolism of folate and purines but appear to have originated by several independent HGT events from different eukaryotic and prokaryotic donors. Surprisingly, the genes are all intact in E. hellem, but in E. romaleae those involved in de novo synthesis of folate are all pseudogenes. Overall, these data suggest that a recent common ancestor of E. hellem and E. romaleae assembled a complete metabolic pathway from multiple independent HGT events and that one descendent already is dispensing with much of this new functionality, highlighting the transient nature of transferred genes.