Genomic adaptation to eutrophication of ammonia‐oxidizing archaea in the Pearl River estuary

Discovery of Candidatus Nitrosomaritimum as a New Genus of Ammonia-Oxidizing Archaea Widespread in Anoxic Saltmarsh Intertidal Aquifers

Ammonia-oxidizing archaea (AOA) are widely distributed in marine and terrestrial habitats, contributing significantly to global nitrogen and carbon cycles. However, their genomic diversity, ecological niches, and metabolic potentials in the anoxic intertidal aquifers remain poorly understood. Here, we discovered and named a novel AOA genus, Candidatus Nitrosomaritimum, from the intertidal aquifers of Yancheng Wetland, showing close metagenomic abundance to the previously acknowledged dominant Nitrosopumilus AOA. Further construction of ammonia monooxygenase-based phylogeny demonstrated the widespread distribution of Nitrosomaritimum AOA in global estuarine-coastal niches and marine sediment. Niche differentiation among sublineages of this new genus in anoxic intertidal aquifers is driven by salinity and dissolved oxygen gradients. Comparative genomics revealed that Candidatus Nitrosomaritimum has the genetic capacity to utilize urea and possesses high-affinity phosphate transporter systems (phnCDE) for surviving phosphorus-limited conditions. Additionally, it contains putative nosZ genes encoding nitrous-oxide (N2O) reductase for reducing N2O to nitrogen gas. Furthermore, we gained first genomic insights into the archaeal phylum Hydrothermarchaeota populations residing in intertidal aquifers and revealed their potential hydroxylamine-detoxification mutualism with AOA through utilizing the AOA-released extracellular hydroxylamine using hydroxylamine oxidoreductase. Together, this study unravels the overlooked role of priorly unknown but abundant AOA lineages of the newly discovered genus Candidatus Nitrosomaritimum in biological nitrogen transformation and their potential for nitrogen pollution mitigation in coastal environments.

Dynamics and activity of an ammonia-oxidizing archaea bloom in South San Francisco Bay

Transient or recurring blooms of ammonia-oxidizing archaea (AOA) have been reported in several estuarine and coastal environments, including recent observations of AOA blooms in South San Francisco Bay. Here, we measured nitrification rates, quantified AOA abundance, and analyzed both metagenomic and metatranscriptomic data to examine the dynamics and activity of nitrifying microorganisms over the course of an AOA bloom in South San Francisco Bay during the autumn of 2018 and seasonally throughout 2019. Nitrification rates were correlated with AOA abundance in quantitative polymerase chain reaction (PCR) data, and both increased several orders of magnitude between the autumn AOA bloom and spring and summer seasons. From bloom samples, we recovered an extremely abundant, high-quality Candidatus Nitrosomarinus catalina-like AOA metagenome-assembled genome that had high transcript abundance during the bloom and expressed >80% of genes in its genome. We also recovered a putative nitrite-oxidizing bacteria metagenome-assembled genome from within the Nitrospinaceae that was of much lower abundance and had lower transcript abundance than AOA. During the AOA bloom, we observed increased transcript abundance for nitrogen uptake and oxidative stress genes in non-nitrifier metagenome-assembled genomes. This study confirms AOA are not only abundant but also highly active during blooms oxidizing large amounts of ammonia to nitrite-a key intermediate in the microbial nitrogen cycle-and producing reactive compounds that may impact other members of the microbial community.

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.

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.

Multiple metal(loid) contamination reshaped the structure and function of soil archaeal community

Archaea are important participants in biogeochemical cycles of metal(loid)-polluted ecosystems, whereas archaeal structure and function in response to metal(loid) contamination remain poorly understood. Here, the effects of multiple metal(loid) pollution on the structure and function of archaeal communities were investigated in three zones within an abandoned sewage reservoir. We found that the high-contamination zone (Zone I) had higher archaeal diversity but a lower habitat niche breadth, relative to the mid-contamination zone (Zone II) and low-contamination zone (Zone III). Particularly, metal-resistant species represented by potential methanogens were markedly enriched in Zone I (cumulative relative abundance: 32.24%) compared to Zone II (1.93%) and Zone III (0.10%), and closer inter-taxon connections and higher network complexity (based on node number, edge number, and degree) were also observed compared to other zones. Meanwhile, the higher abundances of potential metal-resistant and methanogenic functions in Zone I (0.24% and 9.24%, respectively) than in Zone II (0.08% and 7.52%) and Zone III (0.01% and 1.03%) suggested archaeal functional adaptation to complex metal(loid) contamination. More importantly, six bioavailable metal(loid)s (titanium, tin, nickel, chromium, cobalt, and zinc) were the main contributors to archaeal community variations, and metal(loid) pollution reinforced the role of deterministic processes, particularly homogeneous selection, in the archaeal community assembly. Overall, this study provides the first integrated insight into the survival strategies of archaeal communities under multiple metal(loid) contamination, which will be of significant guidance for future bioremediation and environmental governance of metal(loid)-contaminated environments.

Deterministic Factors Determine the Comammox Community Composition in the Pearl River Estuary Ecosystem

Complete ammonia oxidizers (comammox) have been widely detected in riverine and estuarine ecosystems. However, knowledge about the process of comammox community assembly from freshwater to marine environments is still limited. Here, based on deep sequencing, we investigated the community composition of comammox along a salinity gradient in the Pearl River Estuary (PRE), South China. Our results showed that comammox microorganisms in the PRE sediments were extremely diverse and displayed distinct distributional patterns between upstream and downstream habitats. Quantitative PCR demonstrated that comammox was the dominant ammonia-oxidizing microorganism (AOM) in the PRE upstream sediments, and ammonia-oxidizing archaea (AOA) dominated the PRE downstream sediments, while ammonia-oxidizing bacteria (AOB) were not dominant in any section of the PRE. Neutral modeling revealed that stochastic processes explained a limited part of the variation in the comammox community. The majority of beta nearest-taxon index values were higher than 2, indicating that comammox community assembly in the PRE sediments was better explained through a deterministic process than through a stochastic process. Salinity and total nitrogen were the most important contributing factors that shaped the comammox community. This study expanded the current knowledge of the diversity and niche preference of comammox in the estuarine ecosystem, and further enhances our understanding of the assembly of comammox community from freshwater to marine environments.

IMPORTANCE Microbial communities are shaped by stochastic (emigration, immigration, birth, death, and genetic drift of species) and deterministic (e.g., environmental factors) processes. However, it remains unknown as to which type of process is more important in influencing the comammox community assembly from freshwater to marine environments. In this study, we compared the relative importance of stochastic and deterministic processes in shaping the assembly of the comammox community, which demonstrated that the deterministic process was more important in determining the community assembly patterns in the PRE ecosystem.

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.

Ecological features and global distribution of Asgard archaea

Asgard is a newly proposed archaeal superphylum, which has been suggested to hold the key to decipher the origin of Eukaryotes. However, their ecology remains largely unknown. Here, we conducted a meta-analysis of publicly available Asgard-associated 16S rRNA gene fragments, and found that just three previously proposed clades (Lokiarchaeota, Thorarchaeota, and Asgard clade 4) are widely distributed, whereas the other seven clades (phylum or class level) are restricted to the sediment biosphere. Asgard archaea, especially Loki- and Thorarchaeota, seem to adapt to marine sediments, and water depth (the depth of the sediment below water surface) and salinity might be crucial factors for the proportion of these microorganisms as revealed by multivariate regression analyses. However, the abundance of Asgard archaea exhibited distinct environmental drivers at the clade-level; for instance, the proportion of Asgard clade 4 was higher in less saline environments (salinity <6.35 psu), while higher for Heimdallarchaeota-AAG and Asgard clade 2 in more saline environment (salinity ≥35 psu). Furthermore, co-occurrence analysis allowed us to find a significant non-random association of different Asgard clades with other groups (e.g., Lokiarchaeota with Deltaproteobacteria and Anaerolineae; Odinarchaeota with Bathyarchaeota), suggesting different interaction potentials among these clades. Overall, these findings reveal Asgard archaea as a ubiquitous group worldwide and provide initial insights into their ecological features on a global scale.

Expanding ecological assessment by integrating microorganisms into routine freshwater biomonitoring

Bioindication has become an indispensable part of water quality monitoring in most countries of the world, with the presence and abundance of bioindicator taxa, mostly multicellular eukaryotes, used for biotic indices. In contrast, microbes (bacteria, archaea and protists) are seldom used as bioindicators in routine assessments, although they have been recognized for their importance in environmental processes. Recently, the use of molecular methods has revealed unexpected diversity within known functional groups and novel metabolic pathways that are particularly important in energy and nutrient cycling. In various habitats, microbial communities respond to eutrophication, metals, and natural or anthropogenic organic pollutants through changes in diversity and function. In this review, we evaluated the common trends in these changes, documenting that they have value as bioindicators and can be used not only for monitoring but also for improving our understanding of the major processes in lotic and lentic environments. Current knowledge provides a solid foundation for exploiting microbial taxa, community structures and diversity, as well as functional genes, in novel monitoring programs. These microbial community measures can also be combined into biotic indices, improving the resolution of individual bioindicators. Here, we assess particular molecular approaches complemented by advanced bioinformatic analysis, as these are the most promising with respect to detailed bioindication value. We conclude that microbial community dynamics are a missing link important for our understanding of rapid changes in the structure and function of aquatic ecosystems, and should be addressed in the future environmental monitoring of freshwater ecosystems.

Genomic Characteristics of a Novel Species of Ammonia-Oxidizing Archaea from the Jiulong River Estuary

Ammonia-oxidizing archaea (AOA) are ubiquitous in diverse ecosystems and play a pivotal role in global nitrogen and carbon cycling. Although AOA diversity and distribution are widely studied, mainly based on the amoA (alpha subunit of ammonia monooxygenase) genotypes, only limited investigations have addressed the relationship between AOA genetic adaptation, metabolic features, and ecological niches, especially in estuaries. Here, we describe the AOA communities along the Jiulong River estuary in southern China. Nine high-quality AOA metagenome-assembled genomes (MAGs) were obtained by metagenomics. Five of the MAGs are proposed to constitute a new species, "Candidatus Nitrosopumilus aestuariumsis" sp. nov., based on the phylogenies of the 16S and 23S rRNA genes and concatenated ribosomal proteins, as well as the average amino acid identity. Comparative genomic analysis revealed unique features of the new species, including a high number of genes related to diverse carbohydrate-active enzymes, phosphatases, heavy-metal transport systems, flagellation, and chemotaxis. These genes may be crucial for AOA adaptation to the eutrophic and heavy-metal-contaminated Jiulong River estuary. The uncovered detailed genomic characteristics of the new estuarine AOA species highlight AOA contributions to ammonia oxidation in the Jiulong River estuary.IMPORTANCE In this study, AOA communities along a river in southern China were characterized, and metagenome-assembled genomes (MAGs) of a novel AOA clade were also obtained. Based on the characterization of AOA genomes, the study suggests adaptation of the novel AOAs to estuarine environments, providing new information on the ecology of estuarine AOA and the nitrogen cycle in contaminated estuarine environments.

Community, Distribution, and Ecological Roles of Estuarine Archaea

Archaea are diverse and ubiquitous prokaryotes present in both extreme and moderate environments. Estuaries, serving as links between the land and ocean, harbor numerous microbes that are relatively highly active because of massive terrigenous input of nutrients. Archaea account for a considerable portion of the estuarine microbial community. They are diverse and play key roles in the estuarine biogeochemical cycles. Ammonia-oxidizing archaea (AOA) are an abundant aquatic archaeal group in estuaries, greatly contributing estuarine ammonia oxidation. Bathyarchaeota are abundant in sediments, and they may involve in sedimentary organic matter degradation, acetogenesis, and, potentially, methane metabolism, based on genomics. Other archaeal groups are also commonly detected in estuaries worldwide. They include Euryarchaeota, and members of the DPANN and Asgard archaea. Based on biodiversity surveys of the 16S rRNA gene and some functional genes, the distribution and abundance of estuarine archaea are driven by physicochemical factors, such as salinity and oxygen concentration. Currently, increasing amount of genomic information for estuarine archaea is becoming available because of the advances in sequencing technologies, especially for AOA and Bathyarchaeota, leading to a better understanding of their functions and environmental adaptations. Here, we summarized the current knowledge on the community composition and major archaeal groups in estuaries, focusing on AOA and Bathyarchaeota. We also highlighted the unique genomic features and potential adaptation strategies of estuarine archaea, pointing out major unknowns in the field and scope for future research.

Biogeochemical controls on nitrogen transformations in subtropical estuarine wetlands

Changes in soil moisture and salinity are expected to alter gross nitrogen (N) transformations, which can control microbial N dynamics in estuarine wetlands. However, the effects of soil moisture and salinity on microbial N limitation remain poorly understood. To this end, we used a¹⁵N pool dilution approach to characterize the changes in soil gross N transformations with the variations of soil moisture in subtropical estuarine wetlands along a low-level salinity gradient. The results showed that soil gross N mineralization (GNM) increased with the increasing soil moisture. Ammonia immobilization (AIM) and microbial N immobilization (MIM) increased with the increasing soil salinity. High gross nitrification rates were generally found in the wetlands with relatively high ammonia content and low soil moisture. The ratios of MIM to GNM (MIM/GNM), as an indicator of microbial N limitation, decreased in response to the enhancing soil moisture and increased with the increasing soil salinity. Ammonia supply capacity (GNM-AIM) decreased with increasing soil salinity but increased with the increasing soil moisture. These results together indicated that microbial N limitation became stronger in the wetlands characterized with high soil salinity and low soil moisture. Soil moisture and salinity exhibited also indirect effects on microbial N limitation through affecting organic N content, ecological stoichiometry, microbial biomass and enzyme activity. Therefore, soil moisture and salinity levels are crucial for controlling soil N transformations with important implications on microbial N removal and retention in estuarine wetlands.

Hurricane Disturbance Stimulated Nitrification and Altered Ammonia Oxidizer Community Structure in Lake Okeechobee and St. Lucie Estuary (Florida)

Nitrification is an important biological link between oxidized and reduced forms of nitrogen (N). The efficiency of nitrification plays a key role in mitigating excess N in eutrophic systems, including those with cyanobacterial harmful algal blooms (cyanoHABs), since it can be closely coupled with denitrification and removal of excess N. Recent work suggests that competition for ammonium (NH₄⁺) between ammonia oxidizers and cyanoHABs can help determine microbial community structure. Nitrification rates and ammonia-oxidizing archaeal (AOA) and bacterial (AOB) community composition and gene abundances were quantified in Lake Okeechobee and St. Lucie Estuary in southern Florida (United States). We sampled during cyanobacterial (Microcystis) blooms in July 2016 and August 2017 (2 weeks before Hurricane Irma) and 10 days after Hurricane Irma made landfall. Nitrification rates were low during cyanobacterial blooms in Lake Okeechobee and St. Lucie Estuary, while low bloom conditions in St. Lucie Estuary coincided with greater nitrification rates. Nitrification rates in the lake were correlated (R² = 0.94; p = 0.006) with AOA amoA abundance. Following the hurricane, nitrification rates increased by an order of magnitude, suggesting that nitrifiers outcompeted cyanobacteria for NH₄⁺ under turbid, poor light conditions. After Irma, AOA and AOB abundances increased in St. Lucie Estuary, while only AOB increased in Lake Okeechobee. AOA sequences clustered into three major lineages: Nitrosopumilales (NP), Nitrososphaerales (NS), and Nitrosotaleales (NT). Many of the lake OTUs placed within the uncultured and uncharacterized NS δ and NT β clades, suggesting that these taxa are ecologically important along this eutrophic, lacustrine to estuarine continuum. After the hurricane, the AOA community shifted toward dominance by freshwater clades in St. Lucie Estuary and terrestrial genera in Lake Okeechobee, likely due to high rainfall and subsequent increased turbidity and freshwater loading from the lake into the estuary. AOB community structure was not affected by the disturbance. AOA communities were consistently more diverse than AOB, despite fewer sequences recovered, including new, unclassified, eutrophic ecotypes, suggesting a wider ecological biogeography than the oligotrophic niche originally posited. These results and other recent reports contradict the early hypothesis that AOB dominate ammonia oxidation in high-nutrient or terrestrial-influenced systems.

Alternative strategies of nutrient acquisition and energy conservation map to the biogeography of marine ammonia-oxidizing archaea

Ammonia-oxidizing archaea (AOA) are among the most abundant and ubiquitous microorganisms in the ocean, exerting primary control on nitrification and nitrogen oxides emission. Although united by a common physiology of chemoautotrophic growth on ammonia, a corresponding high genomic and habitat variability suggests tremendous adaptive capacity. Here, we compared 44 diverse AOA genomes, 37 from species cultivated from samples collected across diverse geographic locations and seven assembled from metagenomic sequences from the mesopelagic to hadopelagic zones of the deep ocean. Comparative analysis identified seven major marine AOA genotypic groups having gene content correlated with their distinctive biogeographies. Phosphorus and ammonia availabilities as well as hydrostatic pressure were identified as selective forces driving marine AOA genotypic and gene content variability in different oceanic regions. Notably, AOA methylphosphonate biosynthetic genes span diverse oceanic provinces, reinforcing their importance for methane production in the ocean. Together, our combined comparative physiological, genomic, and metagenomic analyses provide a comprehensive view of the biogeography of globally abundant AOA and their adaptive radiation into a vast range of marine and terrestrial habitats.

Metabolic response of prokaryotic microbes to sporadic hypoxia in a eutrophic subtropical estuary

Coastal eutrophication and consequent oxygen depletion (hypoxia) occurs worldwide due to increased human activity. The paucity of genomic information of microbes in hypoxia prone coastal waters have hindered our understanding of microorganism related causation and adaption to the environment. Here, using metagenomic approach, we investigated microbial metabolic capability in heavily polluted Pearl River estuary. Our results highlighted the possible roles of microbial metabolic activity in the formation of bottom water hypoxia by revealing enriched organic degradation related microbial genes in the bottom layer beneath surface phytoplankton bloom. Microbial nitrate reduction in hypoxia layer was low, possibly due to the low pH and fluctuating oxygen level. On contrary, high abundance of sulfate-reducing, and antibiotic and metal resistance related genes were detected in bottom and surface layers, respectively, indicating microbial adaptation to oxygen depletion and pollution. Our study provides gene level information on the interactive relations between microbial functions and environmental stress.