Archaeal and Bacterial Diversity and Distribution Patterns in Mediterranean-Climate Vernal Pools of Mexico and the Western USA.

Biogeographic patterns in microorganisms are poorly understood, despite the importance of microbial communities for a range of ecosystem processes. Our knowledge of microbial ecology and biogeography is particularly deficient in rare and threatened ecosystems. We tested for three ecological patterns in microbial community composition within ephemeral wetlands-vernal pools-located across Baja California (Mexico) and California (USA): (1) habitat filtering; (2) a latitudinal diversity gradient; and (3) distance decay in community composition. Paired water and soil samples were collected along a latitudinal transect of vernal pools, and bacterial and archaeal communities were characterized using 16S rDNA sequencing. We identified two main microbial communities, with one community present in the soil matrix that included archaeal and bacterial soil taxa, and another community present in the overlying water that was dominated by common freshwater bacterial taxa. Aquatic microbial communities were more diverse in the north, and displayed a significant but inverted latitudinal diversity pattern. Aquatic communities also exhibited a significant distance-decay pattern, with geographic proximity, and precipitation explaining part of the community variation. Collectively these results indicate greater sensitivity to spatial and environmental variation in vernal pool aquatic microbial communities than in soil microbial communities. We conclude that vernal pool aquatic microbial communities can display distribution patterns similar to those exhibited by larger organisms, but differ in some key aspects, such as the latitudinal gradient in diversity.

Biogeochemistry and hydrography shape microbial community assembly and activity in the eastern tropical North Pacific Ocean oxygen minimum zone.

Oceanic oxygen minimum zones (OMZs) play a pivotal role in biogeochemical cycles due to extensive microbial activity. How OMZ microbial communities assemble and respond to environmental variation is therefore essential to understanding OMZ functioning and ocean biogeochemistry. Sampling along depth profiles at five stations in the eastern tropical North Pacific Ocean (ETNP), we captured systematic variations in dissolved oxygen (DO) and associated variables (nitrite, chlorophyll, and ammonium) with depth and between stations. We quantitatively analysed relationships between oceanographic gradients and microbial community assembly and activity based on paired 16S rDNA and 16S rRNA sequencing. Overall microbial community composition and diversity were strongly related to regional variations in density, DO, and other variables (regression and redundancy analysis r2  = 0.68-0.82), displaying predictable patterns with depth and between stations. Although similar factors influenced the active community, diversity was substantially lower within the OMZ. We also identified multiple active microbiological networks that tracked specific gradients or features - particularly subsurface ammonium and nitrite maxima. Our findings indicate that overall microbial community assembly is consistently shaped by hydrography and biogeochemistry, while active segments of the community form discrete networks inhabiting distinct portions of the water column, and that both are tightly tuned to environmental conditions in the ETNP.

Climatic, physical, and biogeochemical changes drive rapid oxygen loss and recovery in a marine ecosystem.

Dissolved oxygen (DO) concentrations shape the biogeochemistry and ecological structure of aquatic ecosystems; as a result, understanding how and why DO varies in space and time is of fundamental importance. Using high-resolution, in situ DO time-series collected over the course of a year in a novel marine ecosystem (Jellyfish Lake, Palau), we show that DO declined throughout the marine lake and subsequently recovered in the upper water column. These shifts were accompanied by variations in water temperature and were correlated to changes in wind, precipitation, and especially sea surface height that occurred during the 2015-2016 El Niño-Southern Oscillation event. Multiple approaches used to calculate rates of community respiration, net community production, and gross primary production from DO changes showed that DO consumption and production did not accelerate nor collapse; instead, their variance increased during lake deoxygenation and recovery, and then stabilized. Spatial and temporal variations in rates were significantly related to climatic variability and changes in DO, and causality testing indicated that these relationships were both correlative and causative. Our data indicate that climatic, physical, and biogeochemical properties and processes collectively regulated DO, producing linked feedbacks that drove DO decline and recovery.

Microbes and macro-invertebrates show parallel ß-diversity but contrasting α-diversity patterns in a marine natural experiment.

Documenting ecological patterns across spatially, temporally and taxonomically diverse ecological communities is necessary for a general understanding of the processes shaping biodiversity. A major gap in our understanding remains the comparison of diversity patterns across a broad spectrum of evolutionarily and functionally diverse organisms, particularly in the marine realm. Here, we aim to narrow this gap by comparing the diversity patterns of free-living microbes and macro-invertebrates across a natural experiment provided by the marine lakes of Palau: geographically discrete and environmentally heterogeneous bodies of seawater with comparable geological and climatic history, and a similar regional species pool. We find contrasting patterns of α-diversity but remarkably similar patterns of ß-diversity between microbial and macro-invertebrate communities among lakes. Pairwise dissimilarities in community composition among lakes are positively correlated between microbes and macro-invertebrates, and influenced to a similar degree by marked gradients in oxygen concentration and salinity. Our findings indicate that a shared spatio-temporal and environmental context may result in parallel patterns of ß-diversity in microbes and macro-invertebrates, in spite of key trait differences between these organisms. This raises the possibility that parallel processes also influence transitions among regional biota across the tree of life, at least in the marine realm.

Community ecology across bacteria, archaea and microbial eukaryotes in the sediment and seawater of coastal Puerto Nuevo, Baja California.

Microbial communities control numerous biogeochemical processes critical for ecosystem function and health. Most analyses of coastal microbial communities focus on the characterization of bacteria present in either sediment or seawater, with fewer studies characterizing both sediment and seawater together at a given site, and even fewer studies including information about non-bacterial microbial communities. As a result, knowledge about the ecological patterns of microbial biodiversity across domains and habitats in coastal communities is limited-despite the fact that archaea, bacteria, and microbial eukaryotes are present and known to interact in coastal habitats. To better understand microbial biodiversity patterns in coastal ecosystems, we characterized sediment and seawater microbial communities for three sites along the coastline of Puerto Nuevo, Baja California, Mexico using both 16S and 18S rRNA gene amplicon sequencing. We found that sediment hosted approximately 500-fold more operational taxonomic units (OTUs) for bacteria, archaea, and microbial eukaryotes than seawater (p < 0.001). Distinct phyla were found in sediment versus seawater samples. Of the top ten most abundant classes, Cytophagia (bacterial) and Chromadorea (eukaryal) were specific to the sediment environment, whereas Cyanobacteria and Bacteroidia (bacterial) and Chlorophyceae (eukaryal) were specific to the seawater environment. A total of 47 unique genera were observed to comprise the core taxa community across environment types and sites. No archaeal taxa were observed as part of either the abundant or core taxa. No significant differences were observed for sediment community composition across domains or between sites. For seawater, the bacterial and archaeal community composition was statistically different for the Major Outlet site (p < 0.05), the site closest to a residential area, and the eukaryal community composition was statistically different between all sites (p < 0.05). Our findings highlight the distinct patterns and spatial heterogeneity in microbial communities of a coastal region in Baja California, Mexico.

Microbial community networks associated with variations in community respiration rates during upwelling in nearshore Monterey Bay, California.

Respiration of organic material is a central process in the global carbon (C) cycle catalysed by diverse microbial communities. In the coastal ocean, upwelling can drive variation in both community respiration (CR) and the microbial community, but linkages between the two are not well-understood. We measured CR rates and analysed microbial dynamics via 16S rRNA gene sequencing, to assess whether CR correlated with upwelling irrespective of changes in the microbial community, or if the particular microbial community present was a factor in explaining variations in CR. CR varied significantly over time as a function of temperature, dissolved oxygen (DO) and chlorophyll-all of which are altered by upwelling-but also varied with a 'subnetwork' (i.e., a group of microbial taxa that covaried with one another) of the whole community. One subnetwork was associated with higher CR and warmer temperatures, while another was associated with lower CR and DO. Our results suggest that CR in the coastal ocean varies with both environmental variables, and a portion of the microbial community that is not directly correlated with upwelling intensity.

Microbial diversity and community structure along a lake elevation gradient in Yosemite National Park, California, USA.

Microbial communities are key components of lake ecosystems and play central roles in lake biogeochemical cycles. Freshwater lakes, in turn, have a disproportionate influence on global carbon and nitrogen cycling, while also acting as 'sentinels' of environmental change. Determining what factors regulate microbial community dynamics and their relationship to lake biogeochemistry is therefore essential to understanding global change feedbacks. We used Illumina sequencing of >2 million 16S rRNA genes to examine microbial community structure and diversity in relation to spatial, temporal and biogeochemical variation, within and across lakes located along a 871 m elevation gradient in Yosemite National Park, California, USA. We captured a rich microbial community that included many rare operational taxonomic units (OTUs), but was dominated by a few bacterial classes and OTUs frequently detected in other freshwater ecosystems. Neither richness, evenness nor overall diversity was directly related to elevation. However, redundancy analysis showed that changes in microbial community structure were significantly related to elevation. Along with sampling period and dissolved nutrient concentrations, 29% of the variation in community structure could be explained by measured variables - in congruence with studies in other lakes using different techniques. We also found a distance-decay relationship in microbial community structure across lakes, suggesting that both local environmental factors and dispersal play a role in structuring communities.

Soil microbial community structure is unaltered by plant invasion, vegetation clipping, and nitrogen fertilization in experimental semi-arid grasslands.

Global and regional environmental changes often co-occur, creating complex gradients of disturbance on the landscape. Soil microbial communities are an important component of ecosystem response to environmental change, yet little is known about how microbial structure and function respond to multiple disturbances, or whether multiple environmental changes lead to unanticipated interactive effects. Our study used experimental semi-arid grassland plots in a Mediterranean-climate to determine how soil microbial communities in a seasonally variable ecosystem respond to one, two, or three simultaneous environmental changes: exotic plant invasion, plant invasion + vegetation clipping (to simulate common management practices like mowing or livestock grazing), plant invasion + nitrogen (N) fertilization, and plant invasion + clipping + N fertilization. We examined microbial community structure 5-6 years after plot establishment via sequencing of >1 million 16S rRNA genes. Abiotic soil properties (soil moisture, temperature, pH, and inorganic N) and microbial functioning (nitrification and denitrification potentials) were also measured and showed treatment-induced shifts, including altered NO(-) 3 availability, temperature, and nitrification potential. Despite these changes, bacterial and archaeal communities showed little variation in composition and diversity across treatments. Even communities in plots exposed to three interacting environmental changes were similar to those in restored native grassland plots. Historical exposure to large seasonal and inter-annual variations in key soil properties, in addition to prior site cultivation, may select for a functionally plastic or largely dormant microbial community, resulting in a microbial community that is structurally robust to single and multiple environmental changes.

High abundances of potentially active ammonia-oxidizing bacteria and archaea in oligotrophic, high-altitude lakes of the Sierra Nevada, California, USA.

Nitrification plays a central role in the nitrogen cycle by determining the oxidation state of nitrogen and its subsequent bioavailability and cycling. However, relatively little is known about the underlying ecology of the microbial communities that carry out nitrification in freshwater ecosystems--and particularly within high-altitude oligotrophic lakes, where nitrogen is frequently a limiting nutrient. We quantified ammonia-oxidizing archaea (AOA) and bacteria (AOB) in 9 high-altitude lakes (2289-3160 m) in the Sierra Nevada, California, USA, in relation to spatial and biogeochemical data. Based on their ammonia monooxygenase (amoA) genes, AOB and AOA were frequently detected. AOB were present in 88% of samples and were more abundant than AOA in all samples. Both groups showed >100 fold variation in abundance between different lakes, and were also variable through time within individual lakes. Nutrient concentrations (ammonium, nitrite, nitrate, and phosphate) were generally low but also varied across and within lakes, suggestive of active internal nutrient cycling; AOB abundance was significantly correlated with phosphate (r(2) = 0.32, p<0.1), whereas AOA abundance was inversely correlated with lake elevation (r(2) = 0.43, p<0.05). We also measured low rates of ammonia oxidation--indicating that AOB, AOA, or both, may be biogeochemically active in these oligotrophic ecosystems. Our data indicate that dynamic populations of AOB and AOA are found in oligotrophic, high-altitude, freshwater lakes.

Ocean-scale patterns in community respiration rates along continuous transects across the Pacific Ocean.

Community respiration (CR) of organic material to carbon dioxide plays a fundamental role in ecosystems and ocean biogeochemical cycles, as it dictates the amount of production available to higher trophic levels and for export to the deep ocean. Yet how CR varies across large oceanographic gradients is not well-known: CR is measured infrequently and cannot be easily sensed from space. We used continuous oxygen measurements collected by autonomous gliders to quantify surface CR rates across the Pacific Ocean. CR rates were calculated from changes in apparent oxygen utilization and six different estimates of oxygen flux based on wind speed. CR showed substantial spatial variation: rates were lowest in ocean gyres (mean of 6.93 mmol m(-3) d(-1)±8.0 mmol m(-3) d(-1) standard deviation in the North Pacific Subtropical Gyre) and were more rapid and more variable near the equator (8.69 mmol m(-3) d(-1)±7.32 mmol m(-3) d(-1) between 10°N and 10°S) and near shore (e.g., 5.62 mmol m(-3) d(-1)±45.6 mmol m(-3) d(-1) between the coast of California and 124°W, and 17.0 mmol m(-3) d(-1)±13.9 mmol m(-3) d(-1) between 156°E and the Australian coast). We examined how CR varied with coincident measurements of temperature, turbidity, and chlorophyll concentrations (a proxy for phytoplankton biomass), and found that CR was weakly related to different explanatory variables across the Pacific, but more strongly related to particular variables in different biogeographical areas. Our results indicate that CR is not a simple linear function of chlorophyll or temperature, and that at the scale of the Pacific, the coupling between primary production, ocean warming, and CR is complex and variable. We suggest that this stems from substantial spatial variation in CR captured by high-resolution autonomous measurements.

Deoxygenation alters bacterial diversity and community composition in the ocean's largest oxygen minimum zone.

Oceanic oxygen minimum zones (OMZs) have a central role in biogeochemical cycles and are expanding as a consequence of climate change, yet how deoxygenation will affect the microbial communities that control these cycles is unclear. Here we sample across dissolved oxygen gradients in the oceans' largest OMZ and show that bacterial richness displays a unimodal pattern with decreasing dissolved oxygen, reaching maximum values on the edge of the OMZ and decreasing within it. Rare groups on the OMZ margin are abundant at lower dissolved oxygen concentrations, including sulphur-cycling Chromatiales, for which 16S rRNA was amplified from extracted RNA. Microbial species distribution models accurately replicate community patterns based on multivariate environmental data, demonstrate likely changes in distributions and diversity in the eastern tropical North Pacific Ocean, and highlight the sensitivity of key bacterial groups to deoxygenation. Through these mechanisms, OMZ expansion may alter microbial composition, competition, diversity and function, all of which have implications for biogeochemical cycling in OMZs.

Nitrite oxidation in the upper water column and oxygen minimum zone of the eastern tropical North Pacific Ocean.

Nitrogen (N) is an essential nutrient in the sea and its distribution is controlled by microorganisms. Within the N cycle, nitrite (NO2(-)) has a central role because its intermediate redox state allows both oxidation and reduction, and so it may be used by several coupled and/or competing microbial processes. In the upper water column and oxygen minimum zone (OMZ) of the eastern tropical North Pacific Ocean (ETNP), we investigated aerobic NO2(-) oxidation, and its relationship to ammonia (NH3) oxidation, using rate measurements, quantification of NO2(-)-oxidizing bacteria via quantitative PCR (QPCR), and pyrosequencing. (15)NO2(-) oxidation rates typically exhibited two subsurface maxima at six stations sampled: one located below the euphotic zone and beneath NH3 oxidation rate maxima, and another within the OMZ. (15)NO2(-) oxidation rates were highest where dissolved oxygen concentrations were <5 µM, where NO2(-) accumulated, and when nitrate (NO3(-)) reductase genes were expressed; they are likely sustained by NO3(-) reduction at these depths. QPCR and pyrosequencing data were strongly correlated (r(2)=0.79), and indicated that Nitrospina bacteria numbered up to 9.25% of bacterial communities. Different Nitrospina groups were distributed across different depth ranges, suggesting significant ecological diversity within Nitrospina as a whole. Across the data set, (15)NO2(-) oxidation rates were decoupled from (15)NH4(+) oxidation rates, but correlated with Nitrospina (r(2)=0.246, P<0.05) and NO2(-) concentrations (r(2)=0.276, P<0.05). Our findings suggest that Nitrospina have a quantitatively important role in NO2(-) oxidation and N cycling in the ETNP, and provide new insight into their ecology and interactions with other N-cycling processes in this biogeochemically important region of the ocean.

Oceanographic and biological effects of shoaling of the oxygen minimum zone.

Long-term declines in oxygen concentrations are evident throughout much of the ocean interior and are particularly acute in midwater oxygen minimum zones (OMZs). These regions are defined by extremely low oxygen concentrations (<20-45 µmol kg(-1)), cover wide expanses of the ocean, and are associated with productive oceanic and coastal regions. OMZs have expanded over the past 50 years, and this expansion is predicted to continue as the climate warms worldwide. Shoaling of the upper boundaries of the OMZs accompanies OMZ expansion, and decreased oxygen at shallower depths can affect all marine organisms through multiple direct and indirect mechanisms. Effects include altered microbial processes that produce and consume key nutrients and gases, changes in predator-prey dynamics, and shifts in the abundance and accessibility of commercially fished species. Although many species will be negatively affected by these effects, others may expand their range or exploit new niches. OMZ shoaling is thus likely to have major and far-reaching consequences.

Marine bacterial, archaeal and protistan association networks reveal ecological linkages.

Microbes have central roles in ocean food webs and global biogeochemical processes, yet specific ecological relationships among these taxa are largely unknown. This is in part due to the dilute, microscopic nature of the planktonic microbial community, which prevents direct observation of their interactions. Here, we use a holistic (that is, microbial system-wide) approach to investigate time-dependent variations among taxa from all three domains of life in a marine microbial community. We investigated the community composition of bacteria, archaea and protists through cultivation-independent methods, along with total bacterial and viral abundance, and physico-chemical observations. Samples and observations were collected monthly over 3 years at a well-described ocean time-series site of southern California. To find associations among these organisms, we calculated time-dependent rank correlations (that is, local similarity correlations) among relative abundances of bacteria, archaea, protists, total abundance of bacteria and viruses and physico-chemical parameters. We used a network generated from these statistical correlations to visualize and identify time-dependent associations among ecologically important taxa, for example, the SAR11 cluster, stramenopiles, alveolates, cyanobacteria and ammonia-oxidizing archaea. Negative correlations, perhaps suggesting competition or predation, were also common. The analysis revealed a progression of microbial communities through time, and also a group of unknown eukaryotes that were highly correlated with dinoflagellates, indicating possible symbioses or parasitism. Possible 'keystone' species were evident. The network has statistical features similar to previously described ecological networks, and in network parlance has non-random, small world properties (that is, highly interconnected nodes). This approach provides new insights into the natural history of microbes.

Co-occurrence patterns for abundant marine archaeal and bacterial lineages in the deep chlorophyll maximum of coastal California.

Microorganisms remineralize and respire half of marine primary production, yet the niches occupied by specific microbial groups, and how these different groups may interact, are poorly understood. In this study, we identify co-occurrence patterns for marine Archaea and specific bacterial groups in the chlorophyll maximum of the Southern California Bight. Quantitative PCR time series of marine group 1 (MG1) Crenarchaeota 16S rRNA genes varied substantially over time but were well-correlated (r(2)=0.94, P<0.001) with ammonia monooxygenase subunit A (amoA) genes, and were more weakly related to 16S rRNA genes for all Archaea (r(2)=0.39), indicating that other archaeal groups (for example, Euryarchaeota) were numerically important. These data sets were compared with variability in bacterial community composition based on automated ribosomal intergenic spacer analysis (ARISA). We found that archaeal amoA gene copies and a SAR11 (or Pelagibacter) group Ib operational taxonomic unit (OTU) displayed strong co-variation through time (r(2)=0.55, P<0.05), and archaeal amoA and MG1 16S rRNA genes also co-occurred with two SAR86 and two Bacteroidetes OTUs. The relative abundance of these groups increased and decreased in synchrony over the course of the time series, and peaked during periods of seasonal transition. By using a combination of quantitative and relative abundance estimates, our findings show that abundant microbial OTUs-including the marine Crenarchaeota, SAR11, SAR86 and the Bacteroidetes-co-occur non-randomly; they consequently have important implications for our understanding of microbial community ecology in the sea.

Global declines in oceanic nitrification rates as a consequence of ocean acidification.

Ocean acidification produced by dissolution of anthropogenic carbon dioxide (CO(2)) emissions in seawater has profound consequences for marine ecology and biogeochemistry. The oceans have absorbed one-third of CO(2) emissions over the past two centuries, altering ocean chemistry, reducing seawater pH, and affecting marine animals and phytoplankton in multiple ways. Microbially mediated ocean biogeochemical processes will be pivotal in determining how the earth system responds to global environmental change; however, how they may be altered by ocean acidification is largely unknown. We show here that microbial nitrification rates decreased in every instance when pH was experimentally reduced (by 0.05-0.14) at multiple locations in the Atlantic and Pacific Oceans. Nitrification is a central process in the nitrogen cycle that produces both the greenhouse gas nitrous oxide and oxidized forms of nitrogen used by phytoplankton and other microorganisms in the sea; at the Bermuda Atlantic Time Series and Hawaii Ocean Time-series sites, experimental acidification decreased ammonia oxidation rates by 38% and 36%. Ammonia oxidation rates were also strongly and inversely correlated with pH along a gradient produced in the oligotrophic Sargasso Sea (r(2) = 0.87, P < 0.05). Across all experiments, rates declined by 8-38% in low pH treatments, and the greatest absolute decrease occurred where rates were highest off the California coast. Collectively our results suggest that ocean acidification could reduce nitrification rates by 3-44% within the next few decades, affecting oceanic nitrous oxide production, reducing supplies of oxidized nitrogen in the upper layers of the ocean, and fundamentally altering nitrogen cycling in the sea.

Population ecology of nitrifying archaea and bacteria in the Southern California Bight.

Marine Crenarchaeota are among the most abundant microbial groups in the ocean, and although relatively little is currently known about their biogeochemical roles in marine ecosystems, recognition that Crenarchaeota posses ammonia monooxygenase (amoA) genes and may act as ammonia-oxidizing archaea (AOA) offers another means of probing the ecology of these microorganisms. Here we use a time series approach combining quantification of archaeal and bacterial ammonia oxidizers with bacterial community fingerprints and biogeochemistry, to explore the population and community ecology of nitrification. At multiple depths (150, 500 and 890 m) in the Southern California Bight sampled monthly from 2003 to 2006, AOA were enumerated via quantitative PCR of archaeal amoA and marine group 1 Crenarchaeota 16S rRNA genes. Based on amoA genes, AOA were highly variable in time - a consistent feature of marine Crenarchaeota- however, average values were similar at different depths and ranged from 2.20 to 2.76 x 10(4) amoA copies ml(-1). Archaeal amoA genes were correlated with Crenarchaeota 16S rRNA genes (r(2) = 0.79) and the slope of this relationship was 1.02, demonstrating that the majority of marine group 1 Crenarchaeota present over the dates and depths sampled possessed amoA. Two AOA clades were specifically quantified and compared with betaproteobacterial ammonia-oxidizing bacteria (beta-AOB) amoA genes at 150 m; these AOA groups were found to strongly co-vary in time (r(2) = 0.70, P < 0.001) whereas AOA : beta-AOB ratios ranged from 13 to 5630. Increases in the AOA : beta-AOB ratio correlated with the accumulation of nitrite (r(2) = 0.87, P < 0.001), and may be indicative of differences in substrate affinities and activities leading to periodic decoupling between ammonia and nitrite oxidation. These data capture a dynamic nitrogen cycle in which multiple microbial groups appear to be active participants.

Molecular and biogeochemical evidence for ammonia oxidation by marine Crenarchaeota in the Gulf of California.

Nitrification plays an important role in marine biogeochemistry, yet efforts to link this process to the microorganisms that mediate it are surprisingly limited. In particular, ammonia oxidation is the first and rate-limiting step of nitrification, yet ammonia oxidation rates and the abundance of ammonia-oxidizing bacteria (AOB) have rarely been measured in tandem. Ammonia oxidation rates have not been directly quantified in conjunction with ammonia-oxidizing archaea (AOA), although mounting evidence indicates that marine Crenarchaeota are capable of ammonia oxidation, and they are among the most abundant microbial groups in the ocean. Here, we have directly quantified ammonia oxidation rates by 15N labeling, and AOA and AOB abundances by quantitative PCR analysis of ammonia monooxygenase subunit A (amoA) genes, in the Gulf of California. Based on markedly different archaeal amoA sequence types in the upper water column (60 m) and oxygen minimum zone (OMZ; 450 m), novel amoA PCR primers were designed to specifically target and quantify 'shallow' (group A) and 'deep' (group B) clades. These primers recovered extensive variability with depth. Within the OMZ, AOA were most abundant where nitrification may be coupled to denitrification. In the upper water column, group A tracked variations in nitrogen biogeochemistry with depth and between basins, whereas AOB were present in relatively low numbers or undetectable. Overall, 15NH4+ oxidation rates were remarkably well correlated with AOA group A amoA gene copies (r2=0.90, P<0.001), and with 16S rRNA gene copies from marine Crenarchaeota (r2=0.85, P<0.005). These findings represent compelling evidence for an archaeal role in oceanic nitrification.