Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution

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.

Comprehensive molecular-isotopic characterization of archaeal lipids in the Black Sea water column and underlying sediments

The Black Sea is a permanently anoxic, marine basin serving as model system for the deposition of organic-rich sediments in a highly stratified ocean. In such systems, archaeal lipids are widely used as paleoceanographic and biogeochemical proxies; however, the diverse planktonic and benthic sources as well as their potentially distinct diagenetic fate may complicate their application. To track the flux of archaeal lipids and to constrain their sources and turnover, we quantitatively examined the distributions and stable carbon isotopic compositions (δ13 C) of intact polar lipids (IPLs) and core lipids (CLs) from the upper oxic water column into the underlying sediments, reaching deposits from the last glacial. The distribution of IPLs responded more sensitively to the geochemical zonation than the CLs, with the latter being governed by the deposition from the chemocline. The isotopic composition of archaeal lipids indicates CLs and IPLs in the deep anoxic water column have negligible influence on the sedimentary pool. Archaeol substitutes tetraether lipids as the most abundant IPL in the deep anoxic water column and the lacustrine methanic zone. Its elevated IPL/CL ratios and negative δ13 C values indicate active methane metabolism. Sedimentary CL- and IPL-crenarchaeol were exclusively derived from the water column, as indicated by non-variable δ13 C values that are identical to those in the chemocline and by the low BIT (branched isoprenoid tetraether index). By contrast, in situ production accounts on average for 22% of the sedimentary IPL-GDGT-0 (glycerol dibiphytanyl glycerol tetraether) based on isotopic mass balance using the fermentation product lactate as an endmember for the dissolved substrate pool. Despite the structural similarity, glycosidic crenarchaeol appears to be more recalcitrant in comparison to its non-cycloalkylated counterpart GDGT-0, as indicated by its consistently higher IPL/CL ratio in sediments. The higher TEX86 , CCaT, and GDGT-2/-3 values in glacial sediments could plausibly result from selective turnover of archaeal lipids and/or an archaeal ecology shift during the transition from the glacial lacustrine to the Holocene marine setting. Our in-depth molecular-isotopic examination of archaeal core and intact polar lipids provided new constraints on the sources and fate of archaeal lipids and their applicability in paleoceanographic and biogeochemical studies.

Comparative genomics reveals the adaptation of ammonia-oxidising Thaumarchaeota to arid soils

Thaumarchaeota are predominant in oligotrophic habitats such as deserts and arid soils, but their adaptations to these arid conditions are not well understood. In this study, we assembled 23 Thaumarchaeota genomes from arid and semi-arid soils collected from the Inner Mongolia Steppe and the Qinghai-Tibet Plateau. Using a comparative genomics approach, integrated with 614 Thaumarchaeota genomes from public databases, we identified the traits and evolutionary forces that contribute to their adaptations to aridity. Our results showed that the newly assembled genomes represent an early diverging group within the lineage of ammonia-oxidising Thaumarchaeota. While the genomic functions previously identified in arid soil lineages were conserved across terrestrial, shallow-ocean and deep-ocean lineages, several traits likely contribute to Thaumarchaeota's adaptation to aridity. These include chlorite dismutase, arsenate reductase, V-type ATPase and genes dealing with oxidative stresses. The acquisition and loss of traits at the last common ancestor of arid soil lineages may have facilitated the specialisation of Thaumarchaeota in arid soils. Additionally, the acquisition of unique adaptive traits, such as a urea transporter, Ca2+ :H+ antiporter, mannosyl-3-phosphoglycerate synthase and phosphatase, DNA end-binding protein Ku and phage shock protein A, further distinguishes arid soil Thaumarchaeota. This study provides evidence for the adaptations of Thaumarchaeota to arid soil, enhancing our understanding of the nitrogen and carbon cycling driven by Thaumarchaeota in drylands.

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.

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.

Diversity of Bathyarchaeia viruses in metagenomes and virus-encoded CRISPR system components

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria, and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family "Fuxiviridae" harbor an atypical Type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family "Chiyouviridae" encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibit host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms.

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.

Influence of hydrological parameters on hydroxylated tetraether lipids in a deep Lake Fuxian, China: Implications for their use as environmental proxies

Hydroxylated isoprenoid glycerol dialkyl glycerol tetraethers (OH-GDGTs) have shown their potential in environmental reconstructions. However, the unclear underlying mechanism challenges their application. To elucidate the effects of water parameters on OH-GDGT-derived indices and understand their environmental implications, we investigated the core OH-GDGTs of suspended particulate matter (SPM) from water columns in a year cycle and surface sediments at different water depths along a nearshore-offshore transect in Lake Fuxian, a deep and large lake in southwestern China. OH-GDGTs were primarily found in the hypolimnion and were produced in situ by Group I.1a Thaumarchaeota. The relative abundance of OH-GDGTs (%OH-GDGTs) and ring indices (RI-OH and RI-OH') in the hypolimnion were significantly influenced by dissolved oxygen (DO) and pH, particularly DO, which regulated the inverse physiological functions of the hydroxyl and cyclopentane moieties of archaea. %OH-GDGTs values in SPM were positively correlated with DO and negatively correlated with pH levels, while RI-OH values exhibited an inverse relationship with DO and positive correlation with pH levels. OH-GDGTs in surface sediments appeared to be homologous to that of water columns, indicating that their inferred proxies could be regulated by the configuration of water parameters. The sedimentary %OH-GDGTs values increased as the RI-OH values decreased with water depth along the transect from the lakeshore to the lake center, suggesting their potential as lake-level proxies.

Habitat Transition in the Evolution of Bacteria and Archaea

Related groups of microbes are widely distributed across Earth's habitats, implying numerous dispersal and adaptation events over evolutionary time. However, relatively little is known about the characteristics and mechanisms of these habitat transitions, particularly for populations that reside in animal microbiomes. Here, we review the literature concerning habitat transitions among a variety of bacterial and archaeal lineages, considering the frequency of migration events, potential environmental barriers, and mechanisms of adaptation to new physicochemical conditions, including the modification of protein inventories and other genomic characteristics. Cells dependent on microbial hosts, particularly bacteria from the Candidate Phyla Radiation, have undergone repeated habitat transitions from environmental sources into animal microbiomes. We compare their trajectories to those of both free-living cells-including the Melainabacteria, Elusimicrobia, and methanogenic archaea-and cellular endosymbionts and bacteriophages, which have made similar transitions. We conclude by highlighting major related topics that may be worthy of future study.

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.

Diversity of Bathyarchaeia viruses in metagenomes and virus-encoded CRISPR system components

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. The virome of Bathyarchaeia so far has not been characterized. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria, and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family 'Fuxiviridae' harbor an atypical type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family 'Chiyouviridae' encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibition of host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms.

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.

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.

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.

Metatranscriptomes reveal the diverse responses of Thaumarchaeota ecotypes to environmental variations in the northern slope of the South China Sea

Thaumarchaeota are among the most abundant prokaryotes in the ocean, playing important roles in carbon and nitrogen cycling. Marine Thaumarchaeota ecotypes exhibit depth-related diversification and seasonal changes. However, transcriptomic activities concerning niche partitioning among thaumarchaeal ecotypes remain unclear. Here, we examined the variations in the distribution and transcriptomic activity of marine Thaumarchaeota ecotypes. Three primary ecotypes were identified: a Nitrosopumilus-like clade; a Nitrosopelagicus-like water column A (WCA) clade, thriving in epipelagic water; and a water column B (WCB) clade, dominant in deep water. Depth-related partitioning of the three ecotypes and the seasonal variability of the WCA and WCB ecotypes were observed. Nutrient concentrations, chlorophyll α and salinity were the primary environmental factors. The relative abundance of the WCA ecotype and its transcript abundance of amoA gene were positively correlated with chlorophyll α and salinity, while the WCB ecotype was positively correlated with nitrate and phosphate. Based on high-quality metagenome-assembled genomes, transcriptomic analysis revealed that the three ecotypes exhibited various co-occurring expression patterns of the elemental cycling genes in the nitrogen, carbon, phosphorus, and sulfur cycles. Our results provide transcriptomic evidence of the niche differentiation of marine Thaumarchaeota ecotypes, highlighting the diverse roles of ecotypes and WCA subclades in biogeochemical cycles.

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.

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.

Distinct Features of Sedimentary Archaeal Communities in Hypoxia and Non-Hypoxia Regions off the Changjiang River Estuary

Water hypoxia (DO < 2 mg/L) is a growing global environmental concern that has the potential to significantly influence not only the aquatic ecosystem but also the benthic sedimentary ecosystem. The Changjiang River Estuary hypoxia, classified as one of the world's largest seasonal hypoxic water basins, has been reported to be expanding rapidly in recent decades. However, the microbial community dynamics and responses to this water hypoxia are still unclear. In this study, we examined the abundance, community composition, and distribution of sedimentary archaea, one important component of microbial communities in the Changjiang River Estuary and the East China Sea (ECS). Our results indicated that Thaumarchaeota and Bathyarchaeota were predominant archaeal groups in these research areas, with their 16S rRNA gene abundance ranged from 8.55 × 106 to 7.51 × 108 and 3.18 × 105 to 1.11 × 108 copies/g, respectively. The sedimentary archaeal community was mainly influenced by DO, together with the concentration of ammonium, nitrate, and sulfide. In addition, distinct differences in the archaeal community's composition, abundance, and driving factors were discovered between samples from hypoxia and non-hypoxia stations. Furtherly, microbial networks suggest various microbes leading the different activities in hypoxic and normoxic environments. Bathyarchaeota and Thermoprofundales were "key stone" archaeal members of the low-DO network, whereas Thaumarchaeota constituted a significant component of the high-DO network. Our results provide a clear picture of the sedimentary archaeal community in coastal hypoxia zones and indicates potential distinctions of archaea in hypoxia and non-hypoxia environments, including ecological niches and metabolic functions. IMPORTANCE In this study, the sedimentary archaeal community composition and abundance were detailed revealed and quantified based on 16S rRNA genes off the Changjiang River Estuary. We found that the community composition was distinct between hypoxia and non-hypoxia regions, while Thaumarchaeota and Bathyarchaeota dominated in non-hypoxia and hypoxia samples, respectively. In hypoxia regions, the sedimentary archaea were mainly affected by salinity, ammonium, and nitrate, whereas total organic carbon, total nitrogen, and sulfide were major influencing factors in non-hypoxia regions. The distinct microbial network may suggest the niche difference of archaeal community under various oxygen level.

An Ancient Respiratory System in the Widespread Sedimentary Archaea Thermoprofundales

Thermoprofundales, formerly Marine Benthic Group D (MBG-D), is a ubiquitous archaeal lineage found in sedimentary environments worldwide. However, its taxonomic classification, metabolic pathways, and evolutionary history are largely unexplored because of its uncultivability and limited number of sequenced genomes. In this study, phylogenomic analysis and average amino acid identity values of a collection of 146 Thermoprofundales genomes revealed five Thermoprofundales subgroups (A-E) with distinct habitat preferences. Most of the microorganisms from Subgroups B and D were thermophiles inhabiting hydrothermal vents and hot spring sediments, whereas those from Subgroup E were adapted to surface environments where sunlight is available. H2 production may be featured in Thermoprofundales as evidenced by a gene cluster encoding the ancient membrane-bound hydrogenase (MBH) complex. Interestingly, a unique structure separating the MBH gene cluster into two modular units was observed exclusively in the genomes of Subgroup E, which included a peripheral arm encoding the [NiFe] hydrogenase domain and a membrane arm encoding the Na+/H+ antiporter domain. These two modular structures were confirmed to function independently by detecting the H2-evolving activity in vitro and salt tolerance to 0.2 M NaCl in vivo, respectively. The peripheral arm of Subgroup E resembles the proposed common ancestral respiratory complex of modern respiratory systems, which plays a key role in the early evolution of life. In addition, molecular dating analysis revealed that Thermoprofundales is an early emerging archaeal lineage among the extant MBH-containing microorganisms, indicating new insights into the evolution of this ubiquitous archaea lineage.

Linking microbiomes with per- and poly-fluoroalkyl substances (PFASs) in soil ecosystems: Microbial community assembly, stability, and trophic phylosymbiosis

Microbiomes are vital in promoting nutrient cycling and plant growth in soil ecosystems. However, microbiomes face adverse effects from multiple persistent pollutants, including per- and poly-fluoroalkyl substances (PFASs). PFASs threaten the fertility and health of soil ecosystems, yet the response of microbial community stability and trophic transfer efficiencies to PFASs is still poorly understood. This study explored the spatial patterns of PFASs in topsoil environments from the West Taihu Lake Basin of China and links their presence to soil microbial community stability at compositional and functional levels. Our results revealed that PFBA (13.87%), PFTrDA (11.63%), PFDoA (11.02%), PFOA (10.99%), and PFOS (10.39%) contributed the most to the spatial occurrence of PFASs. Soil properties, including salinity (14.47%), uniformity (9.68%), dissolved inorganic carbon (8.62%), and clay content (8.18%), affected PFASs distribution the most. In soil microbiomes, eukaryotic taxa had wider niche breadths and stronger community stability than prokaryotes when exposed to PFASs (p < 0.05). The presence of PFBA and PFHpA inhibited the functional stability of archaeal and bacterial communities (p < 0.05). PFBA and PFPeA reduced the structural stability of heterotrophic bacteria and Myxobacteria, respectively (p < 0.05). Based on null modeling, PFPeA significantly regulated the assembly processes of most microbial sub-communities (p < 0.01). The trophic transferring efficiencies of autotrophic bacteria to metazoan organisms were directly stimulated by PFASs (p < 0.05), and the potential trophic transferring efficiencies of methanogenic archaea to protozoa were inhibited by PFASs (p < 0.05). This study highlighted the potential contributions of PFASs to soil microbial community stability and food webs during ecological soil management.

Phylogenomics uncovers evolutionary trajectory of nitrogen fixation in Cyanobacteria

Biological nitrogen fixation (BNF) by cyanobacteria is of significant importance for the Earth’s biogeochemical nitrogen cycle but is restricted to a few genera that do not form monophyletic group. To explore the evolutionary trajectory of BNF and investigate the driving forces of its evolution, we analyze 650 cyanobacterial genomes and compile the database of diazotrophic cyanobacteria based on the presence of nitrogen fixation gene clusters (NFGCs). We report that 266 of 650 examined genomes are NFGC-carrying members, and these potentially diazotrophic cyanobacteria are unevenly distributed across the phylogeny of Cyanobacteria, that multiple independent losses shaped the scattered distribution. Among the diazotrophic cyanobacteria, two types of NFGC exist, with one being ancestral and abundant, which have descended from diazotrophic ancestors, and the other being anaerobe-like and sparse, possibly being acquired from anaerobic microbes through horizontal gene transfer. Interestingly, we illustrate that the origin of BNF in Cyanobacteria coincide with two major evolutionary events. One is the origin of multicellularity of cyanobacteria, and the other is concurrent genetic innovations with massive gene gains and expansions, implicating their key roles in triggering the evolutionary transition from nondiazotrophic to diazotrophic cyanobacteria. Additionally, we reveal that genes involved in accelerating respiratory electron transport (coxABC), anoxygenic photosynthetic electron transport (sqr), as well as anaerobic metabolisms (pfor, hemN, nrdG, adhE) are enriched in diazotrophic cyanobacteria, representing adaptive genetic signatures that underpin the diazotrophic lifestyle. Collectively, our study suggests that multicellularity, together with concurrent genetic adaptations contribute to the evolution of diazotrophic cyanobacteria.

Phylogenomic evidence for the Origin of Obligately Anaerobic Anammox Bacteria around the Great Oxidation Event

The anaerobic ammonium oxidation (anammox) bacteria can transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remain enigmatic. Here, we performed a comprehensive phylogenomic and molecular clock analysis of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which include implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based molecular dating approach, we estimated a consistent origin of anammox bacteria at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32–2.5 Ga) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress adaptation, bioenergetics, and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria.

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.

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.

Phylogenomic Analysis of Metagenome-Assembled Genomes Deciphered Novel Acetogenic Nitrogen-Fixing Bathyarchaeota from Hot Spring Sediments

This study describes the phylogenomic analysis and metabolic insights of metagenome-assembled genomes (MAGs) retrieved from hot spring sediment samples. The metagenome-assembled sequences recovered three near-complete genomes belonging to the archaeal phylum. Analysis of genome-wide core genes and 16S rRNA-based phylogeny placed the ILS200 and ILS300 genomes within the uncultivated and largely understudied bathyarchaeal phylum, whereas ILS100 represented the phylum Thaumarchaeota. The average nucleotide identity (ANI) of the bin ILS100 was 76% with Nitrososphaeria_archaeon_isolate_SpSt-1069. However, the bins ILS200 and ILS300 showed ANI values of 75% and 70% with Candidatus_Bathyarchaeota_archaeon_isolate_DRTY-6_2_bin_115 and Candidatus_Bathyarchaeota_archaeon_BA1_ba1_01, respectively. The genomic potential of Bathyarchaeota bins ILS200 and ILS300 showed genes necessary for the Wood-Ljungdahl pathway, and the gene encoding the methyl coenzyme M reductase (mcr) complex essential for methanogenesis was absent. The metabolic potential of the assembled genomes included genes involved in nitrogen assimilation, including nitrogenase and the genes necessary for the urea cycle. The presence of these genes suggested the metabolic potential of Bathyarchaeota to fix nitrogen under extreme environments. In addition, the ILS200 and ILS300 genomes carried genes involved in the tricarboxylic acid (TCA) cycle, glycolysis, and degradation of organic carbons. Finally, we conclude that the reconstructed Bathyarchaeota bins are autotrophic acetogens and organo-heterotrophs.

IMPORTANCE We describe the Bathyarchaeota bins that are likely to be acetogens with a wide range of metabolic potential. These bins did not exhibit methanogenic machinery, suggesting methane production may not occur by all subgroup lineages of Bathyarchaeota. Phylogenetic analysis support that both ILS200 and ILS300 belonged to the Bathyarchaeota. The discovery of new bathyarchaeotal MAGs provides additional knowledge for understanding global carbon and nitrogen metabolism under extreme conditions.

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.

Microbial Activities and Selection from Surface Ocean to Subseafloor on the Namibian Continental Shelf

Oxygen minimum zones (OMZs) are hot spots for redox-sensitive nitrogen transformations fueled by sinking organic matter. In comparison, the regulating role of sulfur-cycling microbes in marine OMZs, their impact on carbon cycling in pelagic and benthic habitats, and activities below the seafloor remain poorly understood. Using ¹³C DNA stable isotope probing (SIP) and metatranscriptomics, we explored microbial guilds involved in sulfur and carbon cycling from the ocean surface to the subseafloor on the Namibian shelf. There was a clear separation in microbial community structure across the seawater-seafloor boundary, which coincided with a 100-fold-increased concentration of microbial biomass and unique gene expression profiles of the benthic communities. ¹³C-labeled 16S rRNA genes in SIP experiments revealed carbon-assimilating taxa and their distribution across the sediment-water interface. Most of the transcriptionally active taxa among water column communities that assimilated ¹³C from diatom exopolysaccharides (mostly Bacteroidetes, Actinobacteria, Alphaproteobacteria, and Planctomycetes) also assimilated ¹³C-bicarbonate under anoxic conditions in sediment incubations. Moreover, many transcriptionally active taxa from the seafloor community (mostly sulfate-reducing Deltaproteobacteria and sulfide-oxidizing Gammaproteobacteria) that assimilated ¹³C-bicarbonate under sediment anoxic conditions also assimilated ¹³C from diatom exopolysaccharides in the surface ocean and OMZ waters. Despite strong selection at the sediment-water interface, many taxa related to either planktonic or benthic communities were found to be present at low abundance and actively assimilating carbon under both sediment and water column conditions. In austral winter, mixing of shelf waters reduces stratification and suspends sediments from the seafloor into the water column, potentially spreading metabolically versatile microbes across niches.

IMPORTANCE Microbial activities in oxygen minimum zones (OMZs) transform inorganic fixed nitrogen into greenhouse gases, impacting the Earth’s climate and nutrient equilibrium. Coastal OMZs are predicted to expand with global change and increase carbon sedimentation to the seafloor. However, the role of sulfur-cycling microbes in assimilating carbon in marine OMZs and related seabed habitats remain poorly understood. Using ¹³C DNA stable isotope probing and metatranscriptomics, we explore microbial guilds involved in sulfur and carbon cycling from ocean surface to subseafloor on the Namibian shelf. Despite strong selection and differential activities across the sediment-water interface, many active taxa were identified in both planktonic and benthic communities, either fixing inorganic carbon or assimilating organic carbon from algal biomass. Our data show that many planktonic and benthic microbes linked to the sulfur cycle can cross redox boundaries when mixing of the shelf waters reduces stratification and suspends seafloor sediment particles into the water column.

Diverse ecophysiological adaptations of subsurface Thaumarchaeota in floodplain sediments revealed through genome-resolved metagenomics

The terrestrial subsurface microbiome contains vastly underexplored phylogenetic diversity and metabolic novelty, with critical implications for global biogeochemical cycling. Among the key microbial inhabitants of subsurface soils and sediments are Thaumarchaeota, an archaeal phylum that encompasses ammonia-oxidizing archaea (AOA) as well as non-ammonia-oxidizing basal lineages. Thaumarchaeal ecology in terrestrial systems has been extensively characterized, particularly in the case of AOA. However, there is little knowledge on the diversity and ecophysiology of Thaumarchaeota in deeper soils, as most lineages, particularly basal groups, remain uncultivated and underexplored. Here we use genome-resolved metagenomics to examine the phylogenetic and metabolic diversity of Thaumarchaeota along a 234 cm depth profile of hydrologically variable riparian floodplain sediments in the Wind River Basin near Riverton, Wyoming. Phylogenomic analysis of the metagenome-assembled genomes (MAGs) indicates a shift in AOA population structure from the dominance of the terrestrial Nitrososphaerales lineage in the well-drained top ~100 cm of the profile to the typically marine Nitrosopumilales in deeper, moister, more energy-limited sediment layers. We also describe two deeply rooting non-AOA MAGs with numerous unexpected metabolic features, including the reductive acetyl-CoA (Wood-Ljungdahl) pathway, tetrathionate respiration, a form III RuBisCO, and the potential for extracellular electron transfer. These MAGs also harbor tungsten-containing aldehyde:ferredoxin oxidoreductase, group 4f [NiFe]-hydrogenases and a canonical heme catalase, typically not found in Thaumarchaeota. Our results suggest that hydrological variables, particularly proximity to the water table, impart a strong control on the ecophysiology of Thaumarchaeota in alluvial sediments.

High Abundance of Thaumarchaeota Found in Deep Metamorphic Subsurface in Eastern China

Members of the Thaumarchaeota phylum play a key role in nitrogen cycling and are prevalent in a variety of environments including soil, sediment, and seawater. However, few studies have shown the presence of Thaumarchaeota in the terrestrial deep subsurface. Using high-throughput 16S rRNA gene sequencing, this study presents evidence for the high relative abundance of Thaumarchaeota in a biofilm sample collected from the well of Chinese Continental Scientific Drilling at a depth of 2000 m. Phylogenetic analysis showed a close relationship of these thaumarchaeotal sequences with known ammonia-oxidizing archaea (AOA) isolates, suggesting the presence of AOA in the deep metamorphic environment of eastern China which is believed to be oxic. Based on fluid geochemistry and FAProTax functional prediction, a pathway of nitrogen cycling is proposed. Firstly, heterotrophic nitrogen fixation is executed by diazotrophic bacteria coupled with methane oxidation. Then, ammonia is oxidized to nitrite by AOA, and nitrite is further oxidized to nitrate by bacteria within the phylum Nitrospirae. Denitrification and anaerobic ammonia oxidation occur slowly, leading to nitrate accumulation in the subsurface. With respect to biogeochemistry, the reaction between downward diffusing O2 and upward diffusing CH4 potentially fuels the ecosystem with a high relative abundance of Thaumarchaeota.

Genome-Resolved Metagenomic Insights into Massive Seasonal Ammonia-Oxidizing Archaea Blooms in San Francisco Bay

Ammonia-oxidizing archaea (AOA) are key for the transformation of ammonia to oxidized forms of nitrogen in aquatic environments around the globe, including nutrient-rich coastal and estuarine waters such as San Francisco Bay (SFB). Using metagenomics and 16S rRNA gene amplicon libraries, we found that AOA are more abundant than ammonia-oxidizing bacteria (AOB) and nitrite-oxidizing bacteria (NOB), except in the freshwater stations in SFB. In South SFB, we observed recurrent AOA blooms of “Candidatus Nitrosomarinus catalina” SPOT01-like organisms, which account for over 20% of 16S rRNA gene amplicons in both surface and bottom waters and co-occur with weeks of high nitrite concentrations (>10 μM) in the oxic water column. We observed pronounced nitrite peaks occurring in the autumn for 7 of the last 9 years (2012 to 2020), suggesting that seasonal AOA blooms are common in South SFB. We recovered two high-quality AOA metagenome-assembled genomes (MAGs), including a Nitrosomarinus-like genome from the South SFB bloom and another Nitrosopumilus genome originating from Suisun Bay in North SFB. Both MAGs cluster with genomes from other estuarine/coastal sites. Analysis of Nitrosomarinus-like genomes show that they are streamlined, with low GC content and high coding density, and harbor urease genes. Our findings support the unique niche of Nitrosomarinus-like organisms which dominate coastal/estuarine waters and provide insights into recurring AOA blooms in SFB.

IMPORTANCE Ammonia-oxidizing archaea (AOA) carry out key transformations of ammonia in estuarine systems such as San Francisco Bay (SFB)—the largest estuary on the west coast of North America—and play a significant role in both local and global nitrogen cycling. Using metagenomics and 16S rRNA gene amplicon libraries, we document a massive, recurrent AOA bloom in South SFB that co-occurs with months of high nitrite concentrations in the oxic water column. Our study is the first to generate metagenome-assembled genomes (MAGs) from SFB, and through this process we recovered two high-quality AOA MAGs, one of which originated from bloom samples. These AOA MAGs yield new insight into the Nitrosopumilus and Nitrosomarinus-like lineages and their potential niches in coastal and estuarine systems. Nitrosomarinus-like AOA are abundant in coastal regions around the globe, and we highlight the common occurrence of urease genes, low GC content, and range of salinity tolerances within this lineage.

Archaeal nitrification without oxygen

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

Progress and Challenges in Studying the Ecophysiology of Archaea

It has been less than two decades since the study of archaeal ecophysiology has become unshackled from the limitations of cultivation and amplicon sequencing through the advent of metagenomics. As a primer to the guide on producing archaeal genomes from metagenomes, we briefly summarize here how different meta'omics, imaging, and wet lab methods have contributed to progress in understanding the ecophysiology of Archaea. We then peer into the history of how our knowledge on two particularly important lineages was assembled: the anaerobic methane and alkane oxidizers, encountered primarily among Euryarchaeota, and the nanosized, mainly parasitic, members of the DPANN superphylum.

Relative time constraints improve molecular dating

Dating the tree of life is central to understanding the evolution of life on Earth. Molecular clocks calibrated with fossils represent the state of the art for inferring the ages of major groups. Yet, other information on the timing of species diversification can be used to date the tree of life. For example, horizontal gene transfer events and ancient coevolutionary interactions such as (endo)symbioses occur between contemporaneous species and thus can imply temporal relationships between two nodes in a phylogeny. Temporal constraints from these alternative sources can be particularly helpful when the geological record is sparse, for example, for microorganisms, which represent the majority of extant and extinct biodiversity. Here, we present a new method to combine fossil calibrations and relative age constraints to estimate chronograms. We provide an implementation of relative age constraints in RevBayes that can be combined in a modular manner with the wide range of molecular dating methods available in the software. We use both realistic simulations and empirical datasets of 40 Cyanobacteria and 62 Archaea to evaluate our method. We show that the combination of relative age constraints with fossil calibrations significantly improves the estimation of node ages. [Archaea, Bayesian analysis, cyanobacteria, dating, endosymbiosis, lateral gene transfer, MCMC, molecular clock, phylogenetic dating, relaxed molecular clock, revbayes, tree of life.]

Cofactor F420: an expanded view of its distribution, biosynthesis and roles in bacteria and archaea

Many bacteria and archaea produce the redox cofactor F420. F420 is structurally similar to the cofactors FAD and FMN but is catalytically more similar to NAD and NADP. These properties allow F420 to catalyze challenging redox reactions, including key steps in methanogenesis, antibiotic biosynthesis and xenobiotic biodegradation. In the last 5 years, there has been much progress in understanding its distribution, biosynthesis, role and applications. Whereas F420 was previously thought to be confined to Actinobacteria and Euryarchaeota, new evidence indicates it is synthesized across the bacterial and archaeal domains, as a result of extensive horizontal and vertical biosynthetic gene transfer. F420 was thought to be synthesized through one biosynthetic pathway; however, recent advances have revealed variants of this pathway and have resolved their key biosynthetic steps. In parallel, new F420-dependent biosynthetic and metabolic processes have been discovered. These advances have enabled the heterologous production of F420 and identified enantioselective F420H2-dependent reductases for biocatalysis. New research has also helped resolve how microorganisms use F420 to influence human and environmental health, providing opportunities for tuberculosis treatment and methane mitigation. A total of 50 years since its discovery, multiple paradigms associated with F420 have shifted, and new F420-dependent organisms and processes continue to be discovered.

Genomes of Thaumarchaeota from deep sea sediments reveal specific adaptations of three independently evolved lineages

Marine sediments represent a vast habitat for complex microbiomes. Among these, ammonia oxidizing archaea (AOA) of the phylum Thaumarchaeota are one of the most common, yet little explored, inhabitants, which seem extraordinarily well adapted to the harsh conditions of the subsurface biosphere. We present 11 metagenome-assembled genomes of the most abundant AOA clades from sediment cores obtained from the Atlantic Mid-Ocean ridge flanks and Pacific abyssal plains. Their phylogenomic placement reveals three independently evolved clades within the order Nitrosopumilales, of which no cultured representative is known yet. In addition to the gene sets for ammonia oxidation and carbon fixation known from other AOA, all genomes encode an extended capacity for the conversion of fermentation products that can be channeled into the central carbon metabolism, as well as uptake of amino acids probably for protein maintenance or as an ammonia source. Two lineages encode an additional (V-type) ATPase and a large repertoire of DNA repair systems that may allow to overcome the challenges of high hydrostatic pressure. We suggest that the adaptive radiation of AOA into marine sediments occurred more than once in evolution and resulted in three distinct lineages with particular adaptations to this extremely energy-limiting and high-pressure environment.

Fine-scale metabolic discontinuity in a stratified prokaryote microbiome of a Red Sea deep halocline

Deep-sea hypersaline anoxic basins are polyextreme environments in the ocean's interior characterized by the high density of brines that prevents mixing with the overlaying seawater, generating sharp chemoclines and redoxclines up to tens of meters thick that host a high concentration of microbial communities. Yet, a fundamental understanding of how such pycnoclines shape microbial life and the associated biogeochemical processes at a fine scale, remains elusive. Here, we applied high-precision sampling of the brine-seawater transition interface in the Suakin Deep, located at 2770 m in the central Red Sea, to reveal previously undocumented fine-scale community structuring and succession of metabolic groups along a salinity gradient only 1 m thick. Metagenomic profiling at a 10-cm-scale resolution highlighted spatial organization of key metabolic pathways and corresponding microbial functional units, emphasizing the prominent role and significance of salinity and oxygen in shaping their ecology. Nitrogen cycling processes are especially affected by the redoxcline with ammonia oxidation processes being taxa and layers specific, highlighting also the presence of novel microorganisms, such as novel Thaumarchaeota and anammox, adapted to the changing conditions of the chemocline. The findings render the transition zone as a critical niche for nitrogen cycling, with complementary metabolic networks, in turn underscoring the biogeochemical complexity of deep-sea brines.

(S)-3-Hydroxybutyryl-CoA Dehydrogenase From the Autotrophic 3-Hydroxypropionate/4-Hydroxybutyrate Cycle in Nitrosopumilus maritimus

Ammonia-oxidizing archaea of the phylum Thaumarchaeota are among the most abundant organisms that exert primary control of oceanic and soil nitrification and are responsible for a large part of dark ocean primary production. They assimilate inorganic carbon via an energetically efficient version of the 3-hydroxypropionate/4-hydroxybutyrate cycle. In this cycle, acetyl-CoA is carboxylated to succinyl-CoA, which is then converted to two acetyl-CoA molecules with 4-hydroxybutyrate as the key intermediate. This conversion includes the (S)-3-hydroxybutyryl-CoA dehydrogenase reaction. Here, we heterologously produced the protein Nmar_1028 catalyzing this reaction in thaumarchaeon Nitrosopumilus maritimus, characterized it biochemically and performed its phylogenetic analysis. This NAD-dependent dehydrogenase is highly active with its substrate, (S)-3-hydroxybutyryl-CoA, and its low K_m value suggests that the protein is adapted to the functioning in the 3-hydroxypropionate/4-hydroxybutyrate cycle. Nmar_1028 is homologous to the dehydrogenase domain of crotonyl-CoA hydratase/(S)-3-hydroxybutyryl-CoA dehydrogenase that is present in many Archaea. Apparently, the loss of the dehydratase domain of the fusion protein in the course of evolution was accompanied by lateral gene transfer of 3-hydroxypropionyl-CoA dehydratase/crotonyl-CoA hydratase from Bacteria. Although (S)-3-hydroxybutyryl-CoA dehydrogenase studied here is neither unique nor characteristic for the HP/HB cycle, Nmar_1028 appears to be the only (S)-3-hydroxybutyryl-CoA dehydrogenase in N. maritimus and is thus essential for the functioning of the 3-hydroxypropionate/4-hydroxybutyrate cycle and for the biology of this important marine archaeon.

The ups and downs of ectoine: structural enzymology of a major microbial stress protectant and versatile nutrient

Ectoine and its derivative 5-hydroxyectoine are compatible solutes and chemical chaperones widely synthesized by Bacteria and some Archaea as cytoprotectants during osmotic stress and high- or low-growth temperature extremes. The function-preserving attributes of ectoines led to numerous biotechnological and biomedical applications and fostered the development of an industrial scale production process. Synthesis of ectoines requires the expenditure of considerable energetic and biosynthetic resources. Hence, microorganisms have developed ways to exploit ectoines as nutrients when they are no longer needed as stress protectants. Here, we summarize our current knowledge on the phylogenomic distribution of ectoine producing and consuming microorganisms. We emphasize the structural enzymology of the pathways underlying ectoine biosynthesis and consumption, an understanding that has been achieved only recently. The synthesis and degradation pathways critically differ in the isomeric form of the key metabolite N-acetyldiaminobutyric acid (ADABA). γ-ADABA serves as preferred substrate for the ectoine synthase, while the α-ADABA isomer is produced by the ectoine hydrolase as an intermediate in catabolism. It can serve as internal inducer for the genetic control of ectoine catabolic genes via the GabR/MocR-type regulator EnuR. Our review highlights the importance of structural enzymology to inspire the mechanistic understanding of metabolic networks at the biological scale.

Molecular analysis of microbial nitrogen transformation and removal potential in mangrove wetlands under anthropogenic nitrogen input

Mangrove ecosystems are natural nitrogen removal systems that are primarily mediated by nitrogen cycle microorganisms, but their relative contributions to nitrogen transformation and removal in mangrove sediments under anthropogenic nitrogen input needs further resolution and characterization. Here, we investigated the responses and the relative contributions of ammonia-oxidizing archaea (AOA), ammonia-oxidizing bacteria (AOB), anaerobic ammonium oxidizing (anammox) bacteria and denitrifying bacteria after spiking urea into mangrove sediments incubated in a laboratory microcosm experiment for four weeks. During incubation, the diversity, abundances and transcription levels of the hzo genes for anammox bacteria, amoA genes for AOA and AOB, and nirS genes for denitrifying bacteria were monitored using targeted gene clone library analyses and quantitative PCR assays at the DNA and RNA levels. The results showed that mangrove sediments harbour habitat-specific anammox bacteria which related to Candidatus Scalindua and Candidatus Kuenenia clades. Mangrove specific AOA related to deep branched clades within Candidatus Nitrososphaera and Candidatus Nitrosotalea, and AOB related to Nitrosomonas and Nitrosospira were also detected in the collected sediment samples. Growth and activity of AOA were detected at all levels of amendment of nitrogen input, whereas AOB growth was detectable only at the high-level nitrogen input (1.5 mg urea per gram of dry sediment) with no amoA transcripts and lower abundance than AOA. The abundance and transcription levels of the nirS gene were higher (~1000 times) than those of the hzo gene in all groups. Pearson correlation analysis demonstrated that the abundance of both AOA and AOB amoA genes had a significant positive correlation with the nirS gene (p < 0.01). These results indicated that nitrification (primarily mediated by the AOA)-denitrification process played the most important role in nitrogen removal from the amendment of nitrogen short-term input in the mangrove sediments.

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

Primordial nitrification processes have been studied extensively using geochemical approaches, but the biological origination of nitrification remains unclear. Ammonia-oxidizing archaea (AOA) are widely distributed nitrifiers and implement the rate-limiting step in nitrification. They are hypothesized to have been important players in the global nitrogen cycle in Earth’s early history. We performed systematic phylogenomic and marker gene analyses to elucidate the diversification timeline of AOA evolution. Our results suggested that the AOA ancestor experienced terrestrial geothermal environments at ∼1,165 Ma (1,928–880 Ma), and gradually evolved into mesophilic soil at ∼652 Ma (767–554 Ma) before diversifying into marine settings at ∼509 Ma (629–412 Ma) and later into shallow and deep oceans, respectively. Corroborated by geochemical evidence and modeling, the timing of key diversification nodes can be linked to the global magmatism and glaciation associated with the assembly and breakup of the supercontinent Rodinia, and the later oxygenation of the deep ocean. Results of this integrated study shed light on the geological forces that may have shaped the evolutionary pathways of the AOA, which played an important role in the ancient global nitrogen cycle.

Marine ammonia-oxidising archaea and bacteria occupy distinct iron and copper niches

Ammonia oxidation by archaea and bacteria (AOA and AOB), is the first step of nitrification in the oceans. As AOA have an ammonium affinity 200-fold higher than AOB isolates, the chemical niche allowing AOB to persist in the oligotrophic ocean remains unclear. Here we show that marine isolates, Nitrosopumilus maritimus strain SCM1 (AOA) and Nitrosococcus oceani strain C-107 (AOB) have contrasting physiologies in response to the trace metals iron (Fe) and copper (Cu), holding potential implications for their niche separation in the oceans. A greater affinity for unchelated Fe may allow AOB to inhabit shallower, euphotic waters where ammonium supply is high, but competition for Fe is rife. In contrast to AOB, AOA isolates have a greater affinity and toxicity threshold for unchelated Cu providing additional explanation to the greater success of AOA in the marine environment where Cu availability can be highly variable. Using comparative genomics, we predict that the proteomic and metal transport basis giving rise to contrasting physiologies in isolates is widespread across phylogenetically diverse marine AOA and AOB that are not yet available in pure culture. Our results develop the testable hypothesis that ammonia oxidation may be limited by Cu in large tracts of the open ocean and suggest a relatively earlier emergence of AOB than AOA when considered in the context of evolving trace metal availabilities over geologic time.

Subsurface and surface halophile communities of the chaotropic Salar de Uyuni

Salar de Uyuni (SdU) is the biggest athalosaline environment on Earth, holding a high percentage of the known world Li reserves. Due to its hypersalinity, temperature and humidity fluctuations, high exposure to UV radiation, and its elevated concentration of chaotropic agents like MgCl₂ , LiCl and NaBr, SdU is considered a polyextreme environment. Here, we report the prokaryotic abundance and diversity of 46 samples obtained in different seasons and geographical areas. The identified bacterial community was found to be more heterogeneous than the archaeal community, with both communities varying geographically. A seasonal difference has been detected for archaea. Salinibacter, Halonotius and Halorubrum were the most abundant genera in Salar de Uyuni. Different unclassified archaea were also detected. In addition, the diversity of two subsurface samples obtained at 20 and 80 m depth was evaluated and compared with the surface data, generating an evolutionary record of a multilayer hypersaline ecosystem.

Genomic Insights of "Candidatus Nitrosocaldaceae" Based on Nine New Metagenome-Assembled Genomes, Including "Candidatus Nitrosothermus" Gen Nov. and Two New Species of "Candidatus Nitrosocaldus"

“Candidatus Nitrosocaldaceae” are globally distributed in neutral or slightly alkaline hot springs and geothermally heated soils. Despite their essential role in the nitrogen cycle in high-temperature ecosystems, they remain poorly understood because they have never been isolated in pure culture, and very few genomes are available. In the present study, a metagenomics approach was employed to obtain “Ca. Nitrosocaldaceae” metagenomic-assembled genomes (MAGs) from hot spring samples collected from India and China. Phylogenomic analysis placed these MAGs within “Ca. Nitrosocaldaceae.” Average nucleotide identity and average amino acid identity analysis suggested the new MAGs represent two novel species of “Candidatus Nitrosocaldus” and a novel genus, herein proposed as “Candidatus Nitrosothermus.” Key genes responsible for chemolithotrophic ammonia oxidation and a thaumarchaeal 3HP/4HB cycle were detected in all MAGs. Furthermore, genes coding for urea degradation were only present in “Ca. Nitrosocaldus,” while biosynthesis of the vitamins, biotin, cobalamin, and riboflavin were detected in almost all MAGs. Comparison of “Ca. Nitrosocaldales/Nitrosocaldaceae” with other AOA revealed 526 specific orthogroups. This included genes related to thermal adaptation (cyclic 2,3-diphosphoglycerate, and S-adenosylmethionine decarboxylase), indicating their importance for life at high temperature. In addition, these MAGs acquired genes from members from archaea (Crenarchaeota) and bacteria (Firmicutes), mainly involved in metabolism and stress responses, which might play a role to allow this group to adapt to thermal habitats.