Phylogenetic reconciliation: making the most of genomes to understand microbial ecology and evolution.

In recent years, phylogenetic reconciliation has emerged as a promising approach for studying microbial ecology and evolution. The core idea is to model how gene trees evolve along a species tree and to explain differences between them via evolutionary events including gene duplications, transfers, and losses. Here, we describe how phylogenetic reconciliation provides a natural framework for studying genome evolution and highlight recent applications including ancestral gene content inference, the rooting of species trees, and the insights into metabolic evolution and ecological transitions they yield. Reconciliation analyses have elucidated the evolution of diverse microbial lineages, from Chlamydiae to Asgard archaea, shedding light on ecological adaptation, host-microbe interactions, and symbiotic relationships. However, there are many opportunities for broader application of the approach in microbiology. Continuing improvements to make reconciliation models more realistic and scalable, and integration of ecological metadata such as habitat, pH, temperature, and oxygen use offer enormous potential for understanding the rich tapestry of microbial life.

The parasitic lifestyle of an archaeal symbiont.

DPANN archaea are a diverse group of microorganisms characterised by small cells and reduced genomes. To date, all cultivated DPANN archaea are ectosymbionts that require direct cell contact with an archaeal host species for growth and survival. However, these interactions and their impact on the host species are poorly understood. Here, we show that a DPANN archaeon (Candidatus Nanohaloarchaeum antarcticus) engages in parasitic interactions with its host (Halorubrum lacusprofundi) that result in host cell lysis. During these interactions, the nanohaloarchaeon appears to enter, or be engulfed by, the host cell. Our results provide experimental evidence for a predatory-like lifestyle of an archaeon, suggesting that at least some DPANN archaea may have roles in controlling host populations and their ecology.

The nature of the last universal common ancestor and its impact on the early Earth system.

The nature of the last universal common ancestor (LUCA), its age and its impact on the Earth system have been the subject of vigorous debate across diverse disciplines, often based on disparate data and methods. Age estimates for LUCA are usually based on the fossil record, varying with every reinterpretation. The nature of LUCA's metabolism has proven equally contentious, with some attributing all core metabolisms to LUCA, whereas others reconstruct a simpler life form dependent on geochemistry. Here we infer that LUCA lived ~4.2 Ga (4.09-4.33 Ga) through divergence time analysis of pre-LUCA gene duplicates, calibrated using microbial fossils and isotope records under a new cross-bracing implementation. Phylogenetic reconciliation suggests that LUCA had a genome of at least 2.5 Mb (2.49-2.99 Mb), encoding around 2,600 proteins, comparable to modern prokaryotes. Our results suggest LUCA was a prokaryote-grade anaerobic acetogen that possessed an early immune system. Although LUCA is sometimes perceived as living in isolation, we infer LUCA to have been part of an established ecological system. The metabolism of LUCA would have provided a niche for other microbial community members and hydrogen recycling by atmospheric photochemistry could have supported a modestly productive early ecosystem.

Minimal and hybrid hydrogenases are active from archaea.

Microbial hydrogen (H2) cycling underpins the diversity and functionality of diverse anoxic ecosystems. Among the three evolutionarily distinct hydrogenase superfamilies responsible, [FeFe] hydrogenases were thought to be restricted to bacteria and eukaryotes. Here, we show that anaerobic archaea encode diverse, active, and ancient lineages of [FeFe] hydrogenases through combining analysis of existing and new genomes with extensive biochemical experiments. [FeFe] hydrogenases are encoded by genomes of nine archaeal phyla and expressed by H2-producing Asgard archaeon cultures. We report an ultraminimal hydrogenase in DPANN archaea that binds the catalytic H-cluster and produces H2. Moreover, we identify and characterize remarkable hybrid complexes formed through the fusion of [FeFe] and [NiFe] hydrogenases in ten other archaeal orders. Phylogenetic analysis and structural modeling suggest a deep evolutionary history of hybrid hydrogenases. These findings reveal new metabolic adaptations of archaea, streamlined H2 catalysts for biotechnological development, and a surprisingly intertwined evolutionary history between the two major H2-metabolizing enzymes.

Selective lipid recruitment by an archaeal DPANN symbiont from its host.

The symbiont Ca. Nanohaloarchaeum antarcticus is obligately dependent on its host Halorubrum lacusprofundi for lipids and other metabolites due to its lack of certain biosynthetic genes. However, it remains unclear which specific lipids or metabolites are acquired from its host, and how the host responds to infection. Here, we explored the lipidome dynamics of the Ca. Nha. antarcticus - Hrr. lacusprofundi symbiotic relationship during co-cultivation. By using a comprehensive untargeted lipidomic methodology, our study reveals that Ca. Nha. antarcticus selectively recruits 110 lipid species from its host, i.e., nearly two-thirds of the total number of host lipids. Lipid profiles of co-cultures displayed shifts in abundances of bacterioruberins and menaquinones and changes in degree of bilayer-forming glycerolipid unsaturation. This likely results in increased membrane fluidity and improved resistance to membrane disruptions, consistent with compensation for higher metabolic load and mechanical stress on host membranes when in contact with Ca. Nha. antarcticus cells. Notably, our findings differ from previous observations of other DPANN symbiont-host systems, where no differences in lipidome composition were reported. Altogether, our work emphasizes the strength of employing untargeted lipidomics approaches to provide details into the dynamics underlying a DPANN symbiont-host system.

ATP synthase evolution on a cross-braced dated tree of life.

The timing of early cellular evolution, from the divergence of Archaea and Bacteria to the origin of eukaryotes, is poorly constrained. The ATP synthase complex is thought to have originated prior to the Last Universal Common Ancestor (LUCA) and analyses of ATP synthase genes, together with ribosomes, have played a key role in inferring and rooting the tree of life. We reconstruct the evolutionary history of ATP synthases using an expanded taxon sampling set and develop a phylogenetic cross-bracing approach, constraining equivalent speciation nodes to be contemporaneous, based on the phylogenetic imprint of endosymbioses and ancient gene duplications. This approach results in a highly resolved, dated species tree and establishes an absolute timeline for ATP synthase evolution. Our analyses show that the divergence of ATP synthase into F- and A/V-type lineages was a very early event in cellular evolution dating back to more than 4 Ga, potentially predating the diversification of Archaea and Bacteria. Our cross-braced, dated tree of life also provides insight into more recent evolutionary transitions including eukaryogenesis, showing that the eukaryotic nuclear and mitochondrial lineages diverged from their closest archaeal (2.67-2.19 Ga) and bacterial (2.58-2.12 Ga) relatives at approximately the same time, with a slightly longer nuclear stem-lineage.

On the origin of the nucleus: a hypothesis.

SUMMARYIn this hypothesis article, we explore the origin of the eukaryotic nucleus. In doing so, we first look afresh at the nature of this defining feature of the eukaryotic cell and its core functions-emphasizing the utility of seeing the eukaryotic nucleoplasm and cytoplasm as distinct regions of a common compartment. We then discuss recent progress in understanding the evolution of the eukaryotic cell from archaeal and bacterial ancestors, focusing on phylogenetic and experimental data which have revealed that many eukaryotic machines with nuclear activities have archaeal counterparts. In addition, we review the literature describing the cell biology of representatives of the TACK and Asgardarchaeaota - the closest known living archaeal relatives of eukaryotes. Finally, bringing these strands together, we propose a model for the archaeal origin of the nucleus that explains much of the current data, including predictions that can be used to put the model to the test.

Defining eukaryotes to dissect eukaryogenesis.

The origin of eukaryotes is among the most contentious debates in evolutionary biology, attracting multiple seemingly incompatible theories seeking to explain the sequence in which eukaryotic characteristics were acquired. Much of the controversy arises from differing views on the defining characteristics of eukaryotes. We argue that eukaryotes should be defined phylogenetically, and that doing so clarifies where competing hypotheses of eukaryogenesis agree and how we may test among aspects of disagreement. Some hypotheses make predictions about the phylogenetic origins of eukaryotic genes and are distinguishable on that basis. However, other hypotheses differ only in the order of key evolutionary steps, like mitochondrial endosymbiosis and nuclear assembly, which cannot currently be distinguished phylogenetically. Stages within eukaryogenesis may be made identifiable through the absolute dating of gene duplicates that map to eukaryotic traits, such as in genes of host or mitochondrial origin that duplicated and diverged functionally prior to emergence of the last eukaryotic common ancestor. In this way, it may finally be possible to distinguish heat from light in the debate over eukaryogenesis.

Parameter Estimation and Species Tree Rooting Using ALE and GeneRax.

ALE and GeneRax are tools for probabilistic gene tree-species tree reconciliation. Based on a common underlying statistical model of how gene trees evolve along species trees, these methods rely on gene vs. species tree discordance to infer gene duplication, transfer, and loss events, map gene family origins, and root species trees. Published analyses have used these methods to root species trees of Archaea, Bacteria, and several eukaryotic groups, as well as to infer ancestral gene repertoires. However, it was recently suggested that reconciliation-based estimates of duplication and transfer events using the ALE/GeneRax model were unreliable, with potential implications for species tree rooting. Here, we assess these criticisms and find that the methods are accurate when applied to simulated data and in generally good agreement with alternative methodological approaches on empirical data. In particular, ALE recovers variation in gene duplication and transfer frequencies across lineages that is consistent with the known biology of studied clades. In plants and opisthokonts, ALE recovers the consensus species tree root; in Bacteria-where there is less certainty about the root position-ALE agrees with alternative approaches on the most likely root region. Overall, ALE and related approaches are promising tools for studying genome evolution.

On distinguishing between canonical tRNA genes and tRNA gene fragments in prokaryotes.

Automated genome annotation is essential for extracting biological information from sequence data. The identification and annotation of tRNA genes is frequently performed by the software package tRNAscan-SE, the output of which is listed for selected genomes in the Genomic tRNA database (GtRNAdb). Here, we highlight a pervasive error in prokaryotic tRNA gene sets on GtRNAdb: the mis-categorization of partial, non-canonical tRNA genes as standard, canonical tRNA genes. Firstly, we demonstrate the issue using the tRNA gene sets of 20 organisms from the archaeal taxon Thermococcaceae. According to GtRNAdb, these organisms collectively deviate from the expected set of tRNA genes in 15 instances, including the listing of eleven putative canonical tRNA genes. However, after detailed manual annotation, only one of these eleven remains; the others are either partial, non-canonical tRNA genes resulting from the integration of genetic elements or CRISPR-Cas activity (seven instances), or attributable to ambiguities in input sequences (three instances). Secondly, we show that similar examples of the mis-categorization of predicted tRNA sequences occur throughout the prokaryotic sections of GtRNAdb. While both canonical and non-canonical prokaryotic tRNA gene sequences identified by tRNAscan-SE are biologically interesting, the challenge of reliably distinguishing between them remains. We recommend employing a combination of (i) screening input sequences for the genetic elements typically associated with non-canonical tRNA genes, and ambiguities, (ii) activating the tRNAscan-SE automated pseudogene detection function, and (iii) scrutinizing predicted tRNA genes with low isotype scores. These measures greatly reduce manual annotation efforts, and lead to improved prokaryotic tRNA gene set predictions.

Is an archaeon the ancestor of eukaryotes?

The origin of complex cellular life is a key puzzle in evolutionary research, which has broad implications for various neighbouring scientific disciplines. Naturally, views on this topic vary widely depending on the world view and context from which this topic is approached. In the following, I will share my perspective about our current scientific knowledge on the origin of eukaryotic cells, that is, eukaryogenesis, from a biological point of view focusing on the question as to whether an archaeon was the ancestor of eukaryotes.

The importance of biofilm formation for cultivation of a Micrarchaeon and its interactions with its Thermoplasmatales host.

Micrarchaeota is a distinctive lineage assigned to the DPANN archaea, which includes poorly characterised microorganisms with reduced genomes that likely depend on interactions with hosts for growth and survival. Here, we report the enrichment of a stable co-culture of a member of the Micrarchaeota (Ca. Micrarchaeum harzensis) together with its Thermoplasmatales host (Ca. Scheffleriplasma hospitalis), as well as the isolation of the latter. We show that symbiont-host interactions depend on biofilm formation as evidenced by growth experiments, comparative transcriptomic analyses and electron microscopy. In addition, genomic, metabolomic, extracellular polymeric substances and lipid content analyses indicate that the Micrarchaeon symbiont relies on the acquisition of metabolites from its host. Our study of the cell biology and physiology of a Micrarchaeon and its host adds to our limited knowledge of archaeal symbioses.

Evolving Perspective on the Origin and Diversification of Cellular Life and the Virosphere.

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

Origin of eukaryotes: What can be learned from the first successfully isolated Asgard archaeon.

The origin of cellular complexity characterizing eukaryotic cells remains a central unresolved issue in the study of diversification of cellular life on Earth. The isolation by Imachi et al.1 of a member of the Asgard archaea2 - a contemporary relative of organisms thought to have given rise to eukaryotic cells about 2 billion years ago - now promises new insight. The complete genome sequence of the isolated Lokiarchaeum strain confirms that the eukaryotic signature proteins (ESPs) previously identified in the Lokiarchaeota3 and other Asgard archaea2 are indeed encoded by these archaeal genomes and do not represent contamination from eukaryotes. These ESPs encode homologs of eukaryotic actins, small GTPases and the ESCRT complex proteins and are required for the functioning of complex eukaryotic cells. The new, slowly growing, anaerobic laboratory strain allows a first direct look at these organisms and provides key insights into the morphology and metabolism of an Asgard archaeal organism. The work has provided valuable information for other laboratories that aim to isolate and characterize related organisms from other environments.

An estimate of the deepest branches of the tree of life from ancient vertically evolving genes.

Core gene phylogenies provide a window into early evolution, but different gene sets and analytical methods have yielded substantially different views of the tree of life. Trees inferred from a small set of universal core genes have typically supported a long branch separating the archaeal and bacterial domains. By contrast, recent analyses of a broader set of non-ribosomal genes have suggested that Archaea may be less divergent from Bacteria, and that estimates of inter-domain distance are inflated due to accelerated evolution of ribosomal proteins along the inter-domain branch. Resolving this debate is key to determining the diversity of the archaeal and bacterial domains, the shape of the tree of life, and our understanding of the early course of cellular evolution. Here, we investigate the evolutionary history of the marker genes key to the debate. We show that estimates of a reduced Archaea-Bacteria (AB) branch length result from inter-domain gene transfers and hidden paralogy in the expanded marker gene set. By contrast, analysis of a broad range of manually curated marker gene datasets from an evenly sampled set of 700 Archaea and Bacteria reveals that current methods likely underestimate the AB branch length due to substitutional saturation and poor model fit; that the best-performing phylogenetic markers tend to support longer inter-domain branch lengths; and that the AB branch lengths of ribosomal and non-ribosomal marker genes are statistically indistinguishable. Furthermore, our phylogeny inferred from the 27 highest-ranked marker genes recovers a clade of DPANN at the base of the Archaea and places the bacterial Candidate Phyla Radiation (CPR) within Bacteria as the sister group to the Chloroflexota.

A rooted phylogeny resolves early bacterial evolution.

A rooted bacterial tree is necessary to understand early evolution, but the position of the root is contested. Here, we model the evolution of 11,272 gene families to identify the root, extent of horizontal gene transfer (HGT), and the nature of the last bacterial common ancestor (LBCA). Our analyses root the tree between the major clades Terrabacteria and Gracilicutes and suggest that LBCA was a free-living flagellated, rod-shaped double-membraned organism. Contrary to recent proposals, our analyses reject a basal placement of the Candidate Phyla Radiation, which instead branches sister to Chloroflexota within Terrabacteria. While most gene families (92%) have evidence of HGT, overall, two-thirds of gene transmissions have been vertical, suggesting that a rooted tree provides a meaningful frame of reference for interpreting bacterial evolution.

Complex subsurface hydrothermal fluid mixing at a submarine arc volcano supports distinct and highly diverse microbial communities.

Hydrothermally active submarine volcanoes are mineral-rich biological oases contributing significantly to chemical fluxes in the deep sea, yet little is known about the microbial communities inhabiting these systems. Here we investigate the diversity of microbial life in hydrothermal deposits and their metagenomics-inferred physiology in light of the geological history and resulting hydrothermal fluid paths in the subsurface of Brothers submarine volcano north of New Zealand on the southern Kermadec arc. From metagenome-assembled genomes we identified over 90 putative bacterial and archaeal genomic families and nearly 300 previously unknown genera, many potentially endemic to this submarine volcanic environment. While magmatically influenced hydrothermal systems on the volcanic resurgent cones of Brothers volcano harbor communities of thermoacidophiles and diverse members of the superphylum "DPANN," two distinct communities are associated with the caldera wall, likely shaped by two different types of hydrothermal circulation. The communities whose phylogenetic diversity primarily aligns with that of the cone sites and magmatically influenced hydrothermal systems elsewhere are characterized predominately by anaerobic metabolisms. These populations are probably maintained by fluids with greater magmatic inputs that have interacted with different (deeper) previously altered mineral assemblages. However, proximal (a few meters distant) communities with gene-inferred aerobic, microaerophilic, and anaerobic metabolisms are likely supported by shallower seawater-dominated circulation. Furthermore, mixing of fluids from these two distinct hydrothermal circulation systems may have an underlying imprint on the high microbial phylogenomic diversity. Collectively our results highlight the importance of considering geologic evolution and history of subsurface processes in studying microbial colonization and community dynamics in volcanic environments.

Hikarchaeia demonstrate an intermediate stage in the methanogen-to-halophile transition.

Halobacteria (henceforth: Haloarchaea) are predominantly aerobic halophiles that are thought to have evolved from anaerobic methanogens. This remarkable transformation most likely involved an extensive influx of bacterial genes. Whether it entailed a single massive transfer event or a gradual stream of transfers remains a matter of debate. To address this, genomes that descend from methanogen-to-halophile intermediates are necessary. Here, we present five such near-complete genomes of Marine Group IV archaea (Hikarchaeia), the closest known relatives of Haloarchaea. Their inclusion in gene tree-aware ancestral reconstructions reveals an intermediate stage that had already lost a large number of genes, including nearly all of those involved in methanogenesis and the Wood-Ljungdahl pathway. In contrast, the last Haloarchaea common ancestor gained a large number of genes and expanded its aerobic respiration and salt/UV resistance gene repertoire. Our results suggest that complex and gradual patterns of gain and loss shaped the methanogen-to-halophile transition.

Chlamydial contribution to anaerobic metabolism during eukaryotic evolution.

The origin of eukaryotes is a major open question in evolutionary biology. Multiple hypotheses posit that eukaryotes likely evolved from a syntrophic relationship between an archaeon and an alphaproteobacterium based on H2 exchange. However, there are no strong indications that modern eukaryotic H2 metabolism originated from archaea or alphaproteobacteria. Here, we present evidence for the origin of H2 metabolism genes in eukaryotes from an ancestor of the Anoxychlamydiales-a group of anaerobic chlamydiae, newly described here, from marine sediments. Among Chlamydiae, these bacteria uniquely encode genes for H2 metabolism and other anaerobiosis-associated pathways. Phylogenetic analyses of several components of H2 metabolism reveal that Anoxychlamydiales homologs are the closest relatives to eukaryotic sequences. We propose that an ancestor of the Anoxychlamydiales contributed these key genes during the evolution of eukaryotes, supporting a mosaic evolutionary origin of eukaryotic metabolism.