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.

Introduction

Polar oceans are subject to strong seasonal variation in primary production, driven by the seasonal cycle of irradiance, vertical mixing and ice cover. The near absence of photoautotrophic production at high latitudes during winter raises the question of how the marine food web, and especially the marine microbial food web, functions during this period of greatly reduced supply of phytoplankton-derived, labile organic carbon. Heterotrophy based on the consumption of moderately recalcitrant dissolved organic carbon stocks built up during the highly productive summer season has been considered as one possibility for supporting the Antarctic winter food web (Azam ). Chemoautotrophic production has also been suggested to account for prokaryotic production in Antarctic waters during winter (Manganelli ), with nitrification as a likely metabolism (Grzymski ; Williams ).

Ammonia oxidation (AO) rates have been measured in most oceans and range from 0 to 100 nm per day (reviewed in Yool ; Ward, 2008, 2011); however, there are few measurements of AO rates from polar regions. Arctic data are limited to a small number of potential rate measurements (addition of 1 μm 15NH4Cl tracer) made on samples collected off Barrow, AK (Christman ). Southern Ocean data are limited to a handful of sites (Olson 1981a, b; Bianchi ) and to depths ⩽100 m, although thaumarchaeotes are abundant in both Arctic (Kalanetra , Christman ) and Antarctic (DeLong ; Massana ; Murray , 1999; Church ; Kalanetra ; Alonso-Sáez ) waters and geochemical evidence suggests significant nitrification in the Antarctic winter mixed layer (Smart ).

The waters of the Antarctic continental shelf West of the Antarctic Peninsula (WAP hereinafter; Supplementary Figure 1a) are characterized by the strong seasonal and interannual variation in irradiance, sea ice cover and phytoplankton biomass and production (Ducklow , 2012) that are typical of polar regions. The region has been studied intensively by the Palmer LTER (PAL-LTER; http://pal.lternet.edu/) since 1990. Water column stratification varies strongly by season (Hofmann and Klinck, 1998; Ducklow ; Martinson ). Using terminology from Church , the winter water column is stratified into two distinct water masses: (1) Antarctic surface water (AASW); and (2) circumpolar deep water (CDW). The surface mixed layer (AASW) extends from the surface to depths of 100–150 m, is colder (−1.8 °C to −0.3 °C) and less saline (salinity 33.3 to 34.3) than the CDW below it (Supplementary Figure 1b). The CDW is found at depths >150 m year-round, and has temperatures ranging between 0.2 °C to 2 °C and salinities ranging from 34.4 to 34.7 (Supplementary Figure 1b). The AASW stratifies into two layers during summer as a result of warming and meltwater input, resulting in three water masses in the summer, coastal water column: (1) summer surface water (SSW, temperature 0.5 °C to 1.5 °C, salinity <33.8, depth 0–50 m); (2) WW (the remainder of the winter mixed layer at depths of 50–150 m); and (3) CDW (Supplementary Figures 1b and c).

Planktonic Archaea are most abundant in the upper water column of this region (AASW) during winter (May–November) and their population is dominated by Thaumarchaeota (Murray ; Church ; Grzymski ). Archaeal abundance decreases during the transition from winter to summer. These reductions are most significant at the surface, with very few Archaea (1–2% of the population of prokaryotes) remaining in the SSW, as compared with higher abundances in the WW (5–10%) and CDW (13–17%) (Murray ; Church ).

Different microbial communities inhabit the three summer water masses (Massana ; Murray ; Church ) and a survey of summer populations of Thaumarchaeota reported lower abundances of Archaea amoA and lower ratios of Archaea amoA to Archaea 16S rRNA genes (rrs) in WW than in CDW samples (Kalanetra ). This suggested functional differentiation of Thaumarchaeota populations in these two water masses, accompanying (or possibly driving) phylogenetic differentiation of the populations, which is also observed (Kalanetra ; Alonso-Sáez ). We hypothesized that AO rate might be affected by community composition at the level of Thaumarchaeota ecotypes, as well as by their abundance. Thaumarchaeota diversity was characterized using high-throughput sequencing of rrs and amoA genes, which allowed robust comparisons of assemblage composition by season, as well as by location and depth, across nearly the entire WAP region. We sampled during late winter and the following summer to determine if the composition of Thaumarchaeota populations changed during the evolution of ASW into summer surface water and WW.

We measured AO rates at multiple depths in continental shelf and slope waters west of the Antarctic Peninsula. We used the measured rates to estimate chemoautotrophic carbon fixation and the potential contribution of AO to bacterioplankton production during winter, as well as the growth rates of thaumarchaeote populations. We compared the AO rates to measurements of gene and transcript abundance with the aim of determining whether rates could be predicted from molecular data. We measured biogenic silica concentrations and diatom 18S rRNA gene abundance to test a hypothesis that Thaumarchaeota distributions in the WAP are coupled to ammonia regeneration from sinking diatom blooms. We also measured the abundance of Nitrospina—a prominent group of nitrite-oxidizing Bacteria—to investigate coupling between these two guilds of nitrifying microorganisms. We quantified Thaumarchaeota urease (ureC) genes in our samples, as their presence has been proposed (Yakimov ; Alonso-Sáez ; Connelly ; Pedneault ) to indicate a role for urea in nitrification and urea use might be tied to the vertical differentiation of Thaumarchaeota populations reported previously (Kalanetra ; Alonso-Sáez ).

Materials and methods

Sample collection

Samples were collected from the ARSV Laurence M Gould during late winter (17–24 September 2010; LMG 10-06; LTER 600 line only) and summer (7–30 January 2011; LMG 11-01, entire WAP region). Stations were located on the Antarctic continental shelf and slope west of the Antarctic Peninsula (Supplementary Figure 1a and Supplementary Table 1) within the PAL-LTER sampling domain (http://pal.lternet.edu/). Sampling focused on the AASW and CDW in late winter and WW and CDW in summer (depths sampled were <500 m; Supplementary Table 1). We did not sample summer surface water as previous work (Church ; Kalanetra ) revealed very low Thaumarchaeota populations in this water mass. Water samples were collected from Go-Flo bottles (General Oceanics Inc., Miami, FL, USA) into opaque HDPE plastic bottles as soon as the rosette was secured on deck. Subsequent processing took place in an adjacent darkened laboratory. RNA samples were filtered immediately, placed in Whirl-pak bags (Nasco, Fort Atkinson, WI, USA), RNAlater (Ambion, Austin, TX, USA) was added, and then they were immediately frozen at −80 °C (<15 min after sample collection). DNA samples were filtered through 0.22 μm pore size Sterivex GVWP filters (EMD Millipore, Billerica, MA, USA), 1.8 ml of lysis buffer (0.75 m sucrose, 40 mm EDTA, 50 mm Tris, pH 8.3) was added, and then the samples were frozen at −80 °C within 30 min of sample collection. Nutrient samples were collected from the Sterivex filtrate, frozen at −20 °C and subsequently analyzed using standard methods: NH4+ (Solórzano, 1969); NO2− and NO3− (Strickland and Parsons, 1972; Jones, 1984). Biogenic silica samples were collected on 0.6 μm pore size, 47 mm diameter, polycarbonate filters using a polycarbonate filter holder, placed in Whirl-pak bags that were frozen at −20 °C and subsequently analyzed as described in Brzezinski .

DNA and RNA analyses

DNA was recovered from frozen filters using a lysozyme and proteinase K digestion, followed by purification using a phenol–chloroform extraction following Bano and Hollibaugh (2000). Relative abundances of Archaea, Bacteria and diatom genes were determined using either SYBR Green I dye (BioRad, Hercules, CA, USA) or TaqMan (Applied Biosystems, Carlsbad, CA, USA) chemistries following published protocols (Kalanetra ; Tolar ) and the primers and probes listed in Supplementary Table 2. Quantitative PCR (qPCR) used an iCycler iQTM Real-Time qPCR detection system (BioRad). Reactions were set up in triplicate and analyzed against a range of standards (101–107 copies  per μl) as described in Tolar . qPCR conditions for Archaea amoA, ammonia-oxidizing Betaproteobacteria amoA and Bacteria rrs genes have been described previously (Kalanetra ; Tolar ). Thaumarchaeota ureC genes were quantified under the same conditions as amoA, with an annealing temperature of 53 °C (from Alonso-Sáez ). Nitrospina rrs genes were quantified as in Mincer and diatom 18S rRNA genes were quantified as described in Nguyen and Baldi .

RNA was extracted from filters as described in Poretsky and Gifford with the substitution of 200-μm-diameter zirconium beads (OPS Diagnostics, Lebanon, NJ, USA) in the initial bead-beating step and a final extract volume of 100 μl. DNA was removed from 30 μl of the RNA extract using the TURBO DNase-Free Kit (Ambion) following the manufacturer's instructions with a second enzyme treatment at 2 × concentration. amoA transcripts from Archaea and Bacteria were quantified as above using the iScript One-Step RT-qPCR Kit with SYBR Green (BioRad) and an additional 10 min reverse transcription step at 50 °C at the beginning of each run.

Raw abundance data (copies per μl of DNA or RNA extract) were converted to concentrations (copies  per liter) of genes or transcripts using the volume filtered and the extract volume and assuming 100% extraction efficiency as in Tolar . The percent of total prokaryotes represented by Thaumarchaeota was calculated using rrs gene abundance (Bacteria plus Thaumarchaeota) determined by qPCR and corrected using an average of 1.8 Bacteria rrs genes per genome (Biers ) and 1.0 Thaumarchaeota rrs gene per genome, as described previously (Kalanetra ). Prokaryote abundance determined by qPCR correlated well with total prokaryote counts made with a flow cytometer by Palmer LTER personnel during LMG 11-01 (slope=0.87; 95% confidence interval=0.52–1.22; Supplementary Figure 2).

Archaeal diversity

Samples for high-throughput sequencing were amplified in triplicate using primers 517Fa and 1058R (Archaea rrs; VAMPS project—http://vamps.mbl.edu/resources/primers.php) or CamoA-19f and CamoA-616r (Archaea amoA; Pester ). The first of the five 517F primers used by the VAMPS project to sequence the Archaea v6-v4 region was found to amplify 93.3% of Marine Group I Archaea and only 11% of Euryarchaeota (with no mismatches allowed), and so was selected for use (Supplementary Table 2). PCR reactions were pooled, purified and quantified, and then modified with Titanium (Lib-L) adaptors and sample-specific barcodes. Tagged amplicons were pooled in equimolar concentrations and submitted for sequencing on a Roche 454 GS-FLX instrument (454 Life Sciences, Branford, CT, USA) with Titanium chemistry maintained by the Georgia Genomics Facility (University of Georgia, Athens, GA, USA). Postsequencing processing of reads followed the Schloss laboratory 454 SOP (http://www.mothur.org/wiki/454_SOP; Schloss , last accessed January 2014).

AO rates

AO rates were measured in 96 h incubations at 0 °C in the dark using 15N-labeled ammonium (>99 at% 15NH4Cl; Cambridge Isotope Laboratories, Tewksbury, MA, USA) added at 50 nm (Santoro ; Beman ) to 200 ml samples contained in acid-washed, sample-rinsed, polycarbonate bottles that were wrapped in aluminum foil, enclosed in cardboard boxes and placed in a dark Percival incubator (Perry, IA, USA). The 50 nm additions were 0.9–93% (median 4.5%) of in situ ammonium concentrations, with only 9 of 66 samples receiving >10% enrichment and only 2 samples receiving >50% enrichments (Supplementary Table 1). There was no statistically significant correlation between AO rate and ammonium enrichment factor for the complete data set. Nine of the 10 samples where enrichments were >10% were from summer CDW samples. The mean±s.e. of AO rates for these samples was 16.6±3.8 versus 17.6±3.3 nm per day for the rest of the summer CDW samples; however, the difference between these subsets was not statistically significant (two-tailed t-test, P=0.84). Treatments were run in duplicate for each sample with a single blank consisting of either filtered sample seawater (Sterivex filtrate) or a subsample frozen at −80 °C immediately after tracer addition. Incubations were terminated by removing 50 ml subsamples from each bottle that were immediately frozen at −80 °C.

The 15NO produced in the incubations was measured using the ‘denitrifier method' (Sigman ) as described in Popp , Dore and Beman . We calculated AO rates from δ15N values of the NO pool as described previously (Beman ; Christman ). Means of the duplicates are reported in Supplementary Table 1. Rates calculated for each replicate assay from a sample differed by >50% of the mean in 16 of 66 samples. These are flagged in Supplementary Table 1, but were retained in the analyses presented below.

We ran time-course incubations with samples taken from two depths (WW and CDW) at two stations to verify that oxidation rates did not change significantly, for example, due to substrate depletion or growth of ammonia oxidizers, during the 96 h incubations. These experiments were set up in 250 ml bottles as above, with six sets of two treatments plus one filtered seawater blank (18 bottles) per sample. Two treatment bottles and one blank were sacrificed at each of six time points taken over time courses of 139 and 181 h. These data are presented in Supplementary Table 3. AO rates were determined from the slope of linear regressions of the data from a given sample. Although the difference between replicates was >50% for some pairs of replicates, which were discarded for the calculation of AO rates from these time courses, there was no indication that rates changed during the course of these prolonged incubations.

Statistical analysis

Principal components analysis (PCA) and non-metric multidimensional scaling (NMDS) were performed in R (http://www.r-project.org/) using the prcomp (stats package; PCA), metaMDS (vegan package; NMDS) and lm (stats package; pairwise linear regression) commands as in Hollibaugh . PCA analysis was run separately on LMG 11-01 and LMG 10-06 data sets to include variables not collected during LMG 10-06 (turbidity, dissolved oxygen, biogenic silica) and to reduce bias from oversampling along the LTER 600 line. A correlation matrix was used for PCA as the variance for all environmental variables was set to 1 (using the option ‘scale.=TRUE' in R). The metaMDS command in R runs a PCA after NMDS to rectify the axes such that the first axis explains the most variation, followed by the second axis, and so on. We compare each PCA axis to both rectified NMDS axes (thus to the entire NMDS data set), but only found significant relationships with NMDS axis 1. Temperature and salinity were used to calculate sigma-t to identify water masses, and NO2− and NO3− were summed as NO for LMG 11-01. Both qPCR-estimated gene abundances and the distribution of rrs and amoA sequences among 454-generated operational taxonomic units (OTUs) were used in NMDS analysis (separately). An analysis of molecular variance (Excoffier ), run in R, was used to test the significance of differences in community composition between samples. Student's t-tests and Pearson's product-moment correlations were also performed in R using the t.test and cor.test commands (stats package), respectively. Model II regressions were calculated using the lmodel2 package (Legendre and Legendre, 1998) in R with major axis regressions.

Data sets

DNA sequences obtained by high-throughput sequencing have been deposited with Genbank in the SRA archive under study PRJNA268106, accession number SRP050086. Nutrient, qPCR and AO rate data have been deposited with the Antarctic Master Directory at http://gcmd.nasa.gov/getdif.htm?NSF-ANT08-38996 and a summary is presented in Supplementary Table 1.

Results

Environmental characteristics

Environmental data are summarized in Table 1 (all data are presented in Supplementary Table 1). PCA of the characteristics of samples taken during summer explained 52.3% (axes PC1 plus PC2) of the total variation in the data set and supported division of samples into WW and CDW water masses (Supplementary Figure 3a), consistent with their characterization based solely on temperature and salinity (Supplementary Figure 1b). The samples did not segregate by sampling location (onshore versus offshore, north versus south). A PCA run on the subset of samples used for high-throughput sequencing demonstrated the same separation of WW and CDW samples as that run for all summer samples (data not shown).

Table 1

Summary of environmental variables measured in the water masses sampled for this study

Water mass	Depth range (m)		Temp (°C)	Salinity (PSU)	Chl a (mg m−3)	Dissolved oxygen (mg l−1)	NO2− (μm)	NO3− (μm)	NH4+ (μm)	Biogenic silicate (μm)
UAASW (winter)	10	Max	−1.76	33.97	0.1	ND	0.2	36.9	1.0	ND
		Min	−1.84	33.90	0.1	ND	0.1	25.7	0.4	ND
		Mean	−1.80	33.93	0.1	ND	0.2	31.2	0.7	ND
		s.e.	0.02	0.01	0.00	ND	0.02	1.6	0.10	ND
LAASW (winter)	55–75	Max	−0.82	34.13	0.1	ND	0.2	36.5	1.1	ND
		Min	−1.83	33.91	0.1	ND	0.1	30.1	0.6	ND
		Mean	−1.62	33.97	0.1	ND	0.2	32.4	0.7	ND
		s.e.	0.16	0.03	0.01	ND	0.01	0.9	0.08	ND
CDW (winter)	260–300	Max	2.09	34.68	0.1	ND	0.2	38.3	1.1	ND
		Min	1.32	34.61	0.1	ND	0.0	29.6	0.5	ND
		Mean	1.55	34.64	0.1	ND	0.1	34.8	0.9	ND
		s.e.	0.11	0.01	0.001	ND	0.02	1.2	0.09	ND
WW (summer)	70–125	Max	0.27	34.03	1.5	8.4	0.3	32.5	5.6	1.4
		Min	−1.77	33.83	0.0	6.5	0.1	13.8	0.6	0.0
		Mean	−1.05	33.94	0.5	7.6	0.2	27.1	2.1	0.2
		s.e.	0.10	0.01	0.09	0.07	0.01	0.9	0.26	0.09
CDW (summer)	285–500	Max	2.04	34.72	0.1	5.7	0.3	38.3	4.0	6.1
		Min	0.19	34.32	0.0	3.9	0.0	12.6	0.1	0.0
		Mean	1.40	34.65	0.0	4.2	0.1	30.2	1.0	0.4
		s.e.	0.09	0.02	0.005	0.07	0.01	1.5	0.20	0.26

Abbreviations: CDW, circumpolar deep water; LAASW, lower Antarctic surface water (winter); ND, not determined; UAASW, upper Antarctic surface water (winter); WW, winter water.

Environmental conditions encountered on the Antarctic continental shelf west of the Antarctic Peninsula on 17–24 September 2010 (LMG 10-06) and 7–30 January 2011 (LMG 11-01). See Supplementary Table 1 for the complete data set. Values shown are maximum (Max); minimum (Min); mean and standard error of the mean (s.e.) for each variable.

AO rates ranged from <0.1 (limit of detection) to 220 nm per day. Late winter AO rates in AASW were not significantly different from rates in late winter CDW samples (Table 2, means±s.e.: 62±16 versus 21±5.7 nm per day, respectively; Student's t-test P=0.11). AO rates in the LAASW (the precursor to WW) were significantly higher compared with that in the WW (79±30 and 13±2.8 for LAASW versus WW, P=0.0003). However, rates in WW samples were not significantly different from rates in summer CDW samples (Table 2; 13±2.8 versus 24±6.6 nm per day; P=0.14).

Table 2

Summary of ammonia oxidation rates, Thaumarchaeota abundance and the contribution of ammonia oxidation to bacterioplankton production on the Antarctic continental shelf

	Water mass
	AASW (0–150 m)	WW (50–150 m)	CDW (>150 m)
Mean AO ratesa (nm per day)
Late Winter (17–24 Sept 2010)	62±16		21±5.7
Summer (7–30 Jan 2011)		13±2.8	24±6.6
Cells producedb (106 cells per liter per day)
Winter	3.1±0.82		1.1±0.28
Summer		0.66±0.14	1.2±0.33
Biomass producedc (106 cells per liter per day)
Winter	2.3±0.60		0.79±0.21
Summer		0.50±0.10	0.90±0.25
Mean abundanced (106 cells per liter)
Winter
Archaea amoA	1.2±0.17		0.78±0.21
Thaum rrs	5.6±0.90		2.9±0.94
Summer
Archaea amoA		4.0±1.0	12±2.7
Thaum rrs		2.9±0.78	14±1.7
Specific growth ratee (per day)
Winter
Archaea amoA	1.2±0.21		0.85±0.21
Thaum rrs	0.46±0.11		0.45±0.17
Summer
Archaea amoA		0.37±0.12	0.16±0.04
Thaum rrs		0.39±0.10	0.12±0.03
Specific nitrification ratef, λnitrif (per day)
Winter	0.078±0.015		0.026±0.007
Summer		0.009±0.002	0.057±0.016
Ammonium pool turnover timef (days)
Winter	30±11		95±54
Summer		608±189	109±30

Abbreviations: AASW, Antarctic surface water; AO, ammonia oxidation; CDW, circumpolar deep water; WW, winter water.

Means and s.e.m. are calculated using data from Supplementary Table 1.

Average ammonia oxidation rates from this paper (Supplementary Table 1).

Means from Supplementary Table 1, calculated assuming 49.5 × 1012 cells produced per mole of NH4+ oxidized (average calculated from Figures 3 and 1 in Konneke and Martens-Habbena , respectively).

Means from Supplementary Table 1 calculated from cells produced assuming 9 fg C per cell (Berg ).

Calculated from average gene abundances from this paper (Supplementary Table 1) assuming amoA and rrs are single-copy genes.

Means from Supplementary Table 1, calculated from cells produced and abundance assuming exponential growth kinetics.

Means from Supplementary Table 1.

Thaumarchaeota gene and transcript abundance

Thaumarchaeota were detected at all stations and depths sampled on both cruises (Supplementary Table 1), with abundance ranging from 1.3 × 104 to 6.5 × 107 cells per liter (estimated from amoA) and 3.6 × 102 to 3.0 × 107 cells per liter (from rrs). Betaproteobacteria amoA was <2% of the abundance of Archaea amoA in all samples and was below the limit of detection in many samples (Supplementary Table 1). Thaumarchaeota accounted for up to 16% of prokaryotes in the AASW (mean=7%) during late winter; up to 11% (mean=1.4%) in the WW; up to 26% (mean=19%) in CDW samples during late winter; and up to 29% in CDW samples during summer (mean=14%).

The differences in the abundance of Thaumarchaeota amoA or rrs genes in AASW versus CDW during late winter were not statistically significant (Figure 1 and Table 2; Student's t-test; P=0.18 for amoA, P=0.085 for rrs). In contrast, Thaumarchaeota amoA and rrs genes were less abundant in WW (means=4.0–2.9 × 106 copies  per liter, respectively) than in CDW (means=12–14 × 106 copies  per liter, respectively) during summer (P=0.0051 for amoA and P<0.0001 for rrs, Student's t-test; Figure 1 and Table 2), and gene abundance in WW samples varied over a wider range (Figure 1). The ratio of archaeal amoA to rrs genes was 0.62 (r2=0.31, 95% confidence interval=0.42–0.87) for all water masses sampled (Figure 1 and Supplementary Table 1), slightly lower compared with the ratio of 1 expected from copy numbers of these genes in sequenced Thaumarchaeota genomes. This discrepancy was taken to indicate the potential for thaumarchaeote heterotrophy in a previous study (Kalanetra ); however, it may also be due to consistent biases in one or both of the primer sets used in qPCR (Tolar ). The abundance of Thaumarchaeota ureC genes was similar to the abundance of Thaumarchaeota amoA and rrs genes (Supplementary Table 1 and Figure 2). The ratio of the abundance of ureC versus rrs genes was 0.28 (r2=0.69, 95% confidence interval=0.24–0.33), and we found no significant differences in the ratio by water mass or season.

Figure 1

Abundance of Thaumarchaeota amoA versus rrs genes in samples from Antarctic continental shelf waters. Symbols represent water masses and dates sampled as shown in the legend. Late winter samples (filled symbols) were collected in September 2010. Summer samples (open symbols) were collected in January 2011. Upper Antarctic surface water (UAASW, filled circles); lower AASW (LAASW, filled squares); and CDW (filled triangles); WW (squares); and CDW (triangles). A line of slope=1 is shown for reference.

Figure 2

Abundance of Thaumarchaeota ureC genes in samples from Antarctic continental shelf waters compared with the abundance of (a) archaeal amoA and (b) Thaumarchaeota rrs genes in the same sample. Symbols represent water masses and dates sampled as in Figure 1. Lines of slope=1 are shown for reference.

NMDS based on qPCR-estimated gene abundances measured in samples from the summer cruise grouped samples by water mass with CDW samples clustering tightly, driven by abundances of Thaumarchaeota (rrs, amoA and ureC) and Nitrospina (rrs) genes, whereas grouping of WW samples was driven by increased Bacteria rrs and diatom 18S rRNA gene abundances (Supplementary Figure 3b). Similar water mass-associated groupings were observed in PCA and NMDS analyses of data from the late winter cruise (Supplementary Figures 4a and b), except that increased AO rates corresponded to the AASW.

Thaumarchaeota amoA transcripts were more abundant in AASW than in CDW samples during late winter (means±s.e.: 2.4±0.69 versus 0.13±0.048 × 106 copies  per liter, respectively; P=0.035), but there was no significant difference in the abundance of amoA transcripts between the lower AASW and WW (3.1±1.2 versus 4.7±1.4 × 106 copies  per liter; P=0.63; Supplementary Table 1 and Figure 3). The abundance of amoA transcripts in CDW samples was not significantly different between late winter versus summer (0.13±0.048 versus 2.7±0.72 × 106 copies  per liter; P=0.096). The ratio of amoA transcripts per amoA gene was 2.8±1.0 for AASW versus 0.16±0.07 for late winter CDW (P=0.089), and 4.4±2.3 for WW versus 0.42±0.13 for summer CDW samples (P=0.11). Betaproteobacteria amoA transcript abundance was 3% (average for all samples) of Archaeal amoA transcript abundance, and was below the limit of detection in many samples (Supplementary Table 1).

Figure 3

Abundance of Archaea amoA genes versus Archaea amoA transcripts in samples from Antarctic coastal waters. Symbols represent water masses and dates sampled as in Figure 1. A line of slope=1 is shown for reference.

The AO rates we measured only correlated significantly with indicators of Thaumarchaeota abundance (rrs and amoA genes), and only for samples collected during summer (all summer samples and the subset of WW samples; Figure 4 and Table 3). These correlations only explained a small portion of the variance in the AO rates (R-values ranged from 0.36 to 0.57). We found no statistically significant correlations between AO rates and either amoA transcript abundance or the ratio of amoA transcripts to genes (Figure 4 and Table 3). We found no statistically significant correlations between AO rates and Betaproteobacteria amoA gene abundances. Nitrospina rrs abundance was significantly correlated with AO rates in both WW and CDW samples during summer (R=0.52, P=0.0076; R=0.42, P=0.05; respectively).

Figure 4

Relationship between AO rate and gene and transcript abundance. Panels show AO rates plotted against: (a) archaeal rrs gene abundance; (b) archaeal amoA gene abundance; (c) archaeal amoA transcript abundance; and (d) the ratio of archaeal amoA transcripts to genes (mRNA/DNA). Symbols represent water masses and dates sampled as in Figure 1.

Table 3

Correlations (r) between ammonia oxidation rates and Thaumarchaeota gene or transcript abundance, measured in Antarctic coastal waters

Sample subset	n	rrs genes	amoA genes	amoA transcripts	Ratio of amoA mRNA:DNA
All winter	18	+0.25	−0.009	+0.41	+0.40
UAASW	6	+0.24	−0.30	−0.38	+0.30
LAASW	6	+0.099	−0.14	+0.082	+0.17
CDW	6	−0.15	+0.21	+0.56	+0.44
All summer	47	+0.45***	+0.36*	+0.025	−0.16
WW	25	+0.52**	+0.57**	+0.33	−0.13
CDW	22	+0.33	+0.23	−0.11	−0.21

Abbreviations: CDW, circumpolar deep water; LAASW, lower Antarctic surface water; UAASW, upper Antarctic surface water; WW, winter water.

Statistically significant correlations are indicated as follows: ***P<0.001, **P<0.01 and *P<0.05.

Thaumarchaeota community composition

A total of 399 389 (rrs; mean length=353 bp) and 154 037 (amoA; mean length=358 bp) useful sequences were retained after processing (Supplementary Table 1). OTUs were defined as 97% and 98% sequence identity for amoA and rrs sequences, respectively. Thaumarchaeota populations, especially those in AASW and WW samples, were strongly dominated by a few clades (Figure 5 and Supplementary Figure 5). The CDW populations of both genes were richer than populations from AASW and WW (Supplementary Figures 6a and b). AASW and WW samples contained the same OTUs (Supplementary Figures 6c and d), with no significant difference in the composition of the populations of rrs or amoA sequences they contain (analysis of molecular variance; P=0.67 and 0.59, respectively). Our analysis revealed that 85% of amoA sequences retrieved from AASW and WW samples fell into one OTU (Figure 5). Most Thaumarchaeota rrs (99.8%) and amoA (99.1%) sequences were common to both late winter and summer CDW populations (Supplementary Figures 6e and f, the distribution of sequences among OTUs is not uniform). However, Thaumarchaeota populations differed significantly between WW and summer CDW samples (as assessed by both rrs and amoA sequences; analysis of molecular variance; P⩽0.001). In summary, WAP Thaumarchaeota populations were strongly differentiated by water mass (AASW or WW versus CDW), but composition within a water mass did not change through time (AASW versus WW, or late winter CDW versus summer CDW), and AASW and WW populations were strongly dominated by one OTU.

Figure 5

Phylogenetic analysis of Archaea amoA sequences retrieved from the study area. Partial sequences (359 bp) of Thaumarchaeota amoA genes obtained by high-throughput sequencing were aligned against the Pester database and the neighbor-joining tree was constructed in ARB. OTUs were defined at 97% similarity. Numbers following each OTU give the number of sequences and % of total sequences it represents. Additional notations apply to four major clades to indicate the % of sequences in these clades from each water mass sampled (AASW—green; WW—blue; CDW—red). Shading of trapezoids representing sequences assigned to four major clades indicates the relative contribution of sequences from AASW+WW (blue) versus all CDW (red) to these clades. Dashed boxes delineate archaeal amoA groups A (surface water) and B (deep water) as defined in Francis . Only bootstrap values ⩾50% (of 1000 iterations) are shown.

NMDS analysis of Thaumarchaeota community composition determined using high-throughput sequencing data also separated WW versus CDW samples (Supplementary Figure 7). Regressions of: AO rates (r2=0.14, P=0.046); AOA amoA transcript abundance (r2=0.34, P<0.001); the ratio of AOA amoA transcripts to genes (r2=0.15, P=0.039); or Nitrospina rrs abundance (r2=0.18, P=0.024) against values for rrs OTU positions on the primary NMDS axis (NMDS1) were significant. There were no statistically significant correlations between any of these variables and amoA OTU positions on the primary NMDS axis (NMDS1).

We found that Nitrospina rrs abundance covaried with Thaumarchaeota gene abundance, in agreement with distributions reported for the Pacific Ocean (Mincer ). We also observed correlations between Nitrospina rrs abundance and AO rates, suggesting a tight coupling between nitrite production by AO and its oxidation by nitrite-oxidizing Bacteria in this region.

Discussion

Distribution of polar Thaumarchaeota

Archaeal amoA and rrs sequences retrieved from the PAL-LTER region in 2010–2011 clustered with sequences obtained from the same area in 2005–2006 by Kalanetra and the phylotypes had similar distributions by water mass in both studies (Supplementary Table 4). Almost all of the WW and AASW amoA sequences (94% and 98%, respectively) from both studies fell into one cluster containing representative ‘near-surface clade' clade A (Francis ) sequences (subcluster 13; Figure 5). This clade was also dominant in libraries constructed from central Arctic Ocean samples (Kalanetra ) and a closely related clade dominated surface samples from the Amundsen Sea (Alonso-Sáez ). Thaumarchaeota diversity was greater in the CDW (Figure 5 and Supplementary Figures 5 and 6), consistent with previous studies of Antarctic coastal waters (Murray ; Kalanetra ; Alonso-Sáez ). The strong grouping of Thaumarchaeota phylotypes by depth, and significantly lower diversity in near-surface versus deep water samples, has been observed in other marine locations that are quite distinct from the Antarctic; for genes other than amoA: rrs, accA and ureC (Yakimov , 2011; Hu ; Tolar ); and the distinction between near-surface and deep water clades appears to be a genomic characteristic (Luo ).

Seasonal shifts in the relative abundance of Thaumarchaeota in surface waters of polar regions have been noted previously (Murray ; Church ; Alonso-Sáez ; Christman ; Luria ), with higher abundances in winter compared with summer. AO rates in AASW versus WW decreased from late winter to summer, with no corresponding change in AO rates in the deeper CDW. This was surprising because we expected CDW populations to respond to enhanced export of organic matter accompanying melting pack ice and the spring phytoplankton bloom (Karl ; Ducklow ).

Ammonia oxidation

Our data indicate that AOB are much less abundant than AOA (Betaproteobacterial AOB amoA genes and transcripts were <3% of Archaeal amoA genes and transcripts in our samples) and we found no correlation between AO rates and the abundance of Betaproteobacterial AOB amoA genes or transcripts, whereas AO rates were significantly correlated with AOA amoA gene abundance in the same samples. AO rates per AOA cell (Supplementary Table 1) are 14±3.4 for the late winter cruise and 5.6±1.2 fmol per ammonia-oxidizer cell per day for the summer cruise. These are in the range of ammonia oxidation rates for Thaumarchaeota isolates reported in Qin : 2.9–12.7 fmol per ammonia-oxidizer cell per day. Although we cannot definitively exclude the possibility of a significant contribution by Betaproteobacterial AOB to the AO rates we measured, our analysis does not suggest excess AO that needs to be accounted for by something other than the Thaumarchaeota we detected in our samples.

There are no other comparable AO rate measurements from Antarctic coastal waters, and only a handful of observations from the Southern Ocean. Bianchi reported AO rates of 25.6 nm per day (averaged over depths 0–100 m) at their southernmost station in the Indian Ocean (52°S, 62°E; north of the Polar Front) during the fall. Olson (1981b) measured AO rates that averaged 2.8 nm per day (n=23, range 0–13 nm per day) in open ocean waters south of the Polar Front in the Scotia Sea (25°W to 60°W and 55°S to 61°S), from depths <100 m, during late winter. These rates are more than an order of magnitude lower compared with the rates we measured in late winter (62±16; Table 2). Olson (1981a) also measured AO rates in two samples retrieved from 30 and 300 m in the Ross Sea during summer. These rates (6.0 and 8.9 nm per day, respectively) are also lower compared with our summer measurements (mean=18 nm per day). It seems unlikely that the differences in analytical methods used by Olson (Wada and Hattori, 1972; Schell, 1978) versus our study would lead to a systematic bias that could account for the differences in rates.

AO rates in our samples (overall mean=27 nm per day) were also higher than potential rates (added 15NH4Cl=1 μm) measured in shallow water (2 m), nearshore Arctic communities, which ranged from 3.6 nm per day (winter mean of four samples) to 0.14 nm per day (summer mean of four samples; Christman ). The abundance of Archaea in the Arctic samples (Christman ) was much lower compared with that we found in Antarctic coastal waters (means of 0.008% and 1.2% of prokaryotes for summer and winter samples, respectively), consistent with a previous comparison (Kalanetra ).

Correlation of rates with gene and transcript abundance

It is of interest to know whether AO rates can be predicted from measurements of gene or transcript abundance (or ratios), because it is typically easier and faster to measure gene or transcript abundance than AO rates. We only found significant (P<0.05) correlations between AO rates and Thaumarchaeota gene abundance for selected subsets of our data: WW, or all summer samples (Table 3). We did not find a relationship between transcript abundance and rates. Although we processed mRNA samples as rapidly as possible after the water sample was available to us, the lack of a correlation might be a result of rapid changes in the mRNA content of Thaumarchaeota (Feike ) as a consequence of changes that take place during the time required for processing samples (Moran ). Alternatively, transcription and AO activity may not be tightly coupled in these populations. Activity requires translation of the mRNA into a protein, the nascent protein may require further processing, including the insertion of cofactors, and AO may be tied more closely to the activity of proteins other than amoA. Thus, lack of a significant correlation between amoA transcript abundance and AO is not surprising. Santoro found no correlation between AO rates and amoA abundance for samples from the central California Current. They postulated that this might be due to differences in cell-specific AO rates among Thaumarchaeota ecotypes. This was shown in studies of Monterey Bay by Smith , who found significant correlations between AO rates and the abundance of amoA genes from the ‘surface water' clade (WCA) of Thaumarchaeota, but not with the ‘deep water' clade (WCB) or with total amoA gene abundance. This depth distribution of ecotypes is similar to the difference we noted between WW and CDW samples on the summer cruise. Although there seems to be a weak relationship between AO and gene abundance in some sets of samples (this study; Beman ), the inconsistency of this relationship cautions against trying to predict AO rates from amoA gene or transcript abundance or ratios.

Geochemical consequences

The rates we measured in the Antarctic are comparable to rates measured in an upwelling zone off Monterey Bay: 0–210 nm per day, with a mean of ~100 nm per day (Santoro ), or 0–70 nm per day, with a mean of ~20 nm per day (Smith ); and in an oceanic oxygen minimum zone: 0–348 nm per day, with a mean of ~50 nm per day (Beman , 2013). Taking the area of oceanic oxygen minimum zones (8% of the global ocean; Paulmier and Ruiz-Pino, 2009) as a proxy for the area of elevated nitrification rates and assuming no nitrification in the rest of the ocean, the contribution of nitrification on the Antarctic continental shelf to the global oceanic total would be ~9% (Supplementary Table 5). Although rough, these calculations serve to demonstrate the potential importance of nitrification in Antarctic coastal waters to the global nitrogen cycle.

Yool , who were interested in the impact of AO on new production measurements (Dugdale and Goering, 1967; Eppley and Peterson, 1979), reported that specific nitrification rates (λnitrif, calculated by dividing the AO rate by the concentration of ammonia in a sample), ranged from 0.001 to 10 per day (mean=0.55 per day) across the global ocean, with Olson's Antarctic samples (Olson 1981a, b) an order of magnitude lower, on average (mean=0.017 per day). λnitrif values calculated from our data (mean for all data=0.041 per day; range=0.0003–0.31 per day; Supplementary Table 1) were on average 2.4 times greater than those reported by Yool for Olson's data.

The inverse of λnitrif is the turnover time of the ammonium pool as a result of its oxidation to nitrite (this does not include ammonium pool turnover driven by phytoplankton or bacterioplankton assimilation). This ranged from 5.2 to 3400 days with a mean of 269 days for our samples (Supplementary Table 1). There were no statistically significant differences in ammonium pool turnover times between AASW, spring CDW and summer CDW samples (Table 2); however, turnover times were significantly longer in WW samples than in AASW (608 versus 30 days, P=0.032) or summer CDW samples (608 versus 109 days, P=0.013).

Nitrate-based phytoplankton growth supplies labile organic nitrogen that can be regenerated rapidly via zooplankton grazing and the microbial loop (Azam ) to replenishing ammonium consumed by AO. However, phytoplankton growth, and thus nitrate assimilation, is greatly reduced during the polar winter, decreasing the flux of ammonium from this source. The mean turnover time of the AASW ammonia pool as a result of AO (30 days; Table 1) is shorter than the duration of the Antarctic winter, indicating the need for an alternate source of ammonium, presumably excretion by eukaryotes or metazoa or the slow degradation of moderately labile organic nitrogen (Manganelli ), to support winter chemoautotrophy based on AO. The AO rates we measured were higher in winter compared with summer, despite the inferred decrease in ammonia regeneration flux in winter versus summer. This suggests that inhibition (e.g., Murray ) may be more important than ammonium supply in controlling thaumarchaeote production during summer.

Contribution of AO to chemoautotrophy and bacterioplankton productivity

Chemoautotrophic carbon fixation by Thaumarchaeota cells was estimated from AO rates (Supplementary Table 1) using 10 mol of ammonium oxidized to nitrite per mole of C fixed (Berg ). We also estimated the production of new thaumarchaeote cells and Thaumarchaeota biomass from the relationships between ammonium consumed and Nitrosopumilus cells produced from Figure 3 in Konneke and Figure 1 in Martens-Habbena : 46.6 and 52.4 (mean=49.5) cells per nmol ammonium oxidized, respectively. Our estimates ranged from 0.05 to 22 nm of C fixed per liter per day and 0.026–10.9 × 106 cells per liter per day. When combined with estimates of Thaumarchaeota abundance (from amoA or rrs abundance), population growth rates averaged 0.50 per day (all amoA data, range 0.01–2.5 per day) or 0.31 per day (all rrs data, range 0.01 to 2.1 per day). Growth rates differed among water masses (analysis of variance P=0.001 for rates calculated using amoA, P=0.013 for rates calculated using rrs), with higher growth rates during late winter than summer (all depths: amoA means 1.1 versus 0.26 per day, P<0.0001; rrs means 0.45 versus 0.25 per day, P=0.092). Calculated growth rates were higher in AASW and WW versus CDW (Table 2), although only the summer comparison based on rrs gene abundance was statistically significant: 0.39±0.10 versus 0.12±0.03, P=0.016.

Taking the rates we measured during late winter to represent those of fall and winter, and integrating over a mean continental shelf water column depth of 475 m, the mean area-specific rates we measured are equivalent to 10.1 and 14.7 mmol of ammonium-N oxidized per m2 per day in summer versus fall+winter+spring, respectively (Supplementary Table 6a). Chemoautotrophic carbon fixation supported by ammonium oxidation can be calculated to be 0.5 mol C  per m2 per year; or about 3% of the annual average photoautotrophic primary production of 15 mol C per m2 per year estimated for the Antarctic continental shelf by Arrigo (Supplementary Table 6b), or 5.7% of an average of 8.1 mol C per m2 per year (range 0.4–16.3% Supplementary Table 6c) estimated for the PAL-LTER study area (from Ducklow ).

Chemoautotrophy has been proposed as a means of maintaining bacterioplankton productivity, and thus the microbial food web, during the Antarctic winter (Manganelli ; Grzymski ). We estimate that production of biomass supported by AO is 1.5 g C  per m2 over the 8 months of fall, winter and spring (Supplementary Table 6a). This production is consistent with estimates based on incorporation of l-leucine and dissolved inorganic carbon reported by Manganelli ). They concluded that chemoautotrophic production could account for 20% of winter bacterioplankton production, or about 3.5 g C  per m2 when integrated over the same depth range we used (see Supplementary Table 6d) in waters of the southern Drake Passage. Our calculations suggest that production of ammonia oxidizer biomass could account for about 8.8% of winter bacterioplankton production at this site.

Annual integrated bacterioplankton production in the Ross Sea is about 23 g C  per m2 per year (Supplementary Table 6e, calculated from Ducklow ), but this includes elevated rates in summer coupled to phytoplankton production. Our estimates suggest that production of ammonia oxidizer biomass could account for about 9% of this annually integrated total. The winter rates of bacterioplankton production they measured (0.8 mmol C per m2 per day) are comparable to the rates we estimate from AO (0.5 mmol C per m2 per day), as are winter rates of bacterioplankton production referenced in Grzymski . Assuming that Thaumarchaeota take up and incorporate exogenous l-leucine, a significant portion of the winter bacterioplankton production measured on the WAP shelf may be growth of the Thaumarchaeota population, as proposed by Grzymski . Alternatively, if Thaumarchaeota do not incorporate l-leucine to the same extent as Bacteria (i.e., data presented in Alonso-Sáez ), then production supported by AO may be in addition to heterotrophic production measured by l-leucine incorporation. In either case, our data indicate that Thaumarchaeota can account for a significant portion of prokaryotic production during the Antarctic winter.

Thaumarchaeota production is unlikely to be able to meet all of the energy demand, because heterotrophic decomposition of organic matter is required to supply the ammonium to be oxidized. Assuming Redfield stoichiometry for labile organic matter (Redfield ) and that ammonium release is linked to complete oxidation of that organic matter, 6.6 mol of C must be mineralized per mole of ammonium regenerated. For each mole of Thaumarchaeota biomass C produced, ~27 mol of ammonia are oxidized by Thaumarchaeota; thus, 178 mol of labile organic C must be oxidized to support chemoautotrophic production of 1 mol of Thaumarchaeota C, unless AOA are mixotrophic or heterotrophic (Spang ; Qin ; Lin ) or they are able to access other forms of reduced nitrogen (e.g., organic nitrogen in urea or other small molecules) more directly. Note that chemoautotrophy fueled by AO is not ‘new production', sensu Eppley and Peterson (1979), even though carbon is being fixed: since AO depends on a source of reduced N supplied by organic matter mineralization, at an ecosystem level it is part of a net oxidative process (Smith and Hollibaugh, 1989).

Conclusions

AO rates measured in Antarctic continental shelf waters ranged from <0.1 to 220 nm per day and are an order of magnitude greater than rates previously reported for other Southern Ocean waters (Scotia Sea, Ross Sea). AOA populations in the surface mixed layer during winter oxidize ammonia at higher rates than those in the summer ‘WW' layer that is derived from the winter mixed layer. AO rates in CDW samples have not been measured previously. They are comparable to rates measured by others in temperate waters, suggesting that the contribution of the Southern Ocean to global nitrification has been underestimated. Chemoautotrophic production supported by AO can account for a ~9% of winter bacterioplankton production in the study area, but must be supported by ongoing ammonium regeneration. amoA transcript abundance was 0.02–40% of amoA gene abundance across all water masses. Despite this large dynamic range, we only found significant correlations between AO rates and gene abundance, and only for a selected subset of our samples. This suggests caution in using amoA transcript abundance or transcript to gene ratios to predict AO rates in the environment.