Phylogenomics demonstrates that breviate flagellates are related to opisthokonts and apusomonads

Is Over-parameterization a Problem for Profile Mixture Models?

Biochemical constraints on the admissible amino acids at specific sites in proteins lead to heterogeneity of the amino acid substitution process over sites in alignments. It is well known that phylogenetic models of protein sequence evolution that do not account for site heterogeneity are prone to long-branch attraction (LBA) artifacts. Profile mixture models were developed to model heterogeneity of preferred amino acids at sites via a finite distribution of site classes each with a distinct set of equilibrium amino acid frequencies. However, it is unknown whether the large number of parameters in such models associated with the many amino acid frequency vectors can adversely affect tree topology estimates because of over-parameterization. Here, we demonstrate theoretically that for long sequences, over-parameterization does not create problems for estimation with profile mixture models. Under mild conditions, tree, amino acid frequencies, and other model parameters converge to true values as sequence length increases, even when there are large numbers of components in the frequency profile distributions. Because large sample theory does not necessarily imply good behavior for shorter alignments we explore the performance of these models with short alignments simulated with tree topologies that are prone to LBA artifacts. We find that over-parameterization is not a problem for complex profile mixture models even when there are many amino acid frequency vectors. In fact, simple models with few site classes behave poorly. Interestingly, we also found that misspecification of the amino acid frequency vectors does not lead to increased LBA artifacts as long as the estimated cumulative distribution function of the amino acid frequencies at sites adequately approximates the true one. In contrast, misspecification of the amino acid exchangeability rates can severely negatively affect parameter estimation. Finally, we explore the effects of including in the profile mixture model an additional "F-class" representing the overall frequencies of amino acids in the data set. Surprisingly, the F-class does not help parameter estimation significantly and can decrease the probability of correct tree estimation, depending on the scenario, even though it tends to improve likelihood scores.

Create, Analyze, and Visualize Phylogenomic Datasets Using PhyloFisher

PhyloFisher is a software package written primarily in Python3 that can be used for the creation, analysis, and visualization of phylogenomic datasets that consist of protein sequences from eukaryotic organisms. Unlike many existing phylogenomic pipelines, PhyloFisher comes with a manually curated database of 240 protein-coding genes, a subset of a previous phylogenetic dataset sampled from 304 eukaryotic taxa. The software package can also utilize a user-created database of eukaryotic proteins, which may be more appropriate for shallow evolutionary questions. PhyloFisher is also equipped with a set of utilities to aid in running routine analyses, such as the prediction of alternative genetic codes, removal of genes and/or taxa based on occupancy/completeness of the dataset, testing for amino acid compositional heterogeneity among sequences, removal of heterotachious and/or fast-evolving sites, removal of fast-evolving taxa, supermatrix creation from randomly resampled genes, and supermatrix creation from nucleotide sequences. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Constructing a phylogenomic dataset Basic Protocol 2: Performing phylogenomic analyses Support Protocol 1: Installing PhyloFisher Support Protocol 2: Creating a custom phylogenomic database.

Molecular and morphological characterization of four new ancyromonad genera and proposal for an updated taxonomy of the Ancyromonadida

Ancyromonads are small biflagellated protists with a bean-shaped morphology. They are cosmopolitan in marine, freshwater, and soil environments, where they attach to surfaces while feeding on bacteria. These poorly known grazers stand out by their uncertain phylogenetic position in the tree of eukaryotes, forming a deep-branching "orphan" lineage that is considered key to a better understanding of the early evolution of eukaryotes. Despite their ecological and evolutionary interest, only limited knowledge exists about their true diversity. Here, we aimed to characterize ancyromonads better by integrating environmental surveys with behavioral observation and description of cell morphology, for which sample isolation and culturing are indispensable. We studied 18 ancyromonad strains, including 14 new isolates and seven new species. We described three new and genetically divergent genera: Caraotamonas, Nyramonas, and Olneymonas, together encompassing four species. The remaining three new species belong to the already-known genera Fabomonas and Ancyromonas. We also raised Striomonas, formerly a subgenus of Nutomonas, to full genus status, on morphological and phylogenetic grounds. We studied the morphology of diverse ancyromonads under light and electron microscopy and carried out molecular phylogenetic analyses, also including 18S rRNA gene sequences from several environmental surveys. Based on these analyses, we have updated the taxonomy of Ancyromonadida.

Compositionally Constrained Sites Drive Long-Branch Attraction

Accurate phylogenies are fundamental to our understanding of the pattern and process of evolution. Yet, phylogenies at deep evolutionary timescales, with correspondingly long branches, have been fraught with controversy resulting from conflicting estimates from models with varying complexity and goodness of fit. Analyses of historical as well as current empirical datasets, such as alignments including Microsporidia, Nematoda, or Platyhelminthes, have demonstrated that inadequate modeling of across-site compositional heterogeneity, which is the result of biochemical constraints that lead to varying patterns of accepted amino acids along sequences, can lead to erroneous topologies that are strongly supported. Unfortunately, models that adequately account for across-site compositional heterogeneity remain computationally challenging or intractable for an increasing fraction of contemporary datasets. Here, we introduce "compositional constraint analysis," a method to investigate the effect of site-specific constraints on amino acid composition on phylogenetic inference. We show that more constrained sites with lower diversity and less constrained sites with higher diversity exhibit ostensibly conflicting signals under models ignoring across-site compositional heterogeneity that lead to long-branch attraction artifacts and demonstrate that more complex models accounting for across-site compositional heterogeneity can ameliorate this bias. We present CAT-posterior mean site frequencies (PMSF), a pipeline for diagnosing and resolving phylogenetic bias resulting from inadequate modeling of across-site compositional heterogeneity based on the CAT model. CAT-PMSF is robust against long-branch attraction in all alignments we have examined. We suggest using CAT-PMSF when convergence of the CAT model cannot be assured. We find evidence that compositionally constrained sites are driving long-branch attraction in two metazoan datasets and recover evidence for Porifera as the sister group to all other animals. [Animal phylogeny; cross-site heterogeneity; long-branch attraction; phylogenomics.].

Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

Expanding the molecular and morphological diversity of Apusomonadida, a deep-branching group of gliding bacterivorous protists

Apusomonads are cosmopolitan bacterivorous biflagellate protists usually gliding on freshwater and marine sediment or wet soils. These nanoflagellates form a sister lineage to opisthokonts and may have retained ancestral features helpful to understanding the early evolution of this large supergroup. Although molecular environmental analyses indicate that apusomonads are genetically diverse, few species have been described. Here, we morphologically characterize 11 new apusomonad strains. Based on molecular phylogenetic analyses of the rRNA gene operon, we describe four new strains of the known species Multimonas media, Podomonas capensis, Apusomonas proboscidea, and Apusomonas australiensis, and rename Thecamonas oxoniensis as Mylnikovia oxoniensis n. gen., n. comb. Additionally, we describe four new genera and six new species: Catacumbia lutetiensis n. gen. n. sp., Cavaliersmithia chaoae n. gen. n. sp., Singekia montserratensis n. gen. n. sp., Singekia franciliensis n. gen. n. sp., Karpovia croatica n. gen. n. sp., and Chelonemonas dolani n. sp. Our comparative analysis suggests that apusomonad ancestor was a fusiform biflagellate with a dorsal pellicle, a plastic ventral surface, and a sleeve covering the anterior flagellum, that thrived in marine, possibly oxygen-poor, environments. It likely had a complex cell cycle with dormant and multiple fission stages, and sex. Our results extend known apusomonad diversity, allow updating their taxonomy, and provide elements to understand early eukaryotic evolution.

Evolutionary diversification of the autophagy-related ubiquitin-like conjugation systems

Two autophagy-related (ATG) ubiquitin-like conjugation systems, the ATG12 and ATG8 systems, play important roles in macroautophagy. While multiple duplications and losses of the ATG conjugation system proteins are found in different lineages, the extent to which the underlying systems diversified across eukaryotes is not fully understood. Here, in order to understand the evolution of the ATG conjugation systems, we constructed a transcriptome database consisting of 94 eukaryotic species covering major eukaryotic clades and systematically identified ATG conjugation system components. Both ATG10 and the C-terminal glycine of ATG12 are essential for the canonical ubiquitin-like conjugation of ATG12 and ATG5. However, loss of ATG10 or the C-terminal glycine of ATG12 occurred at least 16 times in a wide range of lineages, suggesting that possible covalent-to-non-covalent transition is not limited to the species that we previously reported such as Alveolata and some yeast species. Some species have only the ATG8 system (with conjugation enzymes) or only ATG8 (without conjugation enzymes). More than 10 species have ATG8 homologs without the conserved C-terminal glycine, and Tetrahymena has an ATG8 homolog with a predicted transmembrane domain, which may be able to anchor to the membrane independent of the ATG conjugation systems. We discuss the possibility that the ancestor of the ATG12 and ATG8 systems is more similar to ATG8. Overall, our study offers a whole picture of the evolution and diversity of the ATG conjugation systems among eukaryotes, and provides evidence that functional diversifications of the systems are more common than previously thought.Abbreviations: APEAR: ATG8-PE association region; ATG: autophagy-related; LIR: LC3-interacting region; NEDD8: neural precursor cell expressed, developmentally down-regulated gene 8; PE: phosphatidylethanolamine; SAMP: small archaeal modifier protein; SAR: Stramenopiles, Alveolata, and Rhizaria; SMC: structural maintenance of chromosomes; SUMO: small ubiquitin like modifier; TACK: Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota; UBA: ubiquitin like modifier activating enzyme; UFM: ubiquitin fold modifier; URM: ubiquitin related modifier.

New insights on the evolutionary relationships between the major lineages of Amoebozoa

The supergroup Amoebozoa unites a wide diversity of amoeboid organisms and encompasses enigmatic lineages that have been recalcitrant to modern phylogenetics. Deep divergences, taxonomic placement of some key taxa and character evolution in the group largely remain poorly elucidated or controversial. We surveyed available Amoebozoa genomes and transcriptomes to mine conserved putative single copy genes, which were used to enrich gene sampling and generate the largest supermatrix in the group to date; encompassing 824 genes, including gene sequences not previously analyzed. We recovered a well-resolved and supported tree of Amoebozoa, revealing novel deep level relationships and resolving placement of enigmatic lineages congruent with morphological data. In our analysis the deepest branching group is Tubulinea. A recent proposed major clade Tevosa, uniting Evosea and Tubulinea, is not supported. Based on the new phylogenetic tree, paleoecological and paleontological data as well as data on the biology of presently living amoebozoans, we hypothesize that the evolution of Amoebozoa probably was driven by adaptive responses to a changing environment, where successful survival and predation resulted from a capacity to disrupt and graze on microbial mats-a dominant ecosystem of the mid-Proterozoic period of the Earth history.

Ciliary transition zone evolution and the root of the eukaryote tree: implications for opisthokont origin and classification of kingdoms Protozoa, Plantae, and Fungi

I thoroughly discuss ciliary transition zone (TZ) evolution, highlighting many overlooked evolutionarily significant ultrastructural details. I establish fundamental principles of TZ ultrastructure and evolution throughout eukaryotes, inferring unrecognised ancestral TZ patterns for Fungi, opisthokonts, and Corticata (i.e., kingdoms Plantae and Chromista). Typical TZs have a dense transitional plate (TP), with a previously overlooked complex lattice as skeleton. I show most eukaryotes have centriole/TZ junction acorn-V filaments (whose ancestral function was arguably supporting central pair microtubule-nucleating sites; I discuss their role in centriole growth). Uniquely simple malawimonad TZs (without TP, simpler acorn) pinpoint the eukaryote tree's root between them and TP-bearers, highlighting novel superclades. I integrate TZ/ciliary evolution with the best multiprotein trees, naming newly recognised major eukaryote clades and revise megaclassification of basal kingdom Protozoa. Recent discovery of non-photosynthetic phagotrophic flagellates with genome-free plastids (Rhodelphis), the sister group to phylum Rhodophyta (red algae), illuminates plant and chromist early evolution. I show previously overlooked marked similarities in cell ultrastructure between Rhodelphis and Picomonas, formerly considered an early diverging chromist. In both a nonagonal tube lies between their TP and an annular septum surrounding their 9+2 ciliary axoneme. Mitochondrial dense condensations and mitochondrion-linked smooth endomembrane cytoplasmic partitioning cisternae further support grouping Picomonadea and Rhodelphea as new plant phylum Pararhoda. As Pararhoda/Rhodophyta form a robust clade on site-heterogeneous multiprotein trees, I group Pararhoda and Rhodophyta as new infrakingdom Rhodaria of Plantae within subkingdom Biliphyta, which also includes Glaucophyta with fundamentally similar TZ, uniquely in eukaryotes. I explain how biliphyte TZs generated viridiplant stellate-structures.

The closest lineage of Archaeplastida is revealed by phylogenomics analyses that include Microheliella maris

By clarifying the phylogenetic positions of 'orphan' protists (unicellular micro-eukaryotes with no affinity to extant lineages), we may uncover the novel affiliation between two (or more) major lineages in eukaryotes. Microheliella maris was an orphan protist, which failed to be placed within the previously described lineages by pioneering phylogenetic analyses. In this study, we analysed a 319-gene alignment and demonstrated that M. maris represents a basal lineage of one of the major eukaryotic lineages, Cryptista. We here propose a new clade name 'Pancryptista' for Cryptista plus M. maris. The 319-gene analyses also indicated that M. maris is a key taxon to recover the monophyly of Archaeplastida and the sister relationship between Archaeplastida and Pancryptista, which is collectively called 'CAM clade' here. Significantly, Cryptophyceae tend to be attracted to Rhodophyta depending on the taxon sampling (ex., in the absence of M. maris and Rhodelphidia) and the particular phylogenetic 'signal' most likely hindered the stable recovery of the monophyly of Archaeplastida in previous studies.

Neofunctionalization of Glycolytic Enzymes: An Evolutionary Route to Plant Parasitism in the Oomycete Phytophthora nicotianae

Oomycetes, of the genus Phytophthora, comprise of some of the most devastating plant pathogens. Parasitism of Phytophthora results from evolution from an autotrophic ancestor and adaptation to a wide range of environments, involving metabolic adaptation. Sequence mining showed that Phytophthora spp. display an unusual repertoire of glycolytic enzymes, made of multigene families and enzyme replacements. To investigate the impact of these gene duplications on the biology of Phytophthora and, eventually, identify novel functions associated to gene expansion, we focused our study on the first glycolytic step on P. nicotianae, a broad host range pathogen. We reveal that this step is committed by a set of three glucokinase types that differ by their structure, enzymatic properties, and evolutionary histories. In addition, they are expressed differentially during the P. nicotianae life cycle, including plant infection. Last, we show that there is a strong association between the expression of a glucokinase member in planta and extent of plant infection. Together, these results suggest that metabolic adaptation is a component of the processes underlying evolution of parasitism in Phytophthora, which may possibly involve the neofunctionalization of metabolic enzymes.

Quantitative Insights into the Contribution of Nematocysts to the Adaptive Success of Cnidarians Based on Proteomic Analysis

Nematocysts are secretory organelles in cnidarians that play important roles in predation, defense, locomotion, and host invasion. However, the extent to which nematocysts contribute to adaptation and the mechanisms underlying nematocyst evolution are unclear. Here, we investigated the role of the nematocyst in cnidarian evolution based on eight nematocyst proteomes and 110 cnidarian transcriptomes/genomes. We detected extensive species-specific adaptive mutations in nematocyst proteins (NEMs) and evidence for decentralized evolution, in which most evolutionary events involved non-core NEMs, reflecting the rapid diversification of NEMs in cnidarians. Moreover, there was a 33-55 million year macroevolutionary lag between nematocyst evolution and the main phases of cnidarian diversification, suggesting that the nematocyst can act as a driving force in evolution. Quantitative analysis revealed an excess of adaptive changes in NEMs and enrichment for positively selected conserved NEMs. Together, these findings suggest that nematocysts may be key to the adaptive success of cnidarians and provide a reference for quantitative analyses of the roles of phenotypic novelties in adaptation.

The origin of animals: an ancestral reconstruction of the unicellular-to-multicellular transition

How animals evolved from a single-celled ancestor, transitioning from a unicellular lifestyle to a coordinated multicellular entity, remains a fascinating question. Key events in this transition involved the emergence of processes related to cell adhesion, cell–cell communication and gene regulation. To understand how these capacities evolved, we need to reconstruct the features of both the last common multicellular ancestor of animals and the last unicellular ancestor of animals. In this review, we summarize recent advances in the characterization of these ancestors, inferred by comparative genomic analyses between the earliest branching animals and those radiating later, and between animals and their closest unicellular relatives. We also provide an updated hypothesis regarding the transition to animal multicellularity, which was likely gradual and involved the use of gene regulatory mechanisms in the emergence of early developmental and morphogenetic plans. Finally, we discuss some new avenues of research that will complement these studies in the coming years.

Description of Imasa heleensis, gen. nov., sp. nov. (Imasidae, fam. nov.), a Deep-Branching Marine Malawimonad and Possible Key Taxon in Understanding Early Eukaryotic Evolution

Malawimonadida is a deep-level (arguably "kingdom-scale") lineage of eukaryotes whose phylogenetic affinities are uncertain but of great evolutionary interest, as the group is suspected to branch close to the root of the tree of eukaryotes. Part of the difficulty in placing Malawimonadida phylogenetically is its tiny circumscription: at present, it comprises only two described and one cultured but undescribed species, all of them are freshwater suspension-feeding nanoflagellates. In this study, we cultivated and characterised Imasa heleensis gen. nov., sp. nov. (Imasidae fam. nov.), the first marine malawimonad to be described. Light and electron microscopy observations show that Imasa is largely similar to other malawimonads, but more frequently adheres to the substrate, often by means of a pliable posterior extension. Phylogenetic analyses based on two ribosomal RNA genes and four translated protein-coding genes using three different taxon sets place Imasa as sister to the three freshwater malawimonad strains with strong support. Imasa's mitochondrial genome is circular-mapping and shows a similar gene complement to other known malawimonads. We conclude that Imasa represents an important expansion of the range of taxa available for future evolutionary study.

Progress towards the Tree of Eukaryotes

Developing a detailed understanding of how all known forms of life are related to one another in the tree of life has been a major preoccupation of biology since the idea of tree-like evolution first took hold. Since most life is microbial, our intuitive use of morphological comparisons to infer relatedness only goes so far, and molecular sequence data, most recently from genomes and transcriptomes, has been the primary means to infer these relationships. For prokaryotes this presented new challenges, since the degree of horizontal gene transfer led some to question the tree-like depiction of evolution altogether. Most eukaryotes are also microbial, but in contrast to prokaryotic life, the application of large-scale molecular data to the tree of eukaryotes has largely been a constructive process, leading to a small number of very diverse lineages, or 'supergroups'. The tree is not completely resolved, and contentious problems remain, but many well-established supergroups now encompass much more diversity than the traditional kingdoms. Some of the most exciting recent developments come from the discovery of branches in the tree that we previously had no inkling even existed, many of which are of great ecological or evolutionary interest. These new branches highlight the need for more exploration, by high-throughput molecular surveys, but also more traditional means of observations and cultivation.

Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

The Relative Importance of Modeling Site Pattern Heterogeneity Versus Partition-Wise Heterotachy in Phylogenomic Inference

Large taxa-rich genome-scale data sets are often necessary for resolving ancient phylogenetic relationships. But accurate phylogenetic inference requires that they are analyzed with realistic models that account for the heterogeneity in substitution patterns amongst the sites, genes and lineages. Two kinds of adjustments are frequently used: models that account for heterogeneity in amino acid frequencies at sites in proteins, and partitioned models that accommodate the heterogeneity in rates (branch lengths) among different proteins in different lineages (protein-wise heterotachy). Although partitioned and site-heterogeneous models are both widely used in isolation, their relative importance to the inference of correct phylogenies has not been carefully evaluated. We conducted several empirical analyses and a large set of simulations to compare the relative performances of partitioned models, site-heterogeneous models, and combined partitioned site heterogeneous models. In general, site-homogeneous models (partitioned or not) performed worse than site heterogeneous, except in simulations with extreme protein-wise heterotachy. Furthermore, simulations using empirically-derived realistic parameter settings showed a marked long-branch attraction (LBA) problem for analyses employing protein-wise partitioning even when the generating model included partitioning. This LBA problem results from a small sample bias compounded over many single protein alignments. In some cases, this problem was ameliorated by clustering similarly-evolving proteins together into larger partitions using the PartitionFinder method. Similar results were obtained under simulations with larger numbers of taxa or heterogeneity in simulating topologies over genes. For an empirical Microsporidia test data set, all but one tested site-heterogeneous models (with or without partitioning) obtain the correct Microsporidia+Fungi grouping, whereas site-homogenous models (with or without partitioning) did not. The single exception was the fully partitioned site-heterogeneous analysis that succumbed to the compounded small sample LBA bias. In general unless protein-wise heterotachy effects are extreme, it is more important to model site-heterogeneity than protein-wise heterotachy in phylogenomic analyses. Complete protein-wise partitioning should be avoided as it can lead to a serious LBA bias. In cases of extreme protein-wise heterotachy, approaches that cluster similarly-evolving proteins together and coupled with site-heterogeneous models work well for phylogenetic estimation.

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.

All Eukaryotes Are Sexual, unless Proven Otherwise: Many So-Called Asexuals Present Meiotic Machinery and Might Be Able to Have Sex

Here a wide distribution of meiotic machinery is shown, indicating the occurrence of sexual processes in all major eukaryotic groups, without exceptions, including the putative "asexuals." Meiotic machinery has evolved from archaeal DNA repair machinery by means of ancestral gene duplications. Sex is very conserved and widespread in eukaryotes, even though its evolutionary importance is still a matter of debate. The main processes in sex are plasmogamy, followed by karyogamy and meiosis. Meiosis is fundamentally a chromosomal process, which implies recombination and ploidy reduction. Several eukaryotic lineages are proposed to be asexual because their sexual processes are never observed, but presumed asexuality correlates with lack of study. The authors stress the complete lack of meiotic proteins in nucleomorphs and their almost complete loss in the fungus Malassezia. Inversely, complete sets of meiotic proteins are present in fungal groups Glomeromycotina, Trichophyton, and Cryptococcus. Endosymbiont Perkinsela and endoparasitic Microsporidia also present meiotic proteins.

Asgard archaea capable of anaerobic hydrocarbon cycling

Large reservoirs of natural gas in the oceanic subsurface sustain complex communities of anaerobic microbes, including archaeal lineages with potential to mediate oxidation of hydrocarbons such as methane and butane. Here we describe a previously unknown archaeal phylum, Helarchaeota, belonging to the Asgard superphylum and with the potential for hydrocarbon oxidation. We reconstruct Helarchaeota genomes from metagenomic data derived from hydrothermal deep-sea sediments in the hydrocarbon-rich Guaymas Basin. The genomes encode methyl-CoM reductase-like enzymes that are similar to those found in butane-oxidizing archaea, as well as several enzymes potentially involved in alkyl-CoA oxidation and the Wood-Ljungdahl pathway. We suggest that members of the Helarchaeota have the potential to activate and subsequently anaerobically oxidize hydrothermally generated short-chain hydrocarbons.

Natural gas reservoirs in the oceanic subsurface sustain complex communities of anaerobic microbes. Here, Seitz et al. describe a previously unknown archaeal phylum, Helarchaeota, belonging to the Asgard superphylum and with the potential for oxidation of hydrothermally generated short-chain hydrocarbons.

Accelerated Estimation of Frequency Classes in Site-Heterogeneous Profile Mixture Models

As a consequence of structural and functional constraints, proteins tend to have site-specific preferences for particular amino acids. Failing to adjust for heterogeneity of frequencies over sites can lead to artifacts in phylogenetic estimation. Site-heterogeneous mixture-models have been developed to address this problem. However, due to prohibitive computational times, maximum likelihood implementations utilize fixed component frequency vectors inferred from sequences in a database that are external to the alignment under analysis. Here, we propose a composite likelihood approach to estimation of component frequencies for a mixture model that directly uses the data from the alignment of interest. In the common case that the number of taxa under study is not large, several adjustments to the default composite likelihood are shown to be necessary. In simulations, the approach is shown to provide large improvements over hierarchical clustering. For empirical data, substantial improvements in likelihoods are found over mixtures using fixed components.

Comparative genomic analysis of the 'pseudofungus' Hyphochytrium catenoides

Eukaryotic microbes have three primary mechanisms for obtaining nutrients and energy: phagotrophy, photosynthesis and osmotrophy. Traits associated with the latter two functions arose independently multiple times in the eukaryotes. The Fungi successfully coupled osmotrophy with filamentous growth, and similar traits are also manifested in the Pseudofungi (oomycetes and hyphochytriomycetes). Both the Fungi and the Pseudofungi encompass a diversity of plant and animal parasites. Genome-sequencing efforts have focused on host-associated microbes (mutualistic symbionts or parasites), providing limited comparisons with free-living relatives. Here we report the first draft genome sequence of a hyphochytriomycete ‘pseudofungus’; Hyphochytrium catenoides. Using phylogenomic approaches, we identify genes of recent viral ancestry, with related viral derived genes also present on the genomes of oomycetes, suggesting a complex history of viral coevolution and integration across the Pseudofungi. H. catenoides has a complex life cycle involving diverse filamentous structures and a flagellated zoospore with a single anterior tinselate flagellum. We use genome comparisons, drug sensitivity analysis and high-throughput culture arrays to investigate the ancestry of oomycete/pseudofungal characteristics, demonstrating that many of the genetic features associated with parasitic traits evolved specifically within the oomycete radiation. Comparative genomics also identified differences in the repertoire of genes associated with filamentous growth between the Fungi and the Pseudofungi, including differences in vesicle trafficking systems, cell-wall synthesis pathways and motor protein repertoire, demonstrating that unique cellular systems underpinned the convergent evolution of filamentous osmotrophic growth in these two eukaryotic groups.

Triplet-pore structure of a highly divergent TOM complex of hydrogenosomes in Trichomonas vaginalis

Mitochondria originated from proteobacterial endosymbionts, and their transition to organelles was tightly linked to establishment of the protein import pathways. The initial import of most proteins is mediated by the translocase of the outer membrane (TOM). Although TOM is common to all forms of mitochondria, an unexpected diversity of subunits between eukaryotic lineages has been predicted. However, experimental knowledge is limited to a few organisms, and so far, it remains unsettled whether the triplet-pore or the twin-pore structure is the generic form of TOM complex. Here, we analysed the TOM complex in hydrogenosomes, a metabolically specialised anaerobic form of mitochondria found in the excavate Trichomonas vaginalis. We demonstrate that the highly divergent β-barrel T. vaginalis TOM (TvTom)40-2 forms a translocation channel to conduct hydrogenosomal protein import. TvTom40-2 is present in high molecular weight complexes, and their analysis revealed the presence of four tail-anchored (TA) proteins. Two of them, Tom36 and Tom46, with heat shock protein (Hsp)20 and tetratricopeptide repeat (TPR) domains, can bind hydrogenosomal preproteins and most likely function as receptors. A third subunit, Tom22-like protein, has a short cis domain and a conserved Tom22 transmembrane segment but lacks a trans domain. The fourth protein, hydrogenosomal outer membrane protein 19 (Homp19) has no known homology. Furthermore, our data indicate that TvTOM is associated with sorting and assembly machinery (Sam)50 that is involved in β-barrel assembly. Visualisation of TvTOM by electron microscopy revealed that it forms three pores and has an unconventional skull-like shape. Although TvTOM seems to lack Tom7, our phylogenetic profiling predicted Tom7 in free-living excavates. Collectively, our results suggest that the triplet-pore TOM complex, composed of three conserved subunits, was present in the last common eukaryotic ancestor (LECA), while receptors responsible for substrate binding evolved independently in different eukaryotic lineages.

The highly divergent outer membrane translocase (TOM) from the Trichomonas hydrogenosome (an organelle related to mitochondria) is composed of conserved core and lineage-specific subunits, and has an unconventional skull-like triplet-pore structure.

Mitochondria carry out many vital functions in the eukaryotic cells, from energy metabolism to programmed cell death. These organelles descended from bacterial endosymbionts, and during their evolution, the cell established a mechanism to transport nuclear-encoded proteins into mitochondria. Embedded in the mitochondrial outer membrane is a molecular machine, known as the translocase of the outer membrane (TOM) complex, that plays a key role in protein import and biogenesis of the organelle. Here, we provide evidence that the TOM complex of hydrogenosomes, a metabolically specialised anaerobic form of mitochondria in Trichomonas vaginalis, is composed of highly divergent core subunits and lineage-specific peripheral subunits. Despite the evolutionary distance, the T. vaginalis TOM (TvTOM) complex has a conserved triplet-pore structure but with a unique skull-like shape suggesting that the TOM in the early mitochondrion could have formed three pores. Our results contribute to a better understanding of the evolution and adaptation of protein import machinery in anaerobic forms of mitochondria.

Hemimastigophora is a novel supra-kingdom-level lineage of eukaryotes

Almost all eukaryote life forms have now been placed within one of five to eight supra-kingdom-level groups using molecular phylogenetics^1-4. The 'phylum' Hemimastigophora is probably the most distinctive morphologically defined lineage that still awaits such a phylogenetic assignment. First observed in the nineteenth century, hemimastigotes are free-living predatory protists with two rows of flagella and a unique cell architecture^5-7; to our knowledge, no molecular sequence data or cultures are currently available for this group. Here we report phylogenomic analyses based on high-coverage, cultivation-independent transcriptomics that place Hemimastigophora outside of all established eukaryote supergroups. They instead comprise an independent supra-kingdom-level lineage that most likely forms a sister clade to the 'Diaphoretickes' half of eukaryote diversity (that is, the 'stramenopiles, alveolates and Rhizaria' supergroup (Sar), Archaeplastida and Cryptista, as well as other major groups). The previous ranking of Hemimastigophora as a phylum understates the evolutionary distinctiveness of this group, which has considerable importance for investigations into the deep-level evolutionary history of eukaryotic life-ranging from understanding the origins of fundamental cell systems to placing the root of the tree. We have also established the first culture of a hemimastigote (Hemimastix kukwesjijk sp. nov.), which will facilitate future genomic and cell-biological investigations into eukaryote evolution and the last eukaryotic common ancestor.

Comparative Genomics Supports Sex and Meiosis in Diverse Amoebozoa

Sex and reproduction are often treated as a single phenomenon in animals and plants, as in these organisms reproduction implies mixis and meiosis. In contrast, sex and reproduction are independent biological phenomena that may or may not be linked in the majority of other eukaryotes. Current evidence supports a eukaryotic ancestor bearing a mating type system and meiosis, which is a process exclusive to eukaryotes. Even though sex is ancestral, the literature regarding life cycles of amoeboid lineages depicts them as asexual organisms. Why would loss of sex be common in amoebae, if it is rarely lost, if ever, in plants and animals, as well as in fungi? One way to approach the question of meiosis in the "asexuals" is to evaluate the patterns of occurrence of genes for the proteins involved in syngamy and meiosis. We have applied a comparative genomic approach to study the occurrence of the machinery for plasmogamy, karyogamy, and meiosis in Amoebozoa, a major amoeboid supergroup. Our results support a putative occurrence of syngamy and meiotic processes in all major amoebozoan lineages. We conclude that most amoebozoans may perform mixis, recombination, and ploidy reduction through canonical meiotic processes. The present evidence indicates the possibility of sexual cycles in many lineages traditionally held as asexual.

EukRef: Phylogenetic curation of ribosomal RNA to enhance understanding of eukaryotic diversity and distribution

Environmental sequencing has greatly expanded our knowledge of micro-eukaryotic diversity and ecology by revealing previously unknown lineages and their distribution. However, the value of these data is critically dependent on the quality of the reference databases used to assign an identity to environmental sequences. Existing databases contain errors and struggle to keep pace with rapidly changing eukaryotic taxonomy, the influx of novel diversity, and computational challenges related to assembling the high-quality alignments and trees needed for accurate characterization of lineage diversity. EukRef (eukref.org) is an ongoing community-driven initiative that addresses these challenges by bringing together taxonomists with expertise spanning the eukaryotic tree of life and microbial ecologists, who use environmental sequence data to develop reliable reference databases across the diversity of microbial eukaryotes. EukRef organizes and facilitates rigorous mining and annotation of sequence data by providing protocols, guidelines, and tools. The EukRef pipeline and tools allow users interested in a particular group of microbial eukaryotes to retrieve all sequences belonging to that group from International Nucleotide Sequence Database Collaboration (INSDC) (GenBank, the European Nucleotide Archive [ENA], or the DNA DataBank of Japan [DDBJ]), to place those sequences in a phylogenetic tree, and to curate taxonomic and environmental information for the group. We provide guidelines to facilitate the process and to standardize taxonomic annotations. The final outputs of this process are (1) a reference tree and alignment, (2) a reference sequence database, including taxonomic and environmental information, and (3) a list of putative chimeras and other artifactual sequences. These products will be useful for the broad community as they become publicly available (at eukref.org) and are shared with existing reference databases.

Multigene phylogeny and cell evolution of chromist infrakingdom Rhizaria: contrasting cell organisation of sister phyla Cercozoa and Retaria

Infrakingdom Rhizaria is one of four major subgroups with distinct cell body plans that comprise eukaryotic kingdom Chromista. Unlike other chromists, Rhizaria are mostly heterotrophic flagellates, amoebae or amoeboflagellates, commonly with reticulose (net-like) or filose (thread-like) feeding pseudopodia; uniquely for eukaryotes, cilia have proximal ciliary transition-zone hub-lattices. They comprise predominantly flagellate phylum Cercozoa and reticulopodial phylum Retaria, whose exact phylogenetic relationship has been uncertain. Given even less clear relationships amongst cercozoan classes, we sequenced partial transcriptomes of seven Cercozoa representing five classes and endomyxan retarian Filoreta marina to establish 187-gene multiprotein phylogenies. Ectoreta (retarian infraphyla Foraminifera, Radiozoa) branch within classical Cercozoa as sister to reticulose Endomyxa. This supports recent transfer of subphylum Endomyxa from Cercozoa to Retaria alongside subphylum Ectoreta which embraces classical retarians where capsules or tests subdivide cells into organelle-containing endoplasm and anastomosing pseudopodial net-like ectoplasm. Cercozoa are more homogeneously filose, often with filose pseudopodia and/or posterior ciliary gliding motility: zooflagellate Helkesimastix and amoeboid Guttulinopsis form a strongly supported clade, order Helkesida. Cercomonads are polyphyletic (Cercomonadida sister to glissomonads; Paracercomonadida deeper). Thecofilosea are a clade, whereas Imbricatea may not be; Sarcomonadea may be paraphyletic. Helkesea and Metromonadea are successively deeper outgroups within cercozoan subphylum Monadofilosa; subphylum Reticulofilosa (paraphyletic on site-heterogeneous trees) branches earliest, Granofilosea before Chlorarachnea. Our multiprotein trees confirm that Rhizaria are sisters of infrakingdom Halvaria (Alveolata, Heterokonta) within chromist subkingdom Harosa (= SAR); they further support holophyly of chromist subkingdom Hacrobia, and are consistent with holophyly of Chromista as sister of kingdom Plantae. Site-heterogeneous rDNA trees group Kraken with environmental DNA clade 'eSarcomonad', not Paracercomonadida. Ectoretan fossil dates evidence ultrarapid episodic stem sequence evolution. We discuss early rhizarian cell evolution and multigene tree coevolutionary patterns, gene-paralogue evidence for chromist monophyly, and integrate this with fossil evidence for the age of Rhizaria and eukaryote cells, and revise rhizarian classification.

Evolutionary cell biology traces the rise of the exomer complex in Fungi from an ancient eukaryotic component

Cargo is transported from the trans-Golgi Network to the plasma membrane by adaptor complexes, which are pan-eukaryotic components. However, in yeast, cargo can also be exported by the exomer complex, a heterotetrameric protein complex consisting of two copies of Chs5, and any two members of four paralogous proteins (ChAPs). To understand the larger relevance of exomer, its phylogenetic distribution and function outside of yeast need to be explored. We find that the four ChAP proteins are derived from gene duplications after the divergence of Yarrowia from the remaining Saccharomycotina, with BC8 paralogues (Bch2 and Chs6) being more diverged relative to the BB8 paralogues (Bch1 and Bud7), suggesting neofunctionalization. Outside Ascomycota, a single preduplicate ChAP is present in nearly all Fungi and in diverse eukaryotes, but has been repeatedly lost. Chs5, however, is a fungal specific feature, appearing coincidentally with the loss of AP-4. In contrast, the ChAP protein is a wide-spread, yet uncharacterized, membrane-trafficking component, adding one more piece to the increasingly complex machinery deduced as being present in our ancient eukaryotic ancestor.

Enigmatic Diphyllatea eukaryotes: culturing and targeted PacBio RS amplicon sequencing reveals a higher order taxonomic diversity and global distribution

BACKGROUND

The class Diphyllatea belongs to a group of enigmatic unicellular eukaryotes that play a key role in reconstructing the morphological innovation and diversification of early eukaryotic evolution. Despite its evolutionary significance, very little is known about the phylogeny and species diversity of Diphyllatea. Only three species have described morphology, being taxonomically divided by flagella number, two or four, and cell size. Currently, one 18S rRNA Diphyllatea sequence is available, with environmental sequencing surveys reporting only a single partial sequence from a Diphyllatea-like organism. Accordingly, geographical distribution of Diphyllatea based on molecular data is limited, despite morphological data suggesting the class has a global distribution. We here present a first attempt to understand species distribution, diversity and higher order structure of Diphyllatea.

RESULTS

We cultured 11 new strains, characterised these morphologically and amplified their rRNA for a combined 18S-28S rRNA phylogeny. We sampled environmental DNA from multiple sites and designed new Diphyllatea-specific PCR primers for long-read PacBio RSII technology. Near full-length 18S rRNA sequences from environmental DNA, in addition to supplementary Diphyllatea sequence data mined from public databases, resolved the phylogeny into three deeply branching and distinct clades (Diphy I - III). Of these, the Diphy III clade is entirely novel, and in congruence with Diphy II, composed of species morphologically consistent with the earlier described Collodictyon triciliatum. The phylogenetic split between the Diphy I and Diphy II + III clades corresponds with a morphological division of Diphyllatea into bi- and quadriflagellate cell forms.

CONCLUSIONS

This altered flagella composition must have occurred early in the diversification of Diphyllatea and may represent one of the earliest known morphological transitions among eukaryotes. Further, the substantial increase in molecular data presented here confirms Diphyllatea has a global distribution, seemingly restricted to freshwater habitats. Altogether, the results reveal the advantage of combining a group-specific PCR approach and long-read high-throughput amplicon sequencing in surveying enigmatic eukaryote lineages. Lastly, our study shows the capacity of PacBio RS when targeting a protist class for increasing phylogenetic resolution.

Microbial eukaryotes have adapted to hypoxia by horizontal acquisitions of a gene involved in rhodoquinone biosynthesis

Under hypoxic conditions, some organisms use an electron transport chain consisting of only complex I and II (CII) to generate the proton gradient essential for ATP production. In these cases, CII functions as a fumarate reductase that accepts electrons from a low electron potential quinol, rhodoquinol (RQ). To clarify the origins of RQ-mediated fumarate reduction in eukaryotes, we investigated the origin and function of rquA, a gene encoding an RQ biosynthetic enzyme. RquA is very patchily distributed across eukaryotes and bacteria adapted to hypoxia. Phylogenetic analyses suggest lateral gene transfer (LGT) of rquA from bacteria to eukaryotes occurred at least twice and the gene was transferred multiple times amongst protists. We demonstrate that RquA functions in the mitochondrion-related organelles of the anaerobic protist Pygsuia and is correlated with the presence of RQ. These analyses reveal the role of gene transfer in the evolutionary remodeling of mitochondria in adaptation to hypoxia.

Combined morphological and phylogenomic re-examination of malawimonads, a critical taxon for inferring the evolutionary history of eukaryotes

Modern syntheses of eukaryote diversity assign almost all taxa to one of three groups: Amorphea, Diaphoretickes and Excavata (comprising Discoba and Metamonada). The most glaring exception is Malawimonadidae, a group of small heterotrophic flagellates that resemble Excavata by morphology, but branch with Amorphea in most phylogenomic analyses. However, just one malawimonad, Malawimonas jakobiformis, has been studied with both morphological and molecular-phylogenetic approaches, raising the spectre of interpretation errors and phylogenetic artefacts from low taxon sampling. We report a morphological and phylogenomic study of a new deep-branching malawimonad, Gefionella okellyi n. gen. n. sp. Electron microscopy revealed all canonical features of 'typical excavates', including flagellar vanes (as an opposed pair, unlike M. jakobiformis but like many metamonads) and a composite fibre. Initial phylogenomic analyses grouped malawimonads with the Amorphea-related orphan lineage Collodictyon, separate from a Metamonada+Discoba clade. However, support for this topology weakened when more sophisticated evolutionary models were used, and/or fast-evolving sites and long-branching taxa (FS/LB) were excluded. Analyses of '-FS/LB' datasets instead suggested a relationship between malawimonads and metamonads. The 'malawimonad+metamonad signal' in morphological and molecular data argues against a strict Metamonada+Discoba clade (i.e. the predominant concept of Excavata). A Metamonad+Discoba clade should therefore not be assumed when inferring deep-level evolutionary history in eukaryotes.

Modeling Site Heterogeneity with Posterior Mean Site Frequency Profiles Accelerates Accurate Phylogenomic Estimation

Proteins have distinct structural and functional constraints at different sites that lead to site-specific preferences for particular amino acid residues as the sequences evolve. Heterogeneity in the amino acid substitution process between sites is not modeled by commonly used empirical amino acid exchange matrices. Such model misspecification can lead to artefacts in phylogenetic estimation such as long-branch attraction. Although sophisticated site-heterogeneous mixture models have been developed to address this problem in both Bayesian and maximum likelihood (ML) frameworks, their formidable computational time and memory usage severely limits their use in large phylogenomic analyses. Here we propose a posterior mean site frequency (PMSF) method as a rapid and efficient approximation to full empirical profile mixture models for ML analysis. The PMSF approach assigns a conditional mean amino acid frequency profile to each site calculated based on a mixture model fitted to the data using a preliminary guide tree. These PMSF profiles can then be used for in-depth tree-searching in place of the full mixture model. Compared with widely used empirical mixture models with $k$ classes, our implementation of PMSF in IQ-TREE (http://www.iqtree.org) speeds up the computation by approximately $k$/1.5-fold and requires a small fraction of the RAM. Furthermore, this speedup allows, for the first time, full nonparametric bootstrap analyses to be conducted under complex site-heterogeneous models on large concatenated data matrices. Our simulations and empirical data analyses demonstrate that PMSF can effectively ameliorate long-branch attraction artefacts. In some empirical and simulation settings PMSF provided more accurate estimates of phylogenies than the mixture models from which they derive.

A sophisticated, differentiated Golgi in the ancestor of eukaryotes

BACKGROUND

The Golgi apparatus is a central meeting point for the endocytic and exocytic systems in eukaryotic cells, and the organelle's dysfunction results in human disease. Its characteristic morphology of multiple differentiated compartments organized into stacked flattened cisternae is one of the most recognizable features of modern eukaryotic cells, and yet how this is maintained is not well understood. The Golgi is also an ancient aspect of eukaryotes, but the extent and nature of its complexity in the ancestor of eukaryotes is unclear. Various proteins have roles in organizing the Golgi, chief among them being the golgins.

RESULTS

We address Golgi evolution by analyzing genome sequences from organisms which have lost stacked cisternae as a feature of their Golgi and those that have not. Using genomics and immunomicroscopy, we first identify Golgi in the anaerobic amoeba Mastigamoeba balamuthi. We then searched 87 genomes spanning eukaryotic diversity for presence of the most prominent proteins implicated in Golgi structure, focusing on golgins. We show some candidates as animal specific and others as ancestral to eukaryotes.

CONCLUSIONS

None of the proteins examined show a phyletic distribution that correlates with the morphology of stacked cisternae, suggesting the possibility of stacking as an emergent property. Strikingly, however, the combination of golgins conserved among diverse eukaryotes allows for the most detailed reconstruction of the organelle to date, showing a sophisticated Golgi with differentiated compartments and trafficking pathways in the common eukaryotic ancestor.

Kingdom Chromista and its eight phyla: a new synthesis emphasising periplastid protein targeting, cytoskeletal and periplastid evolution, and ancient divergences

In 1981 I established kingdom Chromista, distinguished from Plantae because of its more complex chloroplast-associated membrane topology and rigid tubular multipartite ciliary hairs. Plantae originated by converting a cyanobacterium to chloroplasts with Toc/Tic translocons; most evolved cell walls early, thereby losing phagotrophy. Chromists originated by enslaving a phagocytosed red alga, surrounding plastids by two extra membranes, placing them within the endomembrane system, necessitating novel protein import machineries. Early chromists retained phagotrophy, remaining naked and repeatedly reverted to heterotrophy by losing chloroplasts. Therefore, Chromista include secondary phagoheterotrophs (notably ciliates, many dinoflagellates, Opalozoa, Rhizaria, heliozoans) or walled osmotrophs (Pseudofungi, Labyrinthulea), formerly considered protozoa or fungi respectively, plus endoparasites (e.g. Sporozoa) and all chromophyte algae (other dinoflagellates, chromeroids, ochrophytes, haptophytes, cryptophytes). I discuss their origin, evolutionary diversification, and reasons for making chromists one kingdom despite highly divergent cytoskeletons and trophic modes, including improved explanations for periplastid/chloroplast protein targeting, derlin evolution, and ciliary/cytoskeletal diversification. I conjecture that transit-peptide-receptor-mediated 'endocytosis' from periplastid membranes generates periplastid vesicles that fuse with the arguably derlin-translocon-containing periplastid reticulum (putative red algal trans-Golgi network homologue; present in all chromophytes except dinoflagellates). I explain chromist origin from ancestral corticates and neokaryotes, reappraising tertiary symbiogenesis; a chromist cytoskeletal synapomorphy, a bypassing microtubule band dextral to both centrioles, favoured multiple axopodial origins. I revise chromist higher classification by transferring rhizarian subphylum Endomyxa from Cercozoa to Retaria; establishing retarian subphylum Ectoreta for Foraminifera plus Radiozoa, apicomonad subclasses, new dinozoan classes Myzodinea (grouping Colpovora gen. n., Psammosa), Endodinea, Sulcodinea, and subclass Karlodinia; and ranking heterokont Gyrista as phylum not superphylum.

Origin of animal multicellularity: precursors, causes, consequences-the choanoflagellate/sponge transition, neurogenesis and the Cambrian explosion

Evolving multicellularity is easy, especially in phototrophs and osmotrophs whose multicells feed like unicells. Evolving animals was much harder and unique; probably only one pathway via benthic ‘zoophytes’ with pelagic ciliated larvae allowed trophic continuity from phagocytic protozoa to gut-endowed animals. Choanoflagellate protozoa produced sponges. Converting sponge flask cells mediating larval settling to synaptically controlled nematocysts arguably made Cnidaria. I replace Haeckel's gastraea theory by a sponge/coelenterate/bilaterian pathway: Placozoa, hydrozoan diploblasty and ctenophores were secondary; stem anthozoan developmental mutations arguably independently generated coelomate bilateria and ctenophores. I emphasize animal origin's conceptual aspects (selective, developmental) related to feeding modes, cell structure, phylogeny of related protozoa, sequence evidence, ecology and palaeontology. Epithelia and connective tissue could evolve only by compensating for dramatically lower feeding efficiency that differentiation into non-choanocytes entails. Consequentially, larger bodies enabled filtering more water for bacterial food and harbouring photosynthetic bacteria, together adding more food than cell differentiation sacrificed. A hypothetical presponge of sessile triploblastic sheets (connective tissue sandwiched between two choanocyte epithelia) evolved oogamy through selection for larger dispersive ciliated larvae to accelerate benthic trophic competence and overgrowing protozoan competitors. Extinct Vendozoa might be elaborations of this organismal grade with choanocyte-bearing epithelia, before poriferan water channels and cnidarian gut/nematocysts/synapses evolved.

This article is part of the themed issue ‘Evo-devo in the genomics era, and the origins of morphological diversity’.

Between a Pod and a Hard Test: The Deep Evolution of Amoebae

Amoebozoa is the eukaryotic supergroup sister to Obazoa, the lineage that contains the animals and Fungi, as well as their protistan relatives, and the breviate and apusomonad flagellates. Amoebozoa is extraordinarily diverse, encompassing important model organisms and significant pathogens. Although amoebozoans are integral to global nutrient cycles and present in nearly all environments, they remain vastly understudied. We present a robust phylogeny of Amoebozoa based on broad representative set of taxa in a phylogenomic framework (325 genes). By sampling 61 taxa using culture-based and single-cell transcriptomics, our analyses show two major clades of Amoebozoa, Discosea, and Tevosa. This phylogeny refutes previous studies in major respects. Our results support the hypothesis that the last common ancestor of Amoebozoa was sexual and flagellated, it also may have had the ability to disperse propagules from a sporocarp-type fruiting body. Overall, the main macroevolutionary patterns in Amoebozoa appear to result from the parallel losses of homologous characters of a multiphase life cycle that included flagella, sex, and sporocarps rather than independent acquisition of convergent features.

The Evolutionary Landscape of Dbl-Like RhoGEF Families: Adapting Eukaryotic Cells to Environmental Signals

The dynamics of cell morphology in eukaryotes is largely controlled by small GTPases of the Rho family. Rho GTPases are activated by guanine nucleotide exchange factors (RhoGEFs), of which diffuse B-cell lymphoma (Dbl)-like members form the largest family. Here, we surveyed Dbl-like sequences from 175 eukaryotic genomes and illuminate how the Dbl family evolved in all eukaryotic supergroups. By combining probabilistic phylogenetic approaches and functional domain analysis, we show that the human Dbl-like family is made of 71 members, structured into 20 subfamilies. The 71 members were already present in ancestral jawed vertebrates, but several members were subsequently lost in specific clades, up to 12% in birds. The jawed vertebrate repertoire was established from two rounds of duplications that occurred between tunicates, cyclostomes, and jawed vertebrates. Duplicated members showed distinct tissue distributions, conserved at least in Amniotes. All 20 subfamilies have members in Deuterostomes and Protostomes. Nineteen subfamilies are present in Porifera, the first phylum that diverged in Metazoa, 14 in Choanoflagellida and Filasterea, single-celled organisms closely related to Metazoa and three in Fungi, the sister clade to Metazoa. Other eukaryotic supergroups show an extraordinary variability of Dbl-like repertoires as a result of repeated and independent gain and loss events. Last, we observed that in Metazoa, the number of Dbl-like RhoGEFs varies in proportion of cell signaling complexity. Overall, our analysis supports the conclusion that Dbl-like RhoGEFs were present at the origin of eukaryotes and evolved as highly adaptive cell signaling mediators.

Microbial Eukaryotes that Lack Sterols

It is widely held that sterols are key cyclic triterpenoid lipids in eukaryotic cell membranes and are synthesized through oxygen-dependent multienzyme pathways. However, there are known exceptions-ciliated protozoans, such as Tetrahymena, along with diverse low-oxygen-adapted eukaryotes produce, instead of sterols, the cyclic triterpenoid lipid tetrahymanol that does not require molecular oxygen for its biosynthesis. Here, we report that a number of anaerobic microbial eukaryotes (protists) utilize neither sterols nor tetrahymanol in their membranes. The lack of detectable sterol-like compounds in their membranes may provide an opportunity to reconsider the physiological function of sterols and sterol-like lipids in eukaryotes.

The origin of Metazoa: a unicellular perspective

The first animals evolved from an unknown single-celled ancestor in the Precambrian period. Recently, the identification and characterization of the genomic and cellular traits of the protists most closely related to animals have shed light on the origin of animals. Comparisons of animals with these unicellular relatives allow us to reconstruct the first evolutionary steps towards animal multicellularity. Here, we review the results of these investigations and discuss their implications for understanding the earliest stages of animal evolution, including the origin of metazoan genes and genome function.

Phylogenetic and ecological diversity of apusomonads, a lineage of deep-branching eukaryotes

Apusomonads are a mysterious group of heterotrophic gliding biflagellates branching deeply in the eukaryotic tree of life as sister group to opisthokonts (including animals, fungi, and a variety of unicellular protists). Despite their evolutionary interest, their diversity and ecology remain largely unknown, with very few described species and environmental sequences in databases. Most environmental 18S rRNA gene-based studies generally fail to identify apusomonad sequences, which might be due to primer bias, low abundance, and/or to the fact that their biotopes remain poorly explored. We have carried out an extensive search of 18S rRNA genes using an apusomonad-specific primer in a wide variety of ecosystems. Our study significantly broadens the diversity of apusomonads showing that, despite being mostly rare protists, they often dwell in freshwater and marine benthic environments, generally associated with low-oxygen concentrations. Apusomonads have been identified in environments across a wide salinity range. Some operational taxonomic units (OTUs), occurring in both marine and freshwater ecosystems, seem truly euryhaline, indicating that members of this deep-branching lineage easily cross such ecological barriers.

Expansion of the molecular and morphological diversity of Acanthamoebidae (Centramoebida, Amoebozoa) and identification of a novel life cycle type within the group

BACKGROUND

Acanthamoebidae is a "family" level amoebozoan group composed of the genera Acanthamoeba, Protacanthamoeba, and very recently Luapeleamoeba. This clade of amoebozoans has received considerable attention from the broader scientific community as Acanthamoeba spp. represent both model organisms and human pathogens. While the classical composition of the group (Acanthamoeba + Protacanthamoeba) has been well accepted due to the morphological and ultrastructural similarities of its members, the Acanthamoebidae has never been highly statistically supported in single gene phylogenetic reconstructions of Amoebozoa either by maximum likelihood (ML) or Bayesian analyses.

RESULTS

Here we show using a phylogenomic approach that the Acanthamoebidae is a fully supported monophyletic group within Amoebozoa with both ML and Bayesian analyses. We also expand the known range of morphological and life cycle diversity found in the Acanthamoebidae by demonstrating that the amoebozoans "Protostelium" arachisporum, Dracoamoeba jormungandri n. g. n. sp., and Vacuolamoeba acanthoformis n.g. n.sp., belong within the group. We also found that "Protostelium" pyriformis is clearly a species of Acanthamoeba making it the first reported sporocarpic member of the genus, that is, an amoeba that individually forms a walled, dormant propagule elevated by a non-cellular stalk. Our phylogenetic analyses recover a fully supported Acanthamoebidae composed of five genera. Two of these genera (Acanthamoeba and Luapeleameoba) have members that are sporocarpic.

CONCLUSIONS

Our results provide high statistical support for an Acanthamoebidae that is composed of five distinct genera. This study increases the known morphological diversity of this group and shows that species of Acanthamoeba can include spore-bearing stages. This further illustrates the widespread nature of spore-bearing stages across the tree of Amoebozoa.

REVIEWERS

This article was reviewed by Drs. Eugene Koonin, Purificacion Lopez-Garcia and Sandra Baldauf. Sandra Baldauf was nominated by Purificacion Lopez-Garcia, an Editorial Board member.

Evolution of the Tetrapyrrole Biosynthetic Pathway in Secondary Algae: Conservation, Redundancy and Replacement

Tetrapyrroles such as chlorophyll and heme are indispensable for life because they are involved in energy fixation and consumption, i.e. photosynthesis and oxidative phosphorylation. In eukaryotes, the tetrapyrrole biosynthetic pathway is shaped by past endosymbioses. We investigated the origins and predicted locations of the enzymes of the heme pathway in the chlorarachniophyte Bigelowiella natans, the cryptophyte Guillardia theta, the “green” dinoflagellate Lepidodinium chlorophorum, and three dinoflagellates with diatom endosymbionts (“dinotoms”): Durinskia baltica, Glenodinium foliaceum and Kryptoperidinium foliaceum. Bigelowiella natans appears to contain two separate heme pathways analogous to those found in Euglena gracilis; one is predicted to be mitochondrial-cytosolic, while the second is predicted to be plastid-located. In the remaining algae, only plastid-type tetrapyrrole synthesis is present, with a single remnant of the mitochondrial-cytosolic pathway, a ferrochelatase of G. theta putatively located in the mitochondrion. The green dinoflagellate contains a single pathway composed of mostly rhodophyte-origin enzymes, and the dinotoms hold two heme pathways of apparently plastidal origin. We suggest that heme pathway enzymes in B. natans and L. chlorophorum share a predominantly rhodophytic origin. This implies the ancient presence of a rhodophyte-derived plastid in the chlorarachniophyte alga, analogous to the green dinoflagellate, or an exceptionally massive horizontal gene transfer.

Evolution of bacterial recombinase A (recA) in eukaryotes explained by addition of genomic data of key microbial lineages

Recombinase enzymes promote DNA repair by homologous recombination. The genes that encode them are ancestral to life, occurring in all known dominions: viruses, Eubacteria, Archaea and Eukaryota. Bacterial recombinases are also present in viruses and eukaryotic groups (supergroups), presumably via ancestral events of lateral gene transfer. The eukaryotic recA genes have two distinct origins (mitochondrial and plastidial), whose acquisition by eukaryotes was possible via primary (bacteria-eukaryote) and/or secondary (eukaryote-eukaryote) endosymbiotic gene transfers (EGTs). Here we present a comprehensive phylogenetic analysis of the recA genealogy, with substantially increased taxonomic sampling in the bacteria, viruses, eukaryotes and a special focus on the key eukaryotic supergroup Amoebozoa, earlier represented only by Dictyostelium We demonstrate that several major eukaryotic lineages have lost the bacterial recombinases (including Opisthokonta and Excavata), whereas others have retained them (Amoebozoa, Archaeplastida and the SAR-supergroups). When absent, the bacterial recA homologues may have been lost entirely (secondary loss of canonical mitochondria) or replaced by other eukaryotic recombinases. RecA proteins have a transit peptide for organellar import, where they act. The reconstruction of the RecA phylogeny with its EGT events presented here retells the intertwined evolutionary history of eukaryotes and bacteria, while further illuminating the events of endosymbiosis in eukaryotes by expanding the collection of widespread genes that provide insight to this deep history.

Oligotrophic lagoons of the South Pacific Ocean are home to a surprising number of novel eukaryotic microorganisms

The diversity of microbial eukaryotes was surveyed by environmental sequencing from tropical lagoon sites of the South Pacific, collected through the American Museum of Natural History (AMNH)'s Explore21 expedition to the Solomon Islands in September 2013. The sampled lagoons presented low nutrient concentrations typical of oligotrophic waters, but contained levels of chlorophyll a, a proxy for phytoplankton biomass, characteristic of meso- to eutrophic waters. Two 18S rDNA hypervariable sites, the V4 and V8-V9 regions, were amplified from the total of eight lagoon samples and sequenced on the MiSeq system. After assembly, clustering at 97% similarity, and removal of singletons and chimeras, a total of 2741 (V4) and 2606 (V8-V9) operational taxonomic units (OTUs) were identified. Taxonomic annotation of these reads, including phylogeny, was based on a combination of automated pipeline and manual inspection. About 18.4% (V4) and 13.8% (V8-V9) of the OTUs could not be assigned to any of the known eukaryotic groups. Of these, we focused on OTUs that were not divergent and possessed multiple sources of evidence for their existence. Phylogenetic analyses of these sequences revealed more than ten branches that might represent new deeply-branching lineages of microbial eukaryotes, currently without any cultured representatives or morphological information.

Evolutionary cell biology: functional insight from "endless forms most beautiful"

In animal and fungal model organisms, the complexities of cell biology have been analyzed in exquisite detail and much is known about how these organisms function at the cellular level. However, the model organisms cell biologists generally use include only a tiny fraction of the true diversity of eukaryotic cellular forms. The divergent cellular processes observed in these more distant lineages are still largely unknown in the general scientific community. Despite the relative obscurity of these organisms, comparative studies of them across eukaryotic diversity have had profound implications for our understanding of fundamental cell biology in all species and have revealed the evolution and origins of previously observed cellular processes. In this Perspective, we will discuss the complexity of cell biology found across the eukaryotic tree, and three specific examples of where studies of divergent cell biology have altered our understanding of key functional aspects of mitochondria, plastids, and membrane trafficking.

Environmental Breviatea harbour mutualistic Arcobacter epibionts

Breviatea form a lineage of free living, unicellular protists, distantly related to animals and fungi. This lineage emerged almost one billion years ago, when the oceanic oxygen content was low, and extant Breviatea have evolved or retained an anaerobic lifestyle. Here we report the cultivation of Lenisia limosa, gen. et sp. nov., a newly discovered breviate colonized by relatives of animal-associated Arcobacter. Physiological experiments show that the association of L. limosa with Arcobacter is driven by the transfer of hydrogen and is mutualistic, providing benefits to both partners. With whole-genome sequencing and differential proteomics, we show that an experimentally observed fitness gain of L. limosa could be explained by the activity of a so far unknown type of NAD(P)H-accepting hydrogenase, which is expressed in the presence, but not in the absence, of Arcobacter. Differential proteomics further reveal that the presence of Lenisia stimulates expression of known 'virulence' factors by Arcobacter. These proteins typically enable colonization of animal cells during infection, but may in the present case act for mutual benefit. Finally, re-investigation of two currently available transcriptomic data sets of other Breviatea reveals the presence and activity of related hydrogen-consuming Arcobacter, indicating that mutualistic interaction between these two groups of microbes might be pervasive. Our results support the notion that molecular mechanisms involved in virulence can also support mutualism, as shown here for Arcobacter and Breviatea.

Novel Hydrogenosomes in the Microaerophilic Jakobid Stygiella incarcerata

Mitochondrion-related organelles (MROs) have arisen independently in a wide range of anaerobic protist lineages. Only a few of these organelles and their functions have been investigated in detail, and most of what is known about MROs comes from studies of parasitic organisms such as the parabasalid Trichomonas vaginalis. Here, we describe the MRO of a free-living anaerobic jakobid excavate, Stygiella incarcerata. We report an RNAseq-based reconstruction of S. incarcerata’s MRO proteome, with an associated biochemical map of the pathways predicted to be present in this organelle. The pyruvate metabolism and oxidative stress response pathways are strikingly similar to those found in the MROs of other anaerobic protists, such as Pygsuia and Trichomonas. This elegant example of convergent evolution is suggestive of an anaerobic biochemical ‘module’ of prokaryotic origins that has been laterally transferred among eukaryotes, enabling them to adapt rapidly to anaerobiosis. We also identified genes corresponding to a variety of mitochondrial processes not found in Trichomonas, including intermembrane space components of the mitochondrial protein import apparatus, and enzymes involved in amino acid metabolism and cardiolipin biosynthesis. In this respect, the MROs of S. incarcerata more closely resemble those of the much more distantly related free-living organisms Pygsuia biforma and Cantina marsupialis, likely reflecting these organisms’ shared lifestyle as free-living anaerobes.