Export of Organic Carbon from Reduced Fine-Grained Zones Governs Biogeochemical Reactivity in a Simulated Aquifer.

Sediment interfaces in alluvial aquifers have a disproportionately large influence on biogeochemical activity and, therefore, on groundwater quality. Previous work showed that exports from fine-grained, organic-rich zones sustain reducing conditions in downstream coarse-grained aquifers beyond the influence of reduced aqueous products alone. Here, we show that sustained anaerobic activity can be attributed to the export of organic carbon, including live microorganisms, from fine-grained zones. We used a dual-domain column system with ferrihydrite-coated sand and embedded reduced, fine-grained lenses from Slate River (Crested Butte, CO) and Wind River (Riverton, WY) floodplains. After 50 d of groundwater flow, 8.8 ± 0.7% and 14.8 ± 3.1% of the total organic carbon exported from the Slate and Wind River lenses, respectively, had accumulated in the sand downstream. Furthermore, higher concentrations of dissolved Fe(II) and lower concentrations of dissolved organic carbon in the sand compared to total aqueous transport from the lenses suggest that Fe(II) was produced in situ by microbial oxidation of organic carbon coupled to iron reduction. This was further supported by an elevated abundance of 16S rRNA and iron-reducing (gltA) gene copies. These findings suggest that organic carbon transport across interfaces contributes to downstream biogeochemical reactions in natural alluvial aquifers.

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.

Oxidation of urea-derived nitrogen by thaumarchaeota-dominated marine nitrifying communities.

Urea nitrogen has been proposed to contribute significantly to nitrification by marine thaumarchaeotes. These inferences are based on distributions of thaumarchaeote urease genes rather than activity measurements. We found that ammonia oxidation rates were always higher than oxidation rates of urea-derived N in samples from coastal Georgia, USA (means ± SEM: 382 ± 35 versus 73 ± 24 nmol L-1  d-1 , Mann-Whitney U-test p < 0.0001), and the South Atlantic Bight (20 ± 8.8 versus 2.2 ± 1.7 nmol L-1  d-1 , p = 0.026) but not the Gulf of Alaska (8.8 ± 4.0 versus 1.5 ± 0.6, p > 0.05). Urea-derived N was relatively more important in samples from Antarctic continental shelf waters, though the difference was not statistically significant (19.4 ± 4.8 versus 12.0 ± 2.7 nmol L-1  d-1 , p > 0.05). We found only weak correlations between oxidation rates of urea-derived N and the abundance or transcription of putative Thaumarchaeota ureC genes. Dependence on urea-derived N does not appear to be directly related to pH or ammonium concentrations. Competition experiments and release of 15 NH3 suggest that urea is hydrolyzed to ammonia intracellularly, then a portion is lost to the dissolved pool. The contribution of urea-derived N to nitrification appears to be minor in temperate coastal waters, but may represent a significant portion of the nitrification flux in Antarctic coastal waters.

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.

Structural insights into bifunctional thaumarchaeal crotonyl-CoA hydratase and 3-hydroxypropionyl-CoA dehydratase from Nitrosopumilus maritimus.

The ammonia-oxidizing thaumarchaeal 3-hydroxypropionate/4-hydroxybutyrate (3HP/4HB) cycle is one of the most energy-efficient CO2 fixation cycles discovered thus far. The protein encoded by Nmar_1308 (from Nitrosopumilus maritimus SCM1) is a promiscuous enzyme that catalyzes two essential reactions within the thaumarchaeal 3HP/4HB cycle, functioning as both a crotonyl-CoA hydratase (CCAH) and 3-hydroxypropionyl-CoA dehydratase (3HPD). In performing both hydratase and dehydratase activities, Nmar_1308 reduces the total number of enzymes necessary for CO2 fixation in Thaumarchaeota, reducing the overall cost for biosynthesis. Here, we present the first high-resolution crystal structure of this bifunctional enzyme with key catalytic residues in the thaumarchaeal 3HP/4HB pathway.

Depth-Differentiation and Seasonality of Planktonic Microbial Assemblages in the Monterey Bay Upwelling System.

Coastal upwelling regions are hotspots of biological productivity, supporting diverse communities of microbial life and metabolisms. Monterey Bay (MB), a coastal ocean embayment in central California, experiences seasonal upwelling of cold, nutrient-rich waters that sustain episodes of high phytoplankton production in surface waters. While productivity in surface waters is intimately linked to metabolisms of diverse communities of Archaea and Bacteria, a comprehensive understanding of the microbial community in MB is missing thus far, particularly in relation to the distinct hydrographic seasons characteristic of the MB system. Here we present the results of a 2-year microbial time-series survey in MB, investigating community composition and structure across spatiotemporal gradients. In deciphering these patterns, we used unique sequence variants (SVs) of the 16S rRNA gene (V4-V5 region), complemented with metagenomes and metatranscriptomes representing multiple depth profiles. We found clear depth-differentiation and recurring seasonal abundance patterns within planktonic communities, particularly when analyzed at finer taxonomic levels. Compositional changes were more pronounced in the upper 0-40 m of the water column, whereas deeper depths were characterized by temporally stable populations. In accordance with the dynamic nutrient profiles, the system appears to change from a Bacteroidetes- and Rhodobacterales-dominated upwelling period to an oceanic season dominated by oligotrophic groups such as SAR11 and picocyanobacteria. The cascade of environmental changes brought about by upwelling and relaxation events thus impacts microbial community structure in the bay, with important implications for the temporal variability of nutrient and energy fluxes within the MB ecosystem. Our observations emphasize the need for continued monitoring of planktonic microbial communities in order to predict and manage the behavior of this sensitive marine sanctuary ecosystem, over projected intensification of upwelling in the region.

Correction: Single-cell genomics shedding light on marine Thaumarchaeota diversification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Differential co-occurrence relationships shaping ecotype diversification within Thaumarchaeota populations in the coastal ocean water column.

Ecological factors contributing to depth-related diversification of marine Thaumarchaeota populations remain largely unresolved. To investigate the role of potential microbial associations in shaping thaumarchaeal ecotype diversification, we examined co-occurrence relationships in a community composition dataset (16S rRNA V4-V5 region) collected as part of a 2-year time series in coastal Monterey Bay. Ecotype groups previously defined based on functional gene diversity-water column A (WCA), water column B (WCB) and Nitrosopumilus-like clusters-were recovered in the thaumarchaeal 16S rRNA gene phylogeny. Networks systematically reflected depth-related patterns in the abundances of ecotype populations, suggesting thaumarchaeal ecotypes as keystone members of the microbial community below the euphotic zone. Differential environmental controls on the ecotype populations were further evident in subnetwork modules showing preferential co-occurrence of OTUs belonging to the same ecotype cluster. Correlated abundances of Thaumarchaeota and heterotrophic bacteria (e.g., Bacteroidetes, Marinimicrobia and Gammaproteobacteria) indicated potential reciprocal interactions via dissolved organic matter transformations. Notably, the networks recovered ecotype-specific associations between thaumarchaeal and Nitrospina OTUs. Even at depths where WCB-like Thaumarchaeota dominated, Nitrospina OTUs were found to preferentially co-occur with WCA-like and Nitrosopumilus-like thaumarchaeal OTUs, highlighting the need to investigate the ecological implications of the composition of nitrifier assemblages in marine waters.

Coastal Ocean Metagenomes and Curated Metagenome-Assembled Genomes from Marsh Landing, Sapelo Island (Georgia, USA).

Microbes play a dominant role in the biogeochemistry of coastal waters, which receive organic matter from diverse sources. We present metagenomes and 45 metagenome-assembled genomes (MAGs) from Sapelo Island, Georgia, to further understand coastal microbial populations. Notably, four MAGs are archaea, with two Thaumarchaeota and two marine group II Euryarchaeota.

Light and temperature control the seasonal distribution of thaumarchaeota in the South Atlantic bight.

Mid-summer peaks in the abundance of Thaumarchaeota and nitrite concentration observed on the Georgia, USA, coast could result from in situ activity or advection of populations from another source. We collected data on the distribution of Thaumarchaeota, ammonia-oxidizing betaproteobacteria (AOB), Nitrospina, environmental variables and rates of ammonia oxidation during six cruises in the South Atlantic Bight (SAB) from April to November 2014. These data were used to examine seasonality of nitrification in offshore waters and to test the hypothesis that the bloom was localized to inshore waters. The abundance of Thaumarchaeota marker genes (16S rRNA and amoA) increased at inshore and nearshore stations starting in July and peaked in August at >107 copies L-1. The bloom did not extend onto the mid-shelf, where Thaumarchaeota genes ranged from 103 to 105 copies L-1. Ammonia oxidation rates (AO) were highest at inshore stations during summer (to 840 nmol L-1 d-1) and were always at the limit of detection at mid-shelf stations. Nitrite concentrations were correlated with AO (R = 0.94) and were never elevated at mid-shelf stations. Gene sequences from samples collected at mid-shelf stations generated using Archaea 16S rRNA primers were dominated by Euryarchaeota; sequences from inshore and nearshore stations were dominated by Thaumarchaeota. Thaumarchaeota were also abundant at depth at the shelf-break; however, this population was phylogenetically distinct from the inshore/nearshore population. Our analysis shows that the bloom is confined to inshore waters during summer and suggests that Thaumarchaeota distributions in the SAB are controlled primarily by photoinhibition and secondarily by water temperature.

Nutrient transport suggests an evolutionary basis for charged archaeal surface layer proteins.

Surface layers (S-layers) are two-dimensional, proteinaceous, porous lattices that form the outermost cell envelope component of virtually all archaea and many bacteria. Despite exceptional sequence diversity, S-layer proteins (SLPs) share important characteristics such as their ability to form crystalline sheets punctuated with nano-scale pores, and their propensity for charged amino acids, leading to acidic or basic isoelectric points. However, the precise function of S-layers, or the role of charged SLPs and how they relate to cellular metabolism is unknown. Nano-scale lattices affect the diffusion behavior of low-concentration solutes, even if they are significantly smaller than the pore size. Here, we offer a rationale for charged S-layer proteins in the context of the structural evolution of S-layers. Using the ammonia-oxidizing archaea (AOA) as a model for S-layer geometry, and a 2D electrodiffusion reaction computational framework to simulate diffusion and consumption of the charged solute ammonium (NH4+), we find that the characteristic length scales of nanoporous S-layers elevate the concentration of NH4+ in the pseudo-periplasmic space. Our simulations suggest an evolutionary, mechanistic basis for S-layer charge and shed light on the unique ability of some AOA to oxidize ammonia in environments with nanomolar NH4+ availability, with broad implications for comparisons of ecologically distinct populations.

Integrated structural biology and molecular ecology of N-cycling enzymes from ammonia-oxidizing archaea.

Knowledge of the molecular ecology and environmental determinants of ammonia-oxidizing organisms is critical to understanding and predicting the global nitrogen (N) and carbon cycles, but an incomplete biochemical picture hinders in vitro studies of N-cycling enzymes. Although an integrative structural and dynamic characterization at the atomic scale would advance our understanding of function tremendously, structural knowledge of key N-cycling enzymes from ecologically relevant ammonia oxidizers is unfortunately extremely limited. Here, we discuss the challenges and opportunities for examining the ecology of ammonia-oxidizing organisms, particularly uncultivated Thaumarchaeota, through (meta)genome-driven structural biology of the enzymes ammonia monooxygenase (AMO) and nitrite reductase (NirK).

Contribution of ammonia oxidation to chemoautotrophy in Antarctic coastal waters.

There are few measurements of nitrification in polar regions, yet geochemical evidence suggests that it is significant, and chemoautotrophy supported by nitrification has been suggested as an important contribution to prokaryotic production during the polar winter. This study reports seasonal ammonia oxidation (AO) rates, gene and transcript abundance in continental shelf waters west of the Antarctic Peninsula, where Thaumarchaeota strongly dominate populations of ammonia-oxidizing organisms. Higher AO rates were observed in the late winter surface mixed layer compared with the same water mass sampled during summer (mean±s.e.: 62±16 versus 13±2.8 nm per day, t-test P<0.0005). AO rates in the circumpolar deep water did not differ between seasons (21±5.7 versus 24±6.6 nm per day; P=0.83), despite 5- to 20-fold greater Thaumarchaeota abundance during summer. AO rates correlated with concentrations of Archaea ammonia monooxygenase (amoA) genes during summer, but not with concentrations of Archaea amoA transcripts, or with ratios of Archaea amoA transcripts per gene, or with concentrations of Betaproteobacterial amoA genes or transcripts. The AO rates we report (<0.1-220 nm per day) are ~10-fold greater than reported previously for Antarctic waters and suggest that inclusion of Antarctic coastal waters in global estimates of oceanic nitrification could increase global rate estimates by ~9%. Chemoautotrophic carbon fixation supported by AO was 3-6% of annualized phytoplankton primary production and production of Thaumarchaeota biomass supported by AO could account for ~9% of the bacterioplankton production measured in winter. Growth rates of thaumarchaeote populations inferred from AO rates averaged 0.3 per day and ranged from 0.01 to 2.1 per day.

Single-cell genomics shedding light on marine Thaumarchaeota diversification.

Previous studies based on analysis of amoA, 16S ribosomal RNA or accA gene sequences have established that marine Thaumarchaeota fall into two phylogenetically distinct groups corresponding to shallow- and deep-water clades, but it is not clear how water depth interacts with other environmental factors, including light, temperature and location, to affect this pattern of diversification. Earlier studies focused on single-gene distributions were not able to link phylogenetic structure to other aspects of functional adaptation. Here, we analyzed the genome content of 46 uncultivated single Thaumarchaeota cells sampled from epi- and mesopelagic waters of subtropical, temperate and polar oceans. Phylogenomic analysis showed that populations diverged by depth, as expected, and that mesopelagic populations from different locations were well mixed. Functional analysis showed that some traits, including putative DNA photolyase and catalase genes that may be related to adaptive mechanisms to reduce light-induced damage, were found exclusively in members of the epipelagic clade. Our analysis of partial genomes has thus confirmed the depth differentiation of Thaumarchaeota populations observed previously, consistent with the distribution of putative mechanisms to reduce light-induced damage in shallow- and deep-water populations.

Seasonal variation in the metatranscriptomes of a Thaumarchaeota population from SE USA coastal waters.

We used a combination of metatranscriptomic analyses and quantitative PCR (qPCR) to study seasonal changes in Thaumarchaeota populations from a salt marsh-dominated estuary. Surface waters (0.5 m depth) were sampled quarterly at Marsh Landing, Sapelo Island, GA, USA over a 3-year period. We found a mid-summer peak in Thaumarchaeota abundance measured by qPCR of either 16S rRNA or amoA genes in each of the 3 years. Thaumarchaeota were 100-1000-fold more abundant during the peak than at other times of the year, whereas the abundance of ammonia- and nitrite-oxidizing Bacteria varied <10-fold over the same period. Analysis of the microdiversity of several highly transcribed genes in 20 metatranscriptomes from a 1-year subset of these samples showed that the transcriptionally active population consisted of 2 or 3 dominant phylotypes that differed between successive summers. This shift appeared to have begun during the preceding winter and spring. Transcripts from the same genes dominated the Thaumarchaeota mRNA pool throughout the year, with genes encoding proteins believed to be involved in nitrogen uptake and oxidation, and two hypothetical proteins being the most abundant transcripts in all libraries. Analysis of individual genes over the seasonal cycle suggested that transcription was tied more closely to variation in growth rates than to seasonal changes in environmental conditions. Day-night differences in the relative abundance of transcripts for ribosomal proteins suggested diurnal variation in Thaumarchaeota growth.

An analysis of thaumarchaeota populations from the northern gulf of Mexico.

We sampled Thaumarchaeota populations in the northern Gulf of Mexico, including shelf waters under the Mississippi River outflow plume that are subject to recurrent hypoxia. Data from this study allowed us to: (1) test the hypothesis that Thaumarchaeota would be abundant in this region; (2) assess phylogenetic composition of these populations for comparison with other regions; (3) compare the efficacy of quantitative PCR (qPCR) based on primers for 16S rRNA genes (rrs) with primers for genes in the ammonia oxidation (amoA) and carbon fixation (accA, hcd) pathways; (4) compare distributions obtained by qPCR with the relative abundance of Thaumarchaeota rrs in pyrosequenced libraries; (5) compare Thaumarchaeota distributions with environmental variables to help us elucidate the factors responsible for the distributions; (6) compare the distribution of Thaumarchaeota with Nitrite-Oxidizing Bacteria (NOB) to gain insight into the coupling between ammonia and nitrite oxidation. We found up to 10(8) copies L(-1) of Thaumarchaeota rrs in our samples (up to 40% of prokaryotes) by qPCR, with maximum abundance in slope waters at 200-800 m. Thaumarchaeota rrs were also abundant in pyrosequenced libraries and their relative abundance correlated well with values determined by qPCR (r (2) = 0.82). Thaumarchaeota populations were strongly stratified by depth. Canonical correspondence analysis using a suite of environmental variables explained 92% of the variance in qPCR-estimated gene abundances. Thaumarchaeota rrs abundance was correlated with salinity and depth, while accA abundance correlated with fluorescence and pH. Correlations of Archaeal amoA abundance with environmental variables were primer-dependent, suggesting differential responses of sub-populations to environmental variables. Bacterial amoA was at the limit of qPCR detection in most samples. NOB and Euryarchaeota rrs were found in the pyrosequenced libraries; NOB distribution was correlated with that of Thaumarchaeota (r (2) = 0.49).

Alkane hydroxylase gene (alkB) phylotype composition and diversity in northern Gulf of Mexico bacterioplankton.

Natural and anthropogenic activities introduce alkanes into marine systems where they are degraded by alkane hydroxylases expressed by phylogenetically diverse bacteria. Partial sequences for alkB, one of the structural genes of alkane hydroxylase, have been used to assess the composition of alkane-degrading communities, and to determine their responses to hydrocarbon inputs. We present here the first spatially extensive analysis of alkB in bacterioplankton of the northern Gulf of Mexico (nGoM), a region that experiences numerous hydrocarbon inputs. We have analyzed 401 partial alkB gene sequences amplified from genomic extracts collected during March 2010 from 17 water column samples that included surface waters and bathypelagic depths. Previous analyses of 16S rRNA gene sequences for these and related samples have shown that nGoM bacterial community composition and structure stratify strongly with depth, with distinctly different communities above and below 100 m. Although we hypothesized that alkB gene sequences would exhibit a similar pattern, PCA analyses of operational protein units (OPU) indicated that community composition did not vary consistently with depth or other major physical-chemical variables. We observed 22 distinct OPUs, one of which was ubiquitous and accounted for 57% of all sequences. This OPU clustered with AlkB sequences from known hydrocarbon oxidizers (e.g., Alcanivorax and Marinobacter). Some OPUs could not be associated with known alkane degraders, however, and perhaps represent novel hydrocarbon-oxidizing populations or genes. These results indicate that the capacity for alkane hydrolysis occurs widely in the nGoM, but that alkane degrader diversity varies substantially among sites and responds differently than bulk communities to physical-chemical variables.