ORIGINAL RESEARCH article
Distinct Microbial Assemblage Structure and Archaeal Diversity in Sediments of Arctic Thermokarst Lakes Differing in Methane Sources
- 1Division of Earth and Ecosystem Sciences, Desert Research Institute, Reno, NV, United States
- 2Department of Microbiology and Ecosystem Science, University of Vienna, Vienna, Austria
- 3Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, United States
- 4Department of Land Resources and Environmental Sciences, Montana State University, Bozeman, MT, United States
Developing a microbial ecological understanding of Arctic thermokarst lake sediments in a geochemical context is an essential first step toward comprehending the contributions of these systems to greenhouse gas emissions, and understanding how they may shift as a result of long term changes in climate. In light of this, we set out to study microbial diversity and structure in sediments from four shallow thermokarst lakes in the Arctic Coastal Plain of Alaska. Sediments from one of these lakes (Sukok) emit methane (CH4) of thermogenic origin, as expected for an area with natural gas reserves. However, sediments from a lake 10 km to the North West (Siqlukaq) produce CH4 of biogenic origin. Sukok and Siqlukaq were chosen among the four lakes surveyed to test the hypothesis that active CH4-producing organisms (methanogens) would reflect the distribution of CH4 gas levels in the sediments. We first examined the structure of the little known microbial community inhabiting the thaw bulb of arctic thermokarst lakes near Barrow, AK. Molecular approaches (PCR-DGGE and iTag sequencing) targeting the SSU rRNA gene and rRNA molecule were used to profile diversity, assemblage structure, and identify potentially active members of the microbial assemblages. Overall, the potentially active (rRNA dominant) fraction included taxa that have also been detected in other permafrost environments (e.g., Bacteroidetes, Actinobacteria, Nitrospirae, Chloroflexi, and others). In addition, Siqlukaq sediments were unique compared to the other sites, in that they harbored CH4-cycling organisms (i.e., methanogenic Archaea and methanotrophic Bacteria), as well as bacteria potentially involved in N cycling (e.g., Nitrospirae) whereas Sukok sediments were dominated by taxa typically involved in photosynthesis and biogeochemical sulfur (S) transformations. This study revealed a high degree of archaeal phylogenetic diversity in addition to CH4-producing archaea, which spanned nearly the phylogenetic extent of currently recognized Archaea phyla (e.g., Euryarchaeota, Bathyarchaeota, Thaumarchaeota, Woesearchaeota, Pacearchaeota, and others). Together these results shed light on expansive bacterial and archaeal diversity in Arctic thermokarst lakes and suggest important differences in biogeochemical potential in contrasting Arctic thermokarst lake sediment ecosystems.
The thawing front (talik or thaw bulb) of permafrost underneath shallow (~2–3 m deep) arctic thermokarst lakes is an active area of carbon (C) cycling, where carbon dioxide (CO2) and methane (CH4) are produced and consumed by microbially driven processes (Bretz and Whalen, 2014; Lofton et al., 2014; Heslop et al., 2015). Knowledge of the organisms and communities driving these processes has only recently begun to reveal the extent of microbial communities diversity and structure in thermokarst lakes or similar saturated permafrost features. Initial findings suggest that CH4 cycling microorganisms are important components of these systems. For example, type I and type II bacterial methanotrophs were found in the sediments of an Alaskan Arctic thermokarst lake (Qalluuraq Lake) that emits ebullient CH4 (He et al., 2012b). In another study, CH4 producers (methanogens) affiliated with Methanomicrobiales and Methanosarcinales were observed in sediments of small, shallow ponds in the Canadian Arctic (Negandhi et al., 2013). Likewise, abundant methanogen SSU rRNA gene sequences and high methanogenesis transcript to gene ratios were detected in thermokarst bog soils in the Tanana Flats in central Alaska (Hultman et al., 2015). Denitrifiers and sulfate () reducers co-occurred in the bog soils, demonstrating the interconnectedness between biogeochemical processes in the soils.
The Coastal Plain of northern Alaska is a large (i.e., 22,000 square kilometers; Wang et al., 2012) matrix of permafrost peppered with thermokarst lakes, some of which produce either ebullient CH4 generated by present-day biogenic processes (e.g., Siqlukaq Lake) or as a result of release by thermogenic processes (e.g., Sukok Lake, which overlies a portion of the national petroleum reserve; Matheus Carnevali et al., 2015). To date, permafrost features in this area rather than thermokarst lakes have received the most attention, given the need for gauging response of permafrost to warming temperatures (Graham et al., 2012; Lipson et al., 2013).
As part of an interdisciplinary study on CH4-related process in this region, we hypothesized that the microbial assemblages in the sediments of the ebullient thermogenic CH4 lakes would be distinct from the biogenic CH4 producing lakes. To address this, we first assessed archaeal diversity across four arctic thermokarst lakes (three dominated with biogenic CH4 sources and one with a thermogenic source; Matheus Carnevali, 2015; Matheus Carnevali et al., 2015; Elvert et al., 2016; Hand, unpublished results), and then based on those results, we targeted our focus to two lakes, Siqlukaq (Siq) and Sukok (Suk). In addition to different CH4 sources, sediments of these two lakes harbor contrasting levels of TOC (between 0.2 and 15 wt. %) (Matheus Carnevali et al., 2015). Our analyses targeted both ribosomal RNA genes and the ribosomal RNA molecule of both bacteria and archaea, provided that DNA from non-viable microbial cells may be present in the environment (Willerslev et al., 2004; Blazewicz et al., 2013). In comparison to mRNA, the more environmentally stable rRNA where detected in appreciable levels, may be indicative of potentially active bacteria and archaea (Juottonen et al., 2008; Stoeva et al., 2014; Crevecoeur et al., 2016). For instance, taxonomic affiliations of SSU rRNA iTags and the mRNA fraction (metatranscriptome) assigned to taxonomic bins have been observed to correspond (Tveit et al., 2012; Hultman et al., 2015). We also took a deeper look into the phylogenetic diversity of the Archaea by sequencing cloned archaeal rRNA genes (partially covered).
Materials and Methods
Sampling Sites and Sediment Core Collection
Shallow (~2–3 m deep) arctic thermokarst lakes on the north slope of Alaska were sampled during three different field campaigns (April and October 2010, and October 2011; Table 1). Three of the lakes sampled (Ikroavik, Sukok, and Siqlukaq) are located southeast of Barrow, and the fourth lake (Qalluuraq) is located further south near Atqasuk. Sediment cores were collected through holes drilled in the frozen surface lake ice (0.68–1.40 m in April 2010, 0.25 m in October 2010, and 0.08–0.15 m in October 2011) using a universal percussion corer (Aquatic Research Instruments). Cores were transported by snow mobile to the Naval Arctic Research Laboratory (NARL) and stored frozen at −80°C until subsampling. In depth analyses here focused on sediment cores collected in October 2011: one site in Siqlukaq (Siq) and two sites in Sukok, a site nearby an open CH4 seep (SukS) and a site ~1 km southwest from the seep (SukB). Sukok is underlain by the Walakpa natural gas field and displays ebullition throughout the year, especially in the SukS area. Matheus Carnevali et al. (2015) conducted the first studies on CH4 source from these lakes. The coring sites, sediment intervals sampled, and associated downstream analyses are listed in Table 1.
Nucleic Acid Extractions
Over the course of this study two DNA extraction protocols and one RNA extraction protocol were executed. The different DNA extraction protocols were used for separate downstream purposes (indicated with numerals i and ii).
DNA extraction (i): DNA from 2011 sediments was used for Bacteria and Archaea SSU rRNA gene iTag sequencing and archaeal SSU rRNA gene clone libraries. We studied the upper sediment layers (0–30 cm) given observed differences in CH4 production in these layers (Matheus Carnevali et al., 2015). DNA was extracted from three sections (U-upper, M-middle, L-lower; ~1 g ea.) of the first 30 cm of each 2011 sediment core (Siq, SukS, and SukB; Table 1) using a Power Soil DNA isolation kit (MoBio). DNA was quantified (Quant-iT; Picogreen Assay Kit; Life Technologies) on a fluorescence reader (Spectromax Gemini, Molecular Devices). All comparisons between sites were made with DNA extracted using this method. Additional subsamples from the same core sections preserved at −80°C in sucrose lysis buffer (SLB, 40 mM EDTA, 50 mM Tris-HCl, 0.75 M sucrose) were used to generate DNA for archaeal rRNA gene clone libraries.
DNA extraction (ii): 2010 sediments were used for archaeal PCR-DGGE profiling and 2010 archaeal SSU rRNA gene clone libraries. Subsamples from 2010 cores (1 g ea.; Table 1) were collected from the top ~0–60 cm of Qal10 (6 samples; April), the top ~0–20 cm of Ikr-10A (3 samples; April), the top ~0–30 cm (2 samples; October) of Ikr10-O, the top ~0–30 cm of SukC (4 samples; April), the top ~10–50 cm of Siq10 (3 samples; October), and along the length (~0–110 cm) of the SukS10 core (10 samples; October). Genomic DNA was extracted using phenol-chloroform purification (Massana et al., 1997) and quantified (Nanodrop, ThermoFisher Scientific).
RNA Extraction: SSU rRNA Bacteria and Archaea iTag Sequencing
To target potentially active members of the microbial assemblage, total RNA was extracted from three sections (U-upper, M-middle, L-lower; 4 g ea.) of the first 30 cm of the 2011 cores (Siq, SukB, and SukS; Table 1). The RNA PowerSoil Total RNA Isolation Kit (MoBio) was used for nucleic acid extraction omitting the purification steps (Supplemental Methods), and included reagent-only extraction controls. Total nucleic acid extracts were treated 2–3 times with RNase-free DNase I (4 units each; Life Technologies) at 37°C in 100 μL reaction volumes for 30 min. Total RNA was purified from DNase enzyme following standard procedures (Sambrook et al., 1989) then quantified (Quant-iT RiboGreen RNA Assay Kit) on a fluorescence reader. Extraction controls were also subjected to DNAse treatment.
Reverse Transcription and RNA Amplification of the SSU rRNA Gene
To obtain cDNA from rRNA, total RNA was reverse transcribed with SuperScript III Reverse Transcriptase (Life Technologies) following manufacturer's instructions using a SSU rRNA-specific directed primer (806R; Caporaso et al., 2010). Siq11-U, -M, and -L samples were diluted 1/10 before the reverse transcription step to prevent enzymatic inhibition observed in previous experiments. cDNA products were verified on an agarose gel and quantified (Quant-iT PicoGreen dsDNA Assay Kit). Reagent controls were also subjected to reverse transcription and PCR in which subsequent amplification products were not detected.
RNA yields from five of the nine samples (Siq-U, Siq-M, SukB-U, SukB-M, and SukB-L) were sufficient to generate cDNA directly from the total RNA to be used in iTag sequencing. We then amplified the total RNA for direct comparison of the sequence-based identities of the potentially active bacteria and archaea across the set of 9 lake sediment samples. The MessageAmp II aRNA Amplification Kit (Life Technologies) was used according to manufacturer's instructions and final RNA products (aRNA) were quantified (Quant-iT RiboGreen RNA Assay Kit, Invitrogen) on a fluorescence plate reader. Amplified RNA was then reverse transcribed to cDNA with primer 806R (Caporaso et al., 2010).
DNA (9 samples) and cDNA from both RNA (5 samples) and aRNA (9 samples) were sent to the Joint Genome Institute (JGI, Walnut Creek, CA) for library preparation and paired-end (2 × 250 bp) MiSeq Illumina sequencing of the variable region 4 (V4) SSU rRNA using primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) following Caporaso et al. (2010). Table 1 lists the samples sequenced. The iTag libraries were submitted to GenBank SRA under accession numbers SRR5313010-SRR5313032.
iTag Sequence Processing and Statistical Analyses
Sequence processing included removal of PhiX contaminants and Illumina adapters, then quality was evaluated using the FastQC software (Babraham Bioinformatics). Trimmomatic (Bolger et al., 2014) was then used in the single end mode to filter the low quality sequences resulting in 3,834,911 (85.8%) out of 4,469,505 sequences that passed the quality-filtering step. Sequence file headers formatted with an in-house script were imported into Mothur (Kozich et al., 2013) where paired end reads were mated into contigs, chimeras were removed, sequences were initially classified into the Greengenes gg_13_8_99 taxonomy (Silva v. 123 was referenced for updated archaeal taxonomy), and sequences were clustered into operational taxonomic units (OTUs) at a distance of 0.03. These sequences were used to create an OTU database containing OTU occurrence, a representative sequence, and the consensus taxonomy (Kozich et al., 2013). The Bacterial taxonomy of OTUs that were identified as “potentially active” was revised later using Silva v. 132 taxonomy.
Read depths between different samples were normalized by bootstrap resampling (500 bootstraps) using R (R Core Development Team, 2015) for β-diversity analyses following the methods and custom scripts developed by Herbold et al. (2014). To reduce the effect of sequencing error, OTUs with 5 sequences or less among all of the samples (DNA, aRNA, and RNA from all sites) were removed prior to resampling, resulting in a total of 1,213,686 sequences among the 23 samples. The sample() command was then used to bootstrap each sample at 75% of the smallest library size using the observed OTU counts as a probability vector. To account for stochastic effects on the observance of rare OTUs, the total size of the sample (observed + hidden) was estimated using the Chao1 index (Chao, 1984). The non-observed proportion was evenly divided between all OTUs (observed and non-observed) so that non-observed OTUs in each sample possessed a non-zero probability of being resampled into a bootstrapped community.
Non-metric multidimensional scaling (NMDS) ordinations of DNA and aRNA sequences were generated using the Bray-Curtis dissimilarity values calculated with 20 re-samplings for each sample. Pairwise differences in Bray-Curtis dissimilarity between sites were tested with the Mann-Whitney test on resampled communities. Significance thresholds were defined by the Benjamini-Hochberg procedure with a false discovery rate of 0.05. These statistical analyses were conducted in R using the vegan package as well as the p.adjust(method = “fdr”) commands included in the base stats package.
aRNA to DNA Ratio as an Indicator of Potential Activity
We used both aRNA and DNA sequence types to understand the variance in composition of potentially active organisms between lake sediment sites. To identify OTUs that may have been potentially active in these lake sediments communities were resampled 100,000 times, at 75% of the smallest library size among DNA and aRNA samples (resample size = 31,339 sequences). During each resampling, the resampled aRNA sequences were compared to resampled DNA sequences for each OTU. P-values reflecting the significance of aRNA sequences being greater than DNA sequences, and thus “activity” was estimated by counting the number of times out of 100,000 the resampled aRNA sequences were less than the resampled DNA sequences. The Benjamini–Hochberg procedure was used to determine a p-value equivalent to false discovery rate (FDR) = 0.05 (p ≤ 0.00019). Thus, if the aRNA values were less than the DNA values less than 19 times out of the 100,000 resamplings, that OTU in that sample was classified as “active.” OTUs that may have been overrepresented as a result of RNA amplification were identified by plotting the five samples for which both RNA and aRNA data sets were available (i.e., OTUs 20, 22, 31, 70, 89, and 170 in SukB; Figure S1). These OTUs were included in downstream analyses, but caution was applied in their interpretation (i.e., OTUs with aRNA sequences >100 and RNA sequences <10 were considered suspicious). Additionally, OTUs whose aRNA sequences were exceptionally high (>1,000; represented as 1,000 on Figure 2) than their DNA sequences (< 2), as well as OTUs with DNA sequences < 5, were not taken into consideration. OTUs with aRNA to DNA ratios >1 were broken into three groups as proxies to enable consistent comparisons of varying aRNA:DNA levels between sites. OTUs with ratios >2 and aRNA sequence occurrence levels over 500 were assigned to Group 1 (G1), those with ratios >2 and aRNA sequence levels between 250 and 500 were assigned to Group 2 (G2), and those with ratios between 1 and 2 with aRNA sequence levels ≥250 and DNA levels ≥250 were assigned to Group 3 (G3). OTUs with aRNA levels ≥1,000 or DNA levels ≤5 were not assigned a G-level given the potential for aRNA over-amplification and RNA to DNA ratio inflation with low DNA signals.
PCR-DGGE of Archaeal SSU rRNA Genes
To survey the structure of the archaeal assemblages in the thermokarst lake sediments, core sections collected from all lakes sampled in 2010 (Table 1) were analyzed by PCR-DGGE analysis of the V3 region of the archaeal SSU rRNA gene. Details are presented in the supplemental methods, though briefly, sample loads were normalized to 600 ng of DNA per sample and run a 8% polyacrylamide, 30–65% denaturing gradient gel for 16 h at 60 volts, then the SYBR-Gold-stained gel was visualized and analyzed using GelCompar II software (Applied Maths).
Archaeal SSU rRNA Gene Clone Library Construction and Phylogenetic Analyses
Archaeal clone libraries were prepared for each of the 2010 and 2011 sets of sediment cores (Table 1) in order to maximize the scope of the phylogenetic investigation. Briefly, sediment DNA extracts of the 2010 and 2011 sediment intervals (5 μL each) were pooled for PCR and non-directional cloning (see Supplemental Methods section for additional details). Sequencing was unidirectional and resulted in sequences that partially overlapped. To obtain a broad coverage of Archaea phyla, the 2010 library was generated using SSU rRNA gene primers 109F (Großkopf et al., 1998) and 915R (Stahl and Amann, 1991) and the 2011 library was generated with primers 20F-958R (Delong, 1992), resulting in 120 high quality Sanger sequences for phylogenetic reconstruction. The cloned SSU rRNA gene sequences were deposited in GenBank under accession numbers KY886013-KY886132.
To study the evolutionary relatedness between thermokarst sediment Archaea and known Archaea lineages, the SSU rRNA gene sequences were aligned (Silva's SINA aligner v 1.2.11 under the Archaea variability profile; Pruesse et al., 2012) with a reference set of archaeal rRNA genes. The reference set was built from the closest matches from GenBank, Silva (123), and cultured organisms for each OTU from this study, as well as one representative from each genus for all genera with a sequenced genome in either the NCBI or JGI-IMG databases (Hug et al., 2016). The alignment was edited as follows: columns containing ≥95% gaps were removed, and the alignment was trimmed at the 5′ and 3′ ends based on this study's sequence lengths, to a final length of 662 bp using Geneious version 8 (Kearse et al., 2012). The alignment contained a total of 267 sequences, including 77 from this study. A maximum likelihood phylogenetic tree was built using the RAxML-HPC (Stamatakis, 2014) BlackBox algorithm with default parameters for nucleotide alignments, available through the CIPRES (Cyberinfracstructure for Phylogentic Research) Science Getaway (https://www.phylo.org/). Given the short amplicon size of the archaeal iTag sequences, they were not used in tree reconstruction. However, we included SSU rRNA gene sequences of the closest neighbors (>90% identity) to the iTag sequences present in the Silva database (123) in the phylogenetic tree.
β-Diversity of the Microbial Assemblage Among Lake Study Sites
Sequences were clustered into operational taxonomic units (OTUs; distance of 0.03) by sequence type [i.e., DNA, amplified RNA (aRNA), RNA] forming an average of 1,397 ± 288 DNA OTUs, 1,397 ± 366 aRNA OTUs, and 1,447 ± 232 RNA OTUs; excluding doubletons and singletons (Table S1). Patterns of dissimilarity within and among sites (β-diversity) were revealed with non-metric multidimensional scaling (NMDS) ordination analysis of Bray-Curtis dissimilarity (20 in-silico re-samplings) of DNA and amplified RNA (aRNA) iTag occurrences from the three study sites (Figure 1). Comparisons of DNA iTags suggested distinctions of the microbial composition—such that differences were significant between Siq and SukS (Mann–Whitney test, p = 0.0001, significant in 1,000/1,000 resamples), compared to the difference between Siq and SukB. Significant differences were also found between SukB and SukS (Mann-Whitney test, p = 0.0053, significant in 1,000/1,000 resamples). The NMDS analysis based on aRNA iTag sequences (Figure 1B) showed that Siq samples clustered together, compared to SukB and SukS data sets. However, the Mann-Whitney test indicated that the aRNA sequences from the three sites were not significantly different over all in-silico re-samplings.
Figure 1. Non-metric multidimensional scaling of community dissimilarity (Bray-Curtis) calculated from (A) extracted DNA; stress (0.08) and (B) amplified RNA (aRNA) stress (0.138).
Microbial Assemblage Composition of the Potentially Active Organisms
We observed a clear, positive relationship between aRNA and DNA iTag levels across the different taxa detected. The integrity of the relationship between the aRNA and RNA data sets showed best coherence for those sequences that were highly represented in the data set, while the relationship fell off at lower levels of detection (Figure S1). 7 G1 OTUs, 23 G2 OTUs, and 4 G3 OTUs with aRNA to DNA ratios >1 were identified as taxa with potentially important functions in these coring sites (Figure 2, Table 2). Eleven to sixteen potentially active OTUs and different taxonomic groups (categorized to the family level where possible) fell into G1-3 at each site, with 9 taxonomic groups identified in Siq, 8 in SukB, and 11 in SukS. A higher number of OTUs unique to a site (7) were found in Siq and SukB (6), while SukS, only had 4. The two Sukok sites shared 6 OTUs, only one of which was also active in Siq. In fact, unclassified OTUs were the only ones found to be common among the three coring sites (Siq G2, SukB G1-3, and SukS G2). The SukS data set was generally skewed toward OTUs with higher ratios, as the DNA sequence occurrence levels were low on average.
Figure 2. (A–C) aRNA resampled iTags that were significantly higher (p < 0.00019) than DNA resampled iTags at each of the three sites (panels: Siq, SukB, and SukS). Taxonomic affiliations to Class and/or Order are indicated by symbols in the inset legend. Three groups signaling OTUs with potentially relevant functions were identified as follows. Group 1 (G1): aRNA to DNA ratios > 2 and aRNA occurrence levels over 500. G2 OTUs: aRNA to DNA ratios > 2 and aRNA levels between 250 and 500. G3 abundant OTUs: aRNA to DNA ratios between 1 and 2 with aRNA levels ≥250 and DNA levels ≥250. OTUs whose aRNA iTags were exceptionally higher (>1,000) than their DNA counts (≤5) were set to 1,000. These corresponded to G1 OTUs including an Actinobacteria OTUs in Siq, two Alphaproteobacteria OTUs (Siq and SukS), and one Cyanobacteria OTU and one Gammaproteobacteria OTU in SukB with aRNA occurrences above 1,000 and DNA < 5. Other G1 OTUs with DNA occurrences <5 included three Actinobacteria OTUs (SukS), two Alphaproteobacteria OTUs (SukB and SukS), one Bacteroidetes OTU (Siq), two Cyanobacteria OTUs (SukB and SukS), three Gammaproteobacteria OTUs (Siq, SukB, SukS); and one Firmicutes OTU (SukS) and an Actinobacteria OTU (Siq) with aRNA occurrences ~900.
Of the OTUs unique to Siq, one G1 OTU was affiliated with the Nitrospirae (OTU 4 in Siq-U), while two OTUs related to Methylomonaceae fell into G2 (OTUs 35) and G3 (OTU 21) (Figure 2A). OTUs affiliated with Syntrophaceae (OTU 5) and Pedosphaeraceae (OTU 77), were also unique to the Siq G2 group. Additional G2 Siq OTUs included Hydrogenophilaceae (OTU 34) and Burkholderiaceae (OTU 9), other Gammaprotoebacteria (OTUs 31 and 71), Bacteroidetes (OTU 68), and unclassified Bacteria (OTU 97).
The SukB data set had G1 OTUs affiliated with Cyanobacteria in the families Xenococcaceae and Leptolynbyaceae (OTUs 40 and 119), Intrasporangiaceae in the Actinobacteria (OTU 12), in addition to an unclassified Deltaproteobacteria (OTU 10) and unclassified Bacteria (OTU 322) (Figure 2B). SukB G2 OTUs were associated with uncultured Actinobacteria (OTUs 24 and 98), Desulfobacteraceae (OTUs 7 and 237), Desulfobulbaceae (OTU 62), the family MidBa8 in the Myxococcales (OTU 222), and an unclassified Bacteria (OTU 267). An uncultured G2 Chloroflexi (OTU 43) affiliated with a new candidate class, was unique to this site. G3 OTUs were affiliated with Intrasporangiaceae in the Actinobacteria (OTU 12), Desulfobacteraceae (OTU 15), Desulfarculaceae (OTU 23), and unclassified Bacteria (OTU 8).
Group 1 SukS OTUs were related to Intrasporangiaceae (OTU 12), Lentimicrobiaceae (OTU 13), and unclassified Deltaproteobacteria (OTU 10) (Figure 4C). G2 OTUs included Actinobacteria in the families Microbacteriaceae and Cellulomonadaceae (OTUs 65 and 79), uncultured Armatimonadetes (OTU 187), Xenococcaceae (OTU 40), Desulfarculaceae (OTU 23), Burkholderiaceae (OTU 9), Hydrogenophilaceae (OTU 34), and unclassified Bacteria (OTUs 148 and 370). An OTU (57) related to Desulfobulbaceae was the sole OTU in G3.
Thermokarst Lake Sediment-Associated Archaeal Assemblage Composition
The first indication of high levels of archaeal gene sequence diversity was revealed by PCR-DGGE of archaeal V3 SSU rRNA gene fragments from samples collected in 2010 from discrete depths of the four lakes surveyed (Table 1; Figures S2a,b). Differences between Siq10 and Suk10 lake sediments (Figure S2c) and a high number of archaeal phylotypes per site (18–43) were observed by DGGE. Subsequent analyses based on SSU rRNA gene (DNA) and amplified SSU rRNA (aRNA) iTag sequences from Siq, SukB, and SukS 2011 samples indicated that archaeal OTUs comprised a diverse, but low proportion of the iTag sequences (~1% OTUs excluding singletons and doubletons). Yet, resampling of the data set reproducibly indicated their presence and suggested a distinct composition between the sediments from all three sites. Two-dimensional hierarchical clustering (average linkage) analysis of iTag OTUs based on non-parametric Spearman's rank correlations demonstrated the patterns clearly (Figure 3). Not only were iTag sequences different between sites, but their phylogenetic affiliations were often different; only OTUs related to the Bathyarchaeota and the order Thermoplasmatales in the Euryarchaeota were present at all three sites. OTUs affiliated with Euryarchaeota methanogens were dominant among the DNA and aRNA sequences at the Siq coring site, while OTUs affiliated with Woesearchaeota were only found in DNA and aRNA sequences at the SukB coring site.
Figure 3. Two-dimensional hierarchical clustering (average linkage) analysis of iTag OTUs based on non-parametric Spearman's rank correlations. Abundances were rank-transformed within each resampling. Average neighbor clustering of average rank correlation was used to order OTUs and samples. Only OTUs that were observed in >95% of re-sampled datasets are shown. Sampling sites are shown in color at the top of the graph along with dendogram showing clustering relationship between the sampled aRNA or DNA where SukS samples are in blue, SukB in green, and Siq in orange. The presence of different Archaeal OTUs are shown by colored bars, with those samples represented by aRNA shown with a gray background. Different archaeal phyla are indicated by different shades of blue to purple colors. Different groups of Euryarchaeota (All methanogen-related OTUs, Thermoplasmatales, and ANME OTUs were also indicated by different colors magenta, dark purple, and pink, respectively. Taxonomic affiliation was confirmed by phylogenetic analysis of their closest neighbors (mostly > 97% identity) in the Silva Database (123).
Noteworthy distinctions among sites included OTUs affiliated with Methanobacterium, Methanosarcina, Methanosaeta, and Methanocellales in Siq, while SukB only had a single methanogen-OTU and they were absent in the SukS data set. The SukB iTag composition differed substantially from the other sites: OTUs related to anaerobic methanotrophs of the ANME-2d group and other Euryarchaeota in the order Thermoplasmatales (AMOS1A-4113-D04) were observed, in addition to representatives of the Lokiarchaeota, Woesearchaeota, Thaumarchaeota and Bathyarchaeota. At SukS, Bathyarchaeota-, Thermoplasmatales-, and Thaumarchaeota-related OTUs were present. Additionally, three archaeal OTUs without close relatives (<90% identical in the Silva database) were present in SukB and SukS.
Phylogeny of Arctic Archaea Inferred from the SSU rRNA Gene
Phylogenetic analysis of archaeal SSU rRNA gene sequence region of overlap between the two primer sets (~662 bp) indicated that thermokarst lake sediment Archaea from 2010 and 2011 samples span the a wide range of Archaea phyla recognized to date (Figure 4). Only a few iTag-related sequences fell into clusters not covered by the clone sequences (e.g., Methanobacteriales in the Euryarchaeota, Lokiarchaeota, and Crenarchaeota). The majority of arctic thermokarst lake Archaea cloned sequences and iTag sequences were related to members of the phylum Bathyarchaeota (Meng et al., 2014; Evans et al., 2015), particularly those originating from Sukok Lake. Less well represented were iTag-related sequences within the Thaumarchaeota, a sister group to the Bathyarchaeota; and only one iTag-related sequence was affiliated with Crenarchaeota. Methanogen representatives, Methanosarcinales, Methanomicrobiales, Methanocellales, and a new order of methanogens related to the Thermoplasmatales (Borrel et al., 2013) were found among cloned and iTag-related sequences. Several Thermoplasmatales-related sequences in the Euryarchaeota were also observed in cloned and iTag sequences from Siq and SukB. Only one cloned sequence was related to anaerobic methanotrophs of the ANME-2d group, though several iTag-related sequences grouped in this part of the tree too. Sequences detected in both clone and iTag libraries affiliated with the Woesearchaeota and the Pacearchaeota clustered with numerous sequences from other arctic environments (i.e., Canadian artic lakes and rivers, and the Beaufort Sea). Sequences assigned to the candidate phylum Lokiarchaeota, which has been suggested to hold common ancestry with the Eukaryotes (Spang et al., 2015; Mariotti et al., 2016), were also detected.
Figure 4. Maximum likelihood phylogenetic tree of SSU rRNA gene (626–660 bp). Clone library sequences from both years (2010 and 2011) are coded by symbol as follows: Siq—circles, SukB—hexagons, SukS—squares, and Ikroavik (Ikr)—stars, followed by their corresponding accession numbers. Note that only one representative sequence was included for each taxon in the phylogenetic tree, and all the other sequences are available (KY886013–KY886132). Silva database sequences with ≥97% identity to iTag sequences are shown as red ovals and sequences with ≤97% identity are shown as orange ovals. Environmental sequences or culture isolates in the GenBank that were ≥97% similar to Sanger sequences, and that were obtained from cold environments (e.g., Canadian Arctic, Laptev Sea, Alaska, Antarctica, and others) are represented in pink and nearest neighboring sequences obtained in GenBank or the SILVA database are in purple. Only bootstrap values ≥50 are shown, and bootstrap values ≥80% are represented by black circles. For a full description of the SSU rRNA gene sequences on the tree, including sequences from the JGI IMG database, see Figure S4.
We set out to study the microbial assemblages that may be actively involved in geochemical processes relevant to thermokarst lake sediments. Siq has been reported to have a biogenic CH4 signature and higher levels of TOC and iron (Fe) compared to Sukok, which has a natural gas seep near SukS with a thermogenic CH4 signature, and higher levels of nitrate () and in pore waters at both sites (Matheus Carnevali et al., 2015). Between the two sites in Sukok, the geochemical differences appeared to be modest, though low levels of biogenic CH4 production potential was demonstrated at SukB in addition to higher Fe levels compared with SukS (Matheus Carnevali et al., 2015). Based on these data, and the remarkable presence of biogenic CH4 in an area with natural gas reserves (thermogenic CH4), we focused our attention on these contrasting coring sites. Thus, we hypothesized that the microbial assemblage composition, the potentially active fraction (indicated by high aRNA:DNA levels) in particular, differed following geochemistry.
Microbial Assemblage Structure and Composition
NMDS patterns based on DNA and aRNA suggested differences in microbial assemblage structure among the three sites that were consistent with geochemical profiling of sediment cores (Matheus Carnevali et al., 2015). The lack of statistically significant differences between lakes in the analysis of aRNA sequences may be explained by a high variation in assemblage structure between depths and sites at Sukok. Another potentially confounding influence enhancing distance between sediment samples could be inactive organisms reflected by the DNA-based iTag surveys (Juottonen et al., 2008), though the aRNA results here were largely in support of differences between coring sites.
At the Phylum or Class levels Siqlukaq and Sukok thermokarst lake microbial assemblages are representative of other arctic permafrost and peatland environments. Actinobacteria, Bacteroidetes, and Proteobacteria, which dominated here, also tend to dominate other permafrost environments (Yergeau et al., 2010; Frank-Fahle et al., 2014; Johnston et al., 2016), although Deltaproteobacteria and not Alphaproteobacteria or Betaproteobacteria were more abundant in this study. Nitrospirae were found to be abundant in different kinds of Northern wetlands (Wilhelm et al., 2011; Urbanová and Bárta, 2014). Other studies have also reported Chloroflexi and Cyanobacteria (Ganzert et al., 2014; Hultman et al., 2015). Verrucomicrobia were also found in a Northern wetland, where Armatimonadetes were considered “rare” members of the community (Serkebaeva et al., 2013). Bacteria and Archaea affiliated with the same phyla observed in this study have also been found in Russian permafrost (Shcherbakova et al., 2016), freshwater lakes of the Yunnan Plateau (Zhang et al., 2015), and most importantly in Toolik Lake on the North Slope of Alaska (Crump et al., 2012). Also, Acidobacteria were less abundant in these mineral lake sediments than in acidic peatlands (Wilhelm et al., 2011; Tveit et al., 2012; Serkebaeva et al., 2013).
Major Inferences of Potentially Active Bacteria and Putative Functions in Lake Sediments
The value of using both rRNA and rRNA coding DNA sequence information in this study increased confidence in the interpretation of the molecular ecological results. We observed correlations between RNA and DNA frequencies among different taxa (Campbell et al., 2011) and used rRNA: DNA ratios >1 as a proxy for activity following Wilhelm et al. (2014) and Hultman et al. (2015). These OTUs (Groups 1–3) were considered as candidates for potentially contributing important ecosystem functions and were identified as those with some form of activity. That is, microorganisms may increase the number of ribosomes not only in preparation for protein synthesis to meet real-time growth and production needs, but also as an adaptive strategy to cope with rapid shifts in environmental conditions, or before entering dormancy (Blazewicz et al., 2013). However, we expect that it is more likely to overestimate false negatives (i.e., classify organisms as dormant) than false positives (i.e., classify organisms as active) in communities comprising populations with variable growth strategies (Steven et al., 2017).
Diverse bacterial assemblages with taxa that may be actively involved in C, S, and N cycling were found at all three sediment coring sites. Differences in biogeochemical functional potential between the sites were revealed using rRNA-derived sequence analyses. Group 1 OTUs demonstrated this pattern, where N- or S-cycling Nitrospirae in Siqlukaq; photosynthetic, N2-fixing Cyanobacteria in SukB; organic matter decomposing Actinobacteria and Bacteroidetes in SukS; as well as Actinobacteria and unclassified Deltaproteobacteria (MBNT15) in both Sukok sites, suggested different processes at each site.
Functional inferences regarding methane oxidation in particular, appear to be unique to Siq. Here, the data suggest the prevalence (Figures S2a, S3b) of G2 and G3 OTUs affiliated with Methylobacter and Crenothrix. Some Crenothrix can be responsible for the oxidation of dissolved iron (Fe2+) (Emerson et al., 2010; Hedrich et al., 2011) while others are also capable of aerobic and anaerobic CH4 oxidation (Stoecker et al., 2006; Oswald et al., 2017). When these sequences were compared to the nr database at NCBI using BLAST, these two OTUs (21 and 71) matched Methylobacter tundripaludum (95–100% similarity, 98–100% coverage). The discrepancy in taxonomic affiliations are currently unresolved—regardless, both Crenothrix and Methylobacter are type I methanotrophs. M. tundripaludum is a psychrophile isolated from High Arctic wetland soil (Wartiainen et al., 2006; Svenning et al., 2011) and it has also been observed in Lake Qalluuraq, 60 miles to the South of Siq (one of our initial lakes surveyed for archaeal diversity; He et al., 2012a). The potential metabolisms of these microorganisms reflect the geochemical resources in Siq sediments, which contain the highest amounts of Fe2+ and CH4 of the three sites (Matheus Carnevali et al., 2015).
In a similar vein, though less clear of the functional role(s), Siq contained G1 Nitrospirae affiliated OTUs, not closely related to cultivated members of this group. Nitrospirae are metabolically versatile bacteria that include chemoautotrophic nitrite () oxidizers (Lücker et al., 2010), complete ammonia (NH3) oxidizers (Daims et al., 2015; Koch et al., 2015), and S- or -reducing bacteria (Sekiguchi et al., 2002, 2008).
OTUs that could participate in complex organic matter degradation processes at the coring sites included Actinobacteria and Bacteroidetes in addition to other less well represented, yet potentially active taxa. Bacteroidetes OTUs appeared to be active in SukS (G1) and Siq (G2), and are known to be associated with aerobic and anaerobic cellulose, chitin, and pectin degradation (Kirchman et al., 2014). G1-G3 Actinobacteria OTUs in SukS and SukB (Table 2) may also be contributing to complex C compound degradation (Barka et al., 2016). These sequence-based observations are consistent with analysis of organic matter composition in Suk, which indicated that in comparison to Siq, there was a greater contribution of free hydrocarbons derived from terrestrial plants in these sediments (Matheus Carnevali et al., 2015). At Siq, Verrucomicrobia G2 OTUs were also unique and could contribute to organic matter decomposition with saccharolytic metabolic potential (Bergmann et al., 2011; Hugenholtz et al., 2016).
Sulfur-cycling Deltaproteobacteria in Suk were affiliated with dissimilatory reducers, while Siq Deltaproteobacteria were affiliated with organisms that can be found in syntrophy with methanogens or other H2-consuming bacteria. These organisms could carry out subsequent steps in organic matter mineralization (e.g., fermentation and anaerobic respiration) in the lake sediments. The potential for reducing , sulfite (), and thiosulfate to hydrogen sulfide (H2S) could be inferred from the presence of groups known to perform these functions, with high aRNA to DNA ratios at both Sukok coring sites. Conversely, at the Siq coring site, Syntrophaceae (G2 OTUs) and Hydrogenophilaceae (G2 OTUs) were the OTUs that could potentially be attributed to S cycling. The Hydrogenophilaceae OTU was affiliated with Thiobacillus, chemoautotrophic bacteria capable of coupling Fe or S oxidation with denitrification (Beller et al., 2006).
Expansive Diversity of Archaea and Their Potential Roles in Arctic Lake Sediments
We investigated the previously unknown structure and phylogeny of Archaea in sediments of four arctic thermokarst lakes, and uncovered high archaeal diversity in these sediments. These findings emphasize the 2011 coring sites Siq, SukS, and SukB, even though Ikr10-O was also included in the phylogenetic analyses. Based on SSU rRNA gene clone library sequencing and SSU rRNA gene and SSU rRNA iTag sequencing, distinct Archaea compositions were observed at the three coring sites. These differences were also reflected in the potentially active Archaea (Figure 3). Overall, Archaeal OTUs comprised a low proportion of the total iTag sequences surveyed, which is comparable to other arctic soil and permafrost environments (Wilhelm et al., 2011; Tveit et al., 2012). Likewise, we found higher diversity and higher numbers of potentially active methanogen OTUs in sediments with higher organic levels, similar to findings in the water column of anoxic sub-arctic thaw ponds in Quebec, Canada (Crevecoeur et al., 2016).
OTUs affiliated with methanogens dominated at the Siq coring site, while OTUs affiliated with other archaeal groups were distributed in the Suk sediments. This partition of methanogens mirrored the potential for biological CH4 production observed in Siq and SukB sediments (Matheus Carnevali et al., 2015). The phylogenetic diversity of methanogens in Siq also suggested the potential for multiple pathways of methanogenesis. For example, methanogens in the orders Methanosarcinales and Methanomicrobiales are capable of acetoclastic or hydrogenotrophic methanogenesis, while members in the Methanocellales only produce CH4 from H2 and CO2. Methanosarcinales and Methanomicrobiales were not only found to be potentially active here, but also in a recent study in the anoxic water column of subarctic thaw ponds in Quebec (Crevecoeur et al., 2016). Methanogens related to Methanomassiliicoccus luminyensis are characterized by methylotrophic methanogenesis (Dridi et al., 2012), similar to Archaea in the phylum Bathyarchaeota, which are the first group identified outside the phylum Euryarchaeota that have the potential for CH4 production. The Bathyarchaeota, originally the Miscellaneous Crenarchaeal Group (MCG), is a diverse (76–87.7% intragroup similarity of the SSU rRNA gene) and ubiquitous group of Archaea (Kubo et al., 2012; Meng et al., 2014) that is metabolically versatile, and includes members capable of acetogenesis (He et al., 2016; Lazar et al., 2016).
Of all three sites, SukB had the greatest diversity of archaeal phyla and SukS the lowest. OTUs in SukB were affiliated with Euryarchaeota (e.g., ANME-2d), Thaumarchaeota (e.g., Candidatus Nitrososphaera), Woesearchaeota, and Lokiarchaeota, while OTUs in SukS were limited to Bathyarchaeota, Thaumarchaeota, and Thermoplasmatales (AMOS1A). The presence of OTUs affiliated with ANME-2d in SukB suggested the potential for anaerobic CH4 oxidation in the sediments, as observed in organisms (Candidatus Methanoperedens nitroreducens) that belong to this group (Haroon et al., 2013), which are able to couple the reduction of Fe3+ and Mn4+ to CH4 oxidation (Ettwig et al., 2016). Moreover, the presence of an OTU affiliated with Candidatus Nitrososphaera in SukB and other Thaumarchaeota in SukS, suggested that ammonium-driven microbial processes may occur (Spang et al., 2012; Zhalnina et al., 2014). These findings, in addition to the abundance of N2-fixing Cyanobacteria in Sukok allow us to infer that N cycling may be active in arctic thermokarst lake sediments.
Archaeal phylogenetic analysis also revealed extensive diversity across several poorly characterized and recently recognized archaeal phyla. Here, we found sequences affiliated with recently proposed phyla Woesearchaeota and Pacearchaeota (Baker et al., 2010; Rinke et al., 2013; Castelle et al., 2015), that are currently recognized to be dominant in environments such as oligotrophic alpine Pyrenean lakes (Ortiz-Alvarez and Casamayor, 2016). Analysis of a reconstructed, almost complete Woesearchaeota genome suggested that these organisms have a saccharolytic or fermentative life style, and the ability to metabolize hydrogen (Castelle et al., 2015). Additionally, some Pacearchaeota may have the ability to fix C (Castelle et al., 2015). Interestingly, these phyla may be common to high latitudes, as they were previously recognized as (i) the Rice Cluster V (RC-V), many of which have been reported from Arctic habitats (Großkopf et al., 1998; Galand et al., 2006, 2008; Høj et al., 2008) and (ii) a second geographically widespread cluster reported from an arctic river, a stamukhi lake, and the Arctic Ocean (Galand et al., 2006, 2008), in addition to lake sediments in Germany (Glissman et al., 2004). This group was also previously regarded as Candidate division WCHD3-30 (Galand et al., 2006; Høj et al., 2008).
In summary, our work revealed distinctive bacterial assemblages in the three sites studied, which are also implicated in different biogeochemical processes. Moreover, the high diversity of Archaea observed argues for further study to establish their roles in CH4 production and consumption, not to mention currently unrecognized roles of some pervasive new phyla now recognized in these sediments. Arctic thermokarst lakes are sensitive to climate warming (Osterkamp, 2005; McGuire et al., 2009; Thienpont et al., 2013), and we predict that deepening of the lake talik and shifts in geomorphology might influence these lake ecosystems according to their geochemical properties. Understanding the geochemistry is only part of the picture that underlies the differences in these lakes. To fully understand how small shifts in temperature, seasonal duration, and snowfall, may cascade into shifts in ecosystem function, the regulation of microbial activities is key to predict greenhouse gas emissions.
PM, CH, KH, JP, and AM contributed to the design of the study, analysis of the data, drafting of the manuscript, and approval of the content. All agree to be accountable for all aspects of the work.
Financial support for this research was provided by the NASA Astrobiology Institute, Astrobiology of Icy Worlds program at JPL. iTag sequencing was supported by an award to AE Murray, JGI-634 by the DOE Joint Genome Institute. PM was supported in part by the Division of Earth and Ecosystem Sciences, DRI. KH acknowledges support through the Jet Propulsion Laboratory (JPL), California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Special thanks to H. Adams, D. Berisford, A. Klesh, J. Leichty, M. Rorhssen, P. Santibañez, and the Barrow Arctic Science Consortium for assistance in the field. We are very grateful to E. Kuhn, V. Peng, O. Rassuchine, and G. Trubl for technical assistance in the lab; and to J. Mejia for computational contributions to data processing. An earlier version of this work originally appeared in P. Matheus Carnevali's dissertation (2015). C. Brown, C. Castelle, and L. Hug are kindly acknowledged for their contributions to the phylogenetic analysis, and we are also grateful to L. Hug and K. Anantharanam for useful comments on the manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01192/full#supplementary-material
Baker, B. J., Comolli, L. R., Dick, G. J., Hauser, L. J., Hyatt, D., Dill, B. D., et al. (2010). Enigmatic, ultrasmall, uncultivated Archaea. Proc. Natl. Acad. Sci. U.S.A. 107, 8806–8811. doi: 10.1073/pnas.0914470107
Barka, E. A., Vatsa, P., Sanchez, L., Gaveau-Vaillant, N., Jacquard, C., Klenk, H. P., et al. (2016). Taxonomy, physiology, and natural products of Actinobacteria. Microbiol. Mol. Biol. Rev. 80, 1–43. doi: 10.1128/MMBR.00019-15
Beller, H. R., Chain, P. S., Letain, T. E., Chakicherla, A., Larimer, F. W., Richardson, P. M., et al. (2006). The genome sequence of the obligately chemolithoautotrophic, facultatively anaerobic bacterium Thiobacillus denitrificans. J. Bacteriol. 188, 1473–1488. doi: 10.1128/JB.188.4.1473-1488.2006
Bergmann, G. T., Bates, S. T., Eilers, K. G., Lauber, C. L., Caporaso, J. G., Walters, W. A., et al. (2011). The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol. Biochem. 43, 1450–1455. doi: 10.1016/j.soilbio.2011.03.012
Blazewicz, S. J., Barnard, R. L., Daly, R. A., and Firestone, M. K. (2013). Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses. ISME J. 7, 2061–2068. doi: 10.1038/ismej.2013.102
Borrel, G., O'Toole, P. W., Harris, H. M. B., Peyret, P., Brugère, J.-F., and Gribaldo, S. (2013). Phylogenomic data support a seventh order of methylotrophic methanogens and provide insights into the evolution of methanogenesis. Genome Biol. Evol. 5, 1769–1780. doi: 10.1093/gbe/evt128
Campbell, B. J., Yu, L., Heidelberg, J. F., and Kirchman, D. L. (2011). Activity of abundant and rare bacteria in a coastal ocean. Proc. Natl. Acad. Sci. U.S.A. 108, 12776–12781. doi: 10.1073/pnas.1101405108
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Castelle, C. J., Wrighton, K. C., Thomas, B. C., Hug, L. A., Brown, C. T., Wilkins, M. J., et al. (2015). Genomic expansion of domain archaea highlights roles for organisms from new phyla in anaerobic carbon cycling. Curr. Biol. 25, 690–701. doi: 10.1016/j.cub.2015.01.014
Crump, B. C., Amaral-Zettler, L. A., and Kling, G. W. (2012). Microbial diversity in arctic freshwaters is structured by inoculation of microbes from soils. ISME J. 6, 1629–1639. doi: 10.1038/ismej.2012.9
Dridi, B., Fardeau, M.-L., Ollivier, B., Raoult, D., and Drancourt, M. (2012). Methanomassiliicoccus luminyensis gen. nov., sp. nov., a methanogenic archaeon isolated from human faeces. Int. J. Syst. Evol. Microbiol. 62, 1902–1907. doi: 10.1099/ijs.0.033712-0
Elvert, M., Pohlman, J. W., Becker, K. W., Gaglioti, B., Hinrichs, K. U., and Wooller, M. J. (2016). Methane turnover and environmental change from Holocene lipid biomarker records in a thermokarst lake in Arctic Alaska. Holocene 26, 1766–1777. doi: 10.1177/0959683616645942
Ettwig, K. F., Zhu, B., Speth, D., Keltjens, J. T., Jetten, M. S., and Kartal, B. (2016). Archaea catalyze iron-dependent anaerobic oxidation of methane. Proc. Natl. Acad. Sci. U.S.A. 113, 12792–12796. doi: 10.1073/pnas.1609534113
Evans, P. N., Parks, D. H., Chadwick, G. L., Robbins, S. J., Orphan, V. J., Golding, S. D., et al. (2015). Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science 350, 434–438. doi: 10.1126/science.aac7745
Frank-Fahle, B. A., Yergeau, E., Greer, C. W., Lantuit, H., and Wagner, D. (2014). Microbial functional potential and community composition in permafrost-affected soils of the NW Canadian Arctic. PLoS ONE 9:e84761. doi: 10.1371/journal.pone.0084761
Galand, P. E., Lovejoy, C., and Vincent, W. F. (2006). Remarkably diverse and contrasting archaeal communities in a large arctic river and the coastal Arctic Ocean. Aquat. Microb. Ecol. 44, 115–126. doi: 10.3354/ame044115
Galand, P. E., Lovejoy, C., Pouliot, J., Garneau, M.-È., and Vincent, W. F. (2008). Microbial community diversity and heterotrophic production in a coastal Arctic ecosystem: a stamukhi lake and its source waters. Limnol. Oceanogr. 53, 813–823. doi: 10.4319/lo.2008.53.2.0813
Ganzert, L., Bajerski, F., and Wagner, D. (2014). Bacterial community composition and diversity of five different permafrost-affected soils of Northeast Greenland. FEMS Microbiol. Ecol. 89, 426–441. doi: 10.1111/1574-6941.12352
Glissman, K., Chin, K.-J., Casper, P., and Conrad, R. (2004). Methanogenic pathway and archaeal community structure in the sediment of eutrophic Lake Dagow: effect of temperature. Microb. Ecol. 48, 389–399. doi: 10.1007/s00248-003-2027-2
Graham, D. E., Wallenstein, M. D., Vishnivetskaya, T. A., Waldrop, M. P., Phelps, T. J., Pfiffner, S. M., et al. (2012). Microbes in thawing permafrost: the unknown variable in the climate change equation. ISME J. 6, 709–712. doi: 10.1038/ismej.2011.163
Großkopf, R., Stubner, S., and Liesack, W. (1998). Novel euryarchaeotal lineages detected on rice roots and in the anoxic bulk soil of flooded rice microcosms. Appl. Environ. Microbiol. 64, 4983–4989.
Haroon, M. F., Hu, S., Shi, Y., Imelfort, M., Keller, J., Hugenholtz, P., et al. (2013). Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage. Nature 500, 567–570. doi: 10.1038/nature12375
He, R., Wooller, M. J., Pohlman, J. W., Catranis, C., Quensen, J., Tiedje, J. M., et al. (2012b). Identification of functionally active aerobic methanotrophs in sediments from an arctic lake using stable isotope probing. Environ. Microbiol. 14, 1403–1419. doi: 10.1111/j.1462-2920.2012.02725.x
He, R., Wooller, M. J., Pohlman, J. W., Quensen, J., Tiedje, J. M., and Leigh, M. B. (2012a). Diversity of active aerobic methanotrophs along depth profiles of arctic and subarctic lake water column and sediments. ISME J. 6, 1937–1948. doi: 10.1038/ismej.2012.34
He, Y., Li, M., Perumal, V., Feng, X., Fang, J., Xie, J., et al. (2016). Genomic and enzymatic evidence for acetogenesis among multiple lineages of the archaeal phylum Bathyarchaeota widespread in marine sediments. Nat. Microbiol. 1:16035. doi: 10.1038/nmicrobiol.2016.35
Herbold, C. W., Lee, C. K., McDonald, I. R., and Cary, S. C. (2014). Evidence of global-scale aeolian dispersal and endemism in isolated geothermal microbial communities of Antarctica. Nat. Commun. 5:3875. doi: 10.1038/ncomms4875
Heslop, J. K., Walter Anthony, K. M., Sepulveda-Jauregui, A., Martinez-Cruz, K., Bondurant, A., Grosse, G., et al. (2015). Thermokarst lake methanogenesis along a complete talik profile. Biogeosciences 12, 4317–4331. doi: 10.5194/bg-12-4317-2015
Høj, L., Olsen, R. A., and Torsvik, V. L. (2008). Effects of temperature on the diversity and community structure of known methanogenic groups and other archaea in high Arctic peat. ISME J. 2, 37–48. doi: 10.1038/ismej.2007.84
Hultman, J., Waldrop, M. P., MacKelprang, R., David, M. M., McFarland, J., Blazewicz, S. J., et al. (2015). Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes. Nature 521, 208–212. doi: 10.1038/nature14238
Johnston, E. R., Rodriguez, R. L., Luo, C., Yuan, M. M., Wu, L., He, Z., et al. (2016). Metagenomics reveals pervasive bacterial populations and reduced community diversity across the Alaska tundra ecosystem. Front. Microbiol. 7:579. doi: 10.3389/fmicb.2016.00579.
Juottonen, H., Tuittila, E. S., Juutinen, S., Fritze, H., and Yrjälä, K. (2008). Seasonality of rDNA- and rRNA-derived archaeal communities and methanogenic potential in a boreal mire. ISME J. 2, 1157–1168. doi: 10.1038/ismej.2008.66
Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649. doi: 10.1093/bioinformatics/bts199
Kirchman, D. L., Hanson, T. E., Cottrell, M. T., and Hamdan, L. J. (2014). Metagenomic analysis of organic matter degradation in methane-rich Arctic Ocean sediments. Limnol. Oceanogr. 59, 548–559. doi: 10.4319/lo.2014.59.2.0548
Koch, H., Lücker, S., Albertsen, M., Kitzinger, K., Herbold, C., Spieck, E., et al. (2015). Expanded metabolic versatility of ubiquitous nitrite-oxidizing bacteria from the genus Nitrospira. Proc. Natl. Acad. Sci. U.S.A. 112, 11371–11376. doi: 10.1073/pnas.1506533112
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13
Kubo, K., Lloyd, K. G. F., Biddle, J., Amann, R., Teske, A., and Knittel, K. (2012). Archaea of the miscellaneous crenarchaeotal group are abundant, diverse and widespread in marine sediments. ISME J. 6, 1949–1965. doi: 10.1038/ismej.2012.37
Lazar, C. S., Baker, B. J., Seitz, K., Hyde, A. S., Dick, G. J., Hinrichs, K. U., et al. (2016). Genomic evidence for distinct carbon substrate preferences and ecological niches of Bathyarchaeota in estuarine sediments. Environ. Microbiol. 18, 1200–1211. doi: 10.1111/1462-2920.13142
Lipson, D. A., Haggerty, J. M., Srinivas, A., Raab, T. K., Sathe, S., and Dinsdale, E. A. (2013). Metagenomic insights into anaerobic metabolism along an Arctic peat soil profile. PLoS ONE 8:e64659. doi: 10.1371/journal.pone.0064659
Lofton, D. D., Whalen, S. C., and Hershey, A. E. (2014). Effect of temperature on methane dynamics and evaluation of methane oxidation kinetics in shallow arctic Alaskan lakes. Hydrobiologia 721, 209–222. doi: 10.1007/s10750-013-1663-x
Lücker, S., Wagner, M., Maixner, F., Pelletier, E., Koch, H., Vacherie, B., et al. (2010). A Nitrospira metagenome illuminates the physiology and evolution of globally important nitrite-oxidizing bacteria. Proc. Natl. Acad. Sci. U.S.A. 107, 13479–13484. doi: 10.1073/pnas.1003860107
Mariotti, M., Lobanov, A. V., Manta, B., Santesmasses, D., Bofill, A., Guigo, R., et al. (2016). Lokiarchaeota marks the transition between the archaeal and eukaryotic selenocysteine encoding systems. Mol. Biol. Evol. 33, 2441–2453. doi: 10.1093/molbev/msw122
Massana, R., Murray, A. E., Preston, C. M., and DeLong, E. F. (1997). Vertical distribution and phylogenetic characterization of marine planktonic Archaea in the Santa Barbara Channel. Appl. Environ. Microbiol. 63, 50–56.
Matheus Carnevali, P. B. (2015). Microbial Ecology, Biogeochemistry, and Signatures of Life at Low Temperature in Arctic Thermokarst Lake Sediments and High Sierra Snow Fields. Ph.D. Dissertation, University of Nevada, Reno.
Matheus Carnevali, P. B., Rohrssen, M., Williams, M. R., Michaud, A. B., Adams, H., Berisford, D., et al. (2015). Methane sources in arctic thermokarst lake sediments on the North Slope of Alaska. Geobiology 13, 181–197. doi: 10.1111/gbi.12124
McGuire, A. D., Anderson, L. G., Christensen, T. R., Dallimore, S., Guo, L., Hayes, D. J., et al. (2009). Sensitivity of the carbon cycle in the Arctic to climate change. Ecol. Monogr. 79, 523–555. doi: 10.1890/08-2025.1
Meng, J., Xu, J., Qin, D., He, Y., Xiao, X., and Wang, F. (2014). Genetic and functional properties of uncultivated MCG archaea assessed by metagenome and gene expression analyses. ISME J. 8, 650–659. doi: 10.1038/ismej.2013.174
Negandhi, K., Laurion, I., Whiticar, M. J., Galand, P. E., Xu, X., and Lovejoy, C. (2013). Small thaw ponds: an unaccounted source of methane in the Canadian High Arctic. PLoS ONE 8:e78204. doi: 10.1371/journal.pone.0078204
Ortiz-Alvarez, R., and Casamayor, E. O. (2016). High occurrence of Pacearchaeota and Woesearchaeota (Archaea superphylum DPANN) in the surface waters of oligotrophic high-altitude lakes. Environ. Microbiol. Rep. 8, 210–217. doi: 10.1111/1758-2229.12370
Pruesse, E., Peplies, J., and Glöckner, F. O. (2012). SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 28, 1823–1829. doi: 10.1093/bioinformatics/bts252
Rinke, C., Schwientek, P., Sczyrba, A., Ivanova, N. N., Anderson, I. J., Cheng, J.-F., et al. (2013). Insights into the phylogeny and coding potential of microbial dark matter. Nature 499, 431–437. doi: 10.1038/nature12352
Sekiguchi, Y., Kamagata, Y., Ohashi, A., and Harada, H. (2002). Molecular and conventional analyses of microbial diversity in mesophilic and thermophilic upflow anaerobic sludge blanket granular sludges. Water Sci. Technol. 45, 19–25. doi: 10.2166/wst.2002.0279
Sekiguchi, Y., Muramatsu, M., Imachi, H., Narihiro, T., Ohashi, A., Harada, H., et al. (2008). Thermodesulfovibrio aggregans sp. nov. and Thermodesulfovibrio thiophilus sp. nov., anaerobic, thermophilic, sulfate-reducing bacteria isolated from thermophilic methanogenic sludge, and emended description of the genus Thermodesulfovibrio. Int. J. Syst. Evol. Microbiol. 58, 2541–2548. doi: 10.1099/ijs.0.2008/000893-0
Serkebaeva, Y. M., Kim, Y., Liesack, W., and Dedysh, S. N. (2013). Pyrosequencing-based assessment of the bacteria diversity in surface and subsurface peat layers of a northern wetland, with focus on poorly studied phyla and candidate divisions. PLoS ONE 8:e63994. doi: 10.1371/journal.pone.0063994
Shcherbakova, V., Yoshimura, Y., Ryzhmanova, Y., Taguchi, Y., Segawa, T., Oshurkova, V., et al. (2016). Archaeal communities of Arctic methane-containing permafrost. FEMS Microbiol. Ecol. 92:fiw135. doi: 10.1093/femsec/fiw135
Spang, A., Poehlein, A., Offre, P., Zumbrägel, S., Haider, S., Rychlik, N., et al. (2012). The genome of the ammonia-oxidizing Candidatus Nitrososphaera gargensis: insights into metabolic versatility and environmental adaptations. Environ. Microbiol. 14, 3122–3145. doi: 10.1111/j.1462-2920.2012.02893.x
Spang, A., Saw, J. H., Jørgensen, S. L., Zaremba-Niedzwiedzka, K., Martijn, J., Lind, A. E., et al. (2015). Complex archaea that bridge the gap between prokaryotes and eukaryotes. Nature 521, 173–179. doi: 10.1038/nature14447
Stahl, D. A., and Amann, R. (1991). “Development and application of nucleic acid probes,” in Nucleic Acid Techniques in Bacterial Systematics, eds E. Stackebrandt and M. Goodfello (Chichester: John Wiley & Sons), 205–248.
Steven, B., Hesse, C., Soghigian, J., and Dunbar, J. (2017). Simulated rRNA/DNA ratios show potential to misclassify active populations as dormant. Appl. Environ. Microb. 83, e00696–e00617. doi: 10.1128/AEM.00696-17
Stoecker, K., Bendinger, B., Schöning, B., Nielsen, P. H., Nielsen, J. L., Baranyi, C., et al. (2006). Cohn's Crenothrix is a filamentous methane oxidizer with an unusual methane monooxygenase. Proc. Natl. Acad. Sci. U.S.A. 103, 2363–2367. doi: 10.1073/pnas.0506361103
Stoeva, M. K., Aris-Brosou, S., Chételat, J., Hintelmann, H., Pelletier, P., and Poulain, A. J. (2014). Microbial community structure in lake and wetland sediments from a high Arctic polar desert revealed by targeted transcriptomics. PLoS ONE 9:e89531. doi: 10.1371/journal.pone.0089531
Svenning, M. M., Hestnes, A. G., Wartiainen, I., Stein, L. Y., Klotz, M. G., Kalyuzhnaya, M. G., et al. (2011). Genome sequence of the Arctic methanotroph Methylobacter tundripaludum SV96. J. Bacteriol. 193, 6418–6419. doi: 10.1128/JB.05380-11
Thienpont, J. R., RÜHland, K. M., Pisaric, M. F. J., Kokelj, S. V., Kimpe, L. E., Blais, J. M., et al. (2013). Biological responses to permafrost thaw slumping in Canadian Arctic lakes. Freshw. Biol. 58, 337–353. doi: 10.1111/fwb.12061
Tveit, A., Schwacke, R., Svenning, M. M., and Urich, T. (2012). Organic carbon transformations in High-Arctic peat soils: key functions and microorganisms. ISME J. 7, 299–311. doi: 10.1038/ismej.2012.99
Urbanová, Z., and Bárta, J. (2014). Microbial community composition and in silico predicted metabolic potential reflect biogeochemical gradients between distinct peatland types. FEMS Microbiol. Ecol. 90, 633–646. doi: 10.1111/1574-6941.12422
Wang, J., Sheng, Y., Hinkel, K. M., and Lyons, E. A. (2012). Drained thaw lake basin recovery on the western Arctic Coastal Plain of Alaska using high-resolution digital elevation models and remote sensing imagery. Remote Sens. Environ. 119, 325–336. doi: 10.1016/j.rse.2011.10.027
Wartiainen, I., Hestnes, A. G., McDonald, I. R., and Svenning, M. M. (2006). Methylobacter tundripaludum sp. nov., a methane-oxidizing bacterium from Arctic wetland soil on the Svalbard islands, Norway (78 degrees N). Int. J. Syst. Evol. Microbiol. 56, 109–113. doi: 10.1099/ijs.0.63728-0
Wilhelm, L., Besemer, K., Fasching, C., Urich, T., Singer, G. A., Quince, C., et al. (2014). Rare but active taxa contribute to community dynamics of benthic biofilms in glacier-fed streams. Environ. Microbiol. 16, 2514–2524. doi: 10.1111/1462-2920.12392
Wilhelm, R. C., Niederberger, T. D., Greer, C., and Whyte, L. G. (2011). Microbial diversity of active layer and permafrost in an acidic wetland from the Canadian High Arctic. Can. J. Microbiol. 57, 303–315. doi: 10.1139/w11-004
Yergeau, E., Hogues, H., Whyte, L. G., and Greer, C. W. (2010). The functional potential of high Arctic permafrost revealed by metagenomic sequencing, qPCR and microarray analyses. ISME J. 4, 1206–1214. doi: 10.1038/ismej.2010.41
Zhalnina, K. V., Dias, R., Leonard, M. T., Dorr de Quadros, P., Camargo, F. A. O., Drew, J. C., et al. (2014). Genome sequence of Candidatus Nitrososphaera evergladensis from Group I 1b enriched from everglades soil reveals novel genomic features of the ammonia-oxidizing archaea. PLoS ONE 9:e101648. doi: 10.1371/journal.pone.0101648
Keywords: methane, thermokarst, lake, sediments, Arctic, Alaska, archaea, diversity
Citation: Matheus Carnevali PB, Herbold CW, Hand KP, Priscu JC and Murray AE (2018) Distinct Microbial Assemblage Structure and Archaeal Diversity in Sediments of Arctic Thermokarst Lakes Differing in Methane Sources. Front. Microbiol. 9:1192. doi: 10.3389/fmicb.2018.01192
Received: 26 October 2017; Accepted: 16 May 2018;
Published: 07 June 2018.
Edited by:Mark Alexander Lever, ETH Zürich, Switzerland
Reviewed by:Trinity L. Hamilton, University of Minnesota Twin Cities, United States
Nils Risgaard-Petersen, Aarhus University, Denmark
Copyright © 2018 Matheus Carnevali, Herbold, Hand, Priscu and Murray. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alison E. Murray, firstname.lastname@example.org
†Present Address: Paula B. Matheus Carnevali, Department Earth and Planetary Science, University of California, Berkeley, Berkeley, CA, United States