In-depth Spatiotemporal Characterization of Planktonic Archaeal and Bacterial Communities in North and South San Francisco Bay.

Despite being the largest estuary on the west coast of North America, no in-depth survey of microbial communities in San Francisco Bay (SFB) waters currently exists. In this study, we analyze bacterioplankton and archaeoplankton communities at several taxonomic levels and spatial extents (i.e., North versus South Bay) to reveal patterns in alpha and beta diversity. We assess communities using high-throughput sequencing of the 16S rRNA gene in 177 water column samples collected along a 150-km transect over a 2-year monthly time-series. In North Bay, the microbial community is strongly structured by spatial salinity changes while in South Bay seasonal variations dominate community dynamics. Along the steep salinity gradient in North Bay, we find that operational taxonomic units (OTUs; 97% identity) have higher site specificity than at coarser taxonomic levels and turnover ("species" replacement) is high, revealing a distinct brackish community (in oligo-, meso-, and polyhaline samples) from fresh and marine end-members. At coarser taxonomic levels (e.g., phylum, class), taxa are broadly distributed across salinity zones (i.e., present/abundant in a large number of samples) and brackish communities appear to be a mix of fresh and marine communities. We also observe variations in brackish communities between samples with similar salinities, likely related to differences in water residence times between North and South Bay. Throughout SFB, suspended particulate matter is positively correlated with richness and influences changes in beta diversity. Within several abundant groups, including the SAR11 clade (comprising up to 30% of reads in a sample), OTUs appear to be specialized to a specific salinity range. Some other organisms also showed pronounced seasonal abundance, including Synechococcus, Ca. Actinomarina, and Nitrosopumilus-like OTUs. Overall, this study represents the first in-depth spatiotemporal survey of SFB microbial communities and provides insight into how planktonic microorganisms have specialized to different niches along the salinity gradient.

Diverse Thaumarchaeota Dominate Subsurface Ammonia-oxidizing Communities in Semi-arid Floodplains in the Western United States.

Subsurface microbial communities mediate biogeochemical transformations that drive both local and ecosystem-level cycling of essential elements, including nitrogen. However, their study has been largely limited to the deep ocean, terrestrial mines, caves, and topsoils (< 30 cm). Here, we present regional insights into the microbial ecology of aerobic ammonia oxidation within the terrestrial subsurface of five semi-arid riparian sites spanning a 900-km N-S transect. We sampled sediments, profiled communities to depths of ≤ 10 m, and compared them to reveal trends regionally within and surrounding the Upper Colorado River Basin (CRB). The diversity and abundance of ammonia-oxidizing microbial communities were evaluated in the context of subsurface geochemistry by applying a combination of amoA (encoding ammonia monooxygenase subunit A) gene sequencing, quantitative PCR, and geochemical techniques. Analysis of 898 amoA sequences from ammonia-oxidizing archaea (AOA) and bacteria (AOB) revealed extensive ecosystem-scale diversity, including archaeal amoA sequences from four of the five major AOA lineages currently found worldwide as well as distinct AOA ecotypes associated with naturally reduced zones (NRZs) and hydrogeochemical zones (unsaturated, capillary fringe, and saturated). Overall, AOA outnumber AOB by 2- to 5000-fold over this regional scale, suggesting that AOA may play a prominent biogeochemical role in nitrification within terrestrial subsurface sediments.

Depth distributions of nitrite reductase (nirK) gene variants reveal spatial dynamics of thaumarchaeal ecotype populations in coastal Monterey Bay.

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are key players in nutrient cycling, yet large gaps remain in our understanding of their ecology and metabolism. Despite multiple lines of evidence pointing to a central role for copper-containing nitrite reductase (NirK) in AOA metabolism, the thaumarchaeal nirK gene is rarely studied in the environment. In this study, we examine the diversity of nirK in the marine pelagic environment, in light of previously described ecological patterns of pelagic thaumarchaeal populations. Phylogenetic analyses show that nirK better resolves diversification patterns of marine Thaumarchaeota, compared to the conventionally used marker gene amoA. Specifically, we demonstrate that the three major phylogenetic clusters of marine nirK correspond to the three 'ecotype' populations of pelagic Thaumarchaeota. In this context, we further examine the relative distributions of the three variant groups in metagenomes and metatranscriptomes representing two depth profiles in coastal Monterey Bay. Our results reveal that nirK effectively tracks the dynamics of thaumarchaeal ecotype populations, particularly finer-scale diversification patterns within major lineages. We also find evidence for multiple copies of nirK per genome in a fraction of thaumarchaeal cells in the water column, which must be taken into account when using it as a molecular marker.

Deep nirS amplicon sequencing of San Francisco Bay sediments enables prediction of geography and environmental conditions from denitrifying community composition.

Denitrification is a dominant nitrogen loss process in the sediments of San Francisco Bay. In this study, we sought to understand the ecology of denitrifying bacteria by using next-generation sequencing (NGS) to survey the diversity of a denitrification functional gene, nirS (encoding cytchrome-cd1 nitrite reductase), along the salinity gradient of San Francisco Bay over the course of a year. We compared our dataset to a library of nirS sequences obtained previously from the same samples by standard PCR cloning and Sanger sequencing, and showed that both methods similarly demonstrated geography, salinity and, to a lesser extent, nitrogen, to be strong determinants of community composition. Furthermore, the depth afforded by NGS enabled novel techniques for measuring the association between environment and community composition. We used Random Forests modelling to demonstrate that the site and salinity of a sample could be predicted from its nirS sequences, and to identify indicator taxa associated with those environmental characteristics. This work contributes significantly to our understanding of the distribution and dynamics of denitrifying communities in San Francisco Bay, and provides valuable tools for the further study of this key N-cycling guild in all estuarine systems.

Controls of nitrogen cycling evaluated along a well-characterized climate gradient.

The supply of nitrogen (N) constrains primary productivity in many ecosystems, raising the question "what controls the availability and cycling of N"? As a step toward answering this question, we evaluated N cycling processes and aspects of their regulation on a climate gradient on Kohala Volcano, Hawaii, USA. The gradient extends from sites receiving <300 mm/yr of rain to those receiving >3,000 mm/yr, and the pedology and dynamics of rock-derived nutrients in soils on the gradient are well understood. In particular, there is a soil process domain at intermediate rainfall within which ongoing weathering and biological uplift have enriched total and available pools of rock-derived nutrients substantially; sites at higher rainfall than this domain are acid and infertile as a consequence of depletion of rock-derived nutrients, while sites at lower rainfall are unproductive and subject to wind erosion. We found elevated rates of potential net N mineralization in the domain where rock-derived nutrients are enriched. Higher-rainfall sites have low rates of potential net N mineralization and high rates of microbial N immobilization, despite relatively high rates of gross N mineralization. Lower-rainfall sites have moderately low potential net N mineralization, relatively low rates of gross N mineralization, and rates of microbial N immobilization sufficient to sequester almost all the mineral N produced. Bulk soil δ15 N also varied along the gradient, from +4‰ at high rainfall sites to +14‰ at low rainfall sites, indicating differences in the sources and dynamics of soil N. Our analysis shows that there is a strong association between N cycling and soil process domains that are defined using soil characteristics independent of N along this gradient, and that short-term controls of N cycling can be understood in terms of the supply of and demand for N.

Spatiotemporal Characterization of San Francisco Bay Denitrifying Communities: a Comparison of nirK and nirS Diversity and Abundance.

Denitrifying bacteria play a critical role in the estuarine nitrogen cycle. Through the transformation of nitrate into nitrogen gas, these organisms contribute to the loss of bioavailable (i.e., fixed) nitrogen from low-oxygen environments such as estuary sediments. Denitrifiers have been shown to vary in abundance and diversity across the spatial environmental gradients that characterize estuaries, such as salinity and nitrogen availability; however, little is known about how their communities change in response to temporal changes in those environmental properties. Here, we present a 1-year survey of sediment denitrifier communities along the estuarine salinity gradient of San Francisco Bay. We used quantitative PCR and sequencing of functional genes coding for a key denitrifying enzyme, dissimilatory nitrite reductase, to compare two groups of denitrifiers: those with nirK (encoding copper-dependent nitrite reductase) and those with nirS (encoding the cytochrome-cd 1-dependent variant). We found that nirS was consistently more abundant and more diverse than nirK in all parts of the estuary. The abundances of the two genes were tightly linked across space but differed temporally, with nirK peaking when temperature was low and nirS peaking when nitrate was high. Likewise, the diversity and composition of nirK- versus nirS-type communities differed in their responses to seasonal variations, though both were strongly determined by site. Furthermore, our sequence libraries detected deeply branching clades with no cultured isolates, evidence of enormous diversity within the denitrifiers that remains to be explored.

Spatiotemporal relationships between the abundance, distribution, and potential activities of ammonia-oxidizing and denitrifying microorganisms in intertidal sediments.

The primary objective of this study was to gain an understanding of how key microbial communities involved in nitrogen cycling in estuarine sediments vary over a 12-month period. Furthermore, we sought to determine whether changes in the size of these communities are related to, or indicative of, seasonal patterns in fixed nitrogen dynamics in Elkhorn Slough--a small, agriculturally impacted estuary with a direct connection to Monterey Bay. We assessed sediment and pore water characteristics, abundance of functional genes for nitrification (bacterial and archaeal amoA, encoding ammonia monooxygenase subunit A) and denitrification (nirS and nirK, encoding nitrite reductase), and measurements of potential nitrification and denitrification activities at six sites. No seasonality in the abundance of denitrifier or ammonia oxidizer genes was observed. A strong association between potential nitrification activity and the size of ammonia-oxidizing bacterial communities was observed across the estuary. In contrast, ammonia-oxidizing archaeal abundances remained relatively constant in space and time. Unlike many other estuaries, salinity does not appear to regulate the distribution of ammonia-oxidizing communities in Elkhorn Slough. Instead, their distributions appear to be governed over two different time scales. Long-term niche characteristics selected for the gross size of archaeal and bacterial ammonia-oxidizing communities, yet covariation in their abundances between monthly samples suggests that they respond in a similar manner to short-term changes in their environment. Abundances of denitrifier and ammonia oxidizer genes also covaried, but site-specific differences in this relationship suggest differing levels of interaction (or coupling) between nitrification and denitrification.

Microbial biogeography across a full-scale wastewater treatment plant transect: evidence for immigration between coupled processes.

Wastewater treatment plants use a variety of bioreactor types and configurations to remove organic matter and nutrients. Little is known regarding the effects of different configurations and within-plant immigration on microbial community dynamics. Previously, we found that the structure of ammonia-oxidizing bacterial (AOB) communities in a full-scale dispersed growth activated sludge bioreactor correlated strongly with levels of NO2 (-) entering the reactor from an upstream trickling filter. Here, to further examine this puzzling association, we profile within-plant microbial biogeography (spatial variation) and test the hypothesis that substantial microbial immigration occurs along a transect (raw influent, trickling filter biofilm, trickling filter effluent, and activated sludge) at the same full-scale wastewater treatment plant. AOB amoA gene abundance increased >30-fold between influent and trickling filter effluent concomitant with NO2 (-) production, indicating unexpected growth and activity of AOB within the trickling filter. Nitrosomonas europaea was the dominant AOB phylotype in trickling filter biofilm and effluent, while a distinct "Nitrosomonas-like" lineage dominated in activated sludge. Prior time series indicated that this "Nitrosomonas-like" lineage was dominant when NO2 (-) levels in the trickling filter effluent (i.e., activated sludge influent) were low, while N. europaea became dominant in the activated sludge when NO2 (-) levels were high. This is consistent with the hypothesis that NO2 (-) production may cooccur with biofilm sloughing, releasing N. europaea from the trickling filter into the activated sludge bioreactor. Phylogenetic microarray (PhyloChip) analyses revealed significant spatial variation in taxonomic diversity, including a large excess of methanogens in the trickling filter relative to activated sludge and attenuation of Enterobacteriaceae across the transect, and demonstrated transport of a highly diverse microbial community via the trickling filter effluent to the activated sludge bioreactor. Our results provide compelling evidence that substantial immigration between coupled process units occurs and may exert significant influence over microbial community dynamics within staged bioreactors.

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 (SFB). 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 SFB during the autumn of 2018 and seasonally throughout 2019. Nitrification rates were correlated with AOA abundance in qPCR 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 Ca. Nitrosomarinus catalina-like AOA metagenome-assembled genome (MAG) that had high transcript abundance during the bloom and expressed >80% of genes in its genome. We also recovered a putative nitrite-oxidizing bacteria (NOB) MAG 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 MAGs. 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.

Diverse and unconventional methanogens, methanotrophs, and methylotrophs in metagenome-assembled genomes from subsurface sediments of the Slate River floodplain, Crested Butte, CO, USA.

We use metagenome-assembled genomes (MAGs) to understand single-carbon (C1) compound-cycling-particularly methane-cycling-microorganisms in montane riparian floodplain sediments. We generated 1,233 MAGs (>50% completeness and <10% contamination) from 50- to 150-cm depth below the sediment surface capturing the transition between oxic, unsaturated sediments and anoxic, saturated sediments in the Slate River (SR) floodplain (Crested Butte, CO, USA). We recovered genomes of putative methanogens, methanotrophs, and methylotrophs (n = 57). Methanogens, found only in deep, anoxic depths at SR, originate from three different clades (Methanoregulaceae, Methanotrichaceae, and Methanomassiliicoccales), each with a different methanogenesis pathway; putative methanotrophic MAGs originate from within the Archaea (Candidatus Methanoperedens) in anoxic depths and uncultured bacteria (Ca. Binatia) in oxic depths. Genomes for canonical aerobic methanotrophs were not recovered. Ca. Methanoperedens were exceptionally abundant (~1,400× coverage, >50% abundance in the MAG library) in one sample that also contained aceticlastic methanogens, indicating a potential C1/methane-cycling hotspot. Ca. Methylomirabilis MAGs from SR encode pathways for methylotrophy but do not harbor methane monooxygenase or nitrogen reduction genes. Comparative genomic analysis supports that one clade within the Ca. Methylomirabilis genus is not methanotrophic. The genetic potential for methylotrophy was widespread, with over 10% and 19% of SR MAGs encoding a methanol dehydrogenase or substrate-specific methyltransferase, respectively. MAGs from uncultured Thermoplasmata archaea in the Ca. Gimiplasmatales (UBA10834) contain pathways that may allow for anaerobic methylotrophic acetogenesis. Overall, MAGs from SR floodplain sediments reveal a potential for methane production and consumption in the system and a robust potential for methylotrophy.IMPORTANCEThe cycling of carbon by microorganisms in subsurface environments is of particular relevance in the face of global climate change. Riparian floodplain sediments contain high organic carbon that can be degraded into C1 compounds such as methane, methanol, and methylamines, the fate of which depends on the microbial metabolisms present as well as the hydrological conditions and availability of oxygen. In the present study, we generated over 1,000 MAGs from subsurface sediments from a montane river floodplain and recovered genomes for microorganisms that are capable of producing and consuming methane and other C1 compounds, highlighting a robust potential for C1 cycling in subsurface sediments both with and without oxygen. Archaea from the Ca. Methanoperedens genus were exceptionally abundant in one sample, indicating a potential C1/methane-cycling hotspot in the Slate River floodplain system.

Adaptation of nitrifying microbial biomass to nickel in batch incubations.

Nitrification-microbial oxidation of ammonia to nitrate-is sensitive to an array of inhibitors. Currently, little is known regarding the ecological processes that enable adaptation to inhibitors and recovery of nitrification. This study evaluated inhibition and recovery of nitrification in batch cultures of activated sludge incubated with different levels of nickel (Ni), a model inhibitor. Incubation with 1 mg/L of added Ni did not adversely affect nitrification, and little inhibition occurred at 5 and 10 mg/L Ni. Incubation with 50 mg/L Ni resulted in significant inhibition, decreased amoA transcript abundance, and delayed recovery of nitrification until amoA transcript abundance rebounded after 24 h. For this dosage, recovery of nitrification occurred without a significant change in ammonia-oxidizing bacteria (AOB) community structure. By contrast, incubation with 150 mg/L of added Ni strongly inhibited nitrification and delayed recovery until a shift in AOB community structure occurred after ∼6 weeks of incubation. The results indicate that inhibitor-resistant nitrifying cultures can be obtained from long-term batch incubations of decaying activated sludge incubated with high levels of added inhibitor.

Combined niche and neutral effects in a microbial wastewater treatment community.

It has long been assumed that differences in the relative abundance of taxa in microbial communities reflect differences in environmental conditions. Here we show that in the economically and environmentally important microbial communities in a wastewater treatment plant, the population dynamics are consistent with neutral community assembly, where chance and random immigration play an important and predictable role in shaping the communities. Using dynamic observations, we demonstrate a straightforward calibration of a purely neutral model and a parsimonious method to incorporate environmental influence on the reproduction (or birth) rate of individual taxa. The calibrated model parameters are biologically plausible, with the population turnover and diversity in the heterotrophic community being higher than for the ammonia oxidizing bacteria (AOB) and immigration into AOB community being relatively higher. When environmental factors were incorporated more of the variance in the observations could be explained but immigration and random reproduction and deaths remained the dominant driver in determining the relative abundance of the common taxa. Consequently we suggest that neutral community models should be the foundation of any description of an open biological system.

Pelagic metagenome-assembled genomes from an estuarine salinity gradient in San Francisco Bay.

San Francisco Bay (SFB) is a large and highly human-impacted estuarine system. We produced 449 metagenome-assembled genomes from SFB waters, collected along the salinity gradient, providing a rich data set to compare the metabolic potential of microorganisms from different salinity zones within SFB and to other estuarine systems.

Ecophysiology and genomics of the brackish water adapted SAR11 subclade IIIa.

The Order Pelagibacterales (SAR11) is the most abundant group of heterotrophic bacterioplankton in global oceans and comprises multiple subclades with unique spatiotemporal distributions. Subclade IIIa is the primary SAR11 group in brackish waters and shares a common ancestor with the dominant freshwater IIIb (LD12) subclade. Despite its dominance in brackish environments, subclade IIIa lacks systematic genomic or ecological studies. Here, we combine closed genomes from new IIIa isolates, new IIIa MAGS from San Francisco Bay (SFB), and 460 highly complete publicly available SAR11 genomes for the most comprehensive pangenomic study of subclade IIIa to date. Subclade IIIa represents a taxonomic family containing three genera (denoted as subgroups IIIa.1, IIIa.2, and IIIa.3) that had distinct ecological distributions related to salinity. The expansion of taxon selection within subclade IIIa also established previously noted metabolic differentiation in subclade IIIa compared to other SAR11 subclades such as glycine/serine prototrophy, mosaic glyoxylate shunt presence, and polyhydroxyalkanoate synthesis potential. Our analysis further shows metabolic flexibility among subgroups within IIIa. Additionally, we find that subclade IIIa.3 bridges the marine and freshwater clades based on its potential for compatible solute transport, iron utilization, and bicarbonate management potential. Pure culture experimentation validated differential salinity ranges in IIIa.1 and IIIa.3 and provided detailed IIIa cell size and volume data. This study is an important step forward for understanding the genomic, ecological, and physiological differentiation of subclade IIIa and the overall evolutionary history of SAR11.

Genome sequence of "Candidatus Nitrosoarchaeum limnia" BG20, a low-salinity ammonia-oxidizing archaeon from the San Francisco Bay estuary.

Here, we present the draft genome sequence of "Candidatus Nitrosoarchaeum limnia" BG20, an ammonia-oxidizing archaeon enriched in culture from low-salinity sediments of the San Francisco Bay estuary. The genome sequence revealed many similarities to the previously sequenced genome of "Ca. Nitrosoarchaeum limnia" SFB1 (enriched from a nearby site in San Francisco Bay) and is representative of a clade of ammonia-oxidizing archaea (AOA) found in low-salinity habitats worldwide.

Genome sequence of "Candidatus Nitrosopumilus salaria" BD31, an ammonia-oxidizing archaeon from the San Francisco Bay estuary.

Ammonia-oxidizing archaea (AOA) play important roles in nitrogen and carbon cycling in marine and terrestrial ecosystems. Here, we present the draft genome sequence for the ammonia-oxidizing archaeon "Candidatus Nitrosopumilus salaria" BD31, which was enriched in culture from sediments of the San Francisco Bay estuary. The genome sequences revealed many similarities to the genome of Nitrosopumilus maritimus.

Ecophysiology of an ammonia-oxidizing archaeon adapted to low-salinity habitats.

Ammonia oxidation in marine and terrestrial ecosystems plays a pivotal role in the cycling of nitrogen and carbon. Recent discoveries have shown that ammonia-oxidizing archaea (AOA) are both abundant and diverse in these systems, yet very little is known about their physiology. Here we report a physiological analysis of a novel low-salinity-type AOA enriched from the San Francisco Bay estuary, Candidatus Nitrosoarchaeum limnia strain SFB1. N. limnia has a slower growth rate than Nitrosopumilus maritimus and Nitrososphaera viennensis EN76, the only pure AOA isolates described to date, but the growth rate is comparable to the growth of marine AOA enrichment cultures. The growth rate only slightly decreased when N. limnia was grown under lower-oxygen conditions (5.5 % oxygen in the headspace). Although N. limnia was capable of growth at 75 % of seawater salinity, there was a longer lag time, incomplete oxidation of ammonia to nitrite, and slower overall growth rate. Allylthiourea (ATU) only partially inhibited growth and ammonia oxidation by N. limnia at concentrations known to completely inhibit bacterial ammonia oxidation. Using electron microscopy, we confirmed the presence of flagella as suggested by various flagellar biosynthesis genes in the N. limnia genome. We demonstrate that N. limnia is representative of a low-salinity estuarine AOA ecotype and that more than 85 % of its proteins have highest identity to other coastal and estuarine metagenomic sequences. Our findings further highlight the physiology of N. limnia and help explain its ecological adaptation to low-salinity niches.

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.

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

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

Spatial variability in nitrification rates and ammonia-oxidizing microbial communities in the agriculturally impacted Elkhorn Slough estuary, California.

Ammonia oxidation-the microbial oxidation of ammonia to nitrite and the first step in nitrification-plays a central role in nitrogen cycling in coastal and estuarine systems. Nevertheless, questions remain regarding the connection between this biogeochemical process and the diversity and abundance of the mediating microbial community. In this study, we measured nutrient fluxes and rates of sediment nitrification in conjunction with the diversity and abundance of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing betaproteobacteria (ß-AOB). Sediments were examined from four sites in Elkhorn Slough, a small agriculturally impacted coastal California estuary that opens into Monterey Bay. Using an intact sediment core flowthrough incubation system, we observed significant correlations among NO(3)(-), NO(2)(-), NH(4)(+), and PO(4)(3+) fluxes, indicating a tight coupling of sediment biogeochemical processes. (15)N-based measurements of nitrification rates revealed higher rates at the less impacted, lower-nutrient sites than at the more heavily impacted, nutrient-rich sites. Quantitative PCR analyses revealed that ß-AOB amoA (encoding ammonia monooxygenase subunit A) gene copies outnumbered AOA amoA gene copies by factors ranging from 2- to 236-fold across the four sites. Sites with high nitrification rates primarily contained marine/estuarine Nitrosospira-like bacterial amoA sequences and phylogenetically diverse archaeal amoA sequences. Sites with low nitrification rates were dominated by estuarine Nitrosomonas-like amoA sequences and archaeal amoA sequences similar to those previously described in soils. This is the first report measuring AOA and ß-AOB amoA abundance in conjunction with (15)N-based nitrification rates in estuary sediments.