Extending the conserved phylogenetic core of archaea disentangles the evolution of the third domain of life

Deep origin of eukaryotes outside Heimdallarchaeia within Asgardarchaeota

Research on the morphology, physiology and genomics of Asgard archaea has provided valuable insights into the evolutionary history of eukaryotes1-3. A previous study suggested that eukaryotes are nested within Heimdallarchaeia4, but their exact phylogenetic placement within Asgard archaea remains controversial4,5. This debate complicates understanding of the metabolic features and timescales of early eukaryotic ancestors. Here we generated 223 metagenome-assembled nearly complete genomes of Asgard archaea that have not previously been documented. We identify 16 new lineages at the genus level or higher, which substantially expands the known phylogenetic diversity of Asgard archaea. Through sophisticated phylogenomic analysis of this expanded genomic dataset involving several marker sets we infer that eukaryotes evolved before the diversification of all sampled Heimdallarchaeia, rather than branching with Hodarchaeales within the Heimdallarchaeia. This difference in the placement of eukaryotes is probably caused by the previously underappreciated chimeric nature of Njordarchaeales genomes, which we find are composed of sequences of both Asgard and TACK archaea (Asgard's sister phylum). Using ancestral reconstruction and molecular dating, we infer that the last Asgard archaea and eukaryote common ancestor emerged before the Great Oxidation Event and was probably an anaerobic H2-dependent acetogen. Our findings support the hydrogen hypothesis of eukaryogenesis, which posits that eukaryotes arose from the fusion of a H2-consuming archaeal host and a H2-producing protomitochondrion.

Extremely acidic proteomes and metabolic flexibility in bacteria and highly diversified archaea thriving in geothermal chaotropic brines

Few described archaeal, and fewer bacterial, lineages thrive under salt-saturating conditions, such as solar saltern crystallizers (salinity above 30% w/v). They accumulate molar K+ cytoplasmic concentrations to maintain osmotic balance ('salt-in' strategy) and have proteins adaptively enriched in negatively charged acidic amino acids. Here we analysed metagenomes and metagenome-assembled genomes from geothermally influenced hypersaline ecosystems with increasing chaotropicity in the Danakil Depression. Normalized abundances of universal single-copy genes confirmed that haloarchaea and Nanohaloarchaeota encompass 99% of microbial communities in the near-life-limiting conditions of the Western-Canyon Lakes. Danakil metagenome- and metagenome-assembled-genome-inferred proteomes, compared with those of freshwater, seawater and solar saltern ponds up to saturation (6-14-32% salinity), showed that Western-Canyon Lake archaea encode the most acidic proteomes ever observed (median protein isoelectric points ≤4.4). We identified previously undescribed haloarchaeal families as well as an Aenigmatarchaeota family and a bacterial phylum independently adapted to extreme halophily. Despite phylum-level diversity decreasing with increasing salinity-chaotropicity, and unlike in solar salterns, adapted archaea exceedingly diversified in Danakil ecosystems, challenging the notion of decreasing diversity under extreme conditions. Metabolic flexibility to utilize multiple energy and carbon resources generated by local hydrothermalism along feast-and-famine strategies seemingly shapes microbial diversity in these ecosystems near life limits.

Expanded phylogeny of extremely halophilic archaea shows multiple independent adaptations to hypersaline environments

Extremely halophilic archaea (Haloarchaea, Nanohaloarchaeota, Methanonatronarchaeia and Halarchaeoplasmatales) thrive in saturating salt concentrations where they must maintain osmotic equilibrium with their environment. The evolutionary history of adaptations enabling salt tolerance remains poorly understood, in particular because the phylogeny of several lineages is conflicting. Here we present a resolved phylogeny of extremely halophilic archaea obtained using improved taxon sampling and state-of-the-art phylogenetic approaches designed to cope with the strong compositional biases of their proteomes. We describe two uncultured lineages, Afararchaeaceae and Asbonarchaeaceae, which break the long branches at the base of Haloarchaea and Nanohaloarchaeota, respectively. We obtained 13 metagenome-assembled genomes (MAGs) of these archaea from metagenomes of hypersaline aquatic systems of the Danakil Depression (Ethiopia). Our phylogenomic analyses including these taxa show that at least four independent adaptations to extreme halophily occurred during archaeal evolution. Gene-tree/species-tree reconciliation suggests that gene duplication and horizontal gene transfer played an important role in this process, for example, by spreading key genes (such as those encoding potassium transporters) across extremely halophilic lineages.

Unraveling the phylogenomic diversity of Methanomassiliicoccales and implications for mitigating ruminant methane emissions

BACKGROUND

Methanomassiliicoccales are a recently identified order of methanogens that are diverse across global environments particularly the gastrointestinal tracts of animals; however, their metabolic capacities are defined via a limited number of cultured strains.

RESULTS

Here, we profile and analyze 243 Methanomassiliicoccales genomes assembled from cultured representatives and uncultured metagenomes recovered from various biomes, including the gastrointestinal tracts of different animal species. Our analyses reveal the presence of numerous undefined genera and genetic variability in metabolic capabilities within Methanomassiliicoccales lineages, which is essential for adaptation to their ecological niches. In particular, gastrointestinal tract Methanomassiliicoccales demonstrate the presence of co-diversified members with their hosts over evolutionary timescales and likely originated in the natural environment. We highlight the presence of diverse clades of vitamin transporter BtuC proteins that distinguish Methanomassiliicoccales from other archaeal orders and likely provide a competitive advantage in efficiently handling B12. Furthermore, genome-centric metatranscriptomic analysis of ruminants with varying methane yields reveal elevated expression of select Methanomassiliicoccales genera in low methane animals and suggest that B12 exchanges could enable them to occupy ecological niches that possibly alter the direction of H2 utilization.

CONCLUSIONS

We provide a comprehensive and updated account of divergent Methanomassiliicoccales lineages, drawing from numerous uncultured genomes obtained from various habitats. We also highlight their unique metabolic capabilities involving B12, which could serve as promising targets for mitigating ruminant methane emissions by altering H2 flow.

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.

Asgard archaeal selenoproteome reveals a roadmap for the archaea-to-eukaryote transition of selenocysteine incorporation machinery

Selenocysteine (Sec) is encoded by the UGA codon that normally functions as a stop signal and is specifically incorporated into selenoproteins via a unique recoding mechanism. The translational recoding of UGA as Sec is directed by an unusual RNA structure, the SECIS element. Although archaea and eukaryotes adopt similar Sec encoding machinery, the SECIS elements have no similarities to each other with regard to sequence and structure. We analyzed >400 Asgard archaeal genomes to examine the occurrence of both Sec encoding system and selenoproteins in this archaeal superphylum, the closest prokaryotic relatives of eukaryotes. A comprehensive map of Sec utilization trait has been generated, providing the most detailed understanding of the use of this nonstandard amino acid in Asgard archaea so far. By characterizing the selenoproteomes of all organisms, several selenoprotein-rich phyla and species were identified. Most Asgard archaeal selenoprotein genes possess eukaryotic SECIS-like structures with varying degrees of diversity. Moreover, euryarchaeal SECIS elements might originate from Asgard archaeal SECIS elements via lateral gene transfer, indicating a complex and dynamic scenario of the evolution of SECIS element within archaea. Finally, a roadmap for the transition of eukaryotic SECIS elements from archaea was proposed, and selenophosphate synthetase may serve as a potential intermediate for the generation of ancestral eukaryotic SECIS element. Our results offer new insights into a deeper understanding of the evolution of Sec insertion machinery.

An Unsupervised Classifier for Whole-Genome Phylogenies, the Maxwell© Tool

The development of phylogenetic trees based on RNA or DNA sequences generally requires a precise and limited choice of important RNAs, e.g., messenger RNAs of essential proteins or ribosomal RNAs (like 16S), but rarely complete genomes, making it possible to explain evolution and speciation. In this article, we propose revisiting a classic phylogeny of archaea from only the information on the succession of nucleotides of their entire genome. For this purpose, we use a new tool, the unsupervised classifier Maxwell, whose principle lies in the Burrows-Wheeler compression transform, and we show its efficiency in clustering whole archaeal genomes.

Metabolic versatility of Caldarchaeales from geothermal features of Hawai'i and Chile as revealed by five metagenome-assembled genomes

Members of the archaeal order Caldarchaeales (previously the phylum Aigarchaeota) are poorly sampled and are represented in public databases by relatively few genomes. Additional representative genomes will help resolve their placement among all known members of Archaea and provide insights into their roles in the environment. In this study, we analyzed 16S rRNA gene amplicons belonging to the Caldarchaeales that are available in public databases, which demonstrated that archaea of the order Caldarchaeales are diverse, widespread, and most abundant in geothermal habitats. We also constructed five metagenome-assembled genomes (MAGs) of Caldarchaeales from two geothermal features to investigate their metabolic potential and phylogenomic position in the domain Archaea. Two of the MAGs were assembled from microbial community DNA extracted from fumarolic lava rocks from Mauna Ulu, Hawai'i, and three were assembled from DNA obtained from hot spring sinters from the El Tatio geothermal field in Chile. MAGs from Hawai'i are high quality bins with completeness >95% and contamination <1%, and one likely belongs to a novel species in a new genus recently discovered at a submarine volcano off New Zealand. MAGs from Chile have lower completeness levels ranging from 27 to 70%. Gene content of the MAGs revealed that these members of Caldarchaeales are likely metabolically versatile and exhibit the potential for both chemoorganotrophic and chemolithotrophic lifestyles. The wide array of metabolic capabilities exhibited by these members of Caldarchaeales might help them thrive under diverse harsh environmental conditions. All the MAGs except one from Chile harbor putative prophage regions encoding several auxiliary metabolic genes (AMGs) that may confer a fitness advantage on their Caldarchaeales hosts by increasing their metabolic potential and make them better adapted to new environmental conditions. Phylogenomic analysis of the five MAGs and over 3,000 representative archaeal genomes showed the order Caldarchaeales forms a monophyletic group that is sister to the clade comprising the orders Geothermarchaeales (previously Candidatus Geothermarchaeota), Conexivisphaerales and Nitrososphaerales (formerly known as Thaumarchaeota), supporting the status of Caldarchaeales members as a clade distinct from the Thaumarchaeota.

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.

The methanogen core and pangenome: conservation and variability across biology's growth temperature extremes

Temperature is a key variable in biological processes. However, a complete understanding of biological temperature adaptation is lacking, in part because of the unique constraints among different evolutionary lineages and physiological groups. Here we compared the genomes of cultivated psychrotolerant and thermotolerant methanogens, which are physiologically related and span growth temperatures from -2.5°C to 122°C. Despite being phylogenetically distributed amongst three phyla in the archaea, the genomic core of cultivated methanogens comprises about one-third of a given genome, while the genome fraction shared by any two organisms decreases with increasing phylogenetic distance between them. Increased methanogenic growth temperature is associated with reduced genome size, and thermotolerant organisms-which are distributed across the archaeal tree-have larger core genome fractions, suggesting that genome size is governed by temperature rather than phylogeny. Thermotolerant methanogens are enriched in metal and other transporters, and psychrotolerant methanogens are enriched in proteins related to structure and motility. Observed amino acid compositional differences between temperature groups include proteome charge, polarity and unfolding entropy. Our results suggest that in the methanogens, shared physiology maintains a large, conserved genomic core even across large phylogenetic distances and biology's temperature extremes.

Genomic remnants of ancestral methanogenesis and hydrogenotrophy in Archaea drive anaerobic carbon cycling

Anaerobic methane metabolism is among the hallmarks of Archaea, originating very early in their evolution. Here, we show that the ancestor of methane metabolizers was an autotrophic CO₂-reducing hydrogenotrophic methanogen that possessed the two main complexes, methyl-CoM reductase (Mcr) and tetrahydromethanopterin-CoM methyltransferase (Mtr), the anaplerotic hydrogenases Eha and Ehb, and a set of other genes collectively called "methanogenesis markers" but could not oxidize alkanes. Overturning recent inferences, we demonstrate that methyl-dependent hydrogenotrophic methanogenesis has emerged multiple times independently, either due to a loss of Mtr while Mcr is inherited vertically or from an ancient lateral acquisition of Mcr. Even if Mcr is lost, Mtr, Eha, Ehb, and the markers can persist, resulting in mixotrophic metabolisms centered around the Wood-Ljungdahl pathway. Through their methanogenesis remnants, Thorarchaeia and two newly reconstructed order-level lineages in Archaeoglobi and Bathyarchaeia act as metabolically versatile players in carbon cycling of anoxic environments across the globe.

Recovery of Lutacidiplasmatales archaeal order genomes suggests convergent evolution in Thermoplasmatota

The Terrestrial Miscellaneous Euryarchaeota Group has been identified in various environments, and the single genome investigated thus far suggests that these archaea are anaerobic sulfite reducers. We assemble 35 new genomes from this group that, based on genome analysis, appear to possess aerobic and facultative anaerobic lifestyles and may oxidise rather than reduce sulfite. We propose naming this order (representing 16 genera) "Lutacidiplasmatales" due to their occurrence in various acidic environments and placement within the phylum Thermoplasmatota. Phylum-level analysis reveals that Thermoplasmatota evolution had been punctuated by several periods of high levels of novel gene family acquisition. Several essential metabolisms, such as aerobic respiration and acid tolerance, were likely acquired independently by divergent lineages through convergent evolution rather than inherited from a common ancestor. Ultimately, this study describes the terrestrially prevalent Lutacidiciplasmatales and highlights convergent evolution as an important driving force in the evolution of archaeal lineages.

Investigating Abiotic and Biotic Mechanisms of Pyrite Reduction

Pyrite (FeS2) has a very low solubility and therefore has historically been considered a sink for iron (Fe) and sulfur (S) and unavailable to biology in the absence of oxygen and oxidative weathering. Anaerobic methanogens were recently shown to reduce FeS2 and assimilate Fe and S reduction products to meet nutrient demands. However, the mechanism of FeS2 mineral reduction and the forms of Fe and S assimilated by methanogens remained unclear. Thermodynamic calculations described herein indicate that H2 at aqueous concentrations as low as 10-10 M favors the reduction of FeS2, with sulfide (HS-) and pyrrhotite (Fe1- x S) as products; abiotic laboratory experiments confirmed the reduction of FeS2 with dissolved H2 concentrations greater than 1.98 × 10-4 M H2. Growth studies of Methanosarcina barkeri provided with FeS2 as the sole source of Fe and S resulted in H2 production but at concentrations too low to drive abiotic FeS2 reduction, based on abiotic laboratory experimental data. A strain of M. barkeri with deletions in all [NiFe]-hydrogenases maintained the ability to reduce FeS2 during growth, providing further evidence that extracellular electron transport (EET) to FeS2 does not involve H2 or [NiFe]-hydrogenases. Physical contact between cells and FeS2 was required for mineral reduction but was not required to obtain Fe and S from dissolution products. The addition of a synthetic electron shuttle, anthraquinone-2,6-disulfonate, allowed for biological reduction of FeS2 when physical contact between cells and FeS2 was prohibited, indicating that exogenous electron shuttles can mediate FeS2 reduction. Transcriptomics experiments revealed upregulation of several cytoplasmic oxidoreductases during growth of M. barkeri on FeS2, which may indicate involvement in provisioning low potential electrons for EET to FeS2. Collectively, the data presented herein indicate that reduction of insoluble FeS2 by M. barkeri occurred via electron transfer from the cell surface to the mineral surface resulting in the generation of soluble HS- and mineral-associated Fe1- x S. Solubilized Fe(II), but not HS-, from mineral-associated Fe1- x S reacts with aqueous HS- yielding aqueous iron sulfur clusters (FeS aq ) that likely serve as the Fe and S source for methanogen growth and activity. FeS aq nucleation and subsequent precipitation on the surface of cells may result in accelerated EET to FeS2, resulting in positive feedback between cell activity and FeS2 reduction.

A divide-and-conquer phylogenomic approach based on character supermatrices resolves early steps in the evolution of the Archaea

Background

The recent rise in cultivation-independent genome sequencing has provided key material to explore uncharted branches of the Tree of Life. This has been particularly spectacular concerning the Archaea, projecting them at the center stage as prominently relevant to understand early stages in evolution and the emergence of fundamental metabolisms as well as the origin of eukaryotes. Yet, resolving deep divergences remains a challenging task due to well-known tree-reconstruction artefacts and biases in extracting robust ancient phylogenetic signal, notably when analyzing data sets including the three Domains of Life. Among the various strategies aimed at mitigating these problems, divide-and-conquer approaches remain poorly explored, and have been primarily based on reconciliation among single gene trees which however notoriously lack ancient phylogenetic signal.

Results

We analyzed sub-sets of full supermatrices covering the whole Tree of Life with specific taxonomic sampling to robustly resolve different parts of the archaeal phylogeny in light of their current diversity. Our results strongly support the existence and early emergence of two main clades, Cluster I and Cluster II, which we name Ouranosarchaea and Gaiarchaea, and we clarify the placement of important novel archaeal lineages within these two clades. However, the monophyly and branching of the fast evolving nanosized DPANN members remains unclear and worth of further study.

Conclusions

We inferred a well resolved rooted phylogeny of the Archaea that includes all recently described phyla of high taxonomic rank. This phylogeny represents a valuable reference to study the evolutionary events associated to the early steps of the diversification of the archaeal domain. Beyond the specifics of archaeal phylogeny, our results demonstrate the power of divide-and-conquer approaches to resolve deep phylogenetic relationships, which should be applied to progressively resolve the entire Tree of Life.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01952-0.

Pathways of Iron and Sulfur Acquisition, Cofactor Assembly, Destination, and Storage in Diverse Archaeal Methanogens and Alkanotrophs

Archaeal methanogens, methanotrophs, and alkanotrophs have a high demand for iron (Fe) and sulfur (S); however, little is known of how they acquire, traffic, deploy, and store these elements. Here, we examined the distribution of homologs of proteins mediating key steps in Fe/S metabolism in model microorganisms, including iron(II) sensing/uptake (FeoAB), sulfide extraction from cysteine (SufS), and the biosynthesis of iron-sulfur [Fe-S] clusters (SufBCDE), siroheme (Pch2 dehydrogenase), protoheme (AhbABCD), cytochrome c (Cyt c) (CcmCF), and iron storage/detoxification (Bfr, FtrA, and IssA), among 326 publicly available, complete or metagenome-assembled genomes of archaeal methanogens/methanotrophs/alkanotrophs. The results indicate several prevalent but nonuniversal features, including FeoB, SufBC, and the biosynthetic apparatus for the basic tetrapyrrole scaffold, as well as its siroheme (and F₄₃₀) derivatives. However, several early-diverging genomes lacked SufS and pathways to synthesize and deploy heme. Genomes encoding complete versus incomplete heme biosynthetic pathways exhibited equivalent prevalences of [Fe-S] cluster binding proteins, suggesting an expansion of catalytic capabilities rather than substitution of heme for [Fe-S] in the former group. Several strains with heme binding proteins lacked heme biosynthesis capabilities, while other strains with siroheme biosynthesis capability lacked homologs of known siroheme binding proteins, indicating heme auxotrophy and unknown siroheme biochemistry, respectively. While ferritin proteins involved in ferric oxide storage were widespread, those involved in storing Fe as thioferrate were unevenly distributed. Collectively, the results suggest that differences in the mechanisms of Fe and S acquisition, deployment, and storage have accompanied the diversification of methanogens/methanotrophs/alkanotrophs, possibly in response to differential availability of these elements as these organisms evolved.

IMPORTANCE Archaeal methanogens, methanotrophs, and alkanotrophs, argued to be among the most ancient forms of life, have a high demand for iron (Fe) and sulfur (S) for cofactor biosynthesis, among other uses. Here, using comparative bioinformatic approaches applied to 326 genomes, we show that major differences in Fe/S acquisition, trafficking, deployment, and storage exist in this group. Variation in these characters was generally congruent with the phylogenetic placement of these genomes, indicating that variation in Fe/S usage and deployment has contributed to the diversification and ecology of these organisms. However, incongruency was observed among the distribution of cofactor biosynthesis pathways and known protein destinations for those cofactors, suggesting auxotrophy or yet-to-be-discovered pathways for cofactor biosynthesis.

The Molecular Determinants of Thermoadaptation: Methanococcales as a Case Study

Previous reports have shown that environmental temperature impacts proteome evolution in Bacteria and Archaea. However, it is unknown whether thermoadaptation mainly occurs via the sequential accumulation of substitutions, massive horizontal gene transfers, or both. Measuring the real contribution of amino acid substitution to thermoadaptation is challenging, because of confounding environmental and genetic factors (e.g., pH, salinity, genomic G + C content) that also affect proteome evolution. Here, using Methanococcales, a major archaeal lineage, as a study model, we show that optimal growth temperature is the major factor affecting variations in amino acid frequencies of proteomes. By combining phylogenomic and ancestral sequence reconstruction approaches, we disclose a sequential substitutional scheme in which lysine plays a central role by fine tuning the pool of arginine, serine, threonine, glutamine, and asparagine, whose frequencies are strongly correlated with optimal growth temperature. Finally, we show that colonization to new thermal niches is not associated with high amounts of horizontal gene transfers. Altogether, although the acquisition of a few key proteins through horizontal gene transfer may have favored thermoadaptation in Methanococcales, our findings support sequential amino acid substitutions as the main factor driving thermoadaptation.

Examining pathways of iron and sulfur acquisition, trafficking, deployment, and storage in mineral-grown methanogen cells

Methanogens have a high demand for iron (Fe) and sulfur (S); however, little is known of how they acquire, deploy, and store these elements and how this, in turn, affects their physiology. Methanogens were recently shown to reduce pyrite (FeS₂) generating aqueous iron-sulfide (FeS_(aq)) clusters that are likely assimilated as a source of Fe and S. Here, we compare the phenotype of Methanococcus voltae when grown with FeS₂ or ferrous iron (Fe(II)) and sulfide (HS^-). FeS₂-grown cells are 33% smaller yet have 193% more Fe than Fe(II)/HS^--grown cells. Whole cell EPR revealed similar distributions of paramagnetic Fe, although FeS₂-grown cells showed a broad spectral feature attributed to intracellular thioferrate-like nanoparticles. Differential proteomic analyses showed similar expression of core methanogenesis enzymes, indicating that Fe and S source does not substantively alter the energy metabolism of cells. However, a homolog of the Fe(II) transporter FeoB and its putative transcriptional regulator DtxR were up-expressed in FeS₂-grown cells, suggesting that cells sense Fe(II) limitation. Two homologs of IssA, a protein putatively involved in coordinating thioferrate nanoparticles, were also up-expressed in FeS₂-grown cells. We interpret these data to indicate that, in FeS₂-grown cells, DtxR cannot sense Fe(II) and therefore cannot down-regulate FeoB. We suggest this is due to the transport of Fe(II) complexed with sulfide (FeS_(aq)) leading to excess Fe that is sequestered by IssA as a thioferrate-like species. This model provides a framework for the design of targeted experiments aimed at further characterizing Fe acquisition and homeostasis in M. voltae and other methanogens. IMPORTANCE FeS₂ is the most abundant sulfide mineral in the Earth's crust and is common in environments inhabited by methanogenic archaea. FeS₂ can be reduced by methanogens, yielding aqueous FeS_(aq) clusters that are thought to be a source of Fe and S. Here, we show that growth of Methanococcus voltae on FeS₂ results in smaller cell size and higher Fe content per cell, with Fe likely stored intracellularly as thioferrate-like nanoparticles. Fe(II) transporters and storage proteins were up-regulated in FeS₂-grown cells. These responses are interpreted to result from cells incorrectly sensing Fe(II) limitation due to assimilation of Fe(II) as FeS_(aq). These findings have implications for our understanding of how Fe/S availability influences methanogen physiology and the biogeochemical cycling of these elements.

A standardized archaeal taxonomy for the Genome Taxonomy Database

The accrual of genomic data from both cultured and uncultured microorganisms provides new opportunities to develop systematic taxonomies based on evolutionary relationships. Previously, we established a bacterial taxonomy through the Genome Taxonomy Database. Here, we propose a standardized archaeal taxonomy that is derived from a 122-concatenated-protein phylogeny that resolves polyphyletic groups and normalizes ranks based on relative evolutionary divergence. The resulting archaeal taxonomy, which forms part of the Genome Taxonomy Database, is stable for a range of phylogenetic variables including marker gene selection, inference methods, corrections for rate heterogeneity and compositional bias, tree rooting scenarios and expansion of the genome database. Rank normalization is shown to robustly correct for substitution rates varying up to 30-fold using simulated datasets. Taxonomic curation follows the rules of the International Code of Nomenclature of Prokaryotes while taking into account proposals to formally recognize the rank of phylum and to use genome sequences as type material. This taxonomy is based on 2,392 archaeal genomes, 93.3% of which required one or more changes to their existing taxonomy, mainly owing to incomplete classification. We identify 16 archaeal phyla and reclassify 3 major monophyletic units from the former Euryarchaeota and one phylum that unites the Thaumarchaeota-Aigarchaeota-Crenarchaeota-Korarchaeota (TACK) superphylum into a single phylum.

Programmed Deviations of Ribosomes From Standard Decoding in Archaea

Genetic code decoding, initially considered to be universal and immutable, is now known to be flexible. In fact, in specific genes, ribosomes deviate from the standard translational rules in a programmed way, a phenomenon globally termed recoding. Translational recoding, which has been found in all domains of life, includes a group of events occurring during gene translation, namely stop codon readthrough, programmed ± 1 frameshifting, and ribosome bypassing. These events regulate protein expression at translational level and their mechanisms are well known and characterized in viruses, bacteria and eukaryotes. In this review we summarize the current state-of-the-art of recoding in the third domain of life. In Archaea, it was demonstrated and extensively studied that translational recoding regulates the decoding of the 21st and the 22nd amino acids selenocysteine and pyrrolysine, respectively, and only one case of programmed -1 frameshifting has been reported so far in Saccharolobus solfataricus P2. However, further putative events of translational recoding have been hypothesized in other archaeal species, but not extensively studied and confirmed yet. Although this phenomenon could have some implication for the physiology and adaptation of life in extreme environments, this field is still underexplored and genes whose expression could be regulated by recoding are still poorly characterized. The study of these recoding episodes in Archaea is urgently needed.

Six-state Amino Acid Recoding is not an Effective Strategy to Offset Compositional Heterogeneity and Saturation in Phylogenetic Analyses

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of six-state recoding approaches by comparing the performance of analyses on recoded and non-recoded data sets that have been simulated under gradients of compositional heterogeneity or saturation. In our simulation analyses, non-recoding approaches consistently outperform six-state recoding approaches. Our results suggest that six-state recoding strategies are not effective in the face of high saturation. Furthermore, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies six-state recoding outweighs its benefits. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these consistently outperform six-state recoding. Our analyses of other recoding schemes suggest that under conditions of very high compositional heterogeneity, it may be advantageous to apply recoding using more than six states, but we caution that applying any recoding should include sufficient justification. Our results have important implications for the more than 90 published papers that have incorporated six-state recoding, many of which have significant bearing on relationships across the tree of life. [Compositional heterogeneity; Dayhoff 6-state recoding; S&R 6-state recoding; six-state amino acid recoding; substitution saturation.]

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.

Phylogenomic analysis of novel Diaforarchaea is consistent with sulfite but not sulfate reduction in volcanic environments on early Earth

The origin(s) of dissimilatory sulfate and/or (bi)sulfite reducing organisms (SRO) remains enigmatic despite their importance in global carbon and sulfur cycling since at least 3.4 Ga. Here, we describe novel, deep-branching archaeal SRO populations distantly related to other Diaforarchaea from two moderately acidic thermal springs. Dissimilatory (bi)sulfite reductase homologs, DsrABC, encoded in metagenome assembled genomes (MAGs) from spring sediments comprise one of the earliest evolving Dsr lineages. DsrA homologs were expressed in situ under moderately acidic conditions. MAGs lacked genes encoding proteins that activate sulfate prior to (bi)sulfite reduction. This is consistent with sulfide production in enrichment cultures provided sulfite but not sulfate. We suggest input of volcanic sulfur dioxide to anoxic spring-water yields (bi)sulfite and moderately acidic conditions that favor its stability and bioavailability. The presence of similar volcanic springs at the time SRO are thought to have originated (>3.4 Ga) may have supplied (bi)sulfite that supported ancestral SRO. These observations coincide with the lack of inferred SO42- reduction capacity in nearly all organisms with early-branching DsrAB and which are near universally found in hydrothermal environments.

Different Ways of Doing the Same: Variations in the Two Last Steps of the Purine Biosynthetic Pathway in Prokaryotes

The last two steps of the purine biosynthetic pathway may be catalyzed by different enzymes in prokaryotes. The genes that encode these enzymes include homologs of purH, purP, purO and those encoding the AICARFT and IMPCH domains of PurH, here named purV and purJ, respectively. In Bacteria, these reactions are mainly catalyzed by the domains AICARFT and IMPCH of PurH. In Archaea, these reactions may be carried out by PurH and also by PurP and PurO, both considered signatures of this domain and analogous to the AICARFT and IMPCH domains of PurH, respectively. These genes were searched for in 1,403 completely sequenced prokaryotic genomes publicly available. Our analyses revealed taxonomic patterns for the distribution of these genes and anticorrelations in their occurrence. The analyses of bacterial genomes revealed the existence of genes coding for PurV, PurJ, and PurO, which may no longer be considered signatures of the domain Archaea. Although highly divergent, the PurOs of Archaea and Bacteria show a high level of conservation in the amino acids of the active sites of the protein, allowing us to infer that these enzymes are analogs. Based on the results, we propose that the gene purO was present in the common ancestor of all living beings, whereas the gene encoding PurP emerged after the divergence of Archaea and Bacteria and their isoforms originated in duplication events in the common ancestor of phyla Crenarchaeota and Euryarchaeota. The results reported here expand our understanding of the diversity and evolution of the last two steps of the purine biosynthetic pathway in prokaryotes.

Taxonomic assignment of uncultured prokaryotes with long range PCR targeting the spectinomycin operon

The taxonomic assignment of uncultured prokaryotes to known taxa is a major challenge in microbial systematics. This relies usually on the phylogenetic analysis of the ribosomal small subunit RNA or a few housekeeping genes. Recent works have disclosed ribosomal proteins as valuable markers for systematics and, due to the boom in complete genome sequencing, their use has become widespread. Yet, in the case of uncultured strains, for which complete genome sequences cannot be easily obtained, sequencing many markers is complicated and time consuming. Taking the advantage of the organization of ribosomal protein coding genes in large gene clusters, we amplified a 32 kb conserved region encompassing the spectinomycin (spc) operon using long range PCR from isolated and from uncultured nodular endophytic Frankia strains. The phylogenetic analysis of the 27 ribosomal protein genes contained in this region provided a robust phylogenetic tree consistent with phylogenies based on larger set of markers, indicating that this subset of ribosomal proteins contains enough phylogenetic signal to address systematic issues. This work shows that using long range PCR could break down the barrier preventing the use of ribosomal proteins as phylogenetic markers when complete genome sequences cannot be easily obtained.

Expanding anaerobic alkane metabolism in the domain of Archaea

Methanogenesis and anaerobic methane oxidation through methyl-coenzyme M reductase (MCR) as a key enzyme have been suggested to be basal pathways of archaea¹. How widespread MCR-based alkane metabolism is among archaea, where it occurs and how it evolved remain elusive. Here, we performed a global survey of MCR-encoding genomes based on metagenomic data from various environments. Eleven high-quality mcr-containing metagenomic-assembled genomes were obtained belonging to the Archaeoglobi in the Euryarchaeota, Hadesarchaeota and different TACK superphylum archaea, including the Nezhaarchaeota, Korarchaeota and Verstraetearchaeota. Archaeoglobi WYZ-LMO1 and WYZ-LMO3 and Korarchaeota WYZ-LMO9 encode both the (reverse) methanogenesis and the dissimilatory sulfate reduction pathway, suggesting that they have the genomic potential to couple both pathways in individual organisms. The Hadesarchaeota WYZ-LMO4-6 and Archaeoglobi JdFR-42 encode highly divergent MCRs, enzymes that may enable them to thrive on non-methane alkanes. The occurrence of mcr genes in different archaeal phyla indicates that MCR-based alkane metabolism is common in the domain of Archaea.

Pyrrolysine in archaea: a 22nd amino acid encoded through a genetic code expansion

The 22nd amino acid discovered to be directly encoded, pyrrolysine, is specified by UAG. Until recently, pyrrolysine was only known to be present in archaea from a methanogenic lineage (Methanosarcinales), where it is important in enzymes catalysing anoxic methylamines metabolism, and a few anaerobic bacteria. Relatively new discoveries have revealed wider presence in archaea, deepened functional understanding, shown remarkable carbon source-dependent expression of expanded decoding and extended exploitation of the pyrrolysine machinery for synthetic code expansion. At the same time, other studies have shown the presence of pyrrolysine-containing archaea in the human gut and this has prompted health considerations. The article reviews our knowledge of this fascinating exception to the 'standard' genetic code.

Extreme Deviations from Expected Evolutionary Rates in Archaeal Protein Families

Origin of new biological functions is a complex phenomenon ranging from single-nucleotide substitutions to the gain of new genes via horizontal gene transfer or duplication. Neofunctionalization and subfunctionalization of proteins is often attributed to the emergence of paralogs that are subject to relaxed purifying selection or positive selection and thus evolve at accelerated rates. Such phenomena potentially could be detected as anomalies in the phylogenies of the respective gene families. We developed a computational pipeline to search for such anomalies in 1,834 orthologous clusters of archaeal genes, focusing on lineage-specific subfamilies that significantly deviate from the expected rate of evolution. Multiple potential cases of neofunctionalization and subfunctionalization were identified, including some ancient, house-keeping gene families, such as ribosomal protein S10, general transcription factor TFIIB and chaperone Hsp20. As expected, many cases of apparent acceleration of evolution are associated with lineage-specific gene duplication. On other occasions, long branches in phylogenetic trees correspond to horizontal gene transfer across long evolutionary distances. Significant deceleration of evolution is less common than acceleration, and the underlying causes are not well understood; functional shifts accompanied by increased constraints could be involved. Many gene families appear to be “highly evolvable,” that is, include both long and short branches. Even in the absence of precise functional predictions, this approach allows one to select targets for experimentation in search of new biology.

Neutral Theory and the Somatic Evolution of Cancer

Kimura's neutral theory argued that positive selection was not responsible for an appreciable fraction of molecular substitutions. Correspondingly, quantitative analysis reveals that the vast majority of substitutions in cancer genomes are not detectably under selection. Insights from the somatic evolution of cancer reveal that beneficial substitutions in cancer constitute a small but important fraction of the molecular variants. The molecular evolution of cancer community will benefit by incorporating the neutral theory of molecular evolution into their understanding and analysis of cancer evolution-and accepting the use of tractable, predictive models, even when there is some evidence that they are not perfect.

Metabolic shift at the class level sheds light on adaptation of methanogens to oxidative environments

Methanogens have long been considered strictly anaerobic and oxygen-sensitive microorganisms, but their ability to survive oxygen stress has also been documented. Indeed, methanogens have been found in oxidative environments, and antioxidant genes have been detected in their genomes. How methanogens adapt to oxidative environments, however, remain poorly understood. Here, we systematically predicted and annotated antioxidant features from representative genomes across six well-established methanogen orders. Based on functional gene content involved in production of reactive oxygen species, Hierarchical Clustering analyses grouped methanogens into two distinct clusters, corresponding to the Class I and II methanogens, respectively. Comparative genomics suggested a systematic shift in metabolisms across the two classes, resulting in an enrichment of antioxidant features in the Class II. Moreover, meta-analysis of 16 S rRNA gene sequences obtained from EnvDB indicated that members of Class II were more frequently recovered from microaerophilic and even oxic environments than the Class I members. Phylogenomic analysis suggested that the Class I and II methanogens might have evolved before and around the Great Oxygenation Event, respectively. The enrichment of antioxidant features in the Class II methanogens may have played a key role in the adaption of this group to oxidative environments today and historically.

The growing tree of Archaea: new perspectives on their diversity, evolution and ecology

The Archaea occupy a key position in the Tree of Life, and are a major fraction of microbial diversity. Abundant in soils, ocean sediments and the water column, they have crucial roles in processes mediating global carbon and nutrient fluxes. Moreover, they represent an important component of the human microbiome, where their role in health and disease is still unclear. The development of culture-independent sequencing techniques has provided unprecedented access to genomic data from a large number of so far inaccessible archaeal lineages. This is revolutionizing our view of the diversity and metabolic potential of the Archaea in a wide variety of environments, an important step toward understanding their ecological role. The archaeal tree is being rapidly filled up with new branches constituting phyla, classes and orders, generating novel challenges for high-rank systematics, and providing key information for dissecting the origin of this domain, the evolutionary trajectories that have shaped its current diversity, and its relationships with Bacteria and Eukarya. The present picture is that of a huge diversity of the Archaea, which we are only starting to explore.

Reconstructed ancestral enzymes reveal that negative selection drove the evolution of substrate specificity in ADP-dependent kinases

One central goal in molecular evolution is to pinpoint the mechanisms and evolutionary forces that cause an enzyme to change its substrate specificity; however, these processes remain largely unexplored. Using the glycolytic ADP-dependent kinases of archaea, including the orders Thermococcales, Methanosarcinales, and Methanococcales, as a model and employing an approach involving paleoenzymology, evolutionary statistics, and protein structural analysis, we could track changes in substrate specificity during ADP-dependent kinase evolution along with the structural determinants of these changes. To do so, we studied five key resurrected ancestral enzymes as well as their extant counterparts. We found that a major shift in function from a bifunctional ancestor that could phosphorylate either glucose or fructose 6-phosphate (fructose-6-P) as a substrate to a fructose 6-P-specific enzyme was started by a single amino acid substitution resulting in negative selection with a ground-state mode against glucose and a subsequent 1,600-fold change in specificity of the ancestral protein. This change rendered the residual phosphorylation of glucose a promiscuous and physiologically irrelevant activity, highlighting how promiscuity may be an evolutionary vestige of ancestral enzyme activities, which have been eliminated over time. We also could reconstruct the evolutionary history of substrate utilization by using an evolutionary model of discrete binary characters, indicating that substrate uses can be discretely lost or acquired during enzyme evolution. These findings exemplify how negative selection and subtle enzyme changes can lead to major evolutionary shifts in function, which can subsequently generate important adaptive advantages, for example, in improving glycolytic efficiency in Thermococcales.

A Complex Endomembrane System in the Archaeon Ignicoccus hospitalis Tapped by Nanoarchaeum equitans

Based on serial sectioning, focused ion beam scanning electron microscopy (FIB/SEM), and electron tomography, we depict in detail the highly unusual anatomy of the marine hyperthermophilic crenarchaeon, Ignicoccus hospitalis. Our data support a complex and dynamic endomembrane system consisting of cytoplasmic protrusions, and with secretory function. Moreover, we reveal that the cytoplasm of the putative archaeal ectoparasite Nanoarchaeum equitans can get in direct contact with this endomembrane system, complementing and explaining recent proteomic, transcriptomic and metabolomic data on this inter-archaeal relationship. In addition, we identified a matrix of filamentous structures and/or tethers in the voluminous inter-membrane compartment (IMC) of I. hospitalis, which might be responsible for membrane dynamics. Overall, this unusual cellular compartmentalization, ultrastructure and dynamics in an archaeon that belongs to the recently proposed TACK superphylum prompts speculation that the eukaryotic endomembrane system might originate from Archaea.

Archaeal histone distribution is associated with archaeal genome base composition

A subpopulation of Archaea possesses histones, which are similar to eukaryotic histones H3 and H4. However, archaeal histones are smaller than H3 and H4, and are not post-translationally modified. In addition, the fundamental unit of archaeal histones might be a dimer. The organization of archaeal nucleosomes, therefore, differs from that of eukaryotic nucleosomes. The base compositions of archaeal genome are much more diversified than those of eukaryote and the archaeal histones have more diversified amino acid sequences, which are reflected by their varied isoelectric points. We hypothesized that the highly diversified archaeal genomic DNA base composition may cause the archaeal histone variation. Phylogenetic analysis revealed that the distribution of archaeal histones is associated with their genomic DNA base composition. This result strongly suggests that archaeal histones have evolved concomitantly with their genomic DNA base composition. Eukaryotic histones are one of the most evolutionarily conserved proteins and would limit the diversification of genomic DNA base composition. In contrast, archaeal histones have diversified and would permit the great diversification of genomic DNA base composition.

Arguments Reinforcing the Three-Domain View of Diversified Cellular Life

The archaeal ancestor scenario (AAS) for the origin of eukaryotes implies the emergence of a new kind of organism from the fusion of ancestral archaeal and bacterial cells. Equipped with this “chimeric” molecular arsenal, the resulting cell would gradually accumulate unique genes and develop the complex molecular machineries and cellular compartments that are hallmarks of modern eukaryotes. In this regard, proteins related to phagocytosis and cell movement should be present in the archaeal ancestor, thus identifying the recently described candidate archaeal phylum “Lokiarchaeota” as resembling a possible candidate ancestor of eukaryotes. Despite its appeal, AAS seems incompatible with the genomic, molecular, and biochemical differences that exist between Archaea and Eukarya. In particular, the distribution of conserved protein domain structures in the proteomes of cellular organisms and viruses appears hard to reconcile with the AAS. In addition, concerns related to taxon and character sampling, presupposing bacterial outgroups in phylogenies, and nonuniform effects of protein domain structure rearrangement and gain/loss in concatenated alignments of protein sequences cast further doubt on AAS-supporting phylogenies. Here, we evaluate AAS against the traditional “three-domain” world of cellular organisms and propose that the discovery of Lokiarchaeota could be better reconciled under the latter view, especially in light of several additional biological and technical considerations.

Structure and function of an ancestral-type β-decarboxylating dehydrogenase from Thermococcus kodakarensis

β-Decarboxylating dehydrogenases, which are involved in central metabolism, are considered to have diverged from a common ancestor with broad substrate specificity. In a molecular phylogenetic analysis of 183 β-decarboxylating dehydrogenase homologs from 84 species, TK0280 from Thermococcus kodakarensis was selected as a candidate for an ancestral-type β-decarboxylating dehydrogenase. The biochemical characterization of recombinant TK0280 revealed that the enzyme exhibited dehydrogenase activities toward homoisocitrate, isocitrate, and 3-isopropylmalate, which correspond to key reactions involved in the lysine biosynthetic pathway, tricarboxylic acid cycle, and leucine biosynthetic pathway, respectively. In T. kodakarensis, the growth characteristics of the KUW1 host strain and a TK0280 deletion strain suggested that TK0280 is involved in lysine biosynthesis in this archaeon. On the other hand, gene complementation analyses using Thermus thermophilus as a host revealed that TK0280 functions as both an isocitrate dehydrogenase and homoisocitrate dehydrogenase in this organism, but not as a 3-isopropylmalate dehydrogenase, most probably reflecting its low catalytic efficiency toward 3-isopropylmalate. A crystallographic study on TK0280 binding each substrate indicated that Thr71 and Ser80 played important roles in the recognition of homoisocitrate and isocitrate while the hydrophobic region consisting of Ile82 and Leu83 was responsible for the recognition of 3-isopropylmalate. These analyses also suggested the importance of a water-mediated hydrogen bond network for the stabilization of the β3-α4 loop, including the Thr71 residue, with respect to the promiscuity of the substrate specificity of TK0280.

A Comprehensive Study of Cyanobacterial Morphological and Ecological Evolutionary Dynamics through Deep Geologic Time

Cyanobacteria have exerted a profound influence on the progressive oxygenation of Earth. As a complementary approach to examining the geologic record—phylogenomic and trait evolutionary analyses of extant species can lead to new insights. We constructed new phylogenomic trees and analyzed phenotypic trait data using novel phylogenetic comparative methods. We elucidated the dynamics of trait evolution in Cyanobacteria over billion-year timescales, and provide evidence that major geologic events in early Earth’s history have shaped—and been shaped by—evolution in Cyanobacteria. We identify a robust core cyanobacterial phylogeny and a smaller set of taxa that exhibit long-branch attraction artifacts. We estimated the age of nodes and reconstruct the ancestral character states of 43 phenotypic characters. We find high levels of phylogenetic signal for nearly all traits, indicating the phylogeny carries substantial predictive power. The earliest cyanobacterial lineages likely lived in freshwater habitats, had small cell diameters, were benthic or sessile, and possibly epilithic/endolithic with a sheath. We jointly analyzed a subset of 25 binary traits to determine whether rates of trait evolution have shifted over time in conjunction with major geologic events. Phylogenetic comparative analysis reveal an overriding signal of decreasing rates of trait evolution through time. Furthermore, the data suggest two major rate shifts in trait evolution associated with bursts of evolutionary innovation. The first rate shift occurs in the aftermath of the Great Oxidation Event and “Snowball Earth” glaciations and is associated with decrease in the evolutionary rates around 1.8–1.6 Ga. This rate shift seems to indicate the end of a major diversification of cyanobacterial phenotypes–particularly related to traits associated with filamentous morphology, heterocysts and motility in freshwater ecosystems. Another burst appears around the time of the Neoproterozoic Oxidation Event in the Neoproterozoic, and is associated with the acquisition of traits involved in planktonic growth in marine habitats. Our results demonstrate how uniting genomic and phenotypic datasets in extant bacterial species can shed light on billion-year old events in Earth’s history.

Methanogenesis and the Wood-Ljungdahl Pathway: An Ancient, Versatile, and Fragile Association

Methanogenesis coupled to the Wood–Ljungdahl pathway is one of the most ancient metabolisms for energy generation and carbon fixation in the Archaea. Recent results are sensibly changing our view on the diversity of methane-cycling capabilities in this Domain of Life. The availability of genomic sequences from uncharted branches of the archaeal tree has highlighted the existence of novel methanogenic lineages phylogenetically distant to previously known ones, such as the Methanomassiliicoccales. At the same time, phylogenomic analyses have suggested a methanogenic ancestor for all Archaea, implying multiple independent losses of this metabolism during archaeal diversification. This prediction has been strengthened by the report of genes involved in methane cycling in members of the Bathyarchaeota (a lineage belonging to the TACK clade), representing the first indication of the presence of methanogenesis outside of the Euryarchaeota. In light of these new data, we discuss how the association between methanogenesis and the Wood–Ljungdahl pathway appears to be much more flexible than previously thought, and might provide information on the processes that led to loss of this metabolism in many archaeal lineages. The combination of environmental microbiology, experimental characterization and phylogenomics opens up exciting avenues of research to unravel the diversity and evolutionary history of fundamental metabolic pathways.

Origin of eukaryotes from within archaea, archaeal eukaryome and bursts of gene gain: eukaryogenesis just made easier?

The origin of eukaryotes is a fundamental, forbidding evolutionary puzzle. Comparative genomic analysis clearly shows that the last eukaryotic common ancestor (LECA) possessed most of the signature complex features of modern eukaryotic cells, in particular the mitochondria, the endomembrane system including the nucleus, an advanced cytoskeleton and the ubiquitin network. Numerous duplications of ancestral genes, e.g. DNA polymerases, RNA polymerases and proteasome subunits, also can be traced back to the LECA. Thus, the LECA was not a primitive organism and its emergence must have resulted from extensive evolution towards cellular complexity. However, the scenario of eukaryogenesis, and in particular the relationship between endosymbiosis and the origin of eukaryotes, is far from being clear. Four recent developments provide new clues to the likely routes of eukaryogenesis. First, evolutionary reconstructions suggest complex ancestors for most of the major groups of archaea, with the subsequent evolution dominated by gene loss. Second, homologues of signature eukaryotic proteins, such as actin and tubulin that form the core of the cytoskeleton or the ubiquitin system, have been detected in diverse archaea. The discovery of this ‘dispersed eukaryome’ implies that the archaeal ancestor of eukaryotes was a complex cell that might have been capable of a primitive form of phagocytosis and thus conducive to endosymbiont capture. Third, phylogenomic analyses converge on the origin of most eukaryotic genes of archaeal descent from within the archaeal evolutionary tree, specifically, the TACK superphylum. Fourth, evidence has been presented that the origin of the major archaeal phyla involved massive acquisition of bacterial genes. Taken together, these findings make the symbiogenetic scenario for the origin of eukaryotes considerably more plausible and the origin of the organizational complexity of eukaryotic cells more readily explainable than they appeared until recently.

Archaeal Distribution in Moonmilk Deposits from Alpine Caves and Their Ecophysiological Potential

(Alpine) caves are, in general, windows into the Earth's subsurface. Frequently occurring structures in caves such as moonmilk (secondary calcite deposits) offer the opportunity to study intraterrestrial microbial communities, adapted to oligotrophic and cold conditions. This is an important research field regarding the dimensions of subsurface systems and cold regions on Earth. On a methodological level, moonmilk deposits from 11 caves in the Austrian Alps were collected aseptically and investigated using a molecular (qPCR and DGGE sequencing-based) methodology in order to study the occurrence, abundance, and diversity of the prevailing native Archaea community. Furthermore, these Archaea were enriched in complex media and studied regarding their physiology, with a media selection targeting different physiological requirements, e.g. methanogenesis and ammonia oxidation. The investigation of the environmental samples showed that all moonmilk deposits were characterized by the presence of the same few habitat-specific archaeal species, showing high abundances and constituting about 50 % of the total microbial communities. The largest fraction of these Archaea was ammonia-oxidizing Thaumarchaeota, while another abundant group was very distantly related to extremophilic Euryarchaeota (Moonmilk Archaea). The archaeal community showed a depth- and oxygen-dependent stratification. Archaea were much more abundant (around 80 %), compared to bacteria, in the actively forming surface part of moonmilk deposits, decreasing to about 5 % down to the bedrock. Via extensive cultivation efforts, it was possible to enrich the enigmatic Moonmilk Archaea and also AOA significantly above the level of bacteria. The most expedient prerequisites for cultivating Moonmilk Archaea were a cold temperature, oligotrophic conditions, short incubation times, a moonmilk surface inoculum, the application of erythromycin, and anaerobic (microaerophilic) conditions. On a physiological level, it seems that methanogenesis is of marginal importance, while ammonia oxidation and a still undiscovered metabolic pathway are vital elements in the (archaeal) moonmilk biome.

Gene Acquisitions from Bacteria at the Origins of Major Archaeal Clades Are Vastly Overestimated

In a recent article, Nelson-Sathi et al. (NS) report that the origins of major archaeal lineages (MAL) correspond to massive group-specific gene acquisitions via HGT from bacteria (Nelson-Sathi et al. 2015. Origins of major archaeal clades correspond to gene acquisitions from bacteria. Nature 517(7532):77-80.). If correct, this would have fundamental implications for the process of diversification in microbes. However, a reexamination of these data and results shows that the methodology used by NS systematically inflates the number of genes acquired at the root of each MAL, and incorrectly assumes bacterial origins for these genes. A reanalysis of their data with appropriate phylogenetic models accounting for the dynamics of gene gain and loss between lineages supports the continuous acquisition of genes over long periods in the evolution of Archaea.

Archaeal ancestors of eukaryotes: not so elusive any more

The origin of eukaryotes is one of the hardest problems in evolutionary biology and sometimes raises the ominous specter of irreducible complexity. Reconstruction of the gene repertoire of the last eukaryotic common ancestor (LECA) has revealed a highly complex organism with a variety of advanced features but no detectable evolutionary intermediates to explain their origin. Recently, however, genome analysis of diverse archaea led to the discovery of apparent ancestral versions of several signature eukaryotic systems, such as the actin cytoskeleton and the ubiquitin network, that are scattered among archaea. These findings inspired the hypothesis that the archaeal ancestor of eukaryotes was an unusually complex form with an elaborate intracellular organization. The latest striking discovery made by deep metagenomic sequencing vindicates this hypothesis by showing that in phylogenetic trees eukaryotes fall within a newly identified archaeal group, the Lokiarchaeota, which combine several eukaryotic signatures previously identified in different archaea. The discovery of complex archaea that are the closest living relatives of eukaryotes is most compatible with the symbiogenetic scenario for eukaryogenesis.

The universal tree of life: an update

Biologists used to draw schematic “universal” trees of life as metaphors illustrating the history of life. It is indeed a priori possible to construct an organismal tree connecting the three major domains of ribosome encoding organisms: Archaea, Bacteria and Eukarya, since they originated by cell division from LUCA. Several universal trees based on ribosomal RNA sequence comparisons proposed at the end of the last century are still widely used, although some of their main features have been challenged by subsequent analyses. Several authors have proposed to replace the traditional universal tree with a ring of life, whereas others have proposed more recently to include viruses as new domains. These proposals are misleading, suggesting that endosymbiosis can modify the shape of a tree or that viruses originated from the last universal common ancestor (LUCA). I propose here an updated version of Woese’s universal tree that includes several rootings for each domain and internal branching within domains that are supported by recent phylogenomic analyses of domain specific proteins. The tree is rooted between Bacteria and Arkarya, a new name proposed for the clade grouping Archaea and Eukarya. A consensus version, in which each of the three domains is unrooted, and a version in which eukaryotes emerged within archaea are also presented. This last scenario assumes the transformation of a modern domain into another, a controversial evolutionary pathway. Viruses are not indicated in these trees but are intrinsically present because they infect the tree from its roots to its leaves. Finally, I present a detailed tree of the domain Archaea, proposing the sub-phylum neo-Euryarchaeota for the monophyletic group of euryarchaeota containing DNA gyrase. These trees, that will be easily updated as new data become available, could be useful to discuss controversial scenarios regarding early life evolution.

Quest for Orthologs Entails Quest for Tree of Life: In Search of the Gene Stream

Quest for Orthologs (QfO) is a community effort with the goal to improve and benchmark orthology predictions. As quality assessment assumes prior knowledge on species phylogenies, we investigated the congruency between existing species trees by comparing the relationships of 147 QfO reference organisms from six Tree of Life (ToL)/species tree projects: The National Center for Biotechnology Information (NCBI) taxonomy, Opentree of Life, the sequenced species/species ToL, the 16S ribosomal RNA (rRNA) database, and trees published by Ciccarelli et al. (Ciccarelli FD, et al. 2006. Toward automatic reconstruction of a highly resolved tree of life. Science 311:1283–1287) and by Huerta-Cepas et al. (Huerta-Cepas J, Marcet-Houben M, Gabaldon T. 2014. A nested phylogenetic reconstruction approach provides scalable resolution in the eukaryotic Tree Of Life. PeerJ PrePrints 2:223) Our study reveals that each species tree suggests a different phylogeny: 87 of the 146 (60%) possible splits of a dichotomous and rooted tree are congruent, while all other splits are incongruent in at least one of the species trees. Topological differences are observed not only at deep speciation events, but also within younger clades, such as Hominidae, Rodentia, Laurasiatheria, or rosids. The evolutionary relationships of 27 archaea and bacteria are highly inconsistent. By assessing 458,108 gene trees from 65 genomes, we show that consistent species topologies are more often supported by gene phylogenies than contradicting ones. The largest concordant species tree includes 77 of the QfO reference organisms at the most. Results are summarized in the form of a consensus ToL (http://swisstree.vital-it.ch/species_tree) that can serve different benchmarking purposes.

Uncultivated thermophiles: current status and spotlight on 'Aigarchaeota'

Meta-analysis of cultivation-independent sequence data shows that geothermal systems host an abundance of novel organisms, representing a vast unexplored phylogenetic and functional diversity among yet-uncultivated thermophiles. A number of thermophiles have recently been interrogated using metagenomic and/or single-cell genomic approaches, including members of taxonomic groups that inhabit both thermal and non-thermal environments, such as 'Acetothermia' (OP1) and 'Atribacteria' (OP9/JS1), as well as the exclusively thermophilic lineages 'Korarchaeota', 'Calescamantes' (EM19), 'Fervidibacteria' (OctSpA1-106), and 'Aigarchaeota' (HWCG-I). The 'Aigarchaeota', a sister lineage to the Thaumarchaeota, likely includes both hyperthermophiles and moderate thermophiles. They inhabit terrestrial, marine, and subsurface thermal environments and comprise at least nine genus-level lineages, several of which are globally distributed.

Capturing prokaryotic dark matter genomes

Prokaryotes are the most diverse and abundant cellular life forms on Earth. Most of them, identified by indirect molecular approaches, belong to microbial dark matter. The advent of metagenomic and single-cell genomic approaches has highlighted the metabolic capabilities of numerous members of this dark matter through genome reconstruction. Thus, linking functions back to the species has revolutionized our understanding of how ecosystem function is sustained by the microbial world. This review will present discoveries acquired through the illumination of prokaryotic dark matter genomes by these innovative approaches.