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

Application of single-cell Raman-deuterium isotope probing to reveal the resistance of marine ammonia-oxidizing archaea SCM1 against common antibiotics

Antimicrobial resistance (AMR) in oceans poses a significant threat to human health through the seafood supply chain. Ammonia-oxidizing archaea (AOA) are important marine microorganisms and play a key role in the biogeochemical nitrogen cycle around the world. However, the AMR of marine AOA to aquicultural antibiotics is poorly explored. Here, Raman-deuterium isotope probing (Raman-DIP), a single-cell tool, was developed to reveal the AMR of a typical marine species of AOA, Nitrosopumilus maritimus (designated SCM1), against six antibiotics, including erythromycin, tetracycline, novobiocin, neomycin, bacitracin, and vancomycin. The D2O concentration (30% v/v) and culture period (9 days) were optimized for the precise detection of metabolic activity in SCM1 cells through Raman-DIP. The relative metabolic activity of SCM1 upon exposure to antibiotics was semi-quantitatively calculated based on single-cell Raman spectra. SCM1 exhibited high resistance to erythromycin, tetracycline, novobiocin, neomycin, and vancomycin, with minimum inhibitory concentration (MIC) values between 100 and 400 mg/L while SCM1 is very sensitive to bacitracin (MIC: 0.8 mg/L). Notably, SCM1 cells were completely inactive under the metabolic activity minimum inhibitory concentration conditions (MA-MIC: 1.6∼800 mg/L) for the six antibiotics. Further genomic analysis revealed the antibiotic resistance genes (ARGs) of SCM1 including 14 types categorized into 33 subtypes. This work increases our knowledge of the AMR of marine AOA by linking the resistant phenome to the genome, contributing to the risk assessment of AMR in the underexplored ocean environment. As antibiotic resistance in marine microorganisms is significantly affected by the concentration of antibiotics in coastal environments, we encourage more studies concentrating on both the phenotypic and genotypic antibiotic resistance of marine archaea. This may facilitate a comprehensive evaluation of the capacity of marine microorganisms to spread AMR and the implementation of suitable control measures to protect environmental safety and human health.

15N-DNA stable isotope probing reveals niche differentiation of ammonia oxidizers in paddy soils

Chemoautotrophic canonical ammonia oxidizers (ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB)) and complete ammonia oxidizers (comammox Nitrospira) are accountable for ammonia oxidation, which is a fundamental process of nitrification in terrestrial ecosystems. However, the relationship between autotrophic nitrification and the active nitrifying populations during 15N-urea incubation has not been totally clarified. The 15N-labeled DNA stable isotope probing (DNA-SIP) technique was utilized in order to study the response from the soil nitrification process and the active nitrifying populations, in both acidic and neutral paddy soils, to the application of urea. The presence of C2H2 almost completely inhibited NO3--N production, indicating that autotrophic ammonia oxidation was dominant in both paddy soils. 15N-DNA-SIP technology could effectively distinguish active nitrifying populations in both soils. The active ammonia oxidation groups in both soils were significantly different, AOA (NS (Nitrososphaerales)-Alpha, NS-Gamma, NS-Beta, NS-Delta, NS-Zeta and NT (Ca. Nitrosotaleales)-Alpha), and AOB (Nitrosospira) were functionally active in the acidic paddy soil, whereas comammox Nitrospira clade A and Nitrosospira AOB were functionally active in the neutral paddy soil. This study highlights the effective discriminative effect of 15N-DNA-SIP and niche differentiation of nitrifying populations in these paddy soils. KEY POINTS: • 15N-DNA-SIP technology could effectively distinguish active ammonia oxidizers. • Comammox Nitrospira clade A plays a lesser role than canonical ammonia oxidizers. • The active groups in the acidic and neutral paddy soils were significantly different.

Combinatorial characterization of bacterial taxa-driven differences in the microbiome of oyster reefs

Oyster reefs are invaluable ecosystems that provide a wide array of critical ecosystem services, including water filtration, coastal protection, and habitat provision for various marine species. However, these essential habitats face escalating threats from climate change and anthropogenic stressors. To combat these challenges, numerous oyster restoration initiatives have been undertaken, representing a global effort to preserve and restore these vital ecosystems. A significant, yet poorly understood, component of oyster reefs is the microbial communities. These communities account for a substantial proportion of marine reefs and are pivotal in driving key biogeochemical processes. Particularly, the environmental microbiome plays a crucial role in supporting the health and resilience of oyster populations. In our study, we sought to shed light on the microbiome within oyster reef ecosystems by characterizing the abundance, and diversity of microorganisms in the soil, biofilm, and oysters in 4 sites using a combinatorial approach to identify differentially abundant microbes by sample type and by sampling location. Our investigation revealed distinct microbial taxa in oysters, sediment and biofilm. The maximum Shannon Index indicated a slightly increased diversity in Heron's Head (5.47), followed by Brickyard park (5.35), Dunphy Park (5.17) and Point Pinole (4.85). This is likely to be driven by significantly higher oyster mortality observed at Point Pinole during routine monitoring and restoration efforts. Interestingly Ruminococcus, Streptococcus, Staphylococcus, Prevotella, Porphyromonas, Parvimonas, Neisseria, Lactococcus, Haemophilus, Fusobacterium, Dorea, Clostridium, Campylobacter, Bacteroides, and Akkermansia were positively associated with the biofilm. Yet we have limited understanding of their beneficial and/or detrimental implications to oyster growth and survival. By unraveling the intricate relationships in microbial composition across an oyster reef, our study contributes to advancing the knowledge needed to support effective oyster reef conservation and restoration efforts.

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

Introduction

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

Methods

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

Results and Discussion

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

Rare subcommunity maintains the stability of ecosystem multifunctionality by deterministic assembly processes in subtropical estuaries

Microorganisms, especially rare microbial species, are crucial in estuarine ecosystems for driving biogeochemical processes and preserving biodiversity. However, the understanding of the links between ecosystem multifunctionality (EMF) and the diversity of rare bacterial taxa in estuary ecosystems remains limited. Employing high-throughput sequencing and a variety of statistical methods, we assessed the diversities and assembly process of abundant and rare bacterioplankton and their contributions to EMF in a subtropical estuary. Taxonomic analysis revealed Proteobacteria as the predominant phylum among both abundant and rare bacterial taxa. Notably, rare taxa demonstrated significantly higher taxonomic diversity and a larger species pool than abundant taxa. Additionally, our findings highlighted that deterministic assembly processes predominantly shape microbial communities, with heterogeneous selection exerting a stronger influence on rare taxa. Further analysis reveals that rare bacterial beta-diversity significantly impacts to EMF, whereas alpha diversity did not. The partial least squares path modeling (PLS-PM) analysis demonstrated that the beta diversity of abundant and rare taxa, as the main biotic factor, directly affected EMF, while temperature and total organic carbon (TOC) were additional key factors to determine the relationship between beta diversity and EMF. These findings advance our understanding of the distribution features and ecological knowledge of the abundant and rare taxa in EMF in subtropical estuaries, and provide a reference for exploring the multifunctionality of different biospheres in aquatic environments.

Metagenome-assembled genomes of two nitrifying microorganisms from the freshwater archaeal ammonia-oxidizing enrichment culture BO1

Two metagenome-assembled genomes (MAGs) were recovered from the ammonia-oxidizing enrichment culture BO1 obtained from the sediment of the freshwater reservoir Lake Burr Oak, Ohio, USA. High quality MAGs were assembled for the archaeal ammonia oxidizer Nitrosarchaeum sp. BO1 and the canonical nitrite oxidizer Nitrospira sp. BO1.

A unique subseafloor microbiosphere in the Mariana Trench driven by episodic sedimentation

Hadal trenches are characterized by enhanced and infrequent high-rate episodic sedimentation events that likely introduce not only labile organic carbon and key nutrients but also new microbes that significantly alter the subseafloor microbiosphere. Currently, the role of high-rate episodic sedimentation in controlling the composition of the hadal subseafloor microbiosphere is unknown. Here, analyses of carbon isotope composition in a ~ 750 cm long sediment core from the Challenger Deep revealed noncontinuous deposition, with anomalous 14C ages likely caused by seismically driven mass transport and the funneling effect of trench geomorphology. Microbial community composition and diverse enzyme activities in the upper ~ 27 cm differed from those at lower depths, probably due to sudden sediment deposition and differences in redox condition and organic matter availability. At lower depths, microbial population numbers, and composition remained relatively constant, except at some discrete depths with altered enzyme activity and microbial phyla abundance, possibly due to additional sudden sedimentation events of different magnitude. Evidence is provided of a unique role for high-rate episodic sedimentation events in controlling the subsurface microbiosphere in Earth's deepest ocean floor and highlight the need to perform thorough analysis over a large depth range to characterize hadal benthic populations. Such depositional processes are likely crucial in shaping deep-water geochemical environments and thereby the deep subseafloor biosphere.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00212-y.

Ammonia-oxidizing bacteria and archaea exhibit differential nitrogen source preferences

Ammonia-oxidizing microorganisms (AOM) contribute to one of the largest nitrogen fluxes in the global nitrogen budget. Four distinct lineages of AOM: ammonia-oxidizing archaea (AOA), beta- and gamma-proteobacterial ammonia-oxidizing bacteria (β-AOB and γ-AOB) and complete ammonia oxidizers (comammox), are thought to compete for ammonia as their primary nitrogen substrate. In addition, many AOM species can utilize urea as an alternative energy and nitrogen source through hydrolysis to ammonia. How the coordination of ammonia and urea metabolism in AOM influences their ecology remains poorly understood. Here we use stable isotope tracing, kinetics and transcriptomics experiments to show that representatives of the AOM lineages employ distinct regulatory strategies for ammonia or urea utilization, thereby minimizing direct substrate competition. The tested AOA and comammox species preferentially used ammonia over urea, while β-AOB favoured urea utilization, repressed ammonia transport in the presence of urea and showed higher affinity for urea than for ammonia. Characterized γ-AOB co-utilized both substrates. These results reveal contrasting niche adaptation and coexistence patterns among the major AOM lineages.

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.

Distribution and genomic variation of ammonia-oxidizing archaea in abyssal and hadal surface sediments

Ammonia-oxidizing archaea of the phylum Thaumarchaeota play a central role in the biogeochemical cycling of nitrogen in benthic sediments, at the interface between pelagic and subsurface ecosystems. However, our understanding of their niche separation and of the processes controlling their population structure in hadal and abyssal surface sediments is still limited. Here, we reconstructed 47 AOA metagenome-assembled genomes (MAGs) from surface sediments of the Atacama and Kermadec trench systems. They formed deep-sea-specific groups within the family Nitrosopumilaceae and were assigned to six amoA gene-based clades. MAGs from different clades had distinct distribution patterns along oxygen-ammonium counter gradients in surface sediments. At the species level, MAGs thus seemed to form different ecotypes and follow deterministic niche-based distributions. In contrast, intraspecific population structure, defined by patterns of Single Nucleotide Variants (SNV), seemed to reflect more complex contributions of both deterministic and stochastic processes. Firstly, the bathymetric range had a strong effect on population structure, with distinct populations in abyssal plains and hadal trenches. Then, hadal populations were clearly separated by trench system, suggesting a strong isolation-by-topography effect, whereas abyssal populations were rather controlled by sediment depth or geographic distances, depending on the clade considered. Interestingly, genetic variability between samples was lowest in sediment layers where the mean MAG coverage was highest, highlighting the importance of selective pressure linked with each AOA clade's ecological niche. Overall, our results show that deep-sea AOA genome distributions seem to follow both deterministic and stochastic processes, depending on the genomic variability scale considered.

Transcriptomic Insights into Archaeal Nitrification in the Amundsen Sea Polynya, Antarctica

Antarctic polynyas have the highest Southern Ocean summer primary productivity, and due to anthropogenic climate change, these areas have formed faster recently. Ammonia-oxidizing archaea (AOA) are among the most ubiquitous and abundant microorganisms in the ocean and play a primary role in the global nitrogen cycle. We utilized metagenomics and metatranscriptomics to gain insights into the physiology and metabolism of AOA in polar oceans, which are associated with ecosystem functioning. A polar-specific ecotype of AOA, from the "Candidatus Nitrosomarinus"-like group, was observed to be dominant in the Amundsen Sea Polynya (ASP), West Antarctica, during a succession of summer phytoplankton blooms. AOA had the highest transcriptional activity among prokaryotes during the bloom decline phase (DC). Metatranscriptomic analysis of key genes involved in ammonia oxidation, carbon fixation, transport, and cell division indicated that this polar AOA ecotype was actively involved in nitrification in the bloom DC in the ASP. This study revealed the physiological and metabolic traits of this key polar-type AOA in response to phytoplankton blooms in the ASP and provided insights into AOA functions in polar oceans.

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.

Distinct patterns of distribution, community assembly and cross-domain co-occurrence of planktonic archaea in four major estuaries of China

BACKGROUND

Archaea are key mediators of estuarine biogeochemical cycles, but comprehensive studies comparing archaeal communities among multiple estuaries with unified experimental protocols during the same sampling periods are scarce. Here, we investigated the distribution, community assembly, and cross-domain microbial co-occurrence of archaea in surface waters across four major estuaries (Yellow River, Yangtze River, Qiantang River, and Pearl River) of China cross climatic zones (~ 1,800 km) during the winter and summer cruises.

RESULTS

The relative abundance of archaea in the prokaryotic community and archaeal community composition varied with estuaries, seasons, and stations (reflecting local environmental changes such as salinity). Archaeal communities in four estuaries were overall predominated by ammonia-oxidizing archaea (AOA) (aka. Marine Group (MG) I; primarily Nitrosopumilus), while the genus Poseidonia of Poseidoniales (aka. MGII) was occasionally predominant in Pearl River estuary. The cross-estuary dispersal of archaea was largely limited and the assembly mechanism of archaea varied with estuaries in the winter cruise, while selection governed archaeal assembly in all estuaries in the summer cruise. Although the majority of archaea taxa in microbial networks were peripherals and/or connectors, extensive and distinct cross-domain associations of archaea with bacteria were found across the estuaries, with AOA as the most crucial archaeal group. Furthermore, the expanded associations of MGII taxa with heterotrophic bacteria were observed, speculatively indicating the endogenous demand for co-processing high amount and diversity of organic matters in the estuarine ecosystem highly impacted by terrestrial/anthropogenic input, which is worthy of further study.

CONCLUSIONS

Our results highlight the lack of common patterns in the dynamics of estuarine archaeal communities along the geographic gradient, expanding the understanding of roles of archaea in microbial networks of this highly dynamic ecosystem.

Nitrosopumilus as main source of isoprenoid glycerol dialkyl glycerol tetraether lipids in the central Baltic Sea

Nitrososphaeria in the phylum Crenarchaeota, is a widespread archaeal class in the oceanic realm, playing an important role in the marine carbon and nitrogen cycle. Nitrososphaeria-derived membrane lipids, i.e., isoprenoid glycerol dialkyl glycerol tetraethers (GDGTs), are commonly employed to reconstruct past water temperatures using the TetraEther indeX of 86 carbon atoms (TEX86). This index is of particular importance for the brackish Baltic Sea as to date it appears to be the only applicable organic temperature proxy. In this study, we investigated the distribution of intact and core GDGTs and their potential source organisms in the water column of three deep basins located in the central Baltic Sea to evaluate the application of TEX86. A lipidomic approach on suspended particulate matter was combined with the molecular techniques 16S rRNA gene amplicon sequencing and CARD-FISH. The archaeal community was dominated by Nitrosopumilus (~83-100% of the total archaeal sequences). As other detected taxa known to produce GDGTs each represented less than 2% of the total archaeal sequences, Nitrosopumilus is likely the most dominant GDGT producer in the central Baltic Sea. However, the occurrence of phosphohexose (PH), instead of hexose-phosphohexose (HPH) headgroups, suggested that Nitrosopumilus in the Baltic Sea may differ physiologically from representatives of marine settings and other marginal seas, such as the Black Sea. In the Baltic Sea, Nitrosopumilus is most abundant in the suboxic zone, where intact cells peak according to both CARD-FISH data and intact polar lipid concentrations. The presented data therefore suggest that TEX86 reflects subsurface rather than surface temperature in the central Baltic Sea.

Archaeal lipids

The major archaeal membrane glycerolipids are distinguished from those of bacteria and eukaryotes by the contrasting stereochemistry of their glycerol backbones, and by the use of ether-linked isoprenoid-based alkyl chains rather than ester-linked fatty acyl chains for their hydrophobic moieties. These fascinating compounds play important roles in the extremophile lifestyles of many species, but are also present in the growing numbers of recently discovered mesophilic archaea. The past decade has witnessed significant advances in our understanding of archaea in general and their lipids in particular. Much of the new information has come from the ability to screen large microbial populations via environmental metagenomics, which has revolutionised our understanding of the extent of archaeal biodiversity that is coupled with a strict conservation of their membrane lipid compositions. Significant additional progress has come from new culturing and analytical techniques that are gradually enabling archaeal physiology and biochemistry to be studied in real time. These studies are beginning to shed light on the much-discussed and still-controversial process of eukaryogenesis, which probably involved both bacterial and archaeal progenitors. Puzzlingly, although eukaryotes retain many attributes of their putative archaeal ancestors, their lipid compositions only reflect their bacterial progenitors. Finally, elucidation of archaeal lipids and their metabolic pathways have revealed potentially interesting applications that have opened up new frontiers for biotechnological exploitation of these organisms. This review is concerned with the analysis, structure, function, evolution and biotechnology of archaeal lipids and their associated metabolic pathways.

Comparative Functional Genome Analysis Reveals the Habitat Adaptation and Biocontrol Characteristics of Plant Growth-Promoting Bacteria in NCBI Databases

Plant growth-promoting bacteria (PGPB) are a group of beneficial microorganisms that include 60 bacterial genera, such as Bacillus, Pseudomonas, and Burkholderia, which widely colonize plant leaves and soil, promote plant growth, and/or inhibit pathogen infection. However, the genetic factors underpinning adaptation of PGPB to plant leaves and soil remain poorly understood. In this study, we performed a comparative functional genome analysis approach to investigate the functional genes of 195 leaf-associated (LA) and 283 soil-associated (SA) PGPB strains and their roles in adapting to their environment, using 95 strains from other-associated (OA) environmental habitats with growth-promoting or antimicrobial functions as negative controls. Comparison analysis of the enrichment of nonredundant (NR) protein sequence databases showed that cytochrome P450, DNA repair, and motor chemotaxis genes were significantly enriched in LA PGPB strains related to environmental adaptation, while cell wall-degrading enzymes, TetR transcriptional regulatory factors, and sporulation-related genes were highly enriched in SA PGPB strains. Additionally, analysis of carbohydrate-active enzymes demonstrated that glycosyltransferases (GTs) and glycoside hydrolases (GHs) were abundant families in all PGPB strains, which is in favor of plant growth, and enriched in SA PGPB strains. Except for most Bacillus strains, SA PGPB genomes contained significantly more secondary metabolism clusters than LA PGPB. Most LA PGPB contained hormone biosynthesis genes, which may contribute to plant growth promotion, while SA PGPB harbored numerous carbohydrate and antibiotic metabolism genes. In summary, this study further deepens our understanding of the habitat adaptation and biocontrol characteristics of LA and SA PGPB strains. IMPORTANCE Plant growth-promoting bacteria (PGPB) are essential for the effectiveness of biocontrol agents in plant phyllosphere and rhizosphere. However, little is known about the ecological adaptation of PGPB to different habitats. In this study, comparative functional genome analysis of leaf-associated (LA), soil-associated (SA), and other-associated (OA) PGPB strains was performed. We found that genes related to the metabolism of hormones were enriched in LA PGPB. Carbohydrate and antibiotic metabolism genes were enriched in SA PGPB, which likely facilitated their adaptation to the plant growth environment. Our findings provide genetic insights on LA and SA PGPB strains' ecological adaptation and biocontrol characteristics.

Phylotype resolved spatial variation and association patterns of planktonic Thaumarchaeota in eastern Chinese marginal seas

The majority of marine ammonia oxidizers belong to Thaumarchaeota, a phylum of Archaea, which is distributed throughout the water column. Marine surface waters contain distinct thaumarchaeotal phylotypes compared to the deeper ocean, but spatial dynamics of the surface-associated lineages are largely unsolved. This study of 120 seawater samples from the eastern Chinese marginal seas identified contrasting distribution and association patterns among thaumarchaeotal phylotypes across different dimensions. Horizontally, Nitrosopumilus-like and Nitrosopelagicus-like phylotypes dominated the surface water (3 m) of the Yellow Sea (YS) and East China Sea (ECS), respectively, along with increased abundance of total free-living Thaumarchaeota in ECS. Similar compositional changes were observed in the surface microlayer. The spatial heterogeneity of particle-attached Thaumarchaeota was less clear in surface microlayers than in surface waters. Vertically, the Nitrosopelagicus-like phylotype increased in abundance from surface to 90 m in ECS, which led to an increase in the proportion of Thaumarchaeota relative to total prokaryotes. This occurred mainly in the free-living fraction. These results indicate a clear size-fractionated niche partitioning, which is more pronounced at lower depths than in the surface water/surface microlayer. In addition, associations of Thaumarchaeota with other microbial taxa varied between phylotypes and size fractions. Our results show that a phylotype-resolved and size-fractionated spatial heterogeneity of the thaumarchaeotal community is present in surface oceanic waters and a vertical variation of the Nitrosopelagicus-like phylotype is present in shallow shelf waters.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00169-y.

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.

Temporal profiling resolves the drivers of microbial nitrogen cycling variability in coastal sediments

Here we describe the potential for sediment microbial nitrogen-cycling gene (DNA) and activity (RNA) abundances to spatially resolve coastal areas impacted by seasonal variability in external nutrient inputs. Three sites were chosen within a nitrogen-limited embayment, Port Phillip Bay (PPB), Australia that reflect variability in both proximity to external nutrient inputs and the dominant form of available nitrogen. At three sediment depths (0-1; 1-5; 5-10 cm) across a 2 year study key genes involved in nitrification (archaeal amoA and bacterial β-amoA), nitrite reduction (clade I nirS and cluster I nirK, archaeal nirK-a), anaerobic oxidation of ammonium (anammox 16S rRNA phylogenetic marker) and nitrogen fixation (nifH) were quantified. Sediments impacted by a dominance of organic nitrogen inputs were characterised at all time-points and to sediment depths of 10 cm by the highest transcript abundances of archaeal amoA and archaeal nirk-a. Proximity to a dominance of external nitrate inputs was associated with the highest transcript abundances of nirS which temporally co-varied with seasonal changes in sediment nitrate. Sediments isolated from external inputs displayed the greatest depth-specific decrease in quantifiable transcript abundances. In these isolated sediments bacterial β-amoA transcripts were temporally associated with increased sediment ammonium levels. Across this nitrogen limited system variability in the abundance of bacterial β-amoA, archaeal amoA, archaeal nirk-a or nirS transcripts from the sediment surface (0-1 and 5 cm) demonstrated a capacity to improve our ability to monitor coastal zones impacted by anthropogenic nitrogen inputs. Specifically, the spatial detection sensitivity of bacterial β-amoA transcripts could be developed as a metric to determine spatiotemporal impacts of large external loading events. This temporal study demonstrates a capacity for microbial activity metrics to facilitate coastal management strategies through greater spatial resolution of areas impacted by external nutrient inputs.

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.

Periphytic biofilms-mediated microbial interactions and their impact on the nitrogen cycle in rice paddies

Rice paddies are unique waterlogged wetlands artificially constructed for agricultural production. Periphytic biofilms (PBs) at the soil-water interface play an important role in rice paddies characterized by high nutrient input but low utilization efficiency. PBs are composed of microbial aggregates, including a wide variety of microorganisms (algae, bacteria, fungi, protozoa, and metazoa), extracellular polymeric substances and minerals (iron, aluminum, and calcium), which form an integrated food web and energy flux within a relatively stable micro-ecosystem. PBs are crucial to regulate and streamline the nitrogen cycle by neutralizing nitrogen losses and improving rice production since PBs can serve as both a sink by capturing surplus nitrogen and a source by slowly re-releasing this nitrogen for reutilization. Here the ecological advantages of PBs in regulating the nitrogen cycle in rice paddies are illustrated. We summarize the key functional importance of PBs, including the intricate and delicate community structure, microbial interactions among individual phylotypes, a wide diversity of self-produced organics, the active adaptation of PBs to constantly changing environments, and the intricate mechanisms by which PBs regulate the nitrogen cycle. We also identify the future challenges of microbial interspecific cooperation in PBs and their quantitative contributions to agricultural sustainability, optimizing nitrogen utilization and crop yields in rice paddies.

Comparative Genomics of Thaumarchaeota From Deep-Sea Sponges Reveal Their Niche Adaptation

Thaumarchaeota account for a large portion of microbial symbionts in deep-sea sponges and are even dominant in some cases. In this study, we investigated three new sponge-associated Thaumarchaeota from the deep West Pacific Ocean. Thaumarchaeota were found to be the most dominant phylum in this sponge by both prokaryotic 16S rRNA amplicons and metagenomic sequencing. Fifty-seven published Thaumarchaeota genomes from sponges and other habitats were included for genomic comparison. Similar to shallow sponge-associated Thaumarchaeota, those Thaumarchaeota in deep-sea sponges have extended genome sizes and lower coding density compared with their free-living lineages. Thaumarchaeota in deep-sea sponges were specifically enriched in genes related to stress adapting, symbiotic adhesion and stability, host–microbe interaction and protein transportation. The genes involved in defense mechanisms, such as the restriction-modification system, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas system, and toxin-antitoxin system were commonly enriched in both shallow and deep sponge-associated Thaumarchaeota. Our study demonstrates the significant effects of both depth and symbiosis on forming genomic characteristics of Thaumarchaeota, and provides novel insights into their niche adaptation in deep-sea sponges.

Exaggerated trans-membrane charge of ammonium transporters in nutrient-poor marine environments

Transporter proteins are a vital interface between cells and their environment. In nutrient-limited environments, microbes with transporters that are effective at bringing substrates into their cells will gain a competitive advantage over variants with reduced transport function. Microbial ammonium transporters (Amt) bring ammonium into the cytoplasm from the surrounding periplasm space, but diagnosing Amt adaptations to low nutrient environments solely from sequence data has been elusive. Here, we report altered Amt sequence amino acid distribution from deep marine samples compared to variants sampled from shallow water in two important microbial lineages of the marine water column community-Marine Group I Archaea (Thermoproteota) and the uncultivated gammaproteobacterial lineage SAR86. This pattern indicates an evolutionary pressure towards an increasing dipole in Amt for these clades in deep ocean environments and is predicted to generate stronger electric fields facilitating ammonium acquisition. This pattern of increasing dipole charge with depth was not observed in lineages capable of accessing alternative nitrogen sources, including the abundant alphaproteobacterial clade SAR11. We speculate that competition for ammonium in the deep ocean drives transporter sequence evolution. The low concentration of ammonium in the deep ocean is therefore likely due to rapid uptake by Amts concurrent with decreasing nutrient flux.

Differentiated Evolutionary Strategies of Genetic Diversification in Atlantic and Pacific Thaumarchaeal Populations

Some marine microbes are seemingly “ubiquitous,” thriving across a wide range of environmental conditions. While the increased depth in metagenomic sequencing has led to a growing body of research on within-population heterogeneity in environmental microbial populations, there have been fewer systematic comparisons and characterizations of population-level genetic diversity over broader expanses of time and space. Here, we investigated the factors that govern the diversification of ubiquitous microbial taxa found within and between ocean basins. Specifically, we use mapped metagenomic paired reads to examine the genetic diversity of ammonia-oxidizing archaeal (“Candidatus Nitrosopelagicus brevis”) populations in the Pacific (Hawaii Ocean Time-series [HOT]) and Atlantic (Bermuda Atlantic Time Series [BATS]) Oceans sampled over 2 years. We observed higher nucleotide diversity in “Ca. N. brevis” at HOT, driven by a higher rate of homologous recombination. In contrast, “Ca. N. brevis” at BATS featured a more open pangenome with a larger set of genes that were specific to BATS, suggesting a history of dynamic gene gain and loss events. Furthermore, we identified highly differentiated genes that were regulatory in function, some of which exhibited evidence of recent selective sweeps. These findings indicate that different modes of genetic diversification likely incur specific adaptive advantages depending on the selective pressures that they are under. Within-population diversity generated by the environment-specific strategies of genetic diversification is likely key to the ecological success of “Ca. N. brevis.”

IMPORTANCE Ammonia-oxidizing archaea (AOA) are one of the most abundant chemolithoautotrophic microbes in the marine water column and are major contributors to global carbon and nitrogen cycling. Despite their ecological importance and geographical pervasiveness, there have been limited systematic comparisons and characterizations of their population-level genetic diversity over time and space. Here, we use metagenomic time series from two ocean observatories to address the fundamental questions of how abiotic and biotic factors shape the population-level genetic diversity and how natural microbial populations adapt across diverse habitats. We show that the marine AOA “Candidatus Nitrosopelagicus brevis” in different ocean basins exhibits distinct modes of genetic diversification in response to their selective regimes shaped by nutrient availability and patterns of environmental fluctuations. Our findings specific to “Ca. N. brevis” have broader implications, particularly in understanding the population-level responses to the changing climate and predicting its impact on biogeochemical cycles.

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.

Niche differentiation of ammonia-oxidizing archaea and related autotrophic carbon fixation potential in the water column of the South China Sea

The significant primary production by ammonia-oxidizing archaea (AOA) in the ocean was reported, but the carbon fixation process of AOA and its community composition along the water depth remain unclear. Here, we investigated the abundance, community composition, and potential carbon fixation of AOA in water columns of the South China Sea. Higher abundances of the amoA and accA genes of AOA were found below the euphotic zone. Similarly, higher carbon fixation potential of AOA, evaluated by the ratios of amoA to accA gene, was also observed below euphotic zone and the ratios increased with increasing water depth. The vertical niche differentiation of AOA was further evidenced, with the dominant genus shifting from Nitrosopelagicus in the epipelagic zone to uncultured genus in the meso- and bathypelagic zones. Our findings highlight the higher carbon fixation potential of AOA in deep water and the significance of AOA to the ocean carbon budget.

Comparative Genomics Unveils the Habitat Adaptation and Metabolic Profiles of Clostridium in an Artificial Ecosystem for Liquor Production

Pit mud is a typical artificial ecosystem for Chinese liquor production. Clostridium inhabiting pit mud plays essential roles in the flavor formation of strong-flavor baijiu. The relative abundance of Clostridium increased with pit mud quality, further influencing the quality of baijiu.

Soils and sediments host Thermoplasmata archaea encoding novel copper membrane monooxygenases (CuMMOs)

Copper membrane monooxygenases (CuMMOs) play critical roles in the global carbon and nitrogen cycles. Organisms harboring these enzymes perform the first, and rate limiting, step in aerobic oxidation of ammonia, methane, or other simple hydrocarbons. Within archaea, only organisms in the order Nitrososphaerales (Thaumarchaeota) encode CuMMOs, which function exclusively as ammonia monooxygenases. From grassland and hillslope soils and aquifer sediments, we identified 20 genomes from distinct archaeal species encoding divergent CuMMO sequences. These archaea are phylogenetically clustered in a previously unnamed Thermoplasmatota order, herein named the Ca. Angelarchaeales. The CuMMO proteins in Ca. Angelarchaeales are more similar in structure to those in Nitrososphaerales than those of bacteria, and contain all functional residues required for general monooxygenase activity. Ca. Angelarchaeales genomes are significantly enriched in blue copper proteins (BCPs) relative to sibling lineages, including plastocyanin-like electron carriers and divergent nitrite reductase-like (nirK) 2-domain cupredoxin proteins co-located with electron transport machinery. Ca. Angelarchaeales also encode significant capacity for peptide/amino acid uptake and degradation and share numerous electron transport mechanisms with the Nitrososphaerales. Ca. Angelarchaeales are detected at high relative abundance in some of the environments where their genomes originated from. While the exact substrate specificities of the novel CuMMOs identified here have yet to be determined, activity on ammonia is possible given their metabolic and ecological context. The identification of an archaeal CuMMO outside of the Nitrososphaerales significantly expands the known diversity of CuMMO enzymes in archaea and suggests previously unaccounted organisms contribute to critical global nitrogen and/or carbon cycling functions.

Determination of 15N/14N of Ammonium, Nitrite, Nitrate, Hydroxylamine, and Hydrazine Using Colorimetric Reagents and Matrix-Assisted Laser Desorption Ionization-Time-of-Flight Mass Spectrometry (MALDI-TOF MS)

In the nitrogen (N) cycle, nitrogenous compounds are chemically and biologically converted to various aqueous and gaseous N species. The ¹⁵N-labeling approach is a powerful culture-dependent technique to obtain insights into the complex nitrogen transformation reactions that occur in cultures. In the ¹⁵N-labeling approach, the fates of supplemented ¹⁵N- and/or unlabeled gaseous and aqueous compounds are tracked by mass spectrometry (MS) analysis, whereas MS analysis of aqueous N species requires laborious sample preparation steps and is performed using isotope-ratio mass spectrometry, which requires an expensive mass spectrometer. We developed a simple and high-throughput MS method for determining the ¹⁵N atoms percent of NH₄⁺, NO₂^-, NO₃^-, NH₂OH, and N₂H₄, where liquid samples (<0.5 mL) were mixed with colorimetric reagents (naphthylethylenediamine for NO₂^-, indophenol for NH₄⁺, and p-aminobenzaldehyde for N₂H₄), and the mass spectra of the formed N complex dyes were obtained by matrix-assisted laser desorption ionization-time-of-flight (MALDI-TOF) MS. NH₂OH and NO₃^- were chemically converted to NO₂^- by iodine oxidation and copper/hydrazine reduction reaction, respectively, prior to the above colorimetric reaction. The intensity of the isotope peak (M + 1 or M + 2) increased when the N complex dye was formed by coupling with a ¹⁵N-labeled compound, and a linear relationship was found between the determined ¹⁵N/¹⁴N peak ratio and ¹⁵N atom% for the tested N species. The developed method was applied to bacterial cultures to examine their N-transformation reactions, enabling us to observe the occurrence of NO₂^- oxidation and NO₃^- reduction in a hypoxic Nitrobacter winogradskyi culture. IMPORTANCE ¹⁵N/¹⁴N analysis for aqueous N species is a powerful tool for obtaining insights into the global N cycle, but the procedure is cumbersome and laborious. The combined use of colorimetric reagents and MALDI-TOF MS, designated color MALDI-TOF MS, enabled us to determine the ¹⁵N atom% of common aqueous N species without laborious sample preparation and chromatographic separation steps; for instance, the ¹⁵N atom% of NO₂^- can be determined from >1,000 liquid samples daily at <$1 (U.S.) per 384 samples for routine analysis. This convenient MS method is a powerful tool that will advance our ability to explore the N-transformation reactions that occur in various environments and biological samples.

Functional Traits Resolve Mechanisms Governing the Assembly and Distribution of Nitrogen-Cycling Microbial Communities in the Global Ocean

Microorganisms drive much of the marine nitrogen (N) cycle, which jointly controls the primary production in the global ocean. However, our understanding of the microbial communities driving the global ocean N cycle remains fragmented. Focusing on “who is doing what, where, and how?”, this study draws a clear picture describing the global biogeography of marine N-cycling microbial communities by utilizing the Tara Oceans shotgun metagenomes. The marine N-cycling communities are highly variable taxonomically but relatively even at the functional trait level, showing clear functional redundancy properties. The functional traits and taxonomic groups are shaped by the same set of geo-environmental factors, among which, depth is the major factor impacting marine N-cycling communities, differentiating mesopelagic from epipelagic communities. Latitudinal diversity gradients and distance-decay relationships are observed for taxonomic groups, but rarely or weakly for functional traits. The composition of functional traits is strongly deterministic as revealed by null model analysis, while a higher degree of stochasticity is observed for taxonomic composition. Integrating multiple lines of evidence, in addition to drawing a biogeographic picture of marine N-cycling communities, this study also demonstrated an essential microbial ecological theory—determinism governs the assembly of microbial communities performing essential biogeochemical processes; the environment selects functional traits rather than taxonomic groups; functional redundancy underlies stochastic taxonomic community assembly.

Communality in microbial stress response and differential metabolic interactions revealed by time-series analysis of sponge symbionts

The diversity and function of sponge-associated symbionts is now increasingly understood, however, we lack an understanding on how they dynamically behave to ensure holobiont stability in the face of environmental variation. Here we performed a metatransciptomics analysis of three microbial symbionts of the sponge Cymbastela concentrica in situ over 14 months and through differential gene expression and correlation analysis to environmental variables uncovered differences that speak to their metabolic activities and level of symbiotic and environmental interactions. The nitrite-oxidising Ca. Porinitrospira cymbastela maintained a seemingly stable metabolism, with the few differentially expressed genes related only to stress responses. The heterotrophic Ca. Porivivens multivorans displayed differential use of holobiont-derived compounds and respiration modes, while the ammonium-oxidising archaeon Ca. Nitrosopumilus cymbastelus differentially expressed genes related to phosphate metabolism and symbiosis effectors. One striking similarity between the symbionts was their similar variation in expression of stress-related genes. Our timeseries study showed that the microbial community of C. concentrica undertakes dynamic gene expression adjustments in response to the surroundings, tuned to deal with general stress and metabolic interactions between holobiont members. The success of these dynamic adjustments likely underpins the stability of the sponge holobiont and may provide resilience against environmental change. This article is protected by copyright. All rights reserved.

Microbial methane emissions from the non-methanogenesis processes: A critical review

Methane, a potent greenhouse gas of global importance, has traditionally been considered as an end product of microbial methanogenesis of organic matter. Paradoxically, growing evidence has shown that some microbes, such as cyanobacteria, algae, fungi, purple non-sulfur bacteria, and cryptogamic covers, produce methane in oxygen-saturated aquatic and terrestrial ecosystems. The non-methanogenesis process could be an important potential contributor to methane emissions. This systematic review summarizes the knowledge of microorganisms involved in the non-methanogenesis process and the possible mechanisms of methane formation. Cyanobacteria-derived methane production may be attributed to either demethylation of methyl phosphonates or linked to light-driven primary productivity, while algae produce methane by utilizing methylated sulfur compounds as possible carbon precursors. In addition, fungi produce methane by utilizing methionine as a possible carbon precursor, and purple non-sulfur bacteria reduce carbon dioxide to methane by nitrogenase. The microbial methane distribution from the non-methanogenesis processes in aquatic and terrestrial environments and its environmental significance to global methane emissions, possible mechanisms of methane production in each open water, water-to-air methane fluxes, and the impact of climate change on microorganisms are also discussed. Finally, future perspectives are highlighted, such as establishing more in-situ experiments, quantifying methane flux through optimizing empirical models, distinguishing individual methane sources, and investigating nitrogenase-like enzyme systems to improve our understanding of microbial methane emission from the non-methanogenesis process.

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.

Transcriptomic evidences for microbial carbon and nitrogen cycles in the deoxygenated seawaters of Bohai Sea

Eutrophication-induced water deoxygenation occurs continually in coastal oceans, and alters community structure, metabolic processes, and the energy shunt, resulting in a major threat to the ecological environment. Seasonal deoxygenation events have occurred in the Bohai Sea (China), however, how these affect the functional activity of microorganisms remains unclear. Here, through the use of absolute quantification of 16S rRNA genes amplicon sequencing and metatranscriptomics approaches, we investigated the structure of the microbial community and the patterns of transcriptional activity in deoxygenated seawaters. The dominant phyla were Proteobacteria (average value, 1.4 × 10⁶ copies ml^-1), Cyanobacteria (3.7 × 10⁵ copies ml^-1), Bacteroidetes (2.7 × 10⁵ copies ml^-1), and the ammonia-oxidizing archaea Thaumarchaeota (1.9 × 10⁵ copies ml^-1). Among the various environmental factors, dissolved oxygen, pH and temperature displayed the most significant correlation with microbial community composition and functional activity. Metatranscriptomic data showed high transcriptional activity of Thaumarchaeota in the deoxygenated waters, with a significant increase in the expression of core genes representing ammonia oxidation, ammonia transport, and carbon fixation (3-hydroxypropionic acid/4-hydroxybutyric acid cycle) pathways. The transcripts of Cyanobacteria involved in photosynthesis and carbon fixation (Calvin-Benson-Bassham cycle) significantly decreased in low oxygen waters. Meanwhile, the transcripts for the ribulose bisphosphate carboxylase-encoding gene shifted from being assigned to photoautotrophic to chemoautotrophic organisms in surface and bottom waters, respectively. Moreover, the transcription profile indicated that heterotrophs play a critical role in transforming low-molecular-weight dissolved organic nitrogen. Elevated abundances of transcripts related to microbial antioxidant activity corresponded to an enhanced aerobic metabolism of Thaumarchaeota in the low oxygen seawater. In general, our transcriptional evidences showed a population increase of Thaumarchaeota, especially the coastal ecotype of ammonia oxidizers, in low oxygen aquatic environments, and indicated an enhanced contribution of chemolithoautotrophic carbon fixation to carbon flow.

Integrating experiments with system-level biogeochemical modeling to understand nitrogen cycling of reservoir sediments at elevated hydrostatic pressure

Impoundment of rivers to construct reservoirs for hydropower and irrigation greatly increase the hydrostatic pressure acting on river sediments with potential repercussions for ecosystem-level microbial activity and metabolism. Understanding the functioning and responses of key biogeochemical cycles such as that of nitrogen cycling to shifting hydrostatic pressure is needed to estimate and predict the systemic nutrient dynamics in deep-water reservoirs. We studied the functioning of bacterial communities involved in nitrogen transformation in bioreactors maintained under contrasting hydrostatic pressures (0.5 MPa-3.0 MPa) and complemented the experimental approach with a functional gene-informed biogeochemical model. The model predictions were broadly consistent with observations from the experiment, suggesting that the rates of N₂O production decreased while the sediment concentration of nitrite increased significantly with increasing pressure, at least when exceeding 1.0 MPa. Changes in nitrite reduction (nirS) and aerobic ammonia oxidation (amoA) genes abundances were in accordance with the observed changes in N₂O production and nitrite levels. Moreover, the model predicted that the higher pressures (P > 1.5 MPa) would intensify the inhibition of N₂ production via denitrification and result in an accumulation of ammonia in the sediment along with a decrease in dissolved oxygen. The results imply that increased hydrostatic pressure caused by dam constructions may have a strong effect on microbial nitrogen conversion, and that this may result in lower nitrogen removal.

Linking 16S rRNA Gene Classification to amoA Gene Taxonomy Reveals Environmental Distribution of Ammonia-Oxidizing Archaeal Clades in Peatland Soils

A highly resolved taxonomy for ammonia-oxidizing archaea (AOA) based on the alpha subunit of ammonia monooxygenase (amoA) was recently established, which uncovered novel environmental patterns of AOA, challenging previous generalizations. However, many microbiome studies target the 16S rRNA gene as a marker; thus, the usage of this novel taxonomy is currently limited. Here, we exploited the phylogenetic congruence of archaeal amoA and 16S rRNA genes to link 16S rRNA gene classification to the novel amoA taxonomy. We screened publicly available archaeal genomes and contigs for the co-occurring amoA and 16S rRNA genes and constructed a 16S rRNA gene database with the corresponding amoA clade taxonomy. Phylogenetic trees of both marker genes confirmed congruence, enabling the identification of clades. We validated this approach with 16S rRNA gene amplicon data from peatland soils. We succeeded in linking 16S rRNA gene amplicon sequence variants belonging to the class Nitrososphaeria to seven different AOA (amoA) clades, including two of the most frequently detected clades (Nitrososphaerales γ and δ clades) for which no pure culture is currently available. Water status significantly impacted the distribution of the AOA clades as well as the whole AOA community structure, which was correlated with pH, nitrate, and ammonium, consistent with previous clade predictions. Our study emphasizes the need to distinguish among AOA clades with distinct ecophysiologies and environmental preferences, for a better understanding of the ecology of the globally abundant AOA. IMPORTANCE The recently established phylogeny of amoA provides a finer resolution than previous studies, allowing clustering of AOA beyond the order level and thus revealing novel clades. While the 16S rRNA gene is mostly appreciated in microbiome studies, this novel phylogeny is in limited use. Here, we provide an alternative path to identifying AOA with this novel and highly resolved amoA taxonomy by using 16S rRNA gene sequencing data. We constructed a 16S rRNA gene database with the associated amoA clade taxonomy based on their phylogenetic congruence. With this database, we were able to assign 16S rRNA gene amplicons from peatland soils to different AOA clades, with a level of resolution provided previously only by amoA phylogeny. As 16S rRNA gene amplicon sequencing is still widely employed in microbiome studies, our database may have a broad application for interpreting the ecology of globally abundant AOA.

Nitrosopumilus zosterae sp. nov., an autotrophic ammonia-oxidizing archaeon of phylum Thaumarchaeota isolated from coastal eelgrass sediments of Japan

A novel mesophilic and aerobic ammonia-oxidizing archaeon of the phylum Thaumarchaeota, strain NM25^T, was isolated from coastal eelgrass zone sediment sampled in Shimoda (Japan). The cells were rod-shaped with an S-layer cell wall. The temperature range for growth was 20-37 °C, with an optimum at 30 °C. The pH range for growth was pH 6.1-7.7, with an optimum at pH 7.1. The salinity range for growth was 5-40 %, with an optimum range of 15-32 %. Cells obtained energy from ammonia oxidation and used bicarbonate as a carbon source. Utilization of urea was not observed for energy generation and growth. Strain NM25^T required a hydrogen peroxide scavenger, such as α-ketoglutarate, pyruvate or catalase, for sustained growth on ammonia. Growth of strain NM25^T was inhibited by addition of low concentrations of some organic compounds and organic mixtures, including complete inhibition by glycerol, peptone and yeast extract. Phylogenetic analysis of four concatenated housekeeping genes (16S rRNA, rpoB, rpsI and atpD) and concatenated AmoA, AmoB, AmoC amino acid sequences indicated that the isolate is similar to members of the genus Nitrosopumilus. The closest relative is Nitrosopumilus ureiphilus PS0^T with sequence similarities of 99.5 % for the 16S rRNA gene and 97.2 % for the amoA gene. Genome relatedness between strain NM25^T and N. ureiphilus PS0^T was assessed by average nucleotide identity and digital DNA-DNA hybridization, giving results of 85.4 and 40.2 %, respectively. On the basis of phenotypic, genotypic and phylogenetic data, strain NM25^T represents a novel species of the genus Nitrosopumilus, for which the name sp. nov, is proposed. The type strain is NM25^T (=NBRC 111181^T=ATCC TSD-147^T).

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.

Estuarine microbial diversity and nitrogen cycling increase along sand-mud gradients independent of salinity and distance

Estuaries are depositional environments prone to terrigenous mud sedimentation. While macrofaunal diversity and nitrogen retention are greatly affected by changes in sedimentary mud content, its impact on prokaryotic diversity and nitrogen cycling activity remains understudied. We characterized the composition of estuarine tidal flat prokaryotic communities spanning a habitat range from sandy to muddy sediments, while controlling for salinity and distance. We also determined the diversity, abundance and expression of ammonia oxidizers and N₂ O-reducers within these communities by amoA and clade I nosZ gene and transcript analysis. Results show that prokaryotic communities and nitrogen cycling fractions were sensitive to changes in sedimentary mud content, and that changes in the overall community were driven by a small number of phyla. Significant changes occurred in prokaryotic communities and N₂ O-reducing fractions with only a 3% increase in mud, while thresholds for ammonia oxidizers were less distinct, suggesting other factors are also important for structuring these guilds. Expression of nitrogen cycling genes was substantially higher in muddier sediments, and results indicate that the potential for coupled nitrification-denitrification became increasingly prevalent as mud content increased. Altogether, results demonstrate that mud content is a strong environmental driver of diversity and N-cycling dynamics in estuarine microbial communities.

Proteomic response of the marine ammonia-oxidising archaeon Nitrosopumilus maritimus to iron limitation reveals strategies to compensate for nutrient scarcity

Dissolved iron (Fe) is vanishingly low in the oceans, with ecological success conferred to microorganisms that can restructure their biochemistry to maintain high growth rates during Fe scarcity. Chemolithoautotrophic ammonia-oxidising archaea (AOA) are highly abundant in the oceans, constituting ~30% of cells below the photic zone. Here we examine the proteomic response of the AOA isolate Nitrosopumilus maritimus to growth-limiting Fe concentrations. Under Fe limitation, we observed a significant reduction in the intensity of Fe-dense ferredoxins associated with respiratory complex I whilst complex III and IV proteins with more central roles in the electron transport chain remain unchanged. We concomitantly observed an increase in the intensity of Fe-free functional alternatives such as flavodoxin and plastocyanin, thioredoxin and alkyl hydroperoxide which are known to mediate electron transport and reactive oxygen species detoxification, respectively. Under Fe limitation, we found a marked increase in the intensity of the ABC phosphonate transport system (Phn), highlighting an intriguing link between Fe and P cycling in N. maritimus. We hypothesise that an elevated uptake of exogenous phosphonates under Fe limitation may either supplement N. maritimus' endogenous methylphosphonate biosynthesis pathway - which requires Fe - or enhance the production of phosphonate-containing exopolysaccharides known to efficiently bind environmental Fe.

Genome-enabled exploration of microbial ecology and evolution in the sea: a rising tide lifts all boats

As a young bacteriologist just launching my career during the early days of the 'microbial revolution' in the 1980s, I was fortunate to participate in some early discoveries, and collaborate in the development of cross-disciplinary methods now commonly referred to as "metagenomics". My early scientific career focused on applying phylogenetic and genomic approaches to characterize 'wild' bacteria, archaea and viruses in their natural habitats, with an emphasis on marine systems. These central interests have not changed very much for me over the past three decades, but knowledge, methodological advances and new theoretical perspectives about the microbial world certainly have. In this invited 'How we did it' perspective, I trace some of the trajectories of my lab's collective efforts over the years, including phylogenetic surveys of microbial assemblages in marine plankton and sediments, development of microbial community gene- and genome-enabled surveys, and application of genome-guided, cultivation-independent functional characterization of novel enzymes, pathways and their relationships to in situ biogeochemistry. Throughout this short review, I attempt to acknowledge, all the mentors, students, postdocs and collaborators who enabled this research. Inevitably, a brief autobiographical review like this cannot be fully comprehensive, so sincere apologies to any of my great colleagues who are not explicitly mentioned herein. I salute you all as well!

A framework based on fundamental biochemical principles to engineer microbial community dynamics

Microbial communities are complex but there are basic principles we can apply to constrain the assumed stochasticity of their activity. By understanding the trade-offs behind the kinetic parameters that define microbial growth, we can explain how local interspecies dependencies arise and shape the emerging properties of a community. If we integrate these theoretical descriptions with experimental 'omics' data and bioenergetics analysis of specific environmental conditions, predictions on activity, assembly and spatial structure can be obtained reducing the a priori unpredictable complexity of microbial communities. This information can be used to define the appropriate selective pressures to engineer bioprocesses and propose new hypotheses which can drive experimental research to accelerate innovation in biotechnology.

Exploring Marine Planktonic Archaea: Then and Now

In 1977, Woese and Fox leveraged molecular phylogenetic analyses of ribosomal RNAs and identified a new microbial domain of life on Earth, the Archaebacteria (now known as Archaea). At the time of their discovery, only one archaebacterial group, the strictly anaerobic methanogens, was known. But soon, other phenotypically unrelated microbial isolates were shown to belong to the Archaea, many originating from extreme habitats, including extreme halophiles, extreme thermophiles, and thermoacidophiles. Since most Archaea seemed to inhabit extreme or strictly anoxic habitats, it came as a surprise in 1992 when two new lineages of archaea were reported to be abundant in oxygen rich, temperate marine coastal waters and the deep ocean. Since that time, studies of marine planktonic archaea have revealed many more surprises, including their unexpected ubiquity, unusual symbiotic associations, unpredicted physiologies and biogeochemistry, and global abundance. In this Perspective, early work conducted on marine planktonic Archaea by my lab group and others is discussed in terms of the relevant historical context, some of the original research motivations, and surprises and discoveries encountered along the way.

Impact of Terrestrial Input on Deep-Sea Benthic Archaeal Community Structure in South China Sea Sediments

Archaea are widespread in marine sediments and play important roles in the cycling of sedimentary organic carbon. However, factors controlling the distribution of archaea in marine sediments are not well understood. Here we investigated benthic archaeal communities over glacial-interglacial cycles in the northern South China Sea and evaluated their responses to sediment organic matter sources and inter-species interactions. Archaea in sediments deposited during the interglacial period Marine Isotope Stage (MIS) 1 (Holocene) were significantly different from those in sediments deposited in MIS 2 and MIS 3 of the Last Glacial Period when terrestrial input to the South China Sea was enhanced based on analysis of the long-chain n-alkane C₃₁. The absolute archaeal 16S rRNA gene abundance in subsurface sediments was highest in MIS 2, coincident with high sedimentation rates and high concentrations of total organic carbon. Soil Crenarchaeotic Group (SCG; Nitrososphaerales) species, the most abundant ammonia-oxidizing archaea in soils, increased dramatically during MIS 2, likely reflecting transport of terrestrial archaea during glacial periods with high sedimentation rates. Co-occurrence network analyses indicated significant association of SCG archaea with benthic deep-sea microbes such as Bathyarchaeota and Thermoprofundales in MIS 2 and MIS 3, suggesting potential interactions among these archaeal groups. Meanwhile, Thermoprofundales abundance was positively correlated with total organic carbon (TOC), along with n-alkane C₃₁ and sedimentation rate, indicating that Thermoprofundales may be particularly important in processing of organic carbon in deep-sea sediments. Collectively, these results demonstrate that the composition of heterotrophic benthic archaea in the South China Sea may be influenced by terrestrial organic input in tune with glacial-interglacial cycles, suggesting a plausible link between global climate change and microbial population dynamics in deep-sea marine sediments.

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.