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.

Shuffled ATG8 interacting motifs form an ancestral bridge between UFMylation and autophagy.

UFMylation involves the covalent modification of substrate proteins with UFM1 (Ubiquitin-fold modifier 1) and is important for maintaining ER homeostasis. Stalled translation triggers the UFMylation of ER-bound ribosomes and activates C53-mediated autophagy to clear toxic polypeptides. C53 contains noncanonical shuffled ATG8-interacting motifs (sAIMs) that are essential for ATG8 interaction and autophagy initiation. However, the mechanistic basis of sAIM-mediated ATG8 interaction remains unknown. Here, we show that C53 and sAIMs are conserved across eukaryotes but secondarily lost in fungi and various algal lineages. Biochemical assays showed that the unicellular alga Chlamydomonas reinhardtii has a functional UFMylation pathway, refuting the assumption that UFMylation is linked to multicellularity. Comparative structural analyses revealed that both UFM1 and ATG8 bind sAIMs in C53, but in a distinct way. Conversion of sAIMs into canonical AIMs impaired binding of C53 to UFM1, while strengthening ATG8 binding. Increased ATG8 binding led to the autoactivation of the C53 pathway and sensitization of Arabidopsis thaliana to ER stress. Altogether, our findings reveal an ancestral role of sAIMs in UFMylation-dependent fine-tuning of C53-mediated autophagy activation.

Symbiosis: In search of a deeper understanding.

How do distinct species cofunction in symbiosis, despite conflicting interests? A new collection of articles explores emerging themes as researchers exploit modern research tools and new models to unravel how symbiotic interactions function and evolve.

Duplication and neofunctionalization of a horizontally transferred xyloglucanase as a facet of the Red Queen coevolutionary dynamic.

Oomycete protists share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a distant region of the tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls, a barrier to pathogen invasion and a rich source of carbohydrates. Using a combination of phylogenomics and functional assays, we investigate the diversification of a horizontally transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 xyloglucanase paralogs retained in P. sojae. Using heterologous expression in yeast, we show consistent evidence that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional variants analyzed subtend a phylogenetic node close to the fungi-to-oomycete transfer, suggesting the horizontally transferred gene was a bona fide xyloglucanase. Expression of three xyloglucanase paralogs in Nicotiana benthamiana triggers high-reactive oxygen species (ROS) generation, while others inhibit ROS responses to bacterial immunogens, demonstrating that the paralogs differentially stimulate pattern-triggered immunity. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze the production of variant breakdown profiles, suggesting that secretion of variant xyloglucanases increases efficiency of xyloglucan breakdown as well as diversifying the damage-associated molecular patterns released. We suggest that this pattern of neofunctionalization and the variant host responses represent an aspect of the Red Queen host-pathogen coevolutionary dynamic.

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.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids.

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

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.

Emergent RNA-RNA interactions can promote stability in a facultative phototrophic endosymbiosis.

Eukaryote-eukaryote endosymbiosis was responsible for the spread of chloroplast (plastid) organelles. Stability is required for the metabolic and genetic integration that drives the establishment of new organelles, yet the mechanisms that act to stabilize emergent endosymbioses-between two fundamentally selfish biological organisms-are unclear. Theory suggests that enforcement mechanisms, which punish misbehavior, may act to stabilize such interactions by resolving conflict. However, how such mechanisms can emerge in a facultative endosymbiosis has yet to be explored. Here, we propose that endosymbiont-host RNA-RNA interactions, arising from digestion of the endosymbiont population, can result in a cost to host growth for breakdown of the endosymbiosis. Using the model facultative endosymbiosis between Paramecium bursaria and Chlorella spp., we demonstrate that this mechanism is dependent on the host RNA-interference (RNAi) system. We reveal through small RNA (sRNA) sequencing that endosymbiont-derived messenger RNA (mRNA) released upon endosymbiont digestion can be processed by the host RNAi system into 23-nt sRNA. We predict multiple regions of shared sequence identity between endosymbiont and host mRNA, and demonstrate through delivery of synthetic endosymbiont sRNA that exposure to these regions can knock down expression of complementary host genes, resulting in a cost to host growth. This process of host gene knockdown in response to endosymbiont-derived RNA processing by host RNAi factors, which we term "RNAi collisions," represents a mechanism that can promote stability in a facultative eukaryote-eukaryote endosymbiosis. Specifically, by imposing a cost for breakdown of the endosymbiosis, endosymbiont-host RNA-RNA interactions may drive maintenance of the symbiosis across fluctuating ecological conditions.

Environment-dependent fitness gains can be driven by horizontal gene transfer of transporter-encoding genes.

Many microbes acquire metabolites in a "feeding" process where complex polymers are broken down in the environment to their subunits. The subsequent uptake of soluble metabolites by a cell, sometimes called osmotrophy, is facilitated by transporter proteins. As such, the diversification of osmotrophic microorganisms is closely tied to the diversification of transporter functions. Horizontal gene transfer (HGT) has been suggested to produce genetic variation that can lead to adaptation, allowing lineages to acquire traits and expand niche ranges. Transporter genes often encode single-gene phenotypes and tend to have low protein-protein interaction complexity and, as such, are potential candidates for HGT. Here we test the idea that HGT has underpinned the expansion of metabolic potential and substrate utilization via transfer of transporter-encoding genes. Using phylogenomics, we identify seven cases of transporter-gene HGT between fungal phyla, and investigate compatibility, localization, function, and fitness consequences when these genes are expressed in Saccharomyces cerevisiae Using this approach, we demonstrate that the transporters identified can alter how fungi utilize a range of metabolites, including peptides, polyols, and sugars. We then show, for one model gene, that transporter gene acquisition by HGT can significantly alter the fitness landscape of S. cerevisiae We therefore provide evidence that transporter HGT occurs between fungi, alters how fungi can acquire metabolites, and can drive gain in fitness. We propose a "transporter-gene acquisition ratchet," where transporter repertoires are continually augmented by duplication, HGT, and differential loss, collectively acting to overwrite, fine-tune, and diversify the complement of transporters present in a genome.

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.

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.

Expanded host and geographic range of tadpole associations with the Severe Perkinsea Infection group.

Severe Perkinsea infection is an emerging disease of amphibians, specifically tadpoles. Disease presentation correlates with liver infections of a subclade of Perkinsea (Alveolata) protists, named Pathogenic Perkinsea Clade (PPC). Tadpole mortality events associated with PPC infections have been reported across North America, from Alaska to Florida. Here, we investigate the geographic and host range of PPC associations in seemingly healthy tadpoles sampled from Panama, a biogeographic provenance critically affected by amphibian decline. To complement this work, we also investigate a mortality event among Hyla arborea tadpoles in captive-bred UK specimens. PPC SSU rDNA was detected in 10 of 81 Panama tadpoles tested, and H. arborea tadpoles from the UK. Phylogenies of the Perkinsea SSU rDNA sequences demonstrate they are highly similar to PPC sequences sampled from mortality events in the USA, and phylogenetic analysis of tadpole mitochondrial SSU rDNA demonstrates, for the first time, PPC associations in diverse hylids. These data provide further understanding of the biogeography and host range of this putative pathogenic group, factors likely to be important for conservation planning.

Host-derived viral transporter protein for nitrogen uptake in infected marine phytoplankton.

Phytoplankton community structure is shaped by both bottom-up factors, such as nutrient availability, and top-down processes, such as predation. Here we show that marine viruses can blur these distinctions, being able to amend how host cells acquire nutrients from their environment while also predating and lysing their algal hosts. Viral genomes often encode genes derived from their host. These genes may allow the virus to manipulate host metabolism to improve viral fitness. We identify in the genome of a phytoplankton virus, which infects the small green alga Ostreococcus tauri, a host-derived ammonium transporter. This gene is transcribed during infection and when expressed in yeast mutants the viral protein is located to the plasma membrane and rescues growth when cultured with ammonium as the sole nitrogen source. We also show that viral infection alters the nature of nitrogen compound uptake of host cells, by both increasing substrate affinity and allowing the host to access diverse nitrogen sources. This is important because the availability of nitrogen often limits phytoplankton growth. Collectively, these data show that a virus can acquire genes encoding nutrient transporters from a host genome and that expression of the viral gene can alter the nutrient uptake behavior of host cells. These results have implications for understanding how viruses manipulate the physiology and ecology of phytoplankton, influence marine nutrient cycles, and act as vectors for horizontal gene transfer.

Ancestral Function and Diversification of a Horizontally Acquired Oomycete Carboxylic Acid Transporter.

Horizontal gene transfer (HGT) can equip organisms with novel genes, expanding the repertoire of genetic material available for evolutionary innovation and allowing recipient lineages to colonize new environments. However, few studies have characterized the functions of HGT genes experimentally or examined postacquisition functional divergence. Here, we report the use of ancestral sequence reconstruction and heterologous expression in Saccharomyces cerevisiae to examine the evolutionary history of an oomycete transporter gene family that was horizontally acquired from fungi. We demonstrate that the inferred ancestral oomycete HGT transporter proteins and their extant descendants transport dicarboxylic acids which are intermediates of the tricarboxylic acid cycle. The substrate specificity profile of the most ancestral protein has largely been retained throughout the radiation of oomycetes, including in both plant and animal pathogens and in a free-living saprotroph, indicating that the ancestral HGT transporter function has been maintained by selection across a range of different lifestyles. No evidence of neofunctionalization in terms of substrate specificity was detected for different HGT transporter paralogues which have different patterns of temporal expression. However, a striking expansion of substrate range was observed for one plant pathogenic oomycete, with a HGT derived paralogue from Pythium aphanidermatum encoding a protein that enables tricarboxylic acid uptake in addition to dicarboxylic acid uptake. This demonstrates that HGT acquisitions can provide functional additions to the recipient proteome as well as the foundation material for the evolution of expanded protein functions.

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

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

Phylogeny, Evidence for a Cryptic Plastid, and Distribution of Chytriodinium Parasites (Dinophyceae) Infecting Copepods.

Spores of the dinoflagellate Chytriodinium are known to infest copepod eggs causing their lethality. Despite the potential to control the population of such an ecologically important host, knowledge about Chytriodinium parasites is limited: we know little about phylogeny, parasitism, abundance, or geographical distribution. We carried out genome sequence surveys on four manually isolated sporocytes from the same sporangium, which seemed to be attached to a copepod nauplius, to analyze the phylogenetic position of Chytriodinium based on SSU and concatenated SSU/LSU rRNA gene sequences, and also characterize two genes related to the plastidial heme pathway, hemL and hemY. The results suggest the presence of a cryptic plastid in Chytriodinium and a photosynthetic ancestral state of the parasitic Chytriodinium/Dissodinium clade. Finally, by mapping Tara Oceans V9 SSU amplicon data to the recovered SSU rRNA gene sequences from the sporocytes, we show that globally, Chytriodinium parasites are most abundant within the pico/nano- and mesoplankton of the surface ocean and almost absent within microplankton, a distribution indicating that they generally exist either as free-living spores or host-associated sporangia.

Identifying protist consumers of photosynthetic picoeukaryotes in the surface ocean using stable isotope probing.

Photosynthetic picoeukaryotes contribute a significant fraction of primary production in the upper ocean. Micromonas pusilla is an ecologically relevant photosynthetic picoeukaryote, abundantly and widely distributed in marine waters. Grazing by protists may control the abundance of picoeukaryotes such as M. pusilla, but the diversity of the responsible grazers is poorly understood. To identify protists consuming photosynthetic picoeukaryotes in a productive North Pacific Ocean region, we amended seawater with living 15 N, 13 C-labelled M. pusilla cells in a 24-h replicated bottle experiment. DNA stable isotope probing, combined with high-throughput sequencing of V4 hypervariable regions from 18S rRNA gene amplicons (Tag-SIP), identified 19 operational taxonomic units (OTUs) of microbial eukaryotes that consumed M. pusilla. These OTUs were distantly related to cultured taxa within the dinoflagellates, ciliates, stramenopiles (MAST-1C and MAST-3 clades) and Telonema flagellates, thus, far known only from their environmental 18S rRNA gene sequences. Our discovery of eukaryotic prey consumption by MAST cells confirms that their trophic role in marine microbial food webs includes grazing upon picoeukaryotes. Our study provides new experimental evidence directly linking the genetic identity of diverse uncultivated microbial eukaryotes to the consumption of picoeukaryotic phytoplankton in the upper ocean.

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.

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.