Eukaryotic Tree Of Life Research Articles

On the possibility of yet a third kinetochore system in the protist phylum Euglenozoa.

A macromolecular assembly called the kinetochore is essential for the segregation of genetic material during eukaryotic cell division. Therefore, characterization of kinetochores across species is essential for understanding the mechanisms involved in this key process across the eukaryotic tree of life. In particular, little is known about kinetochores in divergent protists such as Euglenozoa, a group of unicellular flagellates that includes kinetoplastids, euglenids, and diplonemids, the latter being a highly diverse and abundant component of marine plankton. While kinetoplastids have an unconventional kinetochore system and euglenids have a canonical one similar to traditional model eukaryotes, preliminary searches detected neither unconventional nor canonical kinetochore components in diplonemids. Here, we employed state-of-the-art deep homology detection protocols but still could not detect orthologs for the bulk of kinetoplastid-specific nor canonical kinetochore proteins in diplonemids except for a putative centromere-specific histone H3 variant. Our results suggest that diplonemids evolved kinetochores that do not resemble previously known ones.

Mlf mediates proteotoxic response via formation of cellular foci for protein folding and degradation in Giardia.

Myeloid leukemia factor 1 (Mlf1) was identified as a proto-oncoprotein that affects hematopoietic differentiation in humans. However, its cellular function remains elusive, spanning roles from cell cycle regulation to modulation of protein aggregate formation and participation in ciliogenesis. Given that structurally conserved homologs of Mlf1 can be found across the eukaryotic tree of life, we decided to characterize its cellular role underlying this phenotypic pleiotropy. Using a model of the unicellular eukaryote Giardia intestinalis, we demonstrate that its Mlf1 homolog (GiMlf) mainly localizes to two types of cytosolic foci: microtubular structures, where it interacts with Hsp40, and ubiquitin-rich, membraneless compartments, found adjacent to mitochondrion-related organelles known as mitosomes, containing the 26S proteasome regulatory subunit 4. Upon cellular stress, GiMlf either relocates to the affected compartment or disperses across the cytoplasm, subsequently accumulating into enlarged foci during the recovery phase. In vitro assays suggest that GiMlf can be recruited to membranes through its affinity for signaling phospholipids. Importantly, cytosolic foci diminish in the gimlf knockout strain, which exhibits extensive proteomic changes indicative of compromised proteostasis. Consistent with data from other cellular systems, we propose that Mlf acts in the response to proteotoxic stress by mediating the formation of function-specific foci for protein folding and degradation.

Evolutionary origins of the lysosome-related organelle sorting machinery reveal ancient homology in post-endosome trafficking pathways

The major organelles of the endomembrane system were in place by the time of the last eukaryotic common ancestor (LECA) (~1.5 billion years ago). Their acquisitions were defining milestones during eukaryogenesis. Comparative cell biology and evolutionary analyses show multiple instances of homology in the protein machinery controlling distinct interorganelle trafficking routes. Resolving these homologous relationships allows us to explore processes underlying the emergence of additional, distinct cellular compartments, infer ancestral states predating LECA, and explore the process of eukaryogenesis itself. Here, we undertake a molecular evolutionary analysis (including providing a transcriptome of the jakobid flagellate Reclinomonas americana), exploring the origins of the machinery responsible for the biogenesis of lysosome-related organelles (LROs), the Biogenesis of LRO Complexes (BLOCs 1,2, and 3). This pathway has been studied only in animals and is not considered a feature of the basic eukaryotic cell plan. We show that this machinery is present across the eukaryotic tree of life and was likely in place prior to LECA, making it an underappreciated facet of eukaryotic cellular organisation. Moreover, we resolve multiple points of ancient homology between all three BLOCs and other post-endosomal retrograde trafficking machinery (BORC, CCZ1 and MON1 proteins, and an unexpected relationship with the "homotypic fusion and vacuole protein sorting" (HOPS) and "Class C core vacuole/endosomal tethering" (CORVET) complexes), offering a mechanistic and evolutionary unification of these trafficking pathways. Overall, this study provides a comprehensive account of the rise of the LROs biogenesis machinery from before the LECA to current eukaryotic diversity, integrating it into the larger mechanistic framework describing endomembrane evolution.

Lessons from Extremophiles: Functional Adaptations and Genomic Innovations across the Eukaryotic Tree of Life.

From hydrothermal vents, to glaciers, to deserts, research in extreme environments has reshaped our understanding of how and where life can persist. Contained within the genomes of extremophilic organisms are the blueprints for a toolkit to tackle the multitude of challenges of survival in inhospitable environments. As new sequencing technologies have rapidly developed, so too has our understanding of the molecular and genomic mechanisms that have facilitated the success of extremophiles. Although eukaryotic extremophiles remain relatively understudied compared to bacteria and archaea, an increasing number of studies have begun to leverage 'omics tools to shed light on eukaryotic life in harsh conditions. In this perspective paper, we highlight a diverse breadth of research on extremophilic lineages across the eukaryotic tree of life, from microbes to macrobes, that are collectively reshaping our understanding of molecular innovations at life's extremes. These studies are not only advancing our understanding of evolution and biological processes but are also offering a valuable roadmap on how emerging technologies can be applied to identify cellular mechanisms of adaptation to cope with life in stressful conditions, including high and low temperatures, limited water availability, and heavy metal habitats. We shed light on patterns of molecular and organismal adaptation across the eukaryotic tree of life and discuss a few promising research directions, including investigations into the role of horizontal gene transfer in eukaryotic extremophiles and the importance of increasing phylogenetic diversity of model systems.

Somatic genome architecture and molecular evolution are decoupled in "young" linage-specific gene families in ciliates.

The evolution of lineage-specific gene families remains poorly studied across the eukaryotic tree of life, with most analyses focusing on the recent evolution of de novo genes in model species. Here we explore the origins of lineage-specific genes in ciliates, a ~1 billion year old clade of microeukaryotes that are defined by their division of somatic and germline functions into distinct nuclei. Previous analyses on conserved gene families have shown the effect of ciliates' unusual genome architecture on gene family evolution: extensive genome processing-the generation of thousands of gene-sized somatic chromosomes from canonical germline chromosomes-is associated with larger and more diverse gene families. To further study the relationship between ciliate genome architecture and gene family evolution, we analyzed lineage specific gene families from a set of 46 transcriptomes and 12 genomes representing x species from eight ciliate classes. We assess how the evolution lineage-specific gene families occurs among four groups of ciliates: extensive fragmenters with gene-size somatic chromosomes, non-extensive fragmenters with "large'' multi-gene somatic chromosomes, Heterotrichea with highly polyploid somatic genomes and Karyorelictea with 'paradiploid' somatic genomes. Our analyses demonstrate that: 1) most lineage-specific gene families are found at shallow taxonomic scales; 2) extensive genome processing (i.e., gene unscrambling) during development likely influences the size and number of young lineage-specific gene families; and 3) the influence of somatic genome architecture on molecular evolution is increasingly apparent in older gene families. Altogether, these data highlight the influences of genome architecture on the evolution of lineage-specific gene families in eukaryotes.

Voltage-Gated Proton Channels in the Tree of Life.

With a single gene encoding HV1 channel, proton channel diversity is particularly low in mammals compared to other members of the superfamily of voltage-gated ion channels. Nonetheless, mammalian HV1 channels are expressed in many different tissues and cell types where they exert various functions. In the first part of this review, we regard novel aspects of the functional expression of HV1 channels in mammals by differentially comparing their involvement in (1) close conjunction with the NADPH oxidase complex responsible for the respiratory burst of phagocytes, and (2) in respiratory burst independent functions such as pH homeostasis or acid extrusion. In the second part, we dissect expression of HV channels within the eukaryotic tree of life, revealing the immense diversity of the channel in other phylae, such as mollusks or dinoflagellates, where several genes encoding HV channels can be found within a single species. In the last part, a comprehensive overview of the biophysical properties of a set of twenty different HV channels characterized electrophysiologically, from Mammalia to unicellular protists, is given.

Joint inference of evolutionary transitions to self-fertilization and demographic history using whole-genome sequences.

The evolution from outcrossing to selfing is a transition that occurred recurrently throughout the eukaryote tree of life, in plants, animals, fungi and algae. Despite some short-term advantages, selfing is supposed to be an evolutionary dead-end reproductive strategy on the long-term and its tippy distribution on phylogenies suggests that most selfing species are of recent origin. However, dating such transitions is challenging while it is central for this hypothesis. We build on previous theories to explicit the differential effect of past changes in selfing rate or in population size on the probability of recombination events along the genome. This allows us to develop two methods making use of full genome polymorphism data to 1) test if a transition from outcrossing to selfing occurred, and 2) infer its age. The sequentially Markov coalescent based (teSMC) and the Approximate Bayesian Computation (tsABC) methods use a common framework based on a transition matrix summarizing the distribution of times to the most recent common ancestor along the genome, allowing to estimate changes in the ratio of population recombination and mutation rates in time. We first demonstrate that our methods can disentangle between past change in selfing rate from past changes in demographic history. Second, we assess the accuracy of our methods and show that transitions to selfing as old as approximatively 2.5Ne generations can be identified from polymorphism data. Third, our estimates are robust to the presence of linked negative selection on coding sequences. Finally, as a proof of principle, we apply both methods to three populations from Arabidopsis thaliana, recovering a transition to selfing which occurred approximately 600,000 years ago. Our methods pave the way to study recent transitions to predominant self-fertilization in selfing organisms and to better account for variation in mating systems in demographic inferences.

On the evolution of variation in sexual reproduction through the prism of eukaryotic microbes

Almost all eukaryotes undergo sexual reproduction to generate diversity and select for fitness in their population pools. Interestingly, the systems by which sex is defined are highly diverse and can even differ between evolutionarily closely related species. While the most commonly known form of sex determination involves males and females in animals, eukaryotic microbes can have as many as thousands of different mating types for the same species. Furthermore, some species have found alternatives to sexual reproduction and prefer to grow clonally and yet undergo infrequent facultative sexual reproduction. These organisms are mainly invertebrates and microbes, but several examples are also present among vertebrates suggesting that alternative modes of sexual reproduction evolved multiple times throughout evolution. In this review, we summarize the sex-determination modes and variants of sexual reproduction found across the eukaryotic tree of life and suggest that eukaryotic microbes provide unique opportunities to study these processes in detail. We propose that understanding variations in modes of sexual reproduction can serve as a foundation to study the evolution of sex and why and how it evolved in the first place.

A transcriptional switch controls sex determination in Plasmodium falciparum

Sexual reproduction and meiotic sex are deeply rooted in the eukaryotic tree of life, but mechanisms determining sex or mating types are extremely varied and are only well characterized in a few model organisms1. In malaria parasites, sexual reproduction coincides with transmission to the vector host. Sex determination is non-genetic, with each haploid parasite capable of producing either a male or a female gametocyte in the human host2. The hierarchy of events and molecular mechanisms that trigger sex determination and maintenance of sexual identity are yet to be elucidated. Here we show that the male development 1 (md1) gene is both necessary and sufficient for male fate determination in the human malaria parasite Plasmodium falciparum. We show that Md1 has a dual function stemming from two separate domains: in sex determination through its N terminus and in male development from its conserved C-terminal LOTUS/OST-HTH domain. We further identify a bistable switch at the md1 locus, which is coupled with sex determination and ensures that the male-determining gene is not expressed in the female lineage. We describe one of only a few known non-genetic mechanisms of sex determination in a eukaryote and highlight Md1 as a potential target for interventions that block malaria transmission.

Transposable elements drive intron gain in diverse eukaryotes

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Euglena International Network (EIN): Driving euglenoid biotechnology for the benefit of a challenged world.

Euglenoids (Euglenida) are unicellular flagellates possessing exceptionally wide geographical and ecological distribution. Euglenoids combine a biotechnological potential with a unique position in the eukaryotic tree of life. In large part these microbes owe this success to diverse genetics including secondary endosymbiosis and likely additional sources of genes. Multiple euglenoid species have translational applications and show great promise in production of biofuels, nutraceuticals, bioremediation, cancer treatments and more exotically as robotics design simulators. An absence of reference genomes currently limits these applications, including development of efficient tools for identification of critical factors in regulation, growth or optimization of metabolic pathways. The Euglena International Network (EIN) seeks to provide a forum to overcome these challenges. EIN has agreed specific goals, mobilized scientists, established a clear roadmap (Grand Challenges), connected academic and industry stakeholders and is currently formulating policy and partnership principles to propel these efforts in a coordinated and efficient manner.

Molecular diet analysis in mussels and other metazoan filter feeders and an assessment of their utility as natural eDNA samplers.

Molecular gut content analysis is a popular tool to study food web interactions and has recently been suggested as an alternative source for DNA-based biomonitoring. However, the overabundant consumer's DNA often outcompetes that of its diet during PCR. Lineage-specific primers are an efficient means to reduce consumer amplification while retaining broad specificity for dietary taxa. Here, we designed an amplicon sequencing assay to monitor the eukaryotic diet of mussels and other metazoan filter feeders and explore the utility of mussels as natural eDNA samplers to monitor planktonic communities. We designed several lineage-specific rDNA primers with broad taxonomic suitability for eukaryotes. The primers were tested using DNA extracts of different limnic and marine mussel species and the results compared to eDNA water samples collected next to the mussel colonies. In addition, we analysed several 25-year time series samples of mussels from German rivers. Our primer sets efficiently prevent the amplification of mussels and other metazoans. The recovered DNA reflects a broad dietary preference across the eukaryotic tree of life and considerable taxonomic overlap with filtered water samples. We also show the utility of a reversed version of our primers, which prevents amplification of nonmetazoan taxa from complex eukaryote community samples, by enriching fauna associated with the marine brown algae Fucus vesiculosus. Our protocol will enable large-scale dietary analysis in metazoan filter feeders, facilitate aquatic food web analysis and allow surveying of aquacultures for pathogens. Moreover, we show that mussels and other aquatic filter feeders can serve as complementary DNA source for biomonitoring.

Global patterns and rates of habitat transitions across the eukaryotic tree of life

The successful colonization of new habitats has played a fundamental role during the evolution of life. Salinity is one of the strongest barriers for organisms to cross, which has resulted in the evolution of distinct marine and non-marine (including both freshwater and soil) communities. Although microbes represent by far the vast majority of eukaryote diversity, the role of the salt barrier in shaping the diversity across the eukaryotic tree is poorly known. Traditional views suggest rare and ancient marine/non-marine transitions but this view is being challenged by the discovery of several recently transitioned lineages. Here, we investigate habitat evolution across the tree of eukaryotes using a unique set of taxon-rich phylogenies inferred from a combination of long-read and short-read environmental metabarcoding data spanning the ribosomal DNA operon. Our results show that, overall, marine and non-marine microbial communities are phylogenetically distinct but transitions have occurred in both directions in almost all major eukaryotic lineages, with hundreds of transition events detected. Some groups have experienced relatively high rates of transitions, most notably fungi for which crossing the salt barrier has probably been an important aspect of their successful diversification. At the deepest phylogenetic levels, ancestral habitat reconstruction analyses suggest that eukaryotes may have first evolved in non-marine habitats and that the two largest known eukaryotic assemblages (TSAR and Amorphea) arose in different habitats. Overall, our findings indicate that the salt barrier has played an important role during eukaryote evolution and provide a global perspective on habitat transitions in this domain of life.

Conversion of anterograde into retrograde trains is an intrinsic property of intraflagellar transport

SummaryCilia or eukaryotic flagella are microtubule-based organelles found across the eukaryotic tree of life. Their very high aspect ratio and crowded interior are unfavorable to diffusive transport of most components required for their assembly and maintenance. Instead, a system of intraflagellar transport (IFT) trains moves cargo rapidly up and down the cilium (Figure 1A).1, 2, 3 Anterograde IFT, from the cell body to the ciliary tip, is driven by kinesin-II motors, whereas retrograde IFT is powered by cytoplasmic dynein-1b motors.4 Both motors are associated with long chains of IFT protein complexes, known as IFT trains, and their cargoes.5, 6, 7, 8 The conversion from anterograde to retrograde motility at the ciliary tip involves (1) the dissociation of kinesin motors from trains,9 (2) a fundamental restructuring of the train from the anterograde to the retrograde architecture,8,10,11 (3) the unloading and reloading of cargo,2 and (4) the activation of the dynein motors.8,12 A prominent hypothesis is that there is dedicated calcium-dependent protein-based machinery at the ciliary tip to mediate these processes.4,13 However, the mechanisms of IFT turnaround have remained elusive. In this study, we use mechanical and chemical methods to block IFT at intermediate positions along the cilia of the green algae Chlamydomonas reinhardtii, in normal and calcium-depleted conditions. We show that IFT turnaround, kinesin dissociation, and dynein-1b activation can consistently be induced at arbitrary distances from the ciliary tip, with no stationary tip machinery being required. Instead, we demonstrate that the anterograde-to-retrograde conversion is a calcium-independent intrinsic ability of IFT.

The advantage of sex: Reinserting fluctuating selection in the pluralist approach.

The advantage of sex, and its fixation in some clades and species all over the eukaryote tree of life, is considered an evolutionary enigma, especially regarding its assumed two-fold cost. Several likely hypotheses have been proposed such as (1) a better response to the negative frequency-dependent selection imposed by the “Red Queen” hypothesis; (2) the competition between siblings induced by the Tangled Bank hypothesis; (3) the existence of genetic and of (4) ecological factors that can diminish the cost of sex to less than the standard assumed two-fold; and (5) a better maintenance of genetic diversity and its resulting phenotypic variation, providing a selective advantage in randomly fluctuating environments. While these hypotheses have mostly been studied separately, they can also act simultaneously. This was advocated by several studies which presented a pluralist point of view. Only three among the five causes cited above were considered yet in such a framework: the Red Queen hypothesis, the Tangled Bank and the genetic factors lowering the cost of sex. We thus simulated the evolution of a finite mutating population undergoing negative frequency-dependent selection on phenotypes and a two-fold (or less) cost of sexuality, experiencing randomly fluctuating selection along generations. The individuals inherited their reproductive modes, either clonal or sexual. We found that exclusive sexuality begins to fix in populations exposed to environmental variation that exceeds the width of one ecological niche (twice the standard deviation of a Gaussian response to environment). This threshold was lowered by increasing negative frequency-dependent selection and when reducing the two-fold cost of sex. It contributes advocating that the different processes involved in a short-term advantage of sex and recombination can act in combination to favor the fixation of sexual reproduction in populations.

Mixotrophy in the bloom-forming genus Phaeocystis and other haptophytes

Phaeocystis is a globally widespread marine phytoplankton genus, best known for its colony-forming species that can form large blooms and odorous foam during bloom decline. In the North Sea, Phaeocystis globosa typically becomes abundant towards the end of the spring bloom, when nutrients are depleted and the share of mixotrophic protists increases. Although mixotrophy is widespread across the eukaryotic tree of life and is also found amongst haptophytes, a mixotrophic nutrition has not yet been demonstrated in Phaeocystis. Here, we sampled two consecutive Phaeocystis globosa spring blooms in the coastal North Sea. In both years, bacterial cells were observed inside 0.6 – 2% of P. globosa cells using double CARD-FISH hybridizations in combination with laser scanning confocal microscopy. Incubation experiments manipulating light and nutrient availability showed a trend towards higher occurrence of intracellular bacteria under P-deplete conditions. Based on counts of bacteria inside P. globosa cells in combination with theoretical values of prey digestion times, maximum ingestion rates of up to 0.08 bacteria cell−1 h−1 were estimated. In addition, a gene-based predictive model was applied to the transcriptome assemblies of seven Phaeocystis strains and 24 other haptophytes to assess their trophic mode. This model predicted a phago-mixotrophic feeding strategy in several (but not all) strains of P. globosa, P. antarctica and other haptophytes that were previously assumed to be autotrophic. The observation of bacterial cells inside P. globosa and the gene-based model predictions strongly suggest that the phago-mixotrophic feeding strategy is widespread among members of the Phaeocystis genus and other haptophytes, and might contribute to their remarkable success to form nuisance blooms under nutrient-limiting conditions.

Cross‐kingdom transcriptomic trends in the evolution of hybrid gene expression

Hybridization is a route to speciation that occurs widely across the eukaryote tree of life. The success of allopolyploids (hybrid species with increased ploidy) and homoploid hybrids (with unchanged ploidy) is well documented. However, their formation and establishment is not straightforward, with a suite of near‐instantaneous and longer term biological repercussions faced by the new species. Central to these challenges is the rewiring of gene regulatory networks following the merger of distinct genomes inherited from both parental species. Research on the evolution of hybrid gene expression has largely involved studies on a single hybrid species or a few gene families. Here, we present the first standardized transcriptome‐wide study exploring the fates of genes following hybridization across three kingdoms: animals, plants and fungi. Within each kingdom, we pair an allopolyploid system with a closely related homoploid hybrid to decouple the influence of increased ploidy from genome merger. Genome merger, not changes in ploidy, has the greatest effect on posthybridization expression patterns across all study systems. Strikingly, we find that differentially expressed genes in parent species preferentially switch to more similar expression in hybrids across all kingdoms, likely as a consequence of regulatory trans‐acting cross‐talk within the hybrid nucleus. We also highlight the prevalence of gene loss or silencing among extremely differentially expressed genes in hybrid species across all kingdoms. These shared patterns suggest that the evolutionary process of hybridization leads to common high‐level expression outcomes, regardless of the particular species or kingdom.

Convergent evolution of polyploid genomes from across the eukaryotic tree of life.

By modeling the homoeologous gene losses that occurred in 50 genomes deriving from ten distinct polyploidy events, we show that the evolutionary forces acting on polyploids are remarkably similar, regardless of whether they occur in flowering plants, ciliates, fishes, or yeasts. We show that many of the events show a relative rate of duplicate gene loss before the first postpolyploidy speciation that is significantly higher than in later phases of their evolution. The relatively weak selective constraint experienced by the single-copy genes these losses produced leads us to suggest that most of the purely selectively neutral duplicate gene losses occur in the immediate postpolyploid period. Nearly all of the events show strong evidence of biases in the duplicate losses, consistent with them being allopolyploidies, with 2 distinct progenitors contributing to the modern species. We also find ongoing and extensive reciprocal gene losses (alternative losses of duplicated ancestral genes) between these genomes. With the exception of a handful of closely related taxa, all of these polyploid organisms are separated from each other by tens to thousands of reciprocal gene losses. As a result, it is very unlikely that viable diploid hybrid species could form between these taxa, since matings between such hybrids would tend to produce offspring lacking essential genes. It is, therefore, possible that the relatively high frequency of recurrent polyploidies in some lineages may be due to the ability of new polyploidies to bypass reciprocal gene loss barriers.

Conflict over the Eukaryote Root Resides in Strong Outliers, Mosaics and Missing Data Sensitivity of Site-Specific (CAT) Mixture Models.

Phylogenetic reconstruction using concatenated loci ("phylogenomics" or "supermatrix phylogeny") is a powerful tool for solving evolutionary splits that are poorly resolved in single gene/protein trees. However, recent phylogenomic attempts to resolve the eukaryote root have yielded conflicting results, along with claims of various artifacts hidden in the data. We have investigated these conflicts using two new methods for assessing phylogenetic conflict. ConJak uses whole marker (gene or protein) jackknifing to assess deviation from a central mean for each individual sequence, whereas ConWin uses a sliding window to screen for incongruent protein fragments (mosaics). Both methods allow selective masking of individual sequences or sequence fragments in order to minimize missing data, an important consideration for resolving deep splits with limited data. Analyses focused on a set of 76 eukaryotic proteins of bacterial ancestry previously used in various combinations to assess the branching order among the three major divisions of eukaryotes: Amorphea (mainly animals, fungi, and Amoebozoa), Diaphoretickes (most other well-known eukaryotes and nearly all algae) and Excavata, represented here by Discoba (Jakobida, Heterolobosea, and Euglenozoa). ConJak analyses found strong outliers to be concentrated in undersampled lineages, whereas ConWin analyses of Discoba, the most undersampled of the major lineages, detected potentially incongruent fragments scattered throughout. Phylogenetic analyses of the full data using an LG-gamma model support a Discoba sister scenario (neozoan-excavate root), which rises to 99-100% bootstrap support with data masked according to either protocol. However, analyses with two site-specific (CAT) mixture models yielded widely inconsistent results and a striking sensitivity to missing data. The neozoan-excavate root places Amorphea and Diaphoretickes as more closely related to each other than either is to Discoba, a fundamental relationship that should remain unaffected by additional taxa. [CAT-GTR; Discoba; eukaryote tree of life; HGT; jackknife; mixture models; mosaic genes; phylogenomics; sliding window; supermatrix.].

Truly ubiquitous CRESS DNA viruses scattered across the eukaryotic tree of life.

Until recently, most viruses detected and characterized were of economic significance, associated with agricultural and medical diseases. This was certainly true for the eukaryote-infecting circular Rep (replication-associated protein)-encoding single-stranded DNA (CRESS DNA) viruses, which were thought to be a relatively small group of viruses. With the explosion of metagenomic sequencing over the past decade and increasing use of rolling-circle replication for sequence amplification, scientists have identified and annotated copious numbers of novel CRESS DNA viruses - many without known hosts but which have been found in association with eukaryotes. Similar advances in cellular genomics have revealed that many eukaryotes have endogenous sequences homologous to viral Reps, which not only provide 'fossil records' to reconstruct the evolutionary history of CRESS DNA viruses but also reveal potential host species for viruses known by their sequences alone. The Rep protein is a conserved protein that all CRESS DNA viruses use to assist rolling-circle replication that is known to be endogenized in a few eukaryotic species (notably tobacco and water yam). A systematic search for endogenous Rep-like sequences in GenBank's non-redundant eukaryotic database was performed using tBLASTn. We utilized relaxed search criteria for the capture of integrated Rep sequence within eukaryotic genomes, identifying 93 unique species with an endogenized fragment of Rep in their nuclear, plasmid (onespecies), mitochondrial (sixspecies) or chloroplast (eightspecies) genomes. These species come from 19 different phyla, scattered across the eukaryotic tree of life. Exogenous and endogenous CRESS DNA viral Rep tree topology suggested potential hosts for one family of uncharacterized viruses and supports a primarily fungal host range for genomoviruses.