Abstract DNA stable isotope probing (SIP) was used to track the uptake of organic and inorganic carbon sources for TACK archaea (Thaumarchaeota/Aigarchaeota/Crenarchaeota/Korarchaeota) on a cruise of opportunity in the North Atlantic. Due to water limitations, duplicate samples from the deep photic (60–115 m), the mesopelagic zones (local oxygen minimum; 215–835 m) and the bathypelagic zone (2085–2835 m) were amended with various combinations of 12C- or 13C-acetate/urea/bicarbonate to assess cellular carbon acquisition. The SIP results indicated the majority of TACK archaeal operational taxonomic units (OTUs) incorporated 13C from acetate and/or urea into newly synthesized DNA within 48 h. A small fraction (16%) of the OTUs, often representing the most dominant members of the archaeal community, were able to incorporate bicarbonate in addition to organic substrates. Only two TACK archaeal OTUs were found to incorporate bicarbonate but not urea or acetate. These results further demonstrate the utility of SIP to elucidate the metabolic capability of mesothermal archaea in distinct oceanic settings and suggest that TACK archaea play a role in organic carbon recycling in the mid-depth to deep ocean. TACK archaea, SIP, heterotrophy, mixotrophy, mesopelagic, North Atlantic INTRODUCTION Archaea are ubiquitous and abundant in the deep ocean (Delong 1992; Fuhrman, McCallum and Davis 1992; DeLong et al.1994; Stein and Simon 1996; Karner, DeLong and Karl 2001). However, their modes of nutrition remain largely unknown. Much of our current understanding of marine archaeal phyla is derived from DNA sequence data and the few marine members of the domain that are currently in culture. It is therefore necessary to test TACK archaeal (Thaumarchaeota/Aigarchaeota/Crenarchaeota/Korarchaeota; Guy and Ettema 2011) metabolism using field samples with very little knowledge of their metabolic potential. The abundance of archaea throughout marine ecosystems would imply that their activity greatly influences the modes and efficiency of carbon transfer in microbial systems. Additionally, the structure of the marine microbial food web plays a major role in oceanic carbon sequestration and the function of the marine ecosystem as a whole (Limardo and Worden 2015). The archaeal contribution to the microbial loop is, however, still poorly understood. The earliest report of the metabolism of marine thaumarchaea (formerly classified as crenarchaea) from the field was by Ouverney and Fuhrman (2000). Using a microradiographic/fluorescent in situ hybridization method (STARFISH, aka MAR-FISH), the researchers demonstrated that archaea represent 14%–43% of the total microbial community at depth in the Mediterranean and Pacific, and ∼50% of the archaea at 200 m accumulated 3H-amino acids. Later studies demonstrated the uptake of amino acids by archaea using MAR-FISH in the North Atlantic (Teira et al.2004, 2006) and Antarctic Circumpolar Current (Herndl et al.2005); in Arctic waters, the addition of radiolabeled amino acids, whole proteins, glucose and diatom extracellular polymeric substances all produced labeled thaumarchaeal cells (Kirchman et al.2007). Unfortunately, MAR-FISH approaches have limited phylogenetic resolution and the method does not unambiguously demonstrate that a radiolabeled substrate enters central metabolism. Furthermore, it is unclear whether the positive MAR-FISH signal represents a few operational taxonomic units (OTUs) dominating the community with a specific metabolic capability, or if there are many different community members which share a broader metabolic repertoire. In contrast to MAR-FISH, metagenomic analyses have shown that many thaumarchaea possess an ammonium monooxygenase (amoA) gene (Venter et al.2004; Treusch et al.2005), indicating nitrification plays an important role in archaeal metabolism in the water column (Francis et al.2005; Francis, Beman and Kuypers 2007). Initially, an obligately autotrophic culture of Nitrosopumilus maritimus, an ammonia-oxidizing archaeum, was isolated from an aquarium (Könneke et al.2005). Subsequently, other nitrifying archaeal isolates/co-cultures have been obtained from marine sediments, decaying leaf litter and deep agricultural soils (Park et al.2010; Jung et al.2014, 2016). These isolates and other enriched cultures of ammonia-oxidizing archaea (e.g. Nitrosoarchaeum) have demonstrated incorporations of 13CO2 into lipids (Park et al.2010; Jung et al.2014, 2016; Berg et al.2015). Several other studies support the chemoautotrophic fixation of carbon by marine archaea in situ (e.g. Wuchter et al.2003; Ingalls et al.2006). Alternatively, genomic results suggest that marine thaumarchaea may be capable of mixotrophy. Specifically, fosmid sequences from the sponge symbiont Cenarchaeum symbiosum revealed the presence of genes associated with autotrophic carbon assimilation, as well as an oxidative tricarboxylic acid cycle (Hallam et al.2006). Thaumarchaea have also been found to contain genes for 3-hydroxypropionate carbon fixation as well as oligopeptide importers, supporting the concept of mixotrophy (Martin-Cuadrado et al.2008). Two cultures of Nitrosopumilus (strains HCA1 and PS0) were found to require α-ketoglutaric acid to grow in culture (Qin et al.2014; Kim et al.2016). However, the role of α-ketoglutarate as a hydrogen peroxide scavenger rather than a metabolic cofactor has been demonstrated for an ammonium-oxidizing archaeal seawater isolate (DDS1; Kim et al. 2016). These metagenomic studies suggest metabolic potential via gene annotation. However, gene annotations can often be tenuous (low similarity) and it is difficult linking genomic information from TACK archaea in marine samples with specific phylogenetic markers, such as 16S rRNA genes. Furthermore, the presence of a gene does not confirm activity, only potential. Therefore, an empirical testing of potential food sources is essential to definitively prove that a TACK archaeum is capable of metabolizing a particular substrate (Seyler, McGuinness and Kerkhof 2014). Because most TACK archaea are not in culture, we opted for an SIP approach to determine if the TACK archaeal community is composed of varying levels of heterotrophs, autotrophs or mixotrophs. DNA stable isotope probing (SIP) was used on replicate samples from the North Atlantic Basin to determine the uptake of organic and inorganic 13C-carbon for specific OTUs (as defined by terminal restriction fragment length polymorphism (TRFLP) peaks of 16S rRNA genes; Avaniss-Aghajani et al.1994). SIP methods can elucidate active microorganisms capable of incorporating 13C (or 15N) from specific compounds (Radajewski et al.2000; Buckley et al.2007). In SIP, a 13C-labeled substrate is used to generate 13C-labeled DNA from replicating microbes prior to separation from 12C-DNA in cesium gradients, and represents an intermediate approach between metagenomics and MAR-FISH. This method is constrained by the use of specific 13C-labeled substrates to generate 13C-labeled DNA for separation. Although the approach does not provide a broad representation of potential carbon sources, as does metagenomic analysis of metabolic and carbon assimilation pathways, SIP analysis can demonstrate unambiguous 13C incorporation into newly synthesized DNA. 13C-substrate concentrations are not similar to in situ conditions as with MAR-FISH because SIP requires an isotope atom % enrichment above ambient (generally >25%) to physically separate the newly synthesized 13C-DNA obtained during the incubations. SIP has much finer phylogenetic resolution than MAR-FISH, but does not measure in situ incorporation rates. This study focused on carbon uptake in marine thaumarchaea and crenarchaea, referred throughout the text as the TACK archaeal superphylum. MATERIALS AND METHODS In April 2012, samples were taken onboard the University-National Oceanographic Laboratory System (UNOLS) ship R/V Atlantis, as part of the US CLIVAR/CO2 Repeat Hydrography program (Wanninkhof et al. 2013). This study was a cruise of opportunity, replicating the World Ocean Circulation Experiment line A20 (longitude 52°W), previously conducted in 2003 (Fig. S1, Supporting Information). As such, water samples were very limited (often <5 L per hydrocast). Therefore, SIP microcosms were established at six stations (indicated by yellow dots) from the epipelagic (60–115 m), the mesopelagic (215–835 m) and the bathypelagic zones (2085–2835 m) at adjacent stations. Duplicate samples (1 L; at stations 30, 45 and 60) were amended with either 13C sodium bicarbonate (5 mM)/12C sodium acetate (20 μM) or 12C sodium bicarbonate (5 mM)/13C sodium acetate (20 μM) as indicated in Fig. 1. Alternatively, duplicate samples (1 L; at stations 33, 48, 63) were amended with either 13C sodium acetate (20 μM)/12C urea (20 μM) or 13C urea (20 μM)/12C sodium acetate (20 μM). These 13C-substrates and concentrations were chosen based on previous findings in salt marsh sediments where the lower concentrations of 13C-organic substrates produced greater diversity of active thaumarchaea and bathyarchaea (formerly Miscellaneous Crenarchaeal Group) (Seyler, McGuinness and Kerkhof 2014). A complete list of SIP incubations is provided in Table S1 (Supporting Information). Figure 1. View largeDownload slide Classification scheme for determining if 13C-labeled OTUs are autotrophic, heterotrophic or mixotrophic. Duplicate samples at each station/depth are supplemented with the 12C or 13C substrates indicated. Unambiguous growth on the 13C is determined by negative amplification in 12C-amended samples for a particular OTU in the 13C carrier band. Figure 1. View largeDownload slide Classification scheme for determining if 13C-labeled OTUs are autotrophic, heterotrophic or mixotrophic. Duplicate samples at each station/depth are supplemented with the 12C or 13C substrates indicated. Unambiguous growth on the 13C is determined by negative amplification in 12C-amended samples for a particular OTU in the 13C carrier band. The intention of this SIP dual amendment scheme was to probe nearly the same populations of TACK archaea in microcosms from similar stations/depths, while tracing the uptake of specific 13C- labeled substrates. It was assumed that TACK communities at adjacent stations (e.g. stations 45 and 48) were nearly identical. Furthermore, it was assumed that a 48-h time point would be sufficient to assess 13C-growth, based on previous data from the culturing of ammonia-oxidizing archaea (Könneke et al.2005) and personal observations in salt marsh SIP incubations. Because of the limited sample volume, the reciprocal 12C/13C amendment experimental design maximized the possibility of obtaining a ‘positive’ signal from the SIP incubations while delineating ‘negative’ signals. For example, if an OTU was observed in the heavy band for 13C-acetate amendments, that OTU should not be observed in the heavy band with 12C-acetate amendments if the TACK archaeum can only utilize acetate (regardless of additional 13C-amendments with different substrates). Likewise, a similar pattern should be observed with 13C-urea or 13C-bicarbonate if that particular substrate is a sole carbon source. However, this assumption will break down if both 13C-acetate and 13C-urea can be utilized or 13C-acetate and 13C-bicarbonate. This dual amendment approach was also implemented because if traditional 12C negative controls were employed for each station/sample, only a single 13C-carbon source could be tested on the cruise, given the limited amount of water available at each hydrocast station. However, it was possible early in the cruise to establish duplicate 12C-controls and 13C-amendments for acetate at station 2. These station 2 12C and 13C-SIP experiments were incubated for 12, 24, 36, 48 and 72 h to verify the minimum time necessary to obtain a positive 13C-signal for at least a single station. The negative SIP control using 12C sodium acetate was incubated for 72 h. Most negative SIP controls by the criteria outlined above did not amplify under the conditions used in this study. Finally, bicarbonate was only employed in SIP incubations from the epipelagic and mesopelagic zones based on a prior report that found archaeal ammonia monooxygenase [amoA] gene copy numbers were not strongly associated with CO2 fixation rates (e.g. Agogué et al.2008). 13C-Bicarbonate was added at ∼2.5 times the average concentration in seawater as described above so the 13C-uptake signal would not be obscured by high ambient levels of 12C bicarbonate (Seyler, McGuinness and Kerkhof 2014). 12C-Bicarbonate was also added to the 13C-acetate amendment so any small change in pH did not suppress microbial activity. After the amendments, the SIP microcosms were incubated in covered containers on deck for 48 h to promote 13C-incorporation, with a constant inflow of surface sea water to maintain stable temperature (on average, 22°C). Biomass was collected from all SIP microcosms on a 0.2-μm filter using vacuum filtration. DNA extractions were performed using modified phenol-chloroform methods as described previously (Sakano and Kerkhof 1998). Frozen filters were suspended in 100 μl buffer (50 mM glucose–10 mM EDTA–25 mM Tris [pH 8.0]) and rapidly frozen with liquid nitrogen and thawed in a 55°C bath four times. Buffer (125 μl), 75 μl EDTA and 100 μl lysozyme solution (10 mg/ml) were then added, and the samples were agitated for 5 min. Fifty microliters of a 50% sodium dodecyl sulfate solution were then added to lyse the cells. The lysate was extracted twice with 800 μl of a phenol-chloroform-isoamyl alcohol mixture (25:24:1). DNA was precipitated by adding 50 μl of 3.0 M sodium acetate, 2 μl of glycogen (20 mg/ml) and 1 ml of 100% ethanol. The DNA was pelleted by centrifugation (16 000 × g) at 4°C for 20 min, and resuspended in 100 μl sterile ultrapure water. DNA (300 ng) from each sample was fractionated into 12C and 13C bands by isopycnic cesium chloride gradient ultracentrifugation at 80 000 rpm (200 000 × g) for 36 h, in 600 μl of 0.9 g/ml CsCl with 300 ng 13C-labeled E. coli DNA as a carrier (modified from Gallagher et al.2005). The use of carrier DNA produces a visible 13C-lower band in the gradient, increases sensitivity and precludes the need to fractionate the gradient (Gallagher et al.2010; Kerkhof et al.2011; Seyler, McGuinness and Kerkhof 2014; Tuorto et al.2014). Ten microliters were collected by pipette from the 12C (upper) and 13C (lower) bands, dialyzed against 10 mM Tris pH 8, and 2 μl of each sample were amplified by nested PCR with 5΄-fluorescently-labeled, 16S rRNA archaea-specific 21F (5΄-TCCGGTTGATCCYGCCGG) and 958R (5΄-YCCGGCGTTGAMTCCAATT) and thaum/crenarchaeota-specific 7F (5΄-TTCCGGTTGATCCYGCCGGACC) and 518R (5΄-GCTGGTWTTACCGCGGCGGCTGA) (Perevalova et al.2003; Ionescu et al.2009; Hong, Kim and Cho 2014; Seyler, McGuinness and Kerkhof 2014) primers. Amplicons (15 ng) were digested with MnlI in 20 μl volumes for 6 h at 37°C, and then precipitated using sodium acetate and 95% ethanol (McGuinness et al.2006). Precipitated DNA was dried and resuspended in 19.7 μl de-ionized formamide with 0.3 μl ROX 500 size standard (Applied Biosystems (ABI); Foster City, CA). TRFLP fingerprinting was carried out on an ABI 310 genetic analyzer using Genescan software. Peak detection was set at 25 arbitrary fluorescent units. A constant 13C-amplicon mass for SIP/TRFLP profiling was used for direct comparison of the number of active (i.e. 13C label-incorporating) thaumarchaeal OTUs within the various samples. To determine the phylogenetic affiliation of the various 13C-archaeal peaks, an archaeal 13C-amplicon clonal library was constructed from microcosms taken from the mesopelagic at station 63, using the Topo TA cloning kit as per the manufacturer's instruction (Invitrogen, Carlsbad, CA, USA), with 16S rRNA archaea-specific 21F/958R primers. These primers were used for cloning because they overlap with the region for next generation 16S tag sequencing (see below), allowing us to draw connections between next-generation sequencing (NGS) and SIP results. Two hundred recombinant clones were screened in a multiplex format as in Babcock et al. (2007) to find clones (∼900 bp) matching the 13C-SIP peaks. Fourteen clones were identified and sequenced via Sanger methods using M13 primers (Genewiz Inc., South Plainfield, NJ, USA), yielding 11 unique clone sequences. Seven of these unique clones were matched to OTUs obtained via 16S Illumina sequencing using BLAST, and a maximum likelihood phylogenetic tree was reconstructed using 140 unambiguously aligned bases among 95 taxa (including 23 OTUs identified via 16S amplicon sequencing, see below) with Geneious analysis software (Guindon and Gascuel. 2003; Drummond et al.2009). Genbank accession numbers for the cloned sequences are KT372904–KT372917. DNA extractions for NGS were performed on a total of 55 filtered samples taken along the Atlantic transect at 16 stations from depths ranging from the surface to 5450 m (Table S2, Supporting Information), using the PowerWater kit from QIAGEN (Germantown, MD), as per the Earth Microbiome Project protocol (http://www.earthmicrobiome.org/emp-standard-protocols/dna-extraction-protocol/). 16S rRNA gene amplicon sequencing was performed using the Illumina MiSeq platform on each sample using 515F (5΄-AATGATACGGCGACCACCG AGATCTACACTATGGTAATTGTGTGCCAGCMGCCGCGGTAA) and 806R barcoded (5΄-CAAGCAGAAGACGG CATACGAGATXXXXXXXXXXXX AGTCAGTCAGCCGG ACTACHVGGGTWTCTAAT) primers (Caporaso et al.2011). Sample processing, sequencing and core amplicon data analysis were performed by the Earth Microbiome Project (www.earthmicrobiome.org; Thompson et al.2017), and all amplicon sequence data and metadata have been made available through the public data portal (qiita.microbio.me/emp). The sequencing data from two samples were discarded via single rarefaction criteria. Data from the remaining 53 samples has been archived at the European Nucleotide Database (Study: ERP016287). Taxa were identified via comparison with the GreenGenes database (http://greengenes.lbl.gov), and beta diversity (weighted and unweighted unifrac) was summarized using QIIME version 1.7 (Caporaso et al.2010). The 23 OTUs included in the phylogenetic tree were identified by an analysis of variance (ANOVA) test over depth regime (with a Bonferroni correction value of <0.05). Quantitative PCR was performed on a subset (22 total) of these amplicon-sequenced samples to verify if the relative abundance of archaea and bacteria was consistent with other reports on archaeal community structure and if similar patterns of increases/decreases in archaeal signals with depth were observed. Agilent Technologies (Santa Clara, CA) Brilliant II SYBR Green qPCR Low ROX Master Mix was used following the manufacturer's protocol, using archaeal 16S A915F (5΄-AGGAATTGGCGGGGGAGCAC; DeLong 1992)/A1059R (5΄-GCCATGCACCWCCTCT; Yu et al.2005) and bacterial 16S 1114F (5΄-CGGCAACGAGCGCAACCC)/1200R (5΄-CCATTGTAGCACGTGTGTAGCC; Reysenbach and Pace 1995) primers on a Stratagene MX3500P machine. Standard curves were generated using full-length bacterial or archaeal 16S rRNA genes from estuary surface sediments cloned into a vector and purified using a Zymo Research Zyppy Plasmid Prep Kit (Irvine, CA). Because the bacterial 16S primers amplified both the bacterial and archaeal positive controls, these primers were treated as universal prokaryote primers during analysis of the qPCR data. RESULTS The 16S rRNA gene amplicon sequencing results indicated bacteria comprised the majority of the total prokaryotic community, averaging 94% with a range of 83%–99%. Archaea ranged between 0.5% and 16% of the total community, averaging 5% (Fig. S2, Supporting Information). Higher percentages of archaea were most often found in the mesopelagic region (200–1000 m depths). The 16S rRNA gene NGS analysis indicated that archaeal populations were mostly members of phylum Euryarchaeota in the epipelagic zone (60–115 m), and thaumarchaea/crenarchaea between 215 and 5450 m, consistent with previous findings (Karner, DeLong and Karl 2001; Herndl et al.2005). Parvarchaea or ‘other’ archaea were detected in relatively low numbers (between 0.1% and 8%) in the bathypelagic (2085–2835 m) and benthopelagic (4336–5450 m, 10 m above the seafloor) zones (Fig. S3, Supporting Information). Quantitative PCR suggested that the archaea form a greater percentage of the total prokayotic population below the mixocline than would be inferred from the NGS data (Fig. S4, Supporting Information). In many cases, the archaeal percentage determined by qPCR was two to three times greater than the percentage estimated from the NGS sequencing data, with one sample more than an order of magnitude higher, possibly due to the differences in priming sets used for the assays or primer mismatch (Parada, Needham and Furhman 2015). However, the relative increases and decreases in archaeal community contribution with depth between the two methods were largely in agreement with the exception of the 835 m sample (Fig. S4, Supporting Information). SIP results indicated that a minimum incubation time of 48 h was required to observe archaeal DNA in the 13C (lower) band with 13C-amendments and no amplification in the 13C (lower) band with 12C-controls (Fig. S5, Supporting Information). Those microcosms at station 2 that were incubated for 12, 24 and 36 h did not produce signal in the 13C (lower) band using thaum/crenarchaeal primers. TACK 16S rRNA gene amplicons from the 12C (upper) and 13C (lower) bands were profiled by TRFLP, and the resulting peaks were considered OTUs. A total of 76 TACK archaeal OTUs were identified in the CLIVAR cruise samples by TRFLP methods (data not shown). Only 16 of these OTUs were found to be active (i.e. capable of synthesizing 13C-labeled DNA from the various 13C-substrates being provided) in replicate SIP incubations (Table 1). However, 31 of the archaeal OTUs detected along the transect were identified as active in at least one of the SIP incubations (Table S3, Supporting Information). Generally, the highest amplification signal from active TACK OTUs (total arbitrary fluorescence in the entire TRFLP profiles) was observed in the 13C-acetate additions from the mesopelagic zone along the entire transect (Fig. 2). The majority of this 13C-acetate amplification signal was detected in both SIP duplicates. Table 1. Active OTUs detected in both duplicate SIP incubations for a particular depth/station/carbon source are indicated. The different 13C-carbon sources are 13C acetate (light gray), 13C urea (medium gray), and 13C bicarbonate (black). Active OTUs detected in one or both SIP incubations are depicted in Supp. Table 3. View Large Figure 2. View largeDownload slide Histograms of arbitrary fluorescence from active OTUs in the total TRFLP profile based on location in the water column or carbon source (13C-acetate,13C-urea,13C-bicarbonate). Total AFU signal is presented in profiles where a specific 13C-labeled OTU is detected in duplicate incubations (Dups) and in all profiles containing some OTUs only detected in a single incubation (All). The profiles were normalized for the different number of SIP incubations for comparison (as acetate was tested at twice as many stations as urea and bicarbonate). ND = not determined. Figure 2. View largeDownload slide Histograms of arbitrary fluorescence from active OTUs in the total TRFLP profile based on location in the water column or carbon source (13C-acetate,13C-urea,13C-bicarbonate). Total AFU signal is presented in profiles where a specific 13C-labeled OTU is detected in duplicate incubations (Dups) and in all profiles containing some OTUs only detected in a single incubation (All). The profiles were normalized for the different number of SIP incubations for comparison (as acetate was tested at twice as many stations as urea and bicarbonate). ND = not determined. A more detailed breakdown of the active TACK OTUs showed differences with respect to 13C-carbon utilization, depth and station location. No active TACK OTUs were detected in duplicate epipelagic samples (Table 1). The active TACK OTUs in the mesopelagic were 23% autotrophic, 23% heterotrophic and 54% mixotrophic at station 30/33. Only organic carbon uptake was observed in replicates at stations 45/48/60/63 (Fig. 3); however, bicarbonate uptake was observed in one of the replicates in the epipelagic at station 45 and in one bottle in the mesopelagic at station 60 (Table S3, Supporting Information). For the bathypelagic only heterotrophy was tested and most active TACK OTUs incorporated both 13C-acetate and 13C-urea in the SIP incubations (Table 1). For the OTUs that were detected in both duplicate and single SIP incubations, these patterns remained largely the same. In the epipelagic at stations 45/48, there was one active TACK OTU detected on 13C-acetate, one OTU active on 13C-bicarbonate and five OTUs active on 13C-acetate and 13C-bicarbonate. In the mesopelagic, at stations 30/33, 21% of the 19 active OTUs were heterotrophic, 58% mixotrophic and 21% autotrophic; the 9 active OTUs at stations 45/48 were exclusively heterotrophic; and the 17 active OTUs at stations 60/63 were 41% heterotrophic, 59% mixotrophic and 0% autotrophic. Overall, eight OTUs dominated the active TACK SIP/TRFLP profiles with respect to peak area (TRF 66, 82, 83, 116, 117, 124, 223 and 265). Six of these OTUs were mixotrophic (TRFs 66, 83, 117, 124, 223 and 265). The other dominant OTUs (TRF 82 and 116) were only associated with the uptake of 13C-organic carbon. None of the dominant TACK OTUs appeared to be autotrophic. Figure 3. View largeDownload slide Percentage of active OTUs detected in duplicate SIP incubations from the mesopelagic on the various 13C-labeled substrates as indicated. Figure 3. View largeDownload slide Percentage of active OTUs detected in duplicate SIP incubations from the mesopelagic on the various 13C-labeled substrates as indicated. In order to assess the phylogenetic affiliation of the various OTUs, a clone library was generated from the 13C amplicons using universal archaeal primers to obtain longer sequence reads. The long read sequences retrieved in this study were distributed throughout the TACK archaea, and Marine Group II and III euryarchaea. Eleven unique clones were identified, with seven clones being members of the phylum Thaumarchaeota (corresponding to the dominant TRFs 82, 83, 116, 223 and 265 bp; Fig. 4). In addition to the thaumarchaea, four clones were identified as euryarchaeal members of Marine Group II. All four euryarchaeal clones produced a TRF peak of 127 or 128 bp long, clustering closely together in our phylogenetic tree. Figure 4. View largeDownload slide Maximum-likelihood tree of archaeal clones among 95 taxa, including 23 depth-dependent archaeal OTUs from 16S tag sequencing data identified by ANOVA (Bonferroni correction value of <0.05). OTUs that represent >1% of the total archaeal community are indicated with stars; OTUs correlated with oxygen concentration are indicated with an open circle. The numbers indicate prokMSA ID in the GreenGenes database or NCBI accession numbers. Those OTUs not found in GreenGenes are indicated by the word ‘NEW’. OTUs that were identified in the clone library are indicated by TRF[bp length of TRF peak]_[prok MSA ID (if a 100% match was available)]_ [NCBI accession number]. Figure 4. View largeDownload slide Maximum-likelihood tree of archaeal clones among 95 taxa, including 23 depth-dependent archaeal OTUs from 16S tag sequencing data identified by ANOVA (Bonferroni correction value of <0.05). OTUs that represent >1% of the total archaeal community are indicated with stars; OTUs correlated with oxygen concentration are indicated with an open circle. The numbers indicate prokMSA ID in the GreenGenes database or NCBI accession numbers. Those OTUs not found in GreenGenes are indicated by the word ‘NEW’. OTUs that were identified in the clone library are indicated by TRF[bp length of TRF peak]_[prok MSA ID (if a 100% match was available)]_ [NCBI accession number]. A BLAST search of each clone was also performed against the GreenGenes database using the 140-bp fragment covered by the 515-F/806-R universal prokaryote primers to link with the NGS sequencing data. The results of this search identified all of the thaumarchaeal clones as either Cenarchaeum or Nitrosopumilus spp. (with a 91%–94% match) due to the small size of the region used in the search. Many of the thaumarchaeal OTUs identified by 16S tag sequencing were classified as Nitrosopumilus by QIIME using the GreenGenes database (Fig S3). In our phylogenetic tree, the cloned active OTUs clearly fall outside of the cultured Nitrosopumilus representatives. The results suggest the apparent dominance of the archaeal community by Nitrosopumilus at depth is likely an artifact of the short length of the 16S tag sequencing reads, and the lack of identified representative archaea in the GreenGenes database. Beta diversity analysis suggested a correlation between the microbial communities and the water mass where the sample was collected (Fig. S6, Supporting Information). Analysis of similarity of the basin samples (excluding samples taken on the slope or shelf) confirmed that community variance between samples was significantly lower within sampling zones compared to between sampling zones (R = 0.81, P-value = 0.01). No significant correlation was detected between microbial beta-diversity and latitude, depth in meters or temperature. An ANOVA test was used to compare OTU frequencies between sampling zones and determine which archaea were differentially represented in one sampling zone over another. These depth-dependent OTUs are included in the phylogenetic tree (Fig. 4). Similarly, a Pearson correlation test was also performed to assess whether archaeal OTUs relative abundance was strongly correlated with oxygen concentration (in Fig. 4, those OTUs with a Bonferroni correction value of <0.05 are indicated with open circles). DISCUSSION Much of the recent research focus on marine Thaumarchaeota has concentrated on life strategies revolving around ammonia oxidation and autotrophy. Our SIP results demonstrate that inorganic carbon fixation represents a subset of the full metabolic capability of the TACK superphylum as a whole. Only two of the observed active TACK OTUs exclusively incorporated bicarbonate, while half of the TACK OTUs readily incorporated both organic and inorganic carbon. These mixotrophic OTUs represented the largest and most frequently observed peaks in TRFLP analysis. However, it should be noted that multiple TACK 16S rRNA genes could be attributed to a single TRFLP peak (e.g. TRF 223 and TRF 265) and our results could represent an autotrophic and a heterotrophic microorganism rather than two mixotrophic TACK archaea. In a prior study, strict heterotrophy was found to be the prevalent metabolic strategy among the TACK archaea in salt marsh sediments (Seyler, McGuinness and Kerkhof 2014). Likewise, no incorporation of 13C-bicarbonate was detected in the mesopelagic in the Sargasso Sea (stations 45/48; Fig. 3), suggesting that TACK archaea at this location/depth may also adopt a strictly heterotrophic lifestyle, while at other depths/locations these same TRFs were associated with the uptake of both organic and inorganic carbon (Table S3A and B, Supporting Information). However, it should be noted that we did not supplement our 13C-bicarbonate SIP incubations with NH4+ and we may detect increased autotrophy with an increase in ammonium concentration. Unfortunately, the use of carrier DNA in our SIP gradients precludes the use of NGS with universal primers, as most of the data retrieved from such methods would be carrier sequence. The uptake of urea by Thaumarchaeota, particularly in the mesopelagic and bathypelagic, is consistent with the high abundance of urease genes that have been found in this phylum (Alonso-Sáez et al.2012). In our study, incorporation of 13C-urea and 13C-acetate but not 13C-bicarbonate was detected in the mesopelagic at stations 45/48 (Fig. 3). Roughly half of these OTUs did not incorporate 13C-bicarbonate at any station or depth. As with our prior salt marsh study (Seyler, McGuinness and Kerkhof 2014), the 13C-urea is not being incorporated via 13CO2. If that were the case, we should always observe 13C-labeled DNA when the samples are incubated with 13C-bicarbonate. Furthermore, the release of micromolar 13CO2 from 13C-urea would be significantly diluted by the ambient 12C-millimolar bicarbonate in seawater. We should never observe 13C-DNA from 13C-urea if autotrophy were the only pathway to incorporation. Rather this 13C-incorporation from urea in the North Atlantic is consistent with our salt marsh study in which we concluded that TACK archaea run the ornithine cycle in reverse (Seyler, McGuinness and Kerkhof 2014). This incorporation of urea (and acetate) suggests that the role of Thaumarchaeota in organic carbon cycling and remineralization, particularly at depth, may be vital yet underappreciated. Our SIP results also demonstrated that nearly half of the North Atlantic TACK archaeal community are active (40% of the OTUs detected by SIP/TRFLP) and can double within 24 h (two cellular divisions are required to detect a signal in the 100% labeled 13C carrier band; Seyler, McGuinness and Kerkhof 2014). It is possible that the high percentage of active TACK archaea we observe in our SIP/TRFLP studies is a result of our 13C-incubation conditions, which differed markedly from in situ conditions with respect to substrate concentration, temperature and pressure. Yet, changes in incubation temperature, pressure or available substrate cannot bestow a metabolic potential in the TACK archaea where one did not exist in the first place. It is unreasonable to conclude that we are mutating or transforming the TACK archaea and creating a metabolic potential where one did not exist beforehand in our SIP/TRFLP incubations. However, it seems possible that we may be inducing TACK populations that are inactive in situ by incubating at 1 atm and 22°C for 48 h rather than in situ conditions. Alternatively, it also is possible that we could not assign activity to members of the in situ TACK archaeal community simply because of our SIP incubation conditions, including the requirement of a rapid doubling time. Further SIP-type experiments are needed to definitively show which members of the TACK archaeal community are active under in situ conditions. In conclusion, our SIP results demonstrate that inorganic carbon fixation represents a subset of the full metabolic capability of the TACK superphylum as a whole. Our SIP findings indicate that the mesopelagic zone had the highest diversity of TACK archaea which can utilize the 13C-substrates being tested. Roughly half of these TACK OTUs in the North Atlantic incorporated only 13C-acetate and/or 13C-urea for cellular carbon, while the other half (the most dominant members of the TACK community based on SIP/TRFLP profiles) incorporated both organic (either/both 13C-acetate and 13C-urea) and inorganic (13C-bicarbonate) carbon. Only two TACK archaeal OTUs exclusively incorporated 13C-bicarbonate. Given the diversity of TACK OTUs in the mid- to deep ocean and the metabolic capabilities that can be discerned by these SIP methods, the role of archaea should also be considered in future studies on the recycling of dissolved organic matter and inorganic carbon fixation in the meso- and bathypelagic zones. These findings have direct implications concerning the structure of the marine microbial food web, which is a major determinant in the rate of oceanic carbon sequestration and the type/amount of higher consumers in the ecosystem (Limardo and Worden 2015). Comparable studies to directly measure the uptake of carbon substrates will help predict the fate of dissolved and particulate organic carbon in marine ecosystems. SUPPLEMENTARY DATA Supplementary data are available at FEMSEC online. Acknowledgements The authors wish to thank Dr Michael McCartney for serving as principle investigator onboard the R/V Atlantis during CLIVAR leg A20, Jarrad Hampton-Marcell of Argonne National Laboratory for his assistance with 16S rRNA gene tag sequencing and data processing using QIIME, and Kristin Yoshimura of the University of Delaware for her work involving quantitative-PCR. FUNDING This research was made possible through the support of the International Global Ocean Carbon and Repeat Hydrography Program, the Earth Microbiome Project, and NSF Ocean Technology and Interdisciplinary Program grant (#1131022) to LJK. Conflict of interest. None declared. REFERENCES Agogué H, Brink M, Dinasquet J et al. Major gradients in putatively nitrifying and non-nitrifying Archaea in the deep North Atlantic. Nature 2008; 456: 788– 91. Google Scholar CrossRef Search ADS PubMed Alonso-Sáez L, Waller AS, Mende DR et al. Role for urea in nitrification by polar marine Archaea. P Natl Acad Sci USA 2012; 10: 17989– 94. Google Scholar CrossRef Search ADS Avaniss-Aghajani E, Jones K, Chapman D et al. A molecular technique for identification of bacteria using small subunit ribosomal RNA sequences. BioTechniques 1994; 17: 144– 6. Google Scholar PubMed Babcock DA, Wawrik B, Paul JH et al. Rapid screening of a large insert BAC library for specific 16S rRNA genes using TRFLP. J Microbiol Methods 2007; 71: 156– 61. Google Scholar CrossRef Search ADS PubMed Berg C, Listmann L, Vandieken V et al. Chemoautotrophic growth of ammonia-oxidizing Thaumarchaeota enriched from a pelagic redox gradient in the Baltic Sea. Front Microbiol 2015; 5: 786. Google Scholar CrossRef Search ADS PubMed Buckley D, Huangyutitham V, Hsu SF et al. Stable isotope probing with 15N achieved by disentangling the effects of genome G+C content and isotope enrichment on DNA density. Appl Environ Microb 2007; 73: 3189– 95. Google Scholar CrossRef Search ADS Caporaso JG, Lauber CL, Walters WA et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. P Natl Acad Sci USA 2011; 108: 4516– 22. Google Scholar CrossRef Search ADS Caporaso JG, Kuczynski J, Stombaugh J et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 2010, DOI: https://doi.org/10.1038/nmeth.f.3. DeLong EF. Archaea in coastal marine environments. P Natl Acad Sci USA 1992; 89: 5685– 9. Google Scholar CrossRef Search ADS DeLong EF, Wu KY, Prézelin BB et al. High abundance of Archaea in Antarctic marine picoplankton. Nature 1994; 371: 695– 7. Google Scholar CrossRef Search ADS PubMed Drummond AJ, Ashton B, Cheung M et al. Geneious v4.8. Auckland, New Zealand: Biomatters Ltd., 2009. Francis CA, Beman JM, Kuypers MMM. New processes and players in the nitrogen cycle: the microbial ecology of anaerobic and archaeal ammonia oxidation. ISME J 2007; 1: 19– 27. Google Scholar CrossRef Search ADS PubMed Francis CA, Roberts KJ, Beman JM et al. Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. P Natl Acad Sci USA 2005; 102: 14683– 8. Google Scholar CrossRef Search ADS Fuhrman JA, McCallum K, Davis AA. Novel major archaebacterial group from marine plankton. Nature 1992; 356: 148– 9. Google Scholar CrossRef Search ADS PubMed Gallagher E, McGuinness L, Phelps C et al. 13C-Carrier DNA shortens the incubation time needed to detect benzoate-utilizing denitrifying bacteria by stable-isotope probing. Appl Environ Microb 2005; 71: 5192– 6. Google Scholar CrossRef Search ADS Gallagher EM, Young LY, McGuinness LM et al. Detection of 2,4,6-trinitrotoluene-utilizing anaerobic bacteria by 15N and 13C incorporation. Appl Environ Microb 2010; 76: 1695– 8. Google Scholar CrossRef Search ADS Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003; 52: 696– 704. Google Scholar CrossRef Search ADS PubMed Guy L, Ettema TJG. The archaeal ‘TACK’ superphylum and the origin of eukaryotes. Trends Microbiol 2011; 19: 580– 7. Google Scholar CrossRef Search ADS PubMed Hallam SJ, Mincer TJ, Schleper C et al. Pathways of carbon assimilation and ammonia oxidation suggested by environmental genomic analyses of marine Crenarchaeota. PLoS Biol 2006; 4: 520– 36. Google Scholar CrossRef Search ADS Herndl GJ, Reinthaler T, Teira E et al. Contribution of Archaea to total prokaryotic production in the deep Atlantic Ocean. Appl Environ Microb 2005; 71: 2303– 9. Google Scholar CrossRef Search ADS Hong J-K, Kim J-H, Cho J-C. Novel PCR primers for the archaeal phylum Thaumarchaeota designed based on the comparative analysis of 16S rRNA gene sequences. PLoS One 2014; DOI: https://doi.org/10.1371/journal.pone.0096197. Ingalls AE, Shah SR, Hansman RL et al. Quantifying archaeal community autotrophy in the mesopelagic ocean using natural radiocarbon. P Natl Acad Sci USA 2006; 103: 6442– 7. Google Scholar CrossRef Search ADS Ionescu D, Penn S, Haimovich M et al. Archaea in the Gulf of Aqaba. FEMS Microbiol Ecol 2009; 69: 425– 38. Google Scholar CrossRef Search ADS PubMed Jung MY, Kim JG, Sinninghe Damste JS et al. A hydrophobic ammonia-oxidizing archaeon of the Nitrosocosmicus clade isolated from coal tar-contaminated sediment. Environ Microbiol Rep 2016; 8: 983– 92. Google Scholar CrossRef Search ADS PubMed Jung MY, Park SJ, Kim SJ et al. A mesophilic, autotrophic, ammonia-oxidizing archaeon of thaumarchaeal group I.1a cultivated from a deep oligotrophic soil horizon. Appl Environ Microb 2014; 80: 3645– 55. Google Scholar CrossRef Search ADS Karner MB, DeLong EF, Karl DM. Archaeal dominance in the mesopelagic zone of the Pacific Ocean. Nature 2001; 409: 507– 10. Google Scholar CrossRef Search ADS PubMed Kerkhof LJ, Williams KH, Long PE et al. Phase preference by active, acetate-utilizing bacteria at the Rifle, CO integrated field research challenge site. Environ Sci Technol 2011; 45: 1250– 6. Google Scholar CrossRef Search ADS PubMed Kim J-G, Park S-J, Sinninghe Damsté JS et al. Hydrogen peroxide detoxification is a key mechanism for growth of ammonia-oxidizing archaea. P Natl Acad Sci USA 2016; 113: 7888– 93. Google Scholar CrossRef Search ADS Kirchman DL, Elifantz H, Dittel AI et al. Standing stocks and activity of Archaea and Bacteria in the western Arctic Ocean. Limnol Oceanogr 2007; DOI: https://doi.org/10.4319/lo.2007.52.2.0495. Könneke M, Bernhard AE, de la Torre JR et al. Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature 2005; 437: 543– 6. Google Scholar CrossRef Search ADS PubMed Limardo AJ, Worden AZ. Exclusive networks in the sea. Nature 2015; 522: 36– 37. Google Scholar CrossRef Search ADS PubMed McGuinness LM, Salganik M, Vega L et al. Replicability of bacterial communities in denitrifying bioreactors as measured by PCR/T-RFLP analysis. Environ Sci Technol 2006; 40: 509– 15. Google Scholar CrossRef Search ADS PubMed Martin-Cuadrado AB, Rodriguez-Valera F, Moreira D et al. Hindsight in the relative abundance, metabolic potential and genome dynamics of uncultivated marine archaea from comparative metagenomic analyses of bathypelagic plankton of different oceanic regions. ISME J 2008; 2: 865– 86. Google Scholar CrossRef Search ADS PubMed Ouverney CC, Fuhrman JA. Marine planktonic archaea take up amino acids. Appl Environ Microb 2000; 66: 4829– 33. Google Scholar CrossRef Search ADS Parada A, Needham DM, Furhman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time-series and global field samples. Environ Microbiol 2015, DOI: https://doi.org/10.1111/1462-2920.13023. Park BJ, Park SJ, Yoon DN et al. Cultivation of autotrophic ammonia-oxidizing archaea from marine sediments in coculture with sulfur-oxidizing bacteria. Appl Environ Microb 2010; 76: 7575– 87 Google Scholar CrossRef Search ADS Perevalova AA, Lebedinsky AV, Bonch-Osmolovskaya EA et al. Deteaction of hyperthermophilic Archaea of the genus Desulfurococcus by hybridization with oligonucleotide probes. Microbiology 2003; 72: 340– 6. Google Scholar CrossRef Search ADS Qin W, Amin SA, Martens-Habbena W et al. Marine ammonia-oxidizing archaeal isolates display obligate mixotrophy and wide ecotypic variation. P Natl Acad Sci USA 2014; 111: 12504– 9. Google Scholar CrossRef Search ADS Radajewski S, Ineson P, Parekh NR et al. Stable-isotope probing as a tool in microbial ecology. Nature 2000; 403: 646– 9. Google Scholar CrossRef Search ADS PubMed Reysenbach AL, Pace NR. Reliable amplification of hyperthermophilic Archaeal 16S rRNA genes by polymerase chain reaction. In: Robb FT, Place AR, Sowers KR et al. (eds). Archaea, A Laboratory Manual . New York: Cold Spring Harbor Press, 1995, 101– 6. Sakano Y, Kerkhof L. Assessment of changes in microbial community structure during operation of an ammonia biofilter with molecular tools. Appl Environ Microb 1998; 64: 4877– 82. Seyler LM, McGuinness LM, Kerkhof LJ. Crenarchaeal heterotrophy in salt marsh sediments. ISME J 2014; 8: 1534– 43. Google Scholar CrossRef Search ADS PubMed Stein JL, Simon MI. Archaeal ubiquity. P Natl Acad Sci USA 1996; 93: 6228– 30. Google Scholar CrossRef Search ADS Teira E, Reinthaler T, Pernthaler A et al. Combining catalyzed reporter deposition-fluorescence in situ hybridization and microautoradiography to detect substrate utilization by bacteria and archaea in the deep ocean. Appl Environ Microb 2004; 70: 4411– 4. Google Scholar CrossRef Search ADS Teira E, van Aken H, Veth C et al. Archaeal uptake of enantiomeric amino acids in the meso- and bathypelagic waters of the North Atlantic. Limnol Oceanogr 2006; 51: 60– 69. Google Scholar CrossRef Search ADS Thompson LR, Sanger JG, McDonald D et al. A communal catalogue reveals Earth's multiscale microbial diversity. Nature 2017; 551: 457– 63. Google Scholar CrossRef Search ADS PubMed Treusch AH, Leininger S, Kletzin A et al. Novel genes for nitrite reductase and Amo-related proteins indicate a role of uncultivated mesophilic crenarchaeota in nitrogen cycling. Environ Microbiol 2005; 7: 1985– 95. Google Scholar CrossRef Search ADS PubMed Tuorto SJ, Darias P, McGuinness LR et al. Bacterial genome replication at subzero temperatures in permafrost. ISME J 2014; 8: 139– 49. Google Scholar CrossRef Search ADS PubMed Venter JC, Remington K, Heidelberg JF et al. Environmental genome shotgun sequencing of the Sargasso Sea. Science 2004; 304: 66– 74. Google Scholar CrossRef Search ADS PubMed Wanninkhof R, Feely R, Dickson A et al. Carbon Dioxide, Hydrographic, And Chemical Data Obtained During the R/V Atlantis Cruise in the Atlantic Ocean on CLIVAR Repeat Hydrography Section A20 (Apr. 19 - May 15, 2012). http://cdiac.ornl.gov/ftp/oceans/CLIVAR/A20_2012/. 2013, Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tennessee. Wuchter C, Schouten S, Boschker HTS et al. Bicarbonate uptake by marine Crenarchaeota. FEMS Microbiol Lett 2003; 219: 203– 7. Google Scholar CrossRef Search ADS PubMed Yu Y, Lee C, Kim J et al. Group-specific primer and probe sets to detect methanogenic communities using quantitative real-time polymerase chain reaction. Biotechnol Bioeng 2005; 89: 670– 9. Google Scholar CrossRef Search ADS PubMed © FEMS 2018. All rights reserved. For permissions, please e-mail: firstname.lastname@example.org
FEMS Microbiology Ecology – Oxford University Press
Published: Mar 1, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
All the latest content is available, no embargo periods.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud