The evolution pathway of ammonia-oxidizing archaea shaped by major geological events

Exploring protein N-glycosylation in ammonia-oxidizing Nitrososphaerota archaea through glycoproteomic analysis

Ammonia-oxidizing archaea of the phylum Nitrososphaerota, formerly known as Thaumarchaeota, are globally distributed and play critical roles in the nitrogen and carbon cycles, particularly in environments with low ammonia concentrations. Like most archaea, Nitrososphaerota cells are enveloped by S-layer proteins, implicated in concentrating ammonium ions. These proteins are typically modified post-translationally by N-glycans, which often play significant roles in various biological processes, including protein function regulation, protection from phages, and environmental adaptation. Nevertheless, the glycobiological characteristics of Nitrososphaerota remain largely unexplored. Here, we investigated the glycoproteome of ammonia-oxidizing Nitrososphaerota, specifically focusing on the terrestrial Nitrososphaera viennensis and the marine Nitrosopumilus piranensis. Both species exhibited similar protein arrays throughout their growth phases, including those associated with N-glycosylation. Ns. viennensis consistently exhibited N-glycosylation predominantly on an S-layer protein and multicopper oxidase domain-containing proteins throughout all growth phases, with a marked increase during and after the late exponential phase. The glycan, characterized as a novel hexasaccharide with a chitobiose core, is hypothesized to play a role in nitrogen storage due to its probable nitrogen-rich composition, modifying asparagine residues within the conserved triplet sequence (Asn-X-Ser or -Thr). In contrast, Np. piranensis also showed a high abundance of S-layer protein but displayed no apparent N-glycosylation on any protein, suggesting variability in cell surface physical properties between these archaea. Despite similarities in their proteomes and energy metabolism, these two archaea exhibited significant differences in post-translational modification of proteins, revealing previously unrecognized diversity that may have implications for understanding their adaptive transitions to diverse environments.

IMPORTANCE

Autotrophic ammonia-oxidizing archaea of the phylum Nitrososphaerota, formerly known as Thaumarchaeota, are notoriously difficult to culture yet play important roles in the global nitrogen and carbon cycles. Inhabiting environments with extremely low ammonia concentrations, these archaea are expected to conserve ammonia strictly for energy production. However, using advanced liquid chromatography-tandem mass spectrometry and nuclear magnetic resonance techniques, we discovered that one of these archaea decorates its cell surface proteins with the most nitrogen-rich glycan identified to date, suggesting a previously unrecognized function of protein glycosylation in nitrogen storage. This newly identified N-glycan, with a chitobiose core similar to those in Thermoproteota and eukaryotes, not only deepens our understanding of archaeal evolution but also underscores the molecular adaptations enabling these archaea to thrive in diverse environments.

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.

Novel hot spring Thermoproteota support vertical inheritance of ammonia oxidation and carbon fixation in Nitrososphaeria

Aerobic ammonia oxidation is crucial to the nitrogen cycle and is only known to be performed by a small number of bacterial lineages [ammonia-oxidizing bacteria (AOB)] and a single lineage of archaea belonging to the Nitrososphaeria class of Thermoproteota [ammonia-oxidizing Archaea (AOA)]. Most cultivated AOA originate from marine or soil environments, but this may capture only a limited subset of the full diversity of this clade. Here, we describe several genomes of AOA from metagenomic sequencing of a hot spring microbial mat, representing several poorly characterized basal lineages that may be important for understanding the early evolution of archaeal ammonia oxidation. These genomes include a novel genus most closely related to Nitrososphaera as well as novel species belonging to the genera Nitrosotenuis, Nitrososphaera and Nitrosotalea. Furthermore, the distributions and phylogenetic relationships of key metabolic genes support a history of vertical inheritance of ammonia oxidation and carbon fixation from the last common ancestor of crown group AOA.

Niche Partitioning and Intraspecific Variation of Thaumarchaeota in Deep Ocean Sediments

Deep-sea sediments contain a large number of Thaumarchaeota that are phylogenetically distinct from their pelagic counterparts. However, their ecology and evolutionary adaptations are not well understood. Metagenomic analyses were conducted on samples from various depths of a 750-cm sediment core collected from the Mariana Trench Challenger Deep. The abundance of Thaumarchaeota and archaeal amoA generally decreased with depth, except for an unexpected peak midway through the core. The thaumarchaeotal metagenome-assembled genomes were classified into diverse phylogenetic clusters associated with amoA-NP-γ, amoA-NP-θ, and amoA-NP-δ of ammonia-oxidising Thaumarchaeota and non-ammonia-oxidising lineages. The most abundant group was within amoA-NP-γ, which is usually found in coastal and shallow habitats, indicating potential niche expansion from marine shallow to hadal environments. This benthic group showed within-species genomic variations compared to the previously identified Hadal water group, suggesting microdiversification of hadal Thaumarchaeota along with niche separation between benthic and pelagic environments. Evolutionary adaptations associated with the benthic-to-pelagic transition included reduced genome size, loss of motility/cell adhesion, altered energy metabolism, and different mechanisms for substrate acquisition and regulation (e.g., ammonium). These findings offer new insights into the evolution of hadal Thaumarchaeota and demonstrate, for the first time, intraspecies-level genomic variation in Thaumarchaeota related to the benthic-versus-pelagic niche partitioning in the deep ocean.

Phylogeny and evolution of dissimilatory sulfite reduction in prokaryotes

Sulfate is the second most common nonmetallic ion in modern oceans, as its concentration dramatically increased alongside tectonic activity and atmospheric oxidation in the Proterozoic. Microbial sulfate/sulfite metabolism, involving organic carbon or hydrogen oxidation, is linked to sulfur and carbon biogeochemical cycles. However, the coevolution of microbial sulfate/sulfite metabolism and Earth's history remains unclear. Here, we conducted a comprehensive phylogenetic analysis to explore the evolutionary history of the dissimilatory sulfite reduction (Dsr) pathway. The phylogenies of the Dsr-related genes presented similar branching patterns but also some incongruencies, indicating the complex origin and evolution of Dsr. Among these genes, dsrAB is the hallmark of sulfur-metabolizing prokaryotes. Our detailed analyses suggested that the evolution of dsrAB was shaped by vertical inheritance and multiple horizontal gene transfer events and that selection pressure varied across distinct lineages. Dated phylogenetic trees indicated that key evolutionary events of dissimilatory sulfur-metabolizing prokaryotes were related to the Great Oxygenation Event (2.4-2.0 Ga) and several geological events in the "Boring Billion" (1.8-0.8 Ga), including the fragmentation of the Columbia supercontinent (approximately 1.6 Ga), the rapid increase in marine sulfate (1.3-1.2 Ga), and the Neoproterozoic glaciation event (approximately 1.0 Ga). We also proposed that the voluminous iron formations (approximately 1.88 Ga) might have induced the metabolic innovation of iron reduction. In summary, our study provides new insights into Dsr evolution and a systematic view of the coevolution of dissimilatory sulfur-metabolizing prokaryotes and the Earth's environment.

Thaumarchaeota from deep-sea methane seeps provide novel insights into their evolutionary history and ecological implications

BACKGROUND

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota mediate the rate-limiting step of nitrification and remove the ammonia that inhibits the aerobic metabolism of methanotrophs. However, the AOA that inhabit deep-sea methane-seep surface sediments (DMS) are rarely studied. Here, we used global DMS metagenomics and metagenome-assembled genomes (MAGs) to investigate the metabolic activity, evolutionary history, and ecological contributions of AOA. Expression of AOA-specific ammonia-oxidizing gene (amoA) was examined in the sediments collected from the South China Sea (SCS) to identify their active ammonia metabolism in the DMS.

RESULTS

Our analysis indicated that AOA contribute > 75% to the composition of ammonia-utilization genes within the surface layers (above 30 cm) of global DMS. The AOA-specific ammonia-oxidizing gene was actively expressed in the DMS collected from the SCS. Phylogenomic analysis of medium-/high-quality MAGs from 18 DMS-AOA indicated that they evolved from ancestors in the barren deep-sea sediment and then expanded from the DMS to shallow water forming an amoA-NP-gamma clade-affiliated lineage. Molecular dating suggests that the DMS-AOA origination coincided with the Neoproterozoic oxidation event (NOE), which occurred ~ 800 million years ago (mya), and their expansion to shallow water coincided with the Sturtian glaciation (~ 713 mya). Comparative genomic analysis suggests that DMS-AOA exhibit higher requirement of carbon source for protein synthesis with enhanced genomic capability for osmotic regulation, motility, chemotaxis, and utilization of exogenous organic compounds, suggesting it could be more heterotrophic compared with other lineages.

CONCLUSION

Our findings provide new insights into the evolutionary history of AOA within the Thaumarchaeota, highlighting their critical roles in nitrogen cycling in the global DMS ecosystems. Video Abstract.

Dating Ammonia-Oxidizing Bacteria with Abundant Eukaryotic Fossils

Evolution of a complete nitrogen (N) cycle relies on the onset of ammonia oxidation, which aerobically converts ammonia to nitrogen oxides. However, accurate estimation of the antiquity of ammonia-oxidizing bacteria (AOB) remains challenging because AOB-specific fossils are absent and bacterial fossils amenable to calibrate molecular clocks are rare. Leveraging the ancient endosymbiosis of mitochondria and plastid, as well as using state-of-the-art Bayesian sequential dating approach, we obtained a timeline of AOB evolution calibrated largely by eukaryotic fossils. We show that the first AOB evolved in marine Gammaproteobacteria (Gamma-AOB) and emerged between 2.1 and 1.9 billion years ago (Ga), thus postdating the Great Oxidation Event (GOE; 2.4 to 2.32 Ga). To reconcile the sedimentary N isotopic signatures of ammonia oxidation occurring near the GOE, we propose that ammonia oxidation likely occurred at the common ancestor of Gamma-AOB and Gammaproteobacterial methanotrophs, or the actinobacterial/verrucomicrobial methanotrophs which are known to have ammonia oxidation activities. It is also likely that nitrite was transported from the terrestrial habitats where ammonia oxidation by archaea took place. Further, we show that the Gamma-AOB predated the anaerobic ammonia-oxidizing (anammox) bacteria, implying that the emergence of anammox was constrained by the availability of dedicated ammonia oxidizers which produce nitrite to fuel anammox. Our work supports a new hypothesis that N redox cycle involving nitrogen oxides evolved rather late in the ocean.

Ultrastructural insights into cellular organization, energy storage and ribosomal dynamics of an ammonia-oxidizing archaeon from oligotrophic oceans

Introduction

Nitrososphaeria, formerly known as Thaumarchaeota, constitute a diverse and widespread group of ammonia-oxidizing archaea (AOA) inhabiting ubiquitously in marine and terrestrial environments, playing a pivotal role in global nitrogen cycling. Despite their importance in Earth's ecosystems, the cellular organization of AOA remains largely unexplored, leading to a significant unanswered question of how the machinery of these organisms underpins metabolic functions.

Methods

In this study, we combined spherical-chromatic-aberration-corrected cryo-electron tomography (cryo-ET), scanning transmission electron microscopy (STEM), and energy dispersive X-ray spectroscopy (EDS) to unveil the cellular organization and elemental composition of Nitrosopumilus maritimus SCM1, a representative member of marine Nitrososphaeria.

Results and Discussion

Our tomograms show the native ultrastructural morphology of SCM1 and one to several dense storage granules in the cytoplasm. STEM-EDS analysis identifies two types of storage granules: one type is possibly composed of polyphosphate and the other polyhydroxyalkanoate. With precise measurements using cryo-ET, we observed low quantity and density of ribosomes in SCM1 cells, which are in alignment with the documented slow growth of AOA in laboratory cultures. Collectively, these findings provide visual evidence supporting the resilience of AOA in the vast oligotrophic marine environment.

Mapping Geological Events and Nitrogen Fixation Evolution Onto the Timetree of the Evolution of Nitrogen-Fixation Genes

Nitrogen is essential for all organisms, but biological nitrogen fixation (BNF) occurs only in a small fraction of prokaryotes. Previous studies divided nitrogenase-gene-carrying prokaryotes into Groups I to IV and provided evidence that BNF first evolved in bacteria. This study constructed a timetree of the evolution of nitrogen-fixation genes and estimated that archaea evolved BNF much later than bacteria and that nitrogen-fixing cyanobacteria evolved later than 1,900 MYA, considerably younger than the previous estimate of 2,200 MYA. Moreover, Groups III and II/I diverged ∼2,280 MYA, after the Kenorland supercontinent breakup (∼2,500-2,100 MYA) and the Great Oxidation Event (∼2,400-2,100 MYA); Groups III and Vnf/Anf diverged ∼2,086 MYA, after the Yarrabubba impact (∼2,229 MYA); and Groups II and I diverged ∼1,920 MYA, after the Vredefort impact (∼2,023 MYA). In summary, this study provided a timescale of BNF events and discussed the possible effects of geological events on BNF evolution.

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

Oxygen-driven divergence of marine group II archaea reflected by transitions of superoxide dismutases

IMPORTANCE

Reactive oxygen species (ROS), including superoxide anion, is a series of substances that cause oxidative stress for all organisms. Marine group II (MGII) archaea are mainly live in the surface seawater and exposed to considerable ROS. Therefore, it is important to understand the antioxidant capacity of MGII. Our research found that Fe/Mn- superoxide dismutase (Fe/MnSOD) may be more suitable for MGII to resist oxidative damage, and the changes in oxygen concentrations and SOD metallic cofactors play an important role in the selection of SOD by the 17 clades of MGII, which in turn affects the species differentiation of MGII. Overall, this study provides insight into the co-evolutionary history of these uncultivated marine archaea with the earth system.

Niche breadth specialization impacts ecological and evolutionary adaptation following environmental change

Ecological theory predicts that organismal distribution and abundance depend on the ability to adapt to environmental change. It also predicts that eukaryotic specialists and generalists will dominate in extreme environments or following environmental change, respectively. This theory has attracted little attention in prokaryotes, especially in archaea, which drive major global biogeochemical cycles. We tested this concept in Thaumarchaeota using pH niche breadth as a specialization factor. Responses of archaeal growth and activity to pH disturbance were determined empirically in manipulated, long-term, pH-maintained soil plots. The distribution of specialists and generalists was uneven over the pH range, with specialists being more limited to the extreme range. Nonetheless, adaptation of generalists to environmental change was greater than that of specialists, except for environmental changes leading to more extreme conditions. The balance of generalism and specialism over longer timescales was further investigated across evolutionary history. Specialists and generalists diversified at similar rates, reflecting balanced benefits of each strategy, but a higher transition rate from generalists to specialists than the reverse was demonstrated, suggesting that metabolic specialism is more easily gained than metabolic versatility. This study provides evidence for a crucial ecological concept in prokaryotes, significantly extending our understanding of archaeal adaptation to environmental change.

Novel order-level lineage of ammonia-oxidizing archaea widespread in marine and terrestrial environments

Ammonia-oxidizing archaea (AOA) are among the most ubiquitous and abundant archaea on Earth, widely distributed in marine, terrestrial, and geothermal ecosystems. However, the genomic diversity, biogeography, and evolutionary process of AOA populations in subsurface environments are vastly understudied compared to those in marine and soil systems. Here, we report a novel AOA order Candidatus (Ca.) Nitrosomirales which forms a sister lineage to the thermophilic Ca. Nitrosocaldales. Metagenomic and 16S rRNA gene-read mapping demonstrates the abundant presence of Nitrosomirales AOA in various groundwater environments and their widespread distribution across a range of geothermal, terrestrial, and marine habitats. Terrestrial Nitrosomirales AOA show the genetic capacity of using formate as a source of reductant and using nitrate as an alternative electron acceptor. Nitrosomirales AOA appear to have acquired key metabolic genes and operons from other mesophilic populations via horizontal gene transfer, including genes encoding urease, nitrite reductase, and V-type ATPase. The additional metabolic versatility conferred by acquired functions may have facilitated their radiation into a variety of subsurface, marine, and soil environments. We also provide evidence that each of the four AOA orders spans both marine and terrestrial habitats, which suggests a more complex evolutionary history for major AOA lineages than previously proposed. Together, these findings establish a robust phylogenomic framework of AOA and provide new insights into the ecology and adaptation of this globally abundant functional guild.

Temperature, pH, and oxygen availability contributed to the functional differentiation of ancient Nitrososphaeria

Ammonia-oxidizing Nitrososphaeria are among the most abundant archaea on Earth and have profound impacts on the biogeochemical cycles of carbon and nitrogen. In contrast to these well-studied ammonia-oxidizing archaea (AOA), deep-branching non-AOA within this class remain poorly characterized because of a low number of genome representatives. Here, we reconstructed 128 Nitrososphaeria metagenome-assembled genomes from acid mine drainage and hot spring sediment metagenomes. Comparative genomics revealed that extant non-AOA are functionally diverse, with capacity for carbon fixation, carbon monoxide oxidation, methanogenesis, and respiratory pathways including oxygen, nitrate, sulfur, or sulfate, as potential terminal electron acceptors. Despite their diverse anaerobic pathways, evolutionary history inference suggested that the common ancestor of Nitrososphaeria was likely an aerobic thermophile. We further surmise that the functional differentiation of Nitrososphaeria was primarily shaped by oxygen, pH, and temperature, with the acquisition of pathways for carbon, nitrogen, and sulfur metabolism. Our study provides a more holistic and less biased understanding of the diversity, ecology, and deep evolution of the globally abundant Nitrososphaeria.

Metagenomic characterization of a novel non-ammonia-oxidizing Thaumarchaeota from hadal sediment

BACKGROUND

The hadal sediment, found at an ocean depth of more than 6000 m, is geographically isolated and under extremely high hydrostatic pressure, resulting in a unique ecosystem. Thaumarchaeota are ubiquitous marine microorganisms predominantly present in hadal environments. While there have been several studies on Thaumarchaeota there, most of them have primarily focused on ammonia-oxidizing archaea (AOA). However, systematic metagenomic research specifically targeting heterotrophic non-AOA Thaumarchaeota is lacking.

RESULTS

In this study, we explored the metagenomes of Challenger Deep hadal sediment, focusing on the Thaumarchaeota. Functional analysis of sequence reads revealed the potential contribution of Thaumarchaeota to recalcitrant dissolved organic matter degradation. Metagenome assembly binned one new group of hadal sediment-specific and ubiquitously distributed non-AOA Thaumarchaeota, named Group-3.unk. Pathway reconstruction of this new type of Thaumarchaeota also supports heterotrophic characteristics of Group-3.unk, along with ABC transporters for the uptake of amino acids and carbohydrates and catabolic utilization of these substrates. This new clade of Thaumarchaeota also contains aerobic oxidation of carbon monoxide-related genes. Complete glyoxylate cycle is a distinctive feature of this clade in supplying intermediates of anabolic pathways. The pan-genomic and metabolic analyses of metagenome-assembled genomes belonging to Group-3.unk Thaumarchaeota have highlighted distinctions, including the dihydroxy phthalate decarboxylase gene associated with the degradation of aromatic compounds and the absence of genes related to the synthesis of some types of vitamins compared to AOA. Notably, Group-3.unk shares a common feature with deep ocean AOA, characterized by their high hydrostatic pressure resistance, potentially associated with the presence of V-type ATP and di-myo-inositol phosphate syntheses-related genes. The enrichment of organic matter in hadal sediments might be attributed to the high recruitment of sequence reads of the Group-3.unk clade of heterotrophic Thaumarchaeota in the trench sediment. Evolutionary and genetic dynamic analyses suggest that Group-3 non-AOA consists of mesophilic Thaumarchaeota organisms. These results indicate a potential role in the transition from non-AOA to AOA Thaumarchaeota and from thermophilic to mesophilic Thaumarchaeota, shedding light on recent evolutionary pathways.

CONCLUSIONS

One novel clade of heterotrophic non-AOA Thaumarchaeota was identified through metagenome analysis of sediments from Challenger Deep. Our study provides insight into the ecology and genomic characteristics of the new sub-group of heterotrophic non-AOA Thaumarchaeota, thereby extending the knowledge of the evolution of Thaumarchaeota. Video Abstract.

Genomic insight into strategy, interaction and evolution of nitrifiers in metabolizing key labile-dissolved organic nitrogen in different environmental niches

Ammonia-oxidizing archaea (AOA) and bacteria (AOB), nitrite-oxidizing bacteria (NOB), and complete ammonia oxidizers (comammox) are responsible for nitrification in nature; however, some groups have been reported to utilize labile-dissolved organic nitrogen (LDON) for satisfying nitrogen demands. To understand the universality of their capacity of LDON metabolism, we collected 70 complete genomes of AOA, AOB, NOB, and comammox from typical environments for exploring their potentials in the metabolism of representative LDON (urea, polyamines, cyanate, taurine, glycine betaine, and methylamine). Genomic analyses showed that urea was the most popular LDON used by nitrifiers. Each group harbored unique urea transporter genes (AOA: dur3 and utp, AOB: utp, and NOB and comammox: urtABCDE and utp) accompanied by urease genes ureABC. The differentiation in the substrate affinity of these transporters implied the divergence of urea utilization efficiency in nitrifiers, potentially driving them into different niches. The cyanate transporter (cynABD and focA/nirC) and degradation (cynS) genes were detected mostly in NOB, indicating their preference for a wide range of nitrogen substrates to satisfy high nitrogen demands. The lack of genes involved in the metabolism of polyamines, taurine, glycine betaine, and methylamines in most of nitrifiers suggested that they were not able to serve as a source of ammonium, only if they were degraded or oxidized extracellularly as previously reported. The phylogenetic analyses assisted with comparisons of GC% and the Codon Adaptation Index between target genes and whole genomes of nitrifiers implied that urea metabolic genes dur3 and ureC in AOA evolved independently from bacteria during the transition from Thaumarchaeota to AOA, while utp in terrestrial AOA was acquired from bacteria via lateral gene transfer (LGT). Cyanate transporter genes cynS and focA/nirC detected only in a terrestrial AOA Candidadus Nitrsosphaera gargensis Ga9.2 could be gained synchronously with Nitrospira of NOB by an ancient LGT. Our results indicated that LDON utilization was a common feature in nitrifiers, but metabolic potentials were different among nitrifiers, possibly being intensely interacted with their niches, survival strategies, and evolutions.

A timeline of bacterial and archaeal diversification in the ocean

Microbial plankton play a central role in marine biogeochemical cycles, but the timing in which abundant lineages diversified into ocean environments remains unclear. Here, we reconstructed the timeline in which major clades of bacteria and archaea colonized the ocean using a high-resolution benchmarked phylogenetic tree that allows for simultaneous and direct comparison of the ages of multiple divergent lineages. Our findings show that the diversification of the most prevalent marine clades spans throughout a period of 2.2 Ga, with most clades colonizing the ocean during the last 800 million years. The oldest clades - SAR202, SAR324, Ca. Marinimicrobia, and Marine Group II - diversified around the time of the Great Oxidation Event, during which oxygen concentration increased but remained at microaerophilic levels throughout the Mid-Proterozoic, consistent with the prevalence of some clades within these groups in oxygen minimum zones today. We found the diversification of the prevalent heterotrophic marine clades SAR11, SAR116, SAR92, SAR86, and Roseobacter as well as the Marine Group I to occur near to the Neoproterozoic Oxygenation Event (0.8-0.4 Ga). The diversification of these clades is concomitant with an overall increase of oxygen and nutrients in the ocean at this time, as well as the diversification of eukaryotic algae, consistent with the previous hypothesis that the diversification of heterotrophic bacteria is linked to the emergence of large eukaryotic phytoplankton. The youngest clades correspond to the widespread phototrophic clades Prochlorococcus, Synechococcus, and Crocosphaera, whose diversification happened after the Phanerozoic Oxidation Event (0.45-0.4 Ga), in which oxygen concentrations had already reached their modern levels in the atmosphere and the ocean. Our work clarifies the timing at which abundant lineages of bacteria and archaea colonized the ocean, thereby providing key insights into the evolutionary history of lineages that comprise the majority of prokaryotic biomass in the modern ocean.

Diversity and salinity adaptations of ammonia oxidizing archaea in three estuaries of China

Ammonia-oxidizing archaea (AOA) are ubiquitously found in diverse habitats and play pivotal roles in the nitrogen and carbon cycle, especially in estuarine and coastal environments. Despite the fact that the diversity and distribution of AOA are thought to be tightly linked to habitats, little is known about the relationship that underpins their genomic traits, adaptive potentials, and ecological niches. Here, we have characterized and compared the AOA community in three estuaries of China using metagenomics. AOA were the dominant ammonia oxidizers in the three estuaries. Through phylogenetic analyses, five major AOA groups were identified, including the Nitrosomarinus-like, Nitrosopumilus-like, Aestuariumsis-like, Nitrosarchaeum-like, and Nitrosopelagicus-like groups. Statistical analyses showed that the aquatic and sedimentary AOA communities were mainly influenced by spatial factors (latitude and water depth) and environmental factors (salinity, pH, and dissolved oxygen) in estuaries, respectively. Compared to AOA dwelling in terrestrial and marine habitats, estuarine AOA encoded more genes involved in glucose and amino acid metabolism, transport systems, osmotic control, and cell motility. The low proteome isoelectric points (pI), high content of acidic amino acids, and the presence of potassium ion and mechanosensitive channels suggest a "salt-in" strategy for estuarine AOA to counteract high osmolarity in their surroundings. Our findings have indicated potential adaptation strategies and highlighted their importance in the estuarine nitrogen and carbon cycles. KEY POINTS: • Spatial and environmental factors influence water and sediment AOA respectively. • Estuarine AOA share low proteome isoelectric value and high acid amino acids content. • AOA adaptation to estuaries is likely resulted from their unique genomic features.

Nitrogen and phosphorous acquisition strategies drive coexistence patterns among archaeal lineages in soil

Soil represents the largest reservoir of Archaea on Earth. Present-day archaeal diversity in soils globally is dominated by members of the class Nitrososphaeria. The evolutionary radiation of this class is thought to reflect adaptations to a wide range of temperatures, pH, and other environmental conditions. However, the mechanisms that govern competition and coexistence among Nitrososphaeria lineages in soil remain poorly understood. Here we show that predominant soil Nitrososphaeria lineages compose a patchwork of gene inventory and expression profiles for ammonia, urea, and phosphate utilization. In contrast, carbon fixation, respiration, and ATP synthesis genes are conserved and expressed consistently among predominant phylotypes across 12 major evolutionary lineages commonly found in soil. In situ gene expression profiles closely resemble pure culture reference strains under optimal growth conditions. Together, these results reveal resource-based coexistence patterns among Nitrososphaeria lineages and suggest complementary ecophysiological niches associated with differential nutrient acquisition strategies among globally predominant archaeal lineages in soil.

Thermophilic microorganisms involved in the nitrogen cycle in thermal environments: Advances and prospects

Thermophilic microorganisms mediated significant element cycles and material conversion in the early Earth as well as mediating current thermal environments. Over the past few years, versatile microbial communities that drive the nitrogen cycle have been identified in thermal environments. Understanding the microbial-mediated nitrogen cycling processes in these thermal environments has important implications for the cultivation and application of thermal environment microorganisms as well as for exploring the global nitrogen cycle. This work provides a comprehensive review of different thermophilic nitrogen-cycling microorganisms and processes, which are described in detail according to several categories, including nitrogen fixation, nitrification, denitrification, anaerobic ammonium oxidation, and dissimilatory nitrate reduction to ammonium. In particular, we assess the environmental significance and potential applications of thermophilic nitrogen-cycling microorganisms, and highlight knowledge gaps and future research opportunities.

A moderately thermophilic origin of a novel family of marine group II euryarchaeota from deep ocean

Marine group II (MGII) is the most abundant planktonic heterotrophic archaea in the ocean. The evolutionary history of MGII archaea is elusive. In this study, 13 new MGII metagenome-assembled genomes were recovered from surface to the hadal zone in Challenger Deep of the Mariana Trench; four of them from the deep ocean represent a novel group. The optimal growth temperature (OGT) of the common ancestor of MGII has been estimated to be at about 60°C and OGTs of MGIIc, MGIIb, and MGIIa at 47°C-50ºC, 37°C-44ºC, and 30°C-37ºC, respectively, suggesting the adaptation of these species to different temperatures during evolution. The estimated OGT range of MGIIc was supported by experimental measurements of cloned β-galactosidase that showed optimal enzyme activity around 50°C. These results indicate that MGIIc may have originated from a common ancestor that lived in warm or even hot marine environment, such as hydrothermal vents.

Wood-Ljungdahl pathway found in novel marine Korarchaeota groups illuminates their evolutionary history

Korarchaeota, due to its rarity in common environments, is one of the archaeal phyla that has received the least attention from researchers. It was previously thought to consist solely of strict thermophiles. However, our study provides genetic evidence for the presence of korarchaeal members in temperate subsurface seawater. Furthermore, a systematic reclassification of the Korarchaeota based on 16S rRNA genes and genomes has revealed three novel marine groups (Kor-6 to Kor-8) at the root of the Korarchaeota branch. Kor-6 contains microbes that are present in moderate temperatures. All three novel marine phyla possess genes for the Wood-Ljungdahl pathway, and Kor-7 and Kor-8 possess fewer genes encoding oxygen resistance traits than other korarchaeal groups, suggesting a distinct lifestyle for these novel phyla. Our results, together with estimations of Korarchaeota divergence times, suggest that oxygen availability may be one of the important factors that have influenced the evolution of Korarchaeota. IMPORTANCE Korarchaeota were previously thought to inhabit exclusively high-temperature environments. However, our study provides genetic evidence for their unexpected presence in temperate marine waters. Through analysis of publicly available korarchaeal reference data, we have systematically reclassified Korarchaeota and identified the existence of three previously unknown marine groups (Kor-6, Kor-7, and Kor-8) at the root of the Korarchaeota branch. Comparative analysis of their gene content revealed that these novel groups exhibit a lifestyle distinct from other Korarchaeota. Specifically, they have the ability to fix carbon exclusively via the Wood-Ljungdahl (WL) pathway, and the genomes within Kor-7 and Kor-8 contain few genes encoding antioxidant enzymes, indicating their strictly anaerobic lifestyle. Further studies suggest that the genes related to methane metabolism and the WL pathway may have been inherited from a common ancestor of the Korarchaeota and that oxygen availability may be one of the important evolutionary factors that shaped the diversification of this archaeal phylum.

Taxonomic and carbon metabolic diversification of Bathyarchaeia during its coevolution history with early Earth surface environment

Bathyarchaeia, as one of the most abundant microorganisms on Earth, play vital roles in the global carbon cycle. However, our understanding of their origin, evolution, and ecological functions remains poorly constrained. Here, we present the largest dataset of Bathyarchaeia metagenome assembled genome to date and reclassify Bathyarchaeia into eight order-level units corresponding to the former subgroup system. Highly diversified and versatile carbon metabolisms were found among different orders, particularly atypical C1 metabolic pathways, indicating that Bathyarchaeia represent overlooked important methylotrophs. Molecular dating results indicate that Bathyarchaeia diverged at ~3.3 billion years, followed by three major diversifications at ~3.0, ~2.5, and ~1.8 to 1.7 billion years, likely driven by continental emergence, growth, and intensive submarine volcanism, respectively. The lignin-degrading Bathyarchaeia clade emerged at ~300 million years perhaps contributed to the sharply decreased carbon sequestration rate during the Late Carboniferous period. The evolutionary history of Bathyarchaeia potentially has been shaped by geological forces, which, in turn, affected Earth's surface environment.

Thermophilic Dehalococcoidia with unusual traits shed light on an unexpected past

Although the phylum Chloroflexota is ubiquitous, its biology and evolution are poorly understood due to limited cultivability. Here, we isolated two motile, thermophilic bacteria from hot spring sediments belonging to the genus Tepidiforma and class Dehalococcoidia within the phylum Chloroflexota. A combination of cryo-electron tomography, exometabolomics, and cultivation experiments using stable isotopes of carbon revealed three unusual traits: flagellar motility, a peptidoglycan-containing cell envelope, and heterotrophic activity on aromatics and plant-associated compounds. Outside of this genus, flagellar motility has not been observed in Chloroflexota, and peptidoglycan-containing cell envelopes have not been described in Dehalococcoidia. Although these traits are unusual among cultivated Chloroflexota and Dehalococcoidia, ancestral character state reconstructions showed flagellar motility and peptidoglycan-containing cell envelopes were ancestral within the Dehalococcoidia, and subsequently lost prior to a major adaptive radiation of Dehalococcoidia into marine environments. However, despite the predominantly vertical evolutionary histories of flagellar motility and peptidoglycan biosynthesis, the evolution of enzymes for degradation of aromatics and plant-associated compounds was predominantly horizontal and complex. Together, the presence of these unusual traits in Dehalococcoidia and their evolutionary histories raise new questions about the timing and selective forces driving their successful niche expansion into global oceans.

Fine-scale evaluation of two standard 16S rRNA gene amplicon primer pairs for analysis of total prokaryotes and archaeal nitrifiers in differently managed soils

The advance of high-throughput molecular biology tools allows in-depth profiling of microbial communities in soils, which possess a high diversity of prokaryotic microorganisms. Amplicon-based sequencing of 16S rRNA genes is the most common approach to studying the richness and composition of soil prokaryotes. To reliably detect different taxonomic lineages of microorganisms in a single soil sample, an adequate pipeline including DNA isolation, primer selection, PCR amplification, library preparation, DNA sequencing, and bioinformatic post-processing is required. Besides DNA sequencing quality and depth, the selection of PCR primers and PCR amplification reactions arguably have the largest influence on the results. This study tested the performance and potential bias of two primer pairs, i.e., 515F (Parada)-806R (Apprill) and 515F (Parada)-926R (Quince) in the standard pipelines of 16S rRNA gene Illumina amplicon sequencing protocol developed by the Earth Microbiome Project (EMP), against shotgun metagenome-based 16S rRNA gene reads. The evaluation was conducted using five differently managed soils. We observed a higher richness of soil total prokaryotes by using reverse primer 806R compared to 926R, contradicting to in silico evaluation results. Both primer pairs revealed various degrees of taxon-specific bias compared to metagenome-derived 16S rRNA gene reads. Nonetheless, we found consistent patterns of microbial community variation associated with different land uses, irrespective of primers used. Total microbial communities, as well as ammonia oxidizing archaea (AOA), the predominant ammonia oxidizers in these soils, shifted along with increased soil pH due to agricultural management. In the unmanaged low pH plot abundance of AOA was dominated by the acid-tolerant NS-Gamma clade, whereas limed agricultural plots were dominated by neutral-alkaliphilic NS-Delta/NS-Alpha clades. This study stresses how primer selection influences community composition and highlights the importance of primer selection for comparative and integrative studies, and that conclusions must be drawn with caution if data from different sequencing pipelines are to be compared.

Postglacial adaptations enabled colonization and quasi-clonal dispersal of ammonia-oxidizing archaea in modern European large lakes

Ammonia-oxidizing archaea (AOA) play a key role in the aquatic nitrogen cycle. Their genetic diversity is viewed as the outcome of evolutionary processes that shaped ancestral transition from terrestrial to marine habitats. However, current genome-wide insights into AOA evolution rarely consider brackish and freshwater representatives or provide their divergence timeline in lacustrine systems. An unbiased global assessment of lacustrine AOA diversity is critical for understanding their origins, dispersal mechanisms, and ecosystem roles. Here, we leveraged continental-scale metagenomics to document that AOA species diversity in freshwater systems is remarkably low compared to marine environments. We show that the uncultured freshwater AOA, "Candidatus Nitrosopumilus limneticus," is ubiquitous and genotypically static in various large European lakes where it evolved 13 million years ago. We find that extensive proteome remodeling was a key innovation for freshwater colonization of AOA. These findings reveal the genetic diversity and adaptive mechanisms of a keystone species that has survived clonally in lakes for millennia.

River damming enhances ecological functional stability of planktonic microorganisms

Planktonic microorganisms play an important role in maintaining the ecological functions in aquatic ecosystems, but how their structure and function interrelate and respond to environmental changes is still not very clear. Damming interrupts the river continuum and alters river nutrient biogeochemical cycling and biological succession. Considering that river damming decreases the irregular hydrological fluctuation, we hypothesized that it can enhance the ecological functional stability (EFS) of planktonic microorganisms. Therefore, the community composition of planktonic bacteria and archaea, functional genes related to carbon, nitrogen, sulfur, and phosphorus cycling, and relevant environmental factors of four cascade reservoirs in the Pearl River, Southern China, were investigated to understand the impact of damming on microbial community structure and function and verify the above hypothesis. Here, the ratio of function to taxa (F:T) based on Euclidean distance matrix analysis was first proposed to characterize the microbial EFS; the smaller the ratio, the more stable the ecological functions. The results showed that the reservoirs created by river damming had seasonal thermal and chemical stratifications with an increasing hydraulic retention time, which significantly changed the microbial structure and function. The river microbial F:T was significantly higher than that of the reservoirs, indicating that river damming enhances the EFS of the planktonic microorganisms. Structural equation modeling demonstrated that water temperature was an important factor influencing the relationship between the microbial structure and function and thus affected their EFS. In addition, reservoir hydraulic load was found a main factor regulating the seasonal difference in microbial EFS among the reservoirs. This study will help to deepen the understanding of the relationship between microbial structure and function and provide a theoretical basis of assessing the ecological function change after the construction of river damming.

Archaeal lipids trace ecology and evolution of marine ammonia-oxidizing archaea

Archaeal membrane lipids are widely used for paleotemperature reconstructions, yet these molecular fossils also bear rich information about ecology and evolution of marine ammonia-oxidizing archaea (AOA). Here we identified thermal and nonthermal behaviors of archaeal glycerol dialkyl glycerol tetraethers (GDGTs) by comparing the GDGT-based temperature index (TEX₈₆) to the ratio of GDGTs with two and three cyclopentane rings (GDGT-2/GDGT-3). Thermal-dependent biosynthesis should increase TEX₈₆ and decrease GDGT-2/GDGT-3 when the ambient temperature increases. This presumed temperature-dependent (PTD) trend is observed in GDGTs derived from cultures of thermophilic and mesophilic AOA. The distribution of GDGTs in suspended particulate matter (SPM) and sediments collected from above the pycnocline-shallow water samples-also follows the PTD trend. These similar GDGT distributions between AOA cultures and shallow water environmental samples reflect shallow ecotypes of marine AOA. While there are currently no cultures of deep AOA clades, GDGTs derived from deep water SPM and marine sediment samples exhibit nonthermal behavior deviating from the PTD trend. The presence of deep AOA increases the GDGT-2/GDGT-3 ratio and distorts the temperature-controlled correlation between GDGT-2/GDGT-3 and TEX₈₆. We then used Gaussian mixture models to statistically characterize these diagnostic patterns of modern AOA ecology from paleo-GDGT records to infer the evolution of marine AOA from the Mid-Mesozoic to the present. Long-term GDGT-2/GDGT-3 trends suggest a suppression of today's deep water marine AOA during the Mesozoic-early Cenozoic greenhouse climates. Our analysis provides invaluable insights into the evolutionary timeline and the expansion of AOA niches associated with major oceanographic and climate changes.

Candidatus Nitrosopolaris, a genus of putative ammonia-oxidizing archaea with a polar/alpine distribution

Ammonia-oxidizing archaea (AOA) are key players in the nitrogen cycle of polar soils. Here, we analyzed metagenomic data from tundra soils in Rásttigáisá, Norway, and recovered four metagenome-assembled genomes (MAGs) assigned to the genus 'UBA10452', an uncultured lineage of putative AOA in the order Nitrososphaerales ('terrestrial group I.1b'), phylum Thaumarchaeota. Analysis of other eight previously reported MAGs and publicly available amplicon sequencing data revealed that the UBA10452 lineage is predominantly found in acidic polar and alpine soils. In particular, UBA10452 MAGs were more abundant in highly oligotrophic environments such as mineral permafrost than in more nutrient-rich, vegetated tundra soils. UBA10452 MAGs harbour multiple copies of genes related to cold tolerance, particularly genes involved in DNA replication and repair. Based on the phylogenetic, biogeographic, and ecological characteristics of 12 UBA10452 MAGs, which include a high-quality MAG (90.8% complete, 3.9% redundant) with a nearly complete 16S rRNA gene, we propose a novel Candidatus genus, Ca. Nitrosopolaris, with four species representing clear biogeographic/habitat clusters.

Phylogenetic divergence and adaptation of Nitrososphaeria across lake depths and freshwater ecosystems

Thaumarchaeota (now the class Nitrososphaeria in the phylum Thermoproteota in GTDB taxonomy) are abundant across marine and soil habitats; however, their genomic diversity and evolutionary history in freshwater environments remain elusive. Here, we reconstructed 17 high-quality metagenome-assembled genomes of Nitrososphaeria from a deep lake and two great rivers, and compared all available genomes between freshwater and marine habitats regarding their phylogenetic positions, relative abundance, and genomic content. We found that freshwater Nitrososphaeria were dominated by the family Nitrosopumilaceae and could be grouped into three distinct clades closely related to the genera Nitrosopumilus, Nitrosoarchaeum, and Nitrosotenuis. The Nitrosopumilus-like clade was exclusively from deep lakes, while the Nitrosoarchaeum-like clade was dominated by species from deep lakes and rivers, and the Nitrosotenuis-like clade was mainly from rivers, deep lakes, and estuaries. Interestingly, there was vertical niche separation between two clades in deep lakes, showing that the Nitrosopumilus-like species dominated shallow layers, whereas the relative abundance of the Nitrosoarchaeum-like clade increased toward deep waters. Phylogenetic clustering patterns in the Nitrosopumilaceae supported at least one freshwater-to-marine and two marine-to-freshwater transitions, the former of which refined the potential terrestrial-to-marine evolutionary path as previously proposed. The occurrence of the two marine-to-freshwater transitions were accompanied by horizontal transfer of the genes involved in nutrition regulation, osmoregulation, and cell motility during their colonization to freshwater habitats. Specifically, the Nitrosopumilus-like clade showed losses of genes encoding flagella assembly and ion transport, whereas the Nitrosoarchaeum-like clade had losses of intact genes involved in urea uptake and utilization and gains of genes encoding osmolarity-mediated mechanosensitive channels. Collectively, our results reveal for the first time the high genomic diversity of the class Nitrososphaeria across freshwater ecosystems and provide novel insights into their adaptive mechanisms and evolutionary histories.

Archaeal nitrification without oxygen

The single-cell organism can self-produce oxygen for ammonia oxidation.

A Novel Approach to Characterize the Lipidome of Marine Archaeon Nitrosopumilus maritimus by Ion Mobility Mass Spectrometry

Archaea are differentiated from the other two domains of life by their biomolecular characteristics. One such characteristic is the unique structure and composition of their lipids. Characterization of the whole set of lipids in a biological system (the lipidome) remains technologically challenging. This is because the lipidome is innately complex, and not all lipid species are extractable, separable, or ionizable by a single analytical method. Furthermore, lipids are structurally and chemically diverse. Many lipids are isobaric or isomeric and often indistinguishable by the measurement of mass or even their fragmentation spectra. Here we developed a novel analytical protocol based on liquid chromatography ion mobility mass spectrometry to enhance the coverage of the lipidome and characterize the conformations of archaeal lipids by their collision cross-sections (CCSs). The measurements of ion mobility revealed the gas-phase ion chemistry of representative archaeal lipids and provided further insights into their attributions to the adaptability of archaea to environmental stresses. A comprehensive characterization of the lipidome of mesophilic marine thaumarchaeon, Nitrosopumilus maritimus (strain SCM1) revealed potentially an unreported phosphate- and sulfate-containing lipid candidate by negative ionization analysis. It was the first time that experimentally derived CCS values of archaeal lipids were reported. Discrimination of crenarchaeol and its proposed stereoisomer was, however, not achieved with the resolving power of the SYNAPT G2 ion mobility system, and a high-resolution ion mobility system may be required for future work. Structural and spectral libraries of archaeal lipids were constructed in non-vendor-specific formats and are being made available to the community to promote research of Archaea by lipidomics.