Kristina D. Baker
Colleen T. E. Kellogg
James W. McClelland
Kenneth H. Dunton
Byron C. Crump- 1Department of Microbiology, Oregon State University, Corvallis, OR, United States
- 2Hakai Institute, Heriot Bay, BC, Canada
- 3The University of Texas at Austin Marine Science Institute, Port Aransas, TX, United States
- 4College of Earth, Ocean, and Atmospheric Sciences, Oregon State University, Corvallis, OR, United States
In contrast to temperate systems, Arctic lagoons that span the Alaska Beaufort Sea coast face extreme seasonality. Nine months of ice cover up to ∼1.7 m thick is followed by a spring thaw that introduces an enormous pulse of freshwater, nutrients, and organic matter into these lagoons over a relatively brief 2–3 week period. Prokaryotic communities link these subsidies to lagoon food webs through nutrient uptake, heterotrophic production, and other biogeochemical processes, but little is known about how the genomic capabilities of these communities respond to seasonal variability. Replicate water samples from two lagoons and one coastal site near Kaktovik, AK were collected in April (full ice cover), June (ice break up), and August (open water) to represent winter, spring, and summer, respectively. Samples were size fractionated to distinguish free-living and particle-attached microbial communities. Multivariate analysis of metagenomes indicated that seasonal variability in gene abundances was greater than variability between size fractions and sites, and that June differed significantly from the other months. Spring (June) gene abundances reflected the high input of watershed-sourced nutrients and organic matter via spring thaw, featuring indicator genes for denitrification possibly linked to greater organic carbon availability, and genes for processing phytoplankton-derived organic matter associated with spring blooms. Summer featured fewer indicator genes, but had increased abundances of anoxygenic photosynthesis genes, possibly associated with elevated light availability. Winter (April) gene abundances suggested low energy inputs and autotrophic bacterial metabolism, featuring indicator genes for chemoautotrophic carbon fixation, methane metabolism, and nitrification. Winter indicator genes for nitrification belonged to Thaumarchaeota and Nitrosomonadales, suggesting these organisms play an important role in oxidizing ammonium during the under-ice period. This study shows that high latitude estuarine microbial assemblages shift metabolic capabilities as they change phylogenetic composition between these extreme seasons, providing evidence that these communities may be resilient to large hydrological events in a rapidly changing Arctic.
Introduction
The northeastern coast of Alaska is lined with a series of shallow lagoons that are protected from the Beaufort Sea by barrier islands. Of the 1,957 km of coastline along the Alaskan Beaufort Sea coast, 546 km are classified as lagoons with barrier islands (Jorgenson and Brown, 2005). These lagoons are home to a wide array of wildlife such as anadromous and marine fish, and birds from six continents that nest in the summer. They contain a diverse benthic community including bivalves, isopods, gastropods, and mysids. The lagoon food webs possess a high degree of omnivory and rely on terrestrial carbon as an energy source (Dunton et al., 2006; Harris et al., 2018).
The lagoons experience extreme shifts in seasonal environmental conditions, such as nine months of ice cover, water column temperatures ranging from −2 to 14°C, and salinity ranging from 0 to > 45 PSU (Harris et al., 2017). The spring freshet occurs as snowmelt begins and rivers thaw, resulting in over half of the annual freshwater river discharge occurring within a two-week period, usually in late May to early June (McClelland et al., 2014). This freshening not only changes the salinity of the lagoons but introduces terrestrially derived organic matter in dissolved and particulate forms. Terrestrial particulate organic matter (POM) is an important carbon source for the lagoons throughout the year, however, the composition of POM changes seasonally from relatively labile terrestrial material in the spring to more processed, refractory POM in winter (Connelly et al., 2015). Generally, dissolved organic carbon (DOC) is the most abundant form of terrestrial carbon transported by Arctic rivers (Gordeev et al., 1996; Lobbes et al., 2000; Holmes et al., 2012). In spring, this DOC is leached by snowmelt from frozen surface soils and plant litter and is more labile than DOC exported later in the summer (Holmes et al., 2008), likely because low temperature and rapid transport of this material off the frozen land surface limits terrestrial microbial processing. In summer, river DOC is less labile because warmer waters and longer travel times within river networks support greater microbial processing and photodegradation before this material makes it to the ocean (Holmes et al., 2008; Cory et al., 2014).
Previous work studying the food webs in these lagoons shows that organic matter in sediments and the water column have strong terrestrial signatures (Dunton et al., 2012; Harris et al., 2018). Terrestrially derived organic carbon is an important energy subsidy to members of lagoon food webs, especially some deposit feeders and omnivorous fishes who derive > 40% of their diet from terrestrial carbon (Harris et al., 2018). This subsidy is very high because the most abundant prey organisms in these systems, omnivorous detritivores, strongly rely on microbially processed terrestrial carbon (Dunton et al., 2012). Previous work also suggests that sediment microbial biogeochemical processes play a significant role in metazoan consumers assimilating terrestrial carbon (McTigue and Dunton, 2014). An absence of taxa in trophic positions between carbon sources and low-level consumers in the eastern Beaufort Sea suggest that microbes occupy this gap by processing terrestrial organic matter used by higher trophic levels (Bell et al., 2016). Thus, microbes may be a key component linking lagoon food webs to their major source of organic carbon.
Several studies have characterized Arctic microbial communities in freshwater systems and offshore environments (e.g., Garneau et al., 2006; Alonso-Sáez et al., 2008; Vallières et al., 2008; Galand et al., 2009; Underwood et al., 2019), but few studies have focused on the microbial communities of the nearshore Beaufort Sea where inputs of terrestrial material are very high. Changes in microbial community composition have been linked to seasonal variations in physical conditions and water chemistry in the coastal Beaufort Sea based on 16S and 18S rRNA gene data (Kellogg et al., 2019). The microbial community in the under-ice conditions in April were composed of coastal marine and chemoautotrophic taxa that potentially use methane, nitrogen, iron, and sulfur for energy during times when inputs of new/fresh organic matter are low. June communities were composed of freshwater and estuarine taxa, while August communities were composed of oligotrophic marine bacterial taxa (Kellogg et al., 2019).
Changes in microbial communities over seasons may in part be due to seasonal changes in the quantity and composition of dissolved organic matter (DOM). Sipler et al. (2017) found that additions of runoff from thawed permafrost to coastal water of the Chukchi Sea shifted the microbial community from oligotrophic to copiotrophic taxa. Approximately 7% of DOC and 38% of DON in these incubations were determined to be bioavailable within a 4–6 days time frame, showing that organic nutrients are labile and important sources to microbial communities (Sipler et al., 2017). Similarly, the molecular composition of DOM from melting first-year sea ice in the Canadian Arctic Archipelago has been shown to affect microbial community composition and abundance (Underwood et al., 2019). Microbial communities of varying size fractions, either particle-attached or free-living, can respond to changing environmental conditions differently (DeLong et al., 1993; Crump et al., 1999; Ghiglione et al., 2007; Garneau et al., 2009). Particle-attached communities from the Mackenzie River plume had higher activity in spring/summer than in fall/winter, and differences in particle-attached and free-living community composition and production were greatest during times of high POM input from the river (Garneau et al., 2009).
It is important to understand these microbial communities in the context of a changing Arctic climate. Longer periods of open water lead to longer growing seasons for phytoplankton, increasing the annual net primary production for these systems (Arrigo and van Dijken, 2015). Increased runoff of terrestrial organic matter in a warming Arctic may impact growth rates of heterotrophic bacteria, whose growth may be more significantly impacted by DOM than temperature (Peterson, 2002; McClelland et al., 2006; Déry et al., 2009; Kirchman et al., 2009; Rawlins et al., 2010; Romanovsky et al., 2010; Smith and Riseborough, 2010; Barnhart et al., 2014). Therefore, a clear gap in knowledge exists in understanding how marine microbial communities link terrestrial organic matter to the rest of the food web and how these linkages might be altered by a changing environment.
This study characterizes the metabolic capabilities of microbial communities across three seasons in two lagoons and one open coastal site along the eastern Alaska Beaufort Sea coast. We used metagenomic DNA sequence data of bacterial and archaeal water column communities to identify genes of relevant biogeochemical pathways. Results suggest that spring (June) and winter (April) months had the most strongly contrasting patterns in gene abundance. Spring had an increased abundance of genes involved in denitrification and processing of labile carbon from phytoplankton blooms, while winter had an increased abundance of genes involved in nitrification and chemoautotrophic carbon fixation. Summer (August) had few genes of elevated abundance relative to the other two seasons but featured the greatest abundance of genes for anoxygenic photosynthesis. These results suggest the potential for seasonal shifts in microbial community functions associated with different sources of nutrients and organic matter.
Materials and Methods
Study Site
Three sites located near the village of Kaktovik, Alaska were sampled during 2012. Kaktovik is an Iñupiat village located on Barter Island, along the northern-most part of the Arctic National Wildlife Refuge. The sample areas include two lagoons and one open coastal site (Figure 1). Jago Lagoon (JA; 70.106867, −143.443367) receives direct inputs from the Jago River and has direct connections to the Beaufort Sea through two wide inlets. Kaktovik Lagoon (KA; 70.089050, −143.623317) is a comparably smaller, more isolated system with no river input and has limited exchange with the Beaufort Sea owing to one narrow channel. The last site, Bernard Point (BP; 70.135617, −143.637367), is located along the northern shoreline of Barter Island in the Beaufort Sea.
Figure 1. Map of sampling sites located near Kaktovik, Alaska. Three study sites are Bernard Point (BP), Jago Lagoon (JA), and Kaktovik Lagoon (KA).
Sample Collection
Water samples were collected in April (full ice cover), June (ice break-up), and August (open water) 2012 at several locations within each lagoon and outside barrier islands (Figure 1). In April and June, a peristaltic pump was used to collect water 2 m below the surface or top of the ice. August samples were collected by submerging bottles by hand 0.5 m below the surface. Physical and biogeochemical parameters measured included particulate analysis (particulate organic carbon and nitrogen, chlorophyll a) as described in Connelly et al. (2015) and dissolved parameters (dissolved organic carbon and nitrogen, dissolved inorganic nitrogen) as described by the Beaufort Lagoon Ecosystems LTER and Core Program (2020a, b). Salinity, temperature, pH, and dissolved oxygen were measured using a sonde (YSI Corporation). Prokaryotic cell counts were measured by Kellogg et al. (2019). These environmental measurements are also discussed in Kellogg et al. (2019).
Within 6 h of collection, 4–11 L of water was first filtered through a 142 mm diameter, 3 micron pore size polycarbonate filter to collect particle-associated bacteria and then a 142 mm diameter, 0.22 micron pore size Supor filter (Pall) separated by a filter screen. Filters were placed into 15 mL falcon tubes and preserved with 8 mL of DNA extraction buffer (100 mM Tris, 100 mM NaEDTA, 100 mM phosphate buffer, 1.5 M NaCl, 1% CTAB) and frozen until extraction. Two replicates for each sample were collected from each site for a total of 36 samples (3 months, three sites, two size fractions, and replicates of each sample). One replicate BP June free-living sample was not processed due to contamination during transport.
Environmental Sampling and Processing
Known quantities of Thermus thermophilus HB-8 DNA (purchased from ATCC) were added to samples prior to DNA extraction (10.0 ng T. thermophilus DNA per L of sample filtered) following Satinsky et al. (2013) as an internal control to account for DNA lost during extraction and sequencing. Filters then underwent phenol: chloroform extraction following Crump et al. (2003) but scaled up for larger starting volumes as described in Doherty et al. (2017). Library preparation and sequencing were carried out by the Center for Genomic Research and Biocomputing (CGRB) at Oregon State University, Corvallis, OR. Multiplexed, 150 basepair paired-end sequencing was performed on an Illumina HiSeq 2000.
Metagenome Sequencing and Analysis
A detailed description of metagenome sequence analysis is provided in Supplementary Materials. Briefly, paired-end sequence reads were quality trimmed (BBDuk v.38.84, Bushnell, 2014), and internal T. thermophilus standard sequences were identified in a two-step process using BBSplit (BBTools; Bushnell, 2014) with a reference database of the T. thermophilus HB-8 genome, and with Kaiju (Menzel et al., 2016) with a reference database built using RefSeq (O’Leary et al., 2016). T. thermophilus reads were removed from the dataset, and the remaining reads for each sample were then assembled individually with metaSPAdes using default settings (v 3.11.0; Nurk et al., 2017). T. thermophilus sequences in each sample were enumerated (forward + reverse) and used to determine multipliers to calculate genes per liter according to Satinsky et al. (2013) (Supplementary Dataset 11, metadata tab). Briefly, the multipliers for each sample were calculated as Sa/(Sr * volume filtered) in which Sa is the number of molecules of T. thermophilus genomes added to the sample, and Sr is the number of T. thermophilus genomes recovered. Sr was calculated by dividing the number of T. thermophilus sequences recovered by the number of genes in the T. thermophilus genome (2,173).
Contig files for all samples were concatenated and dereplicated, and contigs less than 200 nucleotides in length were removed (Dedupe; Gregg and Eder, 2019). Concatenated contig sequences were submitted to IMG-MER1 for annotation with the DOE-JGI Metagenome Annotation Pipeline (Huntemann et al., 2016). After annotation, coding sequences (i.e., CDS or sequences coding for genes) were extracted from contigs to produce a Bowtie2 database of CDS sequences (Langmead and Salzberg, 2012). The quality-controlled sequences that were used for assembly (not including T. thermophilus sequences) were mapped to this database, and read counts were divided into separate files based on domain (Bacteria, Archaea, Eukaryote, Virus). For this study, only the Bacteria and Archaea sequences were further analyzed (see Supplementary Figure S1 for general workflow).
Read counts for each CDS were normalized following Wagner et al. (2012) for relative abundance and were also used to calculate absolute gene abundance following Satinsky et al. (2013). Briefly, the read count for each CDS was corrected for gene length and read length. To calculate genes per million (GPM) relative abundance for each CDS, this corrected read count (Tg) was scaled to 1 million [GPM = Tg * (1 × 109/ΣTg)]. To calculate genes per liter absolute abundance, Tg was scaled to the original sum of mapped reads [Tg_scaled = Tg * (ΣRg/ΣTg)], and then multiplied by the T. thermophilus multiplier for that sample (see above). Then for each of these measures of abundance we summed abundances for CDS assigned to the same KO number and to the same phylodist string (Supplementary Table S1). Metagenome sequences were deposited in NCBI Sequence Read Archive (SRA) bioproject accession number PRJNA642637 under accessions SRR12147740-SRR12147774.
Statistical Analyses
Environmental data were standardized using the range method and then log transformed using the biostats.R package in R (McGarigal, 2015). The transformed data were then analyzed with Principal Component Analysis (PCA) using the ggfortify package in R (Horikoshi and Tang, 2016; Tang et al., 2016), and plotted with the environmental variables shown as vectors.
To visualize microbial beta-diversity, a Bray-Curtis dissimilarity matrix was calculated from the table of gene abundances (GPM) and used in a non-metric multidimensional scaling (NMDS) analysis using the vegan package in R (Oksanen et al., 2019). Permutational multivariate analysis (PERMANOVA) was performed using the adonis function (999 permutations) from the vegan package in R to determine the variance in gene abundance associated with site, month, and fraction.
Percent relative abundances for taxa were plotted using ggplot2 (v2.2.1) with a cutoff at the class taxonomic level (Wickham, 2016). Abundance counts for replicate samples were averaged together. Counts per sample were normalized to percent relative abundance. Taxa with relative abundances less than 0.01 were classified as “other.”
The DESeq2 package in R was used to find significantly differentially abundant genes between the seasons (Love et al., 2014). While DESeq2 is commonly used in metatranscriptomics studies to analyze RNA sequence data, it has been shown to be a suitable method for identifying differentially abundant genes in metagenomic data (Dugat-Bony et al., 2015; Johnston et al., 2016; Jonsson et al., 2016). Gene counts were modeled following a negative binomial distribution. A likelihood ratio test was used as a hypothesis test, using a full and reduced model. The full model included the terms month and fraction, while the reduced model only included month to test if the fraction term increases the likelihood of the data. The likelihood ratio test is preferable when multiple interaction terms are being tested at once. False discovery rate (FDR) was controlled for using the Benjamini-Hochberg procedure with an alpha cutoff of 0.05. We further filtered CDS by using an adjusted p-value cutoff of < 0.05. Significant samples were transformed using variance stabilized transformation for visualization.
The log2fold changes of marker genes for carbon and nitrogen metabolic processes were assessed in order to look for variability in the potential for several nutrient cycling metabolic processes among months and size fractions. Previously identified marker genes, whose abundance can be used to assess the genetic potential of key steps in biogeochemical pathways, were identified in the DESeq2 output (Lauro et al., 2011; Llorens-Marès et al., 2015; Vavourakis et al., 2016; Kieft et al., 2018). Log2fold changes were plotted using the pheatmap package in R for the six different paired month comparisons. Overall, 47 different KO genes were used to assess 13 different metabolic pathways.
The function multipatt (multi-level pattern analysis) in the indicspecies package available through R (De Cáceres and Legendre, 2009) was used on the GPM gene abundances of only genes identified as significant by DESeq2 to identify significant indicator genes for each month. The setting delug = TRUE was used to limit analyses to single month combinations. This function compares association values between a gene and each of the 3 months, assigning a gene to the month with the highest association value. A permutation test was then performed to test the statistical significance of each relationship. The significant genes (p < 0.05) were analyzed with KEGG mapper to enumerate indicator genes present in KEGG pathways (Kanehisa and Sato, 2019). Only pathways featuring 10 or more indicator genes were analyzed, excluding several global and overview pathways. Pathways of interest were further analyzed by visualizing both relative abundance and genes/L abundance for each indicator gene in the pathway.
Results
Environmental Parameters
Principal component analysis of environmental data shows that samples cluster together by month on a horizontal axis where the primary axis explained 46% of the variation (Figure 2A), reflecting the physicochemical characteristics of the samples. Salinity and nitrate concentrations were highest in April, while chlorophyll, particulate and dissolved organic carbon, particulate nitrogen, and prokaryotic cell counts were highest in June. Water column temperature, however, is not clearly described in this analysis. Water temperature at the lagoon sites (KA and JA) was warmest in August (KA: 11.53°C, JA: 8.27°C), and water temperature at the BP site was warmest in June (1.55°C) (Supplementary Table S2). The secondary axis explained 29% of the variation, appearing to be largely influenced by the April KA sample which had higher total dissolved nitrogen and NH4, and lower dissolved oxygen.
Figure 2. (A) Principal component analysis plot of environmental measurements, and (B) non-metric multidimensional scaling plot of gene per million (GPM) data for KEGG-annotated protein coding genes. Stress of NMDS is 0.085.
Metagenome Dimension Reduction Analysis
For abundances of different genes detected in metagenome samples, month was the most significant factor explaining variation (PERMANOVA, R2 = 0.41, p < 0.001). Size fraction also explained a significant fraction of variation (R2 = 0.20, p < 0.001), but site was not significant (R2 = 0.07, P = 0.336). NMDS analysis of relative gene abundances showed a clear separation between June samples and the other 2 months (Figure 2B). The samples clustered together to a lesser extent based on size fraction. The NMDS had a stress of 0.085, and the range of the Bray-Curtis dissimilarity was 0.027–0.456.
Taxonomic Composition
Together, genes annotated to Alphaproteobacteria, Flavobacteria, and Gammaproteobacteria composed over 50% of the reads that mapped to CDS in all samples except for June free-living samples in which genes annotated to Betaproteobacteria also composed a large percentage of CDS (Supplementary Figure S2). In April, the KA samples differed from other samples, containing fewer Alphaproteobacteria genes and more Flavobacteria and Gammaproteobacteria genes. April samples also contained the highest percentages of genes from Euryarchaeota (ranging from 2.0 to 8.7%) and Thaumarchaeota (ranging from 0.9 to 4.4%). June samples differed from other months, with fewer Alphaproteobacteria genes, more Actinobacteria genes in the free-living fraction, more Cyanobacteria genes in the particle-attached fraction, and more Betaproteobacteria genes overall. The JA August particle-attached sample contained a higher percentage of Flavobacteria genes compared to the other August particle-attached samples. In both April and August, the free-living fraction had higher percentages of Alphaproteobacteria genes. Particle-attached samples in August had high percentages of Clostridia genes, up to 9.1% in the BP sample. Note that these gene abundances were not corrected for genome size, so taxa with larger genomes may be overrepresented relative to cell abundance.
Differential Gene Abundance and Indicator Analyses
DESeq2 identified 4,921 significant differentially abundant genes within the samples. Select marker genes, chosen based on previous studies (Lauro et al., 2011; Llorens-Marès et al., 2015; Hamilton et al., 2016; Kieft et al., 2018; Li et al., 2019; Shen et al., 2019), were used to identify seasonal shifts in genes involved in the nitrogen, carbon, and sulfur cycle (Figure 3). These genes formed four main groups. The first group included marker genes for fermentation and nitrification, that were more abundant in April than in June and August. The second group included marker genes for aerobic and anaerobic carbon fixation, aerobic respiration, organic nitrogen mineralization, inorganic nitrogen assimilation, and dissimilatory sulfate reduction/sulfide oxidation that were more abundant in April and August than in June. The third group included marker genes for denitrification, dissimilatory nitrate reduction, and nitrite oxidation and reduction that were more abundant in April and June than in August. The fourth group included marker genes for aerobic carbon fixation, aerobic respiration, carbon monoxide oxidation, fermentation, inorganic nitrogen assimilation, organic nitrogen mineralization, assimilatory sulfate reduction/oxidation, and sulfur mineralization that were more abundant in June and August than in April. Note that some marker genes used to detect the potential for anaerobic carbon fixation via the reductive TCA cycle (K00244, K00174, K00175) are not unique to that pathway and may participate in other metabolic processes.
Figure 3. Heatmap of z-scaled log2fold change determined by DeSeq2 of marker genes for ecological functions. Log2fold change determined as month to month comparisons, i.e., first column “Apr vs. Jun” is April gene abundance as compared to June gene abundances. Red indicates higher abundance and blue indicates lower abundance.
Indicator analysis found 4,113 genes that were indicators of month using a p-value of 0.05 and a test statistic cutoff of 0.8. Of these, 1,495 genes were indicators for April, 1,889 genes were indicators of June, and 729 genes were indicators of August. When classified into KEGG pathways, differences in gene diversity between months could be seen (Figure 4). Pathways were selected for further study based on their ability to provide insight into ecological function and archaeal abundance, using a cut-off of pathways with 10 or more indicator genes associated with a month. Broadly, these pathways categorize into four classes: energy metabolism, carbohydrate and other metabolism, nucleotide and amino acid biosynthesis and metabolism, and environmental and community interactions (Figure 4). This overview shows that the majority of genes that can be categorized to KEGG pathways are indicators for June and belong to processes such as energy and carbohydrate metabolism, as well as environmental interactions. A closer look at individual pathways revealed fine-scale shifts in indicator gene abundances between the months. Changes in individual gene abundances for nucleotide processing, the nitrogen cycle, two-component systems, and the ABC transporter system are described below.
Figure 4. Number of indicator genes for each month grouped by KEGG pathways. Pathways that included more than 10 indicator genes for all 3 months are shown.
Genes involved in nucleotide processing (DNA replication, RNA transport, and RNA polymerase) all had 11 or more indicator genes for April, with few indicator genes for the other months (Supplementary Table S3). Many of these indicator genes are specific to Archaea, such as RNA polymerase and parts of the DNA replication complex (Supplementary Figure S3). These archaeal nucleotide processing genes had abundances up to 1.6 × 109 genes/L in the April free-living samples.
The nitrogen cycle had a similar number of indicator genes for April and June, with June having 13 indicator genes and April having eight (Supplementary Table S3). April free-living samples had higher abundances of genes belonging to the nitrification process, while June samples (both fractions) had higher abundances of genes for denitrification and assimilatory and dissimilatory nitrate reduction (Figure 5).
Figure 5. Nitrogen cycle indicator genes (KEGG pathway map00910). Heatmap of z-scaled GPM gene abundances for each sample with yellow being higher abundance and blue lower. Bubble plot of genes/L abundance for each month and size-fraction pairing.
June had 68 indicator genes for two-component systems, while August had 15 and April had 13 (Supplementary Table S3). The June two-component system indicator genes were for cellular mobility, amino acid uptake and metabolism, nitrate/nitrite uptake, and aerobic respiration (Figure 6). Based on genes/L abundance, these genes were most abundant in the free-living June samples. Genes involved in aerobic respiration and flagellar assembly had the highest abundances, up to 1.2 × 109 to 2.5 × 109 genes/L, respectively. Genes involved in exopolysaccharide (EPS) production, manganese transport, and hyperosmotic stress/heat shock showed higher abundance in August estuarine (Jago and Kaktovik) samples. There were also six pufX photosystem genes that showed increased abundance in the August estuarine samples. All of these genes showed relatively low abundances in August (< 8.1 × 107 genes/L) but were absent from many of the April and June samples. April and August coastal (Bernard Point) samples showed an increased abundance of genes involved in stress, temperature, and misfolded proteins, with free-living April samples having the highest abundances, up to 2.2 × 108 genes/L.
Figure 6. Two-component system indicator genes (KEGG pathway ko02020). Heatmap of z-scaled GPM gene abundances for each sample with yellow being higher abundance and blue lower. Bubble plot of genes/L abundance for each month and size-fraction pairing.
June also had 58 indicator genes for the ABC transporter system, while April had seven and August had four (Supplementary Table S3). June particle-attached samples had higher abundances of capsular polysaccharide, lipopolysaccharide, xylitol, and s-methylcysteine transporter genes (Supplementary Figure S4). June free-living samples had higher abundances of cell-wall related transporter genes (e.g., teichoic acid and lipo-oligosaccharides), sugar transporters (e.g., erythritol, raffinose/stachyose/melibiose, and a-glucoside), glutamate/aspartate and biotin transporters, nitrogen transporters (e.g., nitrate/nitrite/cyanate), and sulfur transporters (e.g., alkanesulfonate, sulfate/thiosulfate). April samples had high abundances of indicator genes for the sugar transporter for α-1,4-digalacturonate, while August samples had one indicator gene each for vitamin B (biotin), lysine/arginine/ornithine/histidine, and manganese transporters.
Discussion
The Arctic lagoons spanning the northern coast of Alaska support diverse food webs that are subject to extreme changes in environmental conditions among seasons (Dunton and Schell, 1987; Churchwell et al., 2015). While higher trophic level interactions have been studied in this region of the coastal Arctic Ocean (von Biela et al., 2011; Harwood et al., 2015), less is known about the role bacteria and archaea play in C and N cycling and the coastal food web (Boeuf et al., 2013; Kirchman et al., 2014; McTigue et al., 2016; Kellogg et al., 2019). Using shotgun metagenomics, this study found that the metabolic potential of coastal bacterial and archaeal communities is strongly linked to season, with June being the most distinct microbial community (Figure 2). June communities included genes involved in nitrate reduction and in acquiring and processing organic matter from phytoplankton and from terrestrial environments. April communities saw an increase in the abundance of Archaea along with ammonium oxidation and chemoautotrophic carbon fixation. August had fewer indicator genes, suggesting that, metabolically, the summertime open-water microbial community is a transitional community between the spring community and the winter community. To identify ecologically relevant variability in the genomic capabilities of lagoon microbial communities, we investigated differential gene abundance in KEGG pathways involving (a) two-component systems, (b) nucleotide biosynthesis, (c) nitrogen cycling, and (d) carbon cycling.
Two-Component Systems
Two-component systems are signal-transduction pathways used by microbes to sense and respond to environmental stimuli (Hoch, 2000). June samples featured two-component system indicator genes associated with flagellar assembly and twitching mobility. These genes were largely from Alpha-, Beta-, and Gammaproteobacteria (Supplementary Table S4). Motility provides a fitness advantage to bacteria in environments that support strong gradients in resources such as those associated with phytoplankton blooms and the spring pulse of riverborne terrestrial material to the lagoons (Stocker and Seymour, 2012; Watteaux et al., 2015; Brumley et al., 2019). This is consistent with results from Tran et al. (2018) who found elevated expression of genes for flagella and chemotactic abilities by bacteria associated with phytoplankton blooms. Moreover, the high abundance of quorum sensing genes also found in June samples (Figure 4) may be linked to this higher abundance of motility genes, because motility and chemotaxis are partially regulated by quorum sensing (Sperandio et al., 2002; Zhou et al., 2016).
August samples featured only a small number of indicator genes, but these included several indicator genes for the pufX photosystem. Puf genes are commonly found in Rhodobacteraceae and code for light harvesting complex subunits used in aerobic anoxygenic photosynthesis (AAP), a photoheterotrophic process that harvests light for energy while still being capable of heterotrophy using organic carbon (Biel and Marrs, 1983; Zhu and Hearst, 1986; Farchaus et al., 1992). August indicator genes belonging to the puf operon were 60% Rhodobacteraceae when classified at the family level (Supplementary Table S5). An earlier study of 16S rRNA gene analysis at this study site found Rhodobacteraceae to be an indicator taxon for August samples (Kellogg et al., 2019). Rhodobacteraceae abundance has been linked to high nutrient conditions and phytoplankton blooms (González et al., 2000; Brown et al., 2005). One study showed that photoheterotrophic microbial communities off the coast of Utqiagvik, AK showed no change in relative abundance of pufM genes from winter to summer, however this study only looked at one gene and only sampled outside of the lagoons (Cottrell and Kirchman, 2009). Puf genes have also been found to be abundant in summer samples off the coast of Antarctica (Grzymski et al., 2012). Light has been found to enhance the abundance of AAPs (Ferrera et al., 2017), suggesting the open waters of August in the Beaufort lagoons allows this group of organisms to be successful.
Nucleotide Processing
Genes from pathways relating to DNA and RNA processing were indicators for the month of April (Supplementary Figure S3), and a large portion of these are Archaea specific genes as identified by KEGG (Supplementary Figure S2). Most (up to 97.9%) of these indicator genes were annotated to Crenarchaeota, Euryarchaeota, and Thaumarchaeota, with Thaumarchaeota having the highest average GPM of 3.6 (Supplementary Table S6). Across all CDS, the percentage of genes belonging to Archaea vs. Bacteria in April was 6.4%, while it was < 1% for the other 2 months. Paralleling this trend, Thaumarchaeota had a > 1% cumulative relative abundance in April but < 1% abundance in the other months (Supplementary Figure S2). In the Eastern Beaufort Sea, Archaea comprised 16% of total cell counts in the winter but became undetectable in the summer (Alonso-Sáez et al., 2008). It is therefore likely that the archaeal indicator genes for DNA and RNA processing in April are a result of increased abundance of Archaea. Coastal archaeal communities in the Beaufort Sea are taxonomically distinct from river and offshore-marine communities and are dominated by Euryarchaeota, Crenarchaeota, and Thaumarchaeota (Galand et al., 2006). These coastal archaeal communities vary in composition depending on if they are particle-attached or free-living (Galand et al., 2008). Similarly, this study finds Thaumarchaeota to be the most abundant archaeal taxa, with higher abundances in the free-living fraction (Supplementary Figure S2). This study is among the first to compare changes in archaeal communities over multiple seasons in the coastal Arctic Ocean (Christman et al., 2011; Müller et al., 2018).
Nitrogen Cycling
Nitrification
Nitrification genes including methane/ammonia monooxygenase (pmo/amoABC) and hydroxylamine dehydrogenase (hao) were indicators for the month of April, but KEGG does not distinguish between genes for pmoABC and amoABC. Across all samples, these genes belonged to nitrifying Nitrosomonadales (35%) and Thaumarchaeota (23%), and to methane oxidizing Methylococcales (41%), and the genes annotated to each of these groups were more abundant in April than the other months (Supplementary Table S7). The April peak in nitrification genes is partially explained by the increased relative abundance of Thaumarchaeota in April compared to other months (Supplementary Figure S2). However, Betaproteobacteria accounted for a large proportion of nitrification genes but showed the lowest relative abundance in April (Supplementary Figure S2). This suggests that there was a shift in Betaproteobacteria diversity during the winter toward organisms capable of nitrification. Higher abundance of the ammonia oxidizing organisms in April could be due to the higher levels of available nitrogen during this time of year. Several previous studies link seasonal increases in Thaumarchaeota abundance to elevated total nitrogen in the environment (Hollibaugh et al., 2011, 2014; Hong and Cho, 2015). Thaumarchaeota, which perform the first step of nitrification, have been found to become significantly enriched in under-ice communities of freshwater lakes (Butler et al., 2019), potentially playing a role in winter microbial biogeochemical cycling. The genus Nitrosopumilus, which accounts for most Archaea indicator genes for nucleotide processing (see above), are chemoautotrophs that release nitrogen-containing organic compounds and other labile DOM that may be limiting to heterotrophs in the community (Bayer et al., 2019).
Nitrification rates measured in the coastal Arctic Ocean range from 3.6 to 16 μmol N L–1 day–1, and nitrification accounts for the majority of ammonium uptake in the Arctic Ocean water column (Souza et al., 2014). The increased relative abundance of nitrification genes in winter is likely due to several controlling factors. First, ammonium concentration is highest in the lagoons in April (Kellogg et al., 2019). Second, phytoplankton abundance is lowest in the winter, thus ammonia oxidizing bacteria and archaea have less competition for ammonium since phytoplankton can typically outcompete ammonia oxidizing bacteria (AOB) and archaea (AOA) (Ward, 2005; Martens-Habbena et al., 2009). Third, light levels are reduced below ice in winter (Dunton, 1990), and nitrification can be inhibited by light (Ward et al., 1984; Souza et al., 2014). Fourth, nitrification rates in the coastal Chukchi Sea appear to be insensitive to changes in temperature (Baer et al., 2014) suggesting that the colder winter waters do not inhibit this process. One study in the coastal Chukchi and Beaufort Seas found that AOB and AOA amoA genes were 30- to 115-fold more abundant in winter compared to summer, mirroring higher potential nitrification rates in the winter as well as higher concentrations of available ammonium (Christman et al., 2011). Several other studies have also found potential nitrification rates in the Arctic to be higher in the winter (Baer et al., 2014; Souza et al., 2014). Nitrification rates have not been published for Arctic lagoons, but Connelly et al. (2014) found that the microbial community in these lagoons shifts from ammonium uptake in the summer to urea uptake in the winter. Sources of urea for the lagoons include riverine input, zooplankton excretion, melting seasonal ice, and sediment-associated bacteria (Connelly et al., 2014). Urea degradation has been found to enhance nitrification, suggesting that this urea is a potential source for nitrification in April (Swensen and Bakken, 1998; Sliekers et al., 2004).
Denitrification and Nitrate Reduction
Genes for nitrogen reduction were indicators for June including genes for nitrate reduction (nasA, narB), nitrite reduction (nirB, nirD), and denitrification (norC, norB, norZ). The genes nasA and narB are involved in nitrogen assimilation, while the genes nirB and nirD have both assimilative and dissimilative roles (Wang et al., 2019) including participation in the dissimilatory nitrate reduction to ammonia metabolism (DNRA; Heo et al., 2020; Figure 5). The typical marker gene for DNRA (nrfA, K03385; Pandey et al., 2020), though not identified as an indicator for June, was most abundant in June (Figure 3). There was also an elevated abundance of nitrate/nitrite ABC transporter genes in June, which could supply nitrate and nitrite for these metabolic processes (Supplementary Figure S4). The genes nirB and nirD were mainly associated with Gammaproteobacteria (55%) and Betaproteobacteria (19%), and the gene nrfA was associated with Flavobacteriaceae (61%), Gammaproteobacteria (22%), and Planctomycetes (17%). Genes for denitrification (norC, norB, nosZ, nirS, nirK) were mainly associated with Flavobacteriaceae (37%), Gammaproteobacteria (17%), and Alphaproteobacteria (17%) (Supplementary Tables S8, S9). Denitrification and DNRA are carried out by heterotrophic bacteria who use nitrate as an electron acceptor, but they require different sets of genes and have different end products: denitrification produces N2 gas while DNRA produces ammonium. Rates of denitrification and DNRA have not been published for the water column of Arctic lagoons, but previous studies in Arctic coastal sediments showed that sediment denitrification rates do not vary with season (Devol et al., 1997; Chang and Devol, 2009; McTigue et al., 2016). Similarly, a study in a Canadian lake found no change in water column denitrification rates between ice-covered and ice-free seasons (Cavaliere and Baulch, 2018). In contrast, our study suggests that the potential for water column denitrification and DNRA in the lagoons shows clear seasonality and peaks in June during the spring thaw.
When nitrate concentrations are high, availability of organic matter has been found to stimulate denitrification (Arango et al., 2007), and denitrification seems to be influenced by primary productivity and particulate organic carbon (Chang and Devol, 2009). In the Beaufort lagoons, the seasonal pulse of primary production by phytoplankton can reach the sediment before being remineralized, thus acting as a source of organic carbon for denitrification in both the water column and sediment (Sakshaug, 2004; Campbell et al., 2009; McTigue et al., 2016). While some previous studies have found evidence that anaerobic nitrogen processes can occur in the microhabitats of particle aggregates (Ganesh et al., 2014; Stief et al., 2016), other studies have found little difference in the contribution to denitrification between particle-attached and free-living microbial fractions (Donato et al., 2000; Chen et al., 2016). Consistent with this observation, we see little difference in abundance of denitrification genes between size fractions in June (within the same order of magnitude), suggesting that microbes capable of denitrification do not prefer one fraction (Figure 5).
Carbon Cycling
The organic carbon supporting lagoon food webs enters the lagoons from multiple sources during June, when river discharge brings in terrestrial organic matter, diatom blooms fix carbon, and melting ice releases organic matter (Connelly et al., 2015; Traving et al., 2017; Kellogg et al., 2019). This large input of relatively labile organic carbon supports the food web through the rest of the year when inputs are lower and more recalcitrant. The composition and quality of organic matter (OM) changes seasonally with June having high proportions of terrestrial material and diatom production, August being dinoflagellate and/or green algae influenced, and April having low concentrations of predominantly refractory OM (Connelly et al., 2015). Microbial community composition and function often reflect seasonal variability in organic matter composition (Massana et al., 2001; Crump et al., 2003; Teeling et al., 2012; Sipler et al., 2017; Traving et al., 2017). We found similar results for Arctic lagoons in which the abundances of carbon-metabolizing genes followed seasonal patterns similar to those hypothesized to occur in Arctic freshwater lakes (Crump et al., 2003). According to this hypothesis, spring thaw introduces relatively labile terrestrial organic matter as well as nutrients to support phytoplankton production, and this OM supports elevated bacterial productivity. In August, the flux of terrestrial OM decreases and is less labile (Holmes et al., 2008), phytoplankton production slows (e.g., Chl a, Supplementary Table S2), and the microbial food web is supported by OM from degraded terrestrial material and from phytoplankton using recycled nutrients (Connelly et al., 2015; Kellogg et al., 2019). By April, labile components of OM are degraded, phytoplankton production is very low, and the microbial community is supported by more recalcitrant organic matter or alternate energy sources such as chemoautotrophy. We investigated genes involved in carbon fixation and the cycling of phytoplankton organic matter, terrestrially derived organic matter, and methane.
Carbon Fixation
Seasonal indicator genes associated with prokaryotic carbon fixation pathways were most abundant in April (Figure 4) when the carbon pool is dominated by relatively recalcitrant organic matter, and photosynthesis is low. Therefore, microbes must find a way to metabolize this material or find alternative energy pathways, such as chemoautotrophic carbon fixation (Connelly et al., 2015). Most of the April indicator genes associated with carbon fixation belong to pathways used by chemoautotrophic organisms (Supplementary Figure S5) including the hydroxypropiontate-hydroxybutyrate cycle (HHC), found in Thaumarcheaota and other autotrophic Archaea (Berg et al., 2007), and the reductive TCA cycle (rTCA), found in a broad diversity of autotrophic bacteria (Berg, 2011). Most of these genes are not unique identifiers for chemoautotrophic carbon fixation. However, six genes coding for HHC (K18602–K18605, K18593, K18594) were solely annotated to the phylum Thaumarchaeota (Supplementary Table 13), and four genes coding for rTCA (K00169–K00172) were mainly annotated to the phyla Nitrospirae and Nitrospinae (52%), suggesting that these genes are associated with chemoautotrophic nitrifiers. The other genes associated with these pathways were annotated to a broader array of taxa, but primarily to Group II Marine Euryarchaea and Thaumarchaeota.
April peaks in genes associated with chemoautotrophic carbon fixation supports Connelly et al. (2014) who hypothesized microbial dark carbon fixation processes occurring off the coast of Utquiavik, AK based on 13C stable isotope probe data. Genes for HHC have also been found to be abundant in winter coastal communities from Antarctica (Grzymski et al., 2012). In contrast, many of the June indicators for carbon fixation belong to the Calvin Cycle (Supplementary Figure S5), which is used by photosynthetic organisms including cyanobacteria. Most of the indicator genes that coded for the Calvin Cycle are not unique markers of this process, but one Calvin Cycle marker gene, ribulose-bisphosphate carboxylase (K01602), peaked in abundance in the particle-attached fraction in June (Supplementary Figure S2). This gene was annotated to a broad diversity of bacterial including Rhizobiales (38%) and Burkholderiales (23%), but very few cyanobacteria (1%) suggesting that this gene is not associated with photoautotrophs in the lagoons.
Phytoplankton Organic Matter
June communities had a greater number of indicator genes for photosynthesis, pentose and glucuronate interconversions, propanoate metabolism, starch and sucrose metabolism, and fructose and mannose metabolism compared to the other 2 months (Figure 4). The photosynthesis indicator genes were photosystems I and II, and cytochrome b6f complex genes belonging to Chroococcales, Nostocales, Oscillatoriales, and Synechococcales cyanobacteria (Supplementary Table S10). Cyanobacteria had a higher abundance in June particle-attached samples (Supplementary Figure S2), although Kellogg et al. (2019) did not identify cyanobacteria as a June indicator species in unfiltered samples. Previous studies found cyanobacterial abundance to be higher in Arctic rivers and to decrease with increasing salinity along the coast, suggesting they are tied to spring thaw and are allochthonous to the lagoons (Garneau et al., 2006; Waleron et al., 2007; Vallières et al., 2008).
June indicator genes for propanoate metabolism may be linked to the availability of dimethylsulfoniopropionate (DMSP), an organosulfur compound produced by phytoplankton as an osmolyte and possible cryoprotectant (Kirst et al., 1991). Bacterial degradation of DMSP leads to acryloyl-CoA (Curson et al., 2008; Todd et al., 2009) or methionine and cysteine (Kiene and Linn, 2000), and these C3 compounds can be further processed in the propanoate pathway. One fate of DMSP-derived C3 compounds is to be transformed into acetyl-coA, which is then taken to the TCA cycle (Vila-Costa et al., 2010). DMSP is in higher concentrations in polar marine environments than temperate and tropical regions (Lana et al., 2011), and is in very high concentrations in association with sea-ice algae (Galindo et al., 2014). In the Canadian Arctic, DMSP concentrations were correlated with chlorophyll-a concentrations (Luce et al., 2011), and DMSP released from thawing sea ice was found to be rapidly consumed by microbes (Galindo et al., 2014). Increased levels of DMSP have been shown to correlate with the transcription of propanoate metabolism pathway genes in marine microbial communities, suggesting that this pathway is important for DMSP degradation (Vila-Costa et al., 2010). These previous studies are consistent with our finding of abundant genes for propanoate metabolism in June when diatoms were abundant (Connelly et al., 2015; Kellogg et al., 2019) and sea ice was thawing and potentially releasing DMSP into the water column.
Besides DMSP, phytoplankton blooms also release a variety of other low molecular weight (LMW) compounds such as organic acids and sugar alcohols toward the beginning of blooms, and later release high molecular weight (HMW) compounds such as polysaccharides and lipids (Buchan et al., 2014). These compounds impact microbial communities because LMW compounds can be used by most heterotrophic bacteria, but HMW compounds are only available to a restricted set of organisms (Logue et al., 2016). June indicator genes in the ABC transporters pathway included genes for transport of many of these HMW and LMW compounds (Supplementary Figure S4). This is consistent with previous studies suggesting that ABC transporter genes are expressed in response to dissolved phytoplankton material and other DOC (Poretsky et al., 2010; Williams et al., 2013). High abundances of these genes in June as well as genes for starch and sucrose metabolism, fructose and mannose metabolism, and glucose and glucuronate interconversions (Figure 4) may be linked to production of these compounds by spring phytoplankton blooms.
Terrestrial Organic Matter
Terrestrial DOM supplied by Arctic rivers is mainly composed of lignin-like compounds, followed by protein, lipid, and unsaturated hydrocarbon-like compounds (Sipler et al., 2017). The process of microbially mediated lignin degradation is still not well understood (Brown and Chang, 2014; Datta et al., 2017). While there have been genes identified in bacterial lignin degradation, few are annotated to the KEGG database (Datta et al., 2017). However, we found that the dye-decolorizing peroxidase (DyP, EC 1.11.1.19, KO:K15733) gene was significantly differentially abundant, and was most abundant in June and August and close to zero abundance in April. This gene is important in the bacterial oxidation of lignin, and also acts on a broad range of substrates including phenols, hydroquinones, dyes, amines, aromatic alcohols, and xenobiotics (van Bloois et al., 2010; de Gonzalo et al., 2016; Datta et al., 2017). This gene belonged to taxa that were not detected in April metagenomes including the Actinobacteria families Microbacteriaceae and Streptomycetaceae (Supplementary Table S11). In contrast with DyP, most degradation genes for aromatic compounds were indicators for either June or April (Figure 4), with the difference being that April had more genes involved in the anaerobic conversions of aromatic rings (four in April vs. 1 in June). These results suggest that spring and summer microbial communities possess elevated potential to degrade lignin and aromatic components of the terrestrial OM delivered to in spring. By April, however, labile compounds likely become scarce, thus making aromatic compounds (commonly considered refractory) more favorable for microbial degradation. Overall, the degradation of large organic molecules is a complex process, and our results suggest taxa capable of different OM degradation pathways vary in abundance throughout the year in response to the available forms of terrestrial organic matter.
Methane Metabolism
Seasonal indicator genes belonging to the methane metabolism KEGG pathway were most abundant in April. Three of these genes code for methane/ammonia monooxygenase subunits (K10944–K10946; see section “Nitrification”) and were annotated to methane oxidizing Methylococcales (52%), and ammonia oxidizing Nitrosomonadales (31%) and Thaumarchaeota (26%) (Supplementary Table S12). Methylococcales abundance was elevated in Arctic lagoons in April (Kellogg et al., 2019), and one study off the coast of Utquiavik, AK found increased abundance of Gammaproteobacteria methane oxidizing bacteria (MOB) in April seawater and a strong correlation between methylotrophs and the methane oxidation rate constant (Uhlig et al., 2018). Methane concentrations along the Beaufort Shelf increase during periods of ice cover, ranging from 0.14 to 0.43 mg CH4 m–2 day–1 in ice-free conditions to 0.28–1.01 mg CH4 m–2 day–1 during ice-cover (Lorenson et al., 2016). Decomposing sedimentary organic matter is the source of this methane, and microbial degradation of the sedimentary organic matter is the main pathway by which methane is released into the water column (Lorenson et al., 2016). This suggests that further study of the sediment microbial community is needed to better understand methane cycling in this system.
Conclusion
This study suggests that the water column microbial community in Arctic coastal lagoons plays a key role in processing both terrestrial and autochthonous organic matter throughout the year, though how it processes this material varies with season. These findings can inform future studies on Arctic lagoon food webs and the roles of microbes. This will be especially important in the future, because changes affecting seasonality such as rising temperature, changes in river discharge, freshening, and increase in open water season are expected to occur due to climate change (McClelland et al., 2006; Serreze and Barry, 2011; Morison et al., 2012; Barnhart et al., 2014). These changes can impact phytoplankton production, although various factors could increase or decrease production, and in turn alter organic matter availability and heterotrophic processes, likely altering the microbial community composition and function (Wassmann and Reigstad, 2011). Future work looking at the sediment microbial community as well as pairing metatranscriptomics with biogeochemical rate measurements would further our understanding of how microbes process organic matter and their connection to higher trophic levels.
Data Availability Statement
The datasets presented in this study can be found in the NCBI Short Read Archive BioProject# PRJNA642637 under accessions SRR12147740–SRR12147774, https://www.ncbi.nlm.nih.gov/.
Author Contributions
KB performed DNA sequence data analysis and lead writing. CK collected and processed samples, and conducted DNA sequencing. BC contributed to writing and interpretation. BC, JM, and KD conceptualized project. JM and KD led field efforts and collected samples. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
Funding was provided by the NSF grant ARC 1023465 and the Beaufort Lagoon Ecosystems Long Term Ecological Research program (BLE LTER) under NSF award 1656026.
Conflict of Interest
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.
Acknowledgments
We thank R. Thompson, S. Linn, S. Smith, C. Harris, T. Dunton, and J. Dunton for field sampling and L. Fine and M. Stuart for DNA extraction, OSU Center for Genome Research and Biocomputing for sequencing, and J. P. Payet, C. Sullivan, T. Sharpton and B. Kieft for bioinformatics and data analysis. Kaktovik Iñupiat Corporation and CH2MHill Polar Services provided logistical support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.601901/full#supplementary-material
Footnotes
References
Alonso-Sáez, L., Sánchez, O., Gasol, J. M., Balagué, V., and Pedrós-Alio, C. (2008). Winter-to-summer changes in the composition and single-cell activity of near-surface Arctic prokaryotes. Environ. Microbiol. 10, 2444–2454. doi: 10.1111/j.1462-2920.2008.01674.x
Arango, C. P., Tank, J. L., Schaller, J. L., Royer, T. V., Bernot, M. J., and David, M. B. (2007). Benthic organic carbon influences denitrification in streams with high nitrate concentration. Freshw. Biol. 52, 1210–1222. doi: 10.1111/j.1365-2427.2007.01758.x
Arrigo, K. R., and van Dijken, G. L. (2015). Continued increases in Arctic Ocean primary production. Prog. Oceanogr. 136, 60–70. doi: 10.1016/j.pocean.2015.05.002
Baer, S. E., Connelly, T. L., Sipler, R. E., Yager, P. L., and Bronk, D. A. (2014). Effect of temperature on rates of ammonium uptake and nitrification in the western coastal Arctic during winter, spring, and summer. Glob. Biogeochem. Cycles 28, 1455–1466. doi: 10.1002/2013GB004765
Barnhart, K. R., Anderson, R. S., Overeem, I., Wobus, C., Clow, G. D., and Urban, F. E. (2014). Modeling erosion of ice-rich permafrost bluffs along the Alaskan Beaufort Sea coast. J. Geophys. Res. Earth Surf. 119, 1155–1179. doi: 10.1002/2013JF002845
Bayer, B., Hansman, R. L., Bittner, M. J., Noriega-Ortega, B. E., Niggemann, J., Dittmar, T., et al. (2019). Ammonia-oxidizing archaea release a suite of organic compounds potentially fueling prokaryotic heterotrophy in the ocean. Environ. Microbiol. 21, 4062–4075. doi: 10.1111/1462-2920.14755
Beaufort Lagoon Ecosystems LTER and Core Program (2020a). Water column and sediment porewater nutrient concentrations from lagoon, river, and ocean sites along the Alaska Beaufort Sea coast, 2018-ongoing ver 3. EDI Data Portal doi: 10.6073/pasta/9acb5a6f3956476be643e88c403404e3
Beaufort Lagoon Ecosystems LTER and Core Program (2020b). Dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) from river, lagoon, and open ocean sites along the Alaska Beaufort Sea coast, 2018-ongoing ver 4. EDI Data Portal doi: 10.6073/pasta/4dcefec5d550301065619ed15b2192cc
Bell, L., Bluhm, B., and Iken, K. (2016). Influence of terrestrial organic matter in marine food webs of the Beaufort Sea shelf and slope. Mar. Ecol. Prog. Ser. 550, 1–24. doi: 10.3354/meps11725
Berg, I. A. (2011). Ecological aspects of the distribution of different autotrophic CO 2 fixation pathways. Appl. Environ. Microbiol. 77, 1925–1936. doi: 10.1128/AEM.02473-10
Berg, I. A., Kockelkorn, D., Buckel, W., and Fuchs, G. (2007). A 3-hydroxypropionate/4-hydroxybutyrate autotrophic carbon dioxide assimilation pathway in archaea. Science 318, 1782–1786. doi: 10.1126/science.1149976
Biel, A. J., and Marrs, B. L. (1983). Transcriptional regulation of several genes for bacteriochlorophyll biosynthesis in Rhodopseudomonas capsulata in response to oxygen. J. Bacteriol. 156, 686–694. doi: 10.1128/jb.156.2.686-694.1983
Boeuf, D., Cottrell, M. T., Kirchman, D. L., Lebaron, P., and Jeanthon, C. (2013). Summer community structure of aerobic anoxygenic phototrophic bacteria in the western Arctic Ocean. FEMS Microbiol. Ecol. 85, 417–432. doi: 10.1111/1574-6941.12130
Brown, M. E., and Chang, M. C. (2014). Exploring bacterial lignin degradation. Curr. Opin. Chem. Biol. 19, 1–7. doi: 10.1016/j.cbpa.2013.11.015
Brown, M. V., Schwalbach, M. S., Hewson, I., and Fuhrman, J. A. (2005). Coupling 16S-ITS rDNA clone libraries and automated ribosomal intergenic spacer analysis to show marine microbial diversity: development and application to a time series. Environ. Microbiol. 7, 1466–1479. doi: 10.1111/j.1462-2920.2005.00835.x
Brumley, D. R., Carrara, F., Hein, A. M., Yawata, Y., Levin, S. A., and Stocker, R. (2019). Bacteria push the limits of chemotactic precision to navigate dynamic chemical gradients. Proc. Natl. Acad. Sci. U.S.A. 116, 10792–10797. doi: 10.1073/pnas.1816621116
Buchan, A., LeCleir, G. R., Gulvik, C. A., and González, J. M. (2014). Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698. doi: 10.1038/nrmicro3326
Bushnell, B. (2014). BBMap: A Fast, Accurate, Splice-Aware Aligner (No. LBNL-7065E). Berkeley, CA: Lawrence Berkeley National Lab (LBNL).
Butler, T. M., Wilhelm, A.-C., Dwyer, A. C., Webb, P. N., Baldwin, A. L., and Techtmann, S. M. (2019). Microbial community dynamics during lake ice freezing. Sci. Rep. 9:6231. doi: 10.1038/s41598-019-42609-9
Campbell, R. G., Sherr, E. B., Ashjian, C. J., Plourde, S., Sherr, B. F., Hill, V., et al. (2009). Mesozooplankton prey preference and grazing impact in the western Arctic Ocean. Deep Sea Res. II Top. Stud. Oceanogr. 56, 1274–1289. doi: 10.1016/j.dsr2.2008.10.027
Cavaliere, E., and Baulch, H. M. (2018). Denitrification under lake ice. Biogeochemistry 137, 285–295. doi: 10.1007/s10533-018-0419-0
Chang, B. X., and Devol, A. H. (2009). Seasonal and spatial patterns of sedimentary denitrification rates in the Chukchi sea. Deep Sea Res. II Top. Stud. Oceanogr. 56, 1339–1350. doi: 10.1016/j.dsr2.2008.10.024
Chen, X., Jiang, H., Sun, X., Zhu, Y., and Yang, L. (2016). Nitrification and denitrification by algae-attached and free-living microorganisms during a cyanobacterial bloom in Lake Taihu, a shallow Eutrophic Lake in China. Biogeochemistry 131, 135–146. doi: 10.1007/s10533-016-0271-z
Christman, G. D., Cottrell, M. T., Popp, B. N., Gier, E., and Kirchman, D. L. (2011). Abundance, diversity, and activity of ammonia-oxidizing prokaryotes in the coastal arctic ocean in summer and winter. Appl. Environ. Microbiol. 77, 2026–2034. doi: 10.1128/AEM.01907-10
Churchwell, R. T., Kendall, S. J., Blanchard, A. L., Dunton, K. H., and Powell, A. N. (2015). Natural disturbance shapes benthic intertidal macroinvertebrate communities of high latitude river deltas. Estuaries Coasts 39, 798–814. doi: 10.1007/s12237-015-0028-2
Connelly, T. L., Baer, S. E., Cooper, J. T., Bronk, D. A., and Wawrik, B. (2014). Urea uptake and carbon fixation by marine pelagic bacteria and archaea during the arctic Summer and Winter Seasons. Appl. Environ. Microbiol. 80, 6013–6022. doi: 10.1128/AEM.01431-14
Connelly, T. L., McClelland, J. W., Crump, B. C., Kellogg, C. T. E., and Dunton, K. H. (2015). Seasonal changes in quantity and composition of suspended particulate organic matter in lagoons of the Alaskan Beaufort Sea. Mar. Ecol. Prog. Ser. 527, 31–45. doi: 10.3354/meps11207
Cory, R. M., Ward, C. P., Crump, B. C., and Kling, G. W. (2014). Sunlight controls water column processing of carbon in arctic fresh waters. Science 345, 925–928. doi: 10.1126/science.1253119
Cottrell, M. T., and Kirchman, D. L. (2009). Photoheterotrophic microbes in the arctic ocean in summer and winter. Appl. Environ. Microbiol. 75, 4958–4966. doi: 10.1128/AEM.00117-09
Crump, B. C., Armbrust, E. V., and Baross, J. A. (1999). Phylogenetic analysis of particle-attached and free-living bacterial communities in the Columbia river, its estuary, and the Adjacent Coastal Ocean. Appl. Environ. Microbiol. 65, 3192–3204. doi: 10.1128/aem.65.7.3192-3204.1999
Crump, B. C., Kling, G. W., Bahr, M., and Hobbie, J. E. (2003). Bacterioplankton community shifts in an arctic lake correlate with seasonal changes in organic matter source. Appl. Environ. Microbiol. 69, 2253–2268. doi: 10.1128/AEM.69.4.2253-2268.2003
Curson, A. R. J., Rogers, R., Todd, J. D., Brearley, C. A., and Johnston, A. W. B. (2008). Molecular genetic analysis of a dimethylsulfoniopropionate lyase that liberates the climate-changing gas dimethylsulfide in several marine α-proteobacteria and Rhodobacter sphaeroides. Environ. Microbiol. 10, 757–767. doi: 10.1111/j.1462-2920.2007.01499.x
Datta, R., Kelkar, A., Baraniya, D., Molaei, A., Moulick, A., Meena, R., et al. (2017). Enzymatic degradation of lignin in soil: a review. Sustainability 9:1163. doi: 10.3390/su9071163
De Cáceres, M., and Legendre, P. (2009). Associations between species and groups of sites: indices and statistical inference. Ecology 90, 3566–3574. doi: 10.1890/08-1823.1
de Gonzalo, G., Colpa, D. I., Habib, M. H. M., and Fraaije, M. W. (2016). Bacterial enzymes involved in lignin degradation. J. Biotechnol. 236, 110–119. doi: 10.1016/j.jbiotec.2016.08.011
DeLong, E. F., Franks, D. G., and Alldredge, A. L. (1993). Phylogenetic diversity of aggregate-attached vs. free-living marine bacterial assemblages. Limnol. Oceanogr. 38, 924–934. doi: 10.4319/lo.1993.38.5.0924
Déry, S. J., Hernández-Henríquez, M. A., Burford, J. E., and Wood, E. F. (2009). Observational evidence of an intensifying hydrological cycle in northern Canada. Geophys. Res. Lett. 36:e038852. doi: 10.1029/2009GL038852
Devol, A. H., Codispoti, L. A., and Christensen, J. P. (1997). Summer and winter denitrification rates in western Arctic shelf sediments. Continent. Shelf Res. 17, 1029–1050. doi: 10.1016/S0278-4343(97)00003-4
Doherty, M., Yager, P. L., Moran, M. A., Coles, V. J., Fortunato, C. S., Krusche, A. V., et al. (2017). Bacterial biogeography across the Amazon River-ocean continuum. Front. Microbiol. 8:882. doi: 10.3389/fmicb.2017.00882
Donato, S. A., Kato, K., Donald, D. A., Terai, H., and Tundisi, J. G. (2000). Contribution of free-living and attached bacteria to denitrification in the hypolimnion of a mesotrophic Japanese lake. Microb. Environ. 15, 93–101. doi: 10.1264/jsme2.2000.93
Dugat-Bony, E., Straub, C., Teissandier, A., Onésime, D., Loux, V., Monnet, C., et al. (2015). Overview of a surface-ripened cheese community functioning by meta-omics analyses. PLoS One 10:e0124360. doi: 10.1371/journal.pone.0124360
Dunton, K. H. (1990). Growth and production inLaminaria solidungula: relation to continuous underwater light levels in the Alaskan High Arctic. Mar. Biol. 106, 297–304. doi: 10.1007/BF01314813
Dunton, K. H., and Schell, D. M. (1987). Dependence of consumers on macroalgal (Laminaria solidungula) carbon in an Arctic kelp community: δ13C evidence. Mar. Biol. 93, 615–625. doi: 10.1007/bf00392799
Dunton, K. H., Schonberg, S. V., and Cooper, L. W. (2012). Food web structure of the alaskan nearshore shelf and estuarine lagoons of the Beaufort Sea. Estuar. Coasts 35, 416–435. doi: 10.1007/s12237-012-9475-1
Dunton, K. H., Weingartner, T., and Carmack, E. C. (2006). The nearshore western Beaufort Sea ecosystem: circulation and importance of terrestrial carbon in Arctic coastal food webs. Prog. Oceanogr. 71, 362–378. doi: 10.1016/j.pocean.2006.09.011
Farchaus, J. W., Barz, W. P., Grünberg, H., and Oesterhelt, D. (1992). Studies on the expression of the pufX polypeptide and its requirement for photoheterotrophic growth in Rhodobacter sphaeroides. EMBO J. 11, 2779–2788. doi: 10.1002/j.1460-2075.1992.tb05345.x
Ferrera, I., Sánchez, O., Kolářová, E., Koblížek, M., and Gasol, J. M. (2017). Light enhances the growth rates of natural populations of aerobic anoxygenic phototrophic bacteria. ISME J. 11, 2391–2393. doi: 10.1038/ismej.2017.79
Galand, P., Lovejoy, C., and Vincent, W. (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., Casamayor, E. O., Kirchman, D. L., Potvin, M., and Lovejoy, C. (2009). Unique archaeal assemblages in the Arctic Ocean unveiled by massively parallel tag sequencing. ISME J. 3, 860–869. doi: 10.1038/ismej.2009.23
Galand, P. E., Lovejoy, C., Pouliot, J., and Vincent, W. F. (2008). Heterogeneous archaeal communities in the particle-rich environment of an Arctic shelf ecosystem. J. Mar. Syst. 74, 774–782. doi: 10.1016/j.jmarsys.2007.12.001
Galindo, V., Levasseur, M., Mundy, C. J., Gosselin, M., Tremblay, J. -É, Scarratt, M., et al. (2014). Biological and physical processes influencing sea ice, under-ice algae, and dimethylsulfoniopropionate during spring in the Canadian Arctic Archipelago. J. Geophys. Res. Oceans 119, 3746–3766. doi: 10.1002/2013JC009497
Ganesh, S., Parris, D. J., DeLong, E. F., and Stewart, F. J. (2014). Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone. ISME J. 8, 187–211. doi: 10.1038/ismej.2013.144
Garneau, M. -È, Vincent, W. F., Alonso-Sáez, L., Gratton, Y., and Lovejoy, C. (2006). Prokaryotic community structure and heterotrophic production in a river-influenced coastal Arctic ecosystem. Aquat. Microb. Ecol. 42, 27–40. doi: 10.3354/ame042027
Garneau, M. -È, Vincent, W. F., Terrado, R., and Lovejoy, C. (2009). Importance of particle-associated bacterial heterotrophy in a coastal Arctic ecosystem. J. Mar. Syst. 75, 185–197. doi: 10.1016/j.jmarsys.2008.09.002
Ghiglione, J. F., Mevel, G., Pujo-Pay, M., Mousseau, L., Lebaron, P., and Goutx, M. (2007). Diel and seasonal variations in abundance, activity, and community structure of particle-attached and free-living bacteria in NW mediterranean sea. Microb. Ecol. 54, 217–231. doi: 10.1007/s00248-006-9189-7
González, J. M., Simó, R., Massana, R., Covert, J. S., Casamayor, E. O., Pedrós-Alió, C., et al. (2000). Bacterial community structure associated with a dimethylsulfoniopropionate-producing North Atlantic algal bloom. Appl. Environ. Microbiol. 66, 4237–4246. doi: 10.1128/AEM.66.10.4237-4246.2000
Gordeev, V. V., Martin, J. M. I, Sidorov, S., and Sidorova, M. V. (1996). A Reassessment of the Eurasian River Input of Water, Sediment, Major Elements, and Nutrients to the Arctic Ocean. Available online at: http://www.ajsonline.org/content/296/6/664.full.pdf (accessed June 18, 2019).
Gregg, F., and Eder, D. (2019). Dedupe. Available online at: https://github.com/dedupeio/dedupe (accessed January 2019).
Grzymski, J. J., Riesenfeld, C. S., Williams, T. J., Dussaq, A. M., Ducklow, H., Erickson, M., et al. (2012). A metagenomic assessment of winter and summer bacterioplankton from Antarctica Peninsula coastal surface waters. ISME J. 6, 1901–1915. doi: 10.1038/ismej.2012.31
Hamilton, T. L., Bovee, R. J., Sattin, S. R., Mohr, W., Gilhooly, W. P., Lyons, T. W., et al. (2016). Carbon and sulfur cycling below the chemocline in a meromictic lake and the identification of a novel taxonomic lineage in the FCB superphylum, Candidatus aegiribacteria. Front. Microbiol. 7:598. doi: 10.3389/fmicb.2016.00598
Harris, C. M., McClelland, J. W., Connelly, T. L., Crump, B. C., and Dunton, K. H. (2017). Salinity and temperature regimes in eastern alaskan beaufort sea lagoons in relation to Source Water Contributions. Estuar. Coasts 40, 50–62. doi: 10.1007/s12237-016-0123-z
Harris, C. M., McTigue, N. D., McClelland, J. W., and Dunton, K. H. (2018). Do high Arctic coastal food webs rely on a terrestrial carbon subsidy? Food Webs 15:e00081. doi: 10.1016/j.fooweb.2018.e00081
Harwood, L. A., Smith, T. G., George, J. C., Sandstrom, S. J., Walkusz, W., and Divoky, G. J. (2015). Change in the Beaufort Sea ecosystem: diverging trends in body condition and/or production in five marine vertebrate species. Prog. Oceanogr. 136, 263–273. doi: 10.1016/j.pocean.2015.05.003
Heo, H., Kwon, M., Song, B., and Yoon, S. (2020). Involvement of NO3- in ecophysiological regulation of dissimilatory nitrate/nitrite reduction to ammonium (dnra) is implied by physiological characterization of soil dnra bacteria isolated via a colorimetric screening method. Appl. Environ. Microbiol. 86:e01054-20. doi: 10.1128/AEM.01054-20
Hoch, J. A. (2000). Two-component and phosphorelay signal transduction. Curr. Opin. Microbiol. 3, 165–170. doi: 10.1016/S1369-5274(00)00070-9
Hollibaugh, J. T., Gifford, S., Sharma, S., Bano, N., and Moran, M. A. (2011). Metatranscriptomic analysis of ammonia-oxidizing organisms in an estuarine bacterioplankton assemblage. ISME J. 5, 866–878. doi: 10.1038/ismej.2010.172
Hollibaugh, J. T., Gifford, S. M., Moran, M. A., Ross, M. J., Sharma, S., and Tolar, B. B. (2014). Seasonal variation in the metratranscriptomes of a Thaumarchaeota population from SE USA coastal waters. ISME J. 8, 685–698. doi: 10.1038/ismej.2013.171
Holmes, R. M., McClelland, J. W., Peterson, B. J., Tank, S. E., Bulygina, E., Eglinton, T. I., et al. (2012). Seasonal and annual fluxes of nutrients and organic matter from large rivers to the Arctic Ocean and Surrounding Seas. Estuaries Coasts 35, 369–382. doi: 10.1007/s12237-011-9386-6
Holmes, R. M., McClelland, J. W., Raymond, P. A., Frazer, B. B., Peterson, B. J., and Stieglitz, M. (2008). Lability of DOC transported by Alaskan rivers to the Arctic Ocean. Geophys. Res. Lett. 35:e032837. doi: 10.1029/2007GL032837
Hong, J.-K., and Cho, J.-C. (2015). Environmental variables shaping the ecological niche of thaumarchaeota in soil: direct and indirect causal effects. PLoS One 10:e0133763. doi: 10.1371/journal.pone.0133763
Horikoshi, M., and Tang, Y. (2016). ggfortify: Data Visualization Tools for Statistical Analysis Results. Available online at: https://CRAN.R-project.org/package=ggfortify (accessed April 2018).
Huntemann, M., Ivanova, N. N., Mavromatis, K., Tripp, H. J., Paez-Espino, D., Tennessen, K., et al. (2016). The standard operating procedure of the DOE-JGI Metagenome Annotation Pipeline (MAP v.4). Stand. Genomic Sci. 11:17. doi: 10.1186/s40793-016-0138-x
Johnston, E. R., Rodriguez-R, L. M., 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
Jonsson, V., Österlund, T., Nerman, O., and Kristiansson, E. (2016). Statistical evaluation of methods for identification of differentially abundant genes in comparative metagenomics. BMC Genomics 17:78. doi: 10.1186/s12864-016-2386-y
Jorgenson, M. T., and Brown, J. (2005). Classification of the Alaskan Beaufort Sea Coast and estimation of carbon and sediment inputs from coastal erosion. Geo Mar. Lett. 25, 69–80. doi: 10.1007/s00367-004-0188-8
Kanehisa, M., and Sato, Y. (2019). KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. 29, 28–35. doi: 10.1002/pro.3711
Kellogg, C. T. E., McClelland, J. W., Dunton, K. H., and Crump, B. C. (2019). Strong seasonality in Arctic estuarine microbial food webs. Front. Microbiol. 10:2628. doi: 10.3389/fmicb.2019.02628
Kieft, B., Li, Z., Bryson, S., Crump, B. C., Hettich, R., Pan, C., et al. (2018). Microbial community structure–function relationships in yaquina bay estuary reveal spatially distinct carbon and nitrogen cycling capacities. Front. Microbiol. 9:1282. doi: 10.3389/fmicb.2018.01282
Kiene, R. P., and Linn, L. J. (2000). The fate of dissolved dimethylsulfoniopropionate (DMSP) in seawater: tracer studies using 35S-DMSP. Geochim. Cosmochim. Acta 64, 2797–2810. doi: 10.1016/S0016-7037(00)00399-9
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
Kirchman, D. L., Morán, X. A. G., and Ducklow, H. (2009). Microbial growth in the polar oceans — role of temperature and potential impact of climate change. Nat. Rev. Microbiol. 7, 451–459. doi: 10.1038/nrmicro2115
Kirst, G. O., Thiel, C., Wolff, H., Nothnagel, J., Wanzek, M., and Ulmke, R. (1991). Dimethylsulfoniopropionate (DMSP) in icealgae and its possible biological role. Mar. Chem. 35.1, 381–388. doi: 10.1016/s0304-4203(09)90030-5
Lana, A., Bell, T. G., Simó, R., Vallina, S. M., Ballabrera-Poy, J., Kettle, A. J., et al. (2011). An updated climatology of surface dimethlysulfide concentrations and emission fluxes in the global ocean. Glob. Biogeochem. Cycles 25:GB1004. doi: 10.1029/2010GB003850
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Lauro, F. M., DeMaere, M. Z., Yau, S., Brown, M. V., Ng, C., Wilkins, D., et al. (2011). An integrative study of a meromictic lake ecosystem in Antarctica. ISME J. 5, 879–895. doi: 10.1038/ismej.2010.185
Li, Y., Cha, Q.-Q., Dang, Y.-R., Chen, X.-L., Wang, M., McMinn, A., et al. (2019). Reconstruction of the functional ecosystem in the high light, low temperature union glacier region, Antarctica. Front. Microbiol. 10:2408. doi: 10.3389/fmicb.2019.02408
Llorens-Marès, T., Yooseph, S., Goll, J., Hoffman, J., Vila-Costa, M., Borrego, C. M., et al. (2015). Connecting biodiversity and potential functional role in modern euxinic environments by microbial metagenomics. ISME J. 9, 1648–1661. doi: 10.1038/ismej.2014.254
Lobbes, J. M., Fitznar, H. P., and Kattner, G. (2000). Biogeochemical characteristics of dissolved and particulate organic matter in Russian rivers entering the Arctic Ocean. Geochim. Cosmochim. Acta 64, 2973–2983. doi: 10.1016/S0016-7037(00)00409-9
Logue, J. B., Stedmon, C. A., Kellerman, A. M., Nielsen, N. J., Andersson, A. F., Laudon, H., et al. (2016). Experimental insights into the importance of aquatic bacterial community composition to the degradation of dissolved organic matter. ISME J. 10, 533–545. doi: 10.1038/ismej.2015.131
Lorenson, T. D., Greinert, J., and Coffin, R. B. (2016). Dissolved methane in the Beaufort Sea and the Arctic Ocean, 1992–2009; sources and atmospheric flux. Limnol. Oceanogr. 61, S300–S323. doi: 10.1002/lno.10457
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Luce, M., Levasseur, M., Scarratt, M. G., Michaud, S., Royer, S.-J., Kiene, R., et al. (2011). Distribution and microbial metabolism of dimethylsulfoniopropionate and dimethylsulfide during the 2007 Arctic ice minimum. J. Geophys. Res. Oceans 116:e006914. doi: 10.1029/2010JC006914
Martens-Habbena, W., Berube, P. M., Urakawa, H., de la Torre, J. R., and Stahl, D. A. (2009). Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature 461, 976–979. doi: 10.1038/nature08465
Massana, R., PedrósߚAlió, C., Casamayor, E. O., and Gasol, J. M. (2001). Changes in marine bacterioplankton phylogenetic composition during incubations designed to measure biogeochemically significant parameters. Limnol. Oceanogr. 46, 1181–1188. doi: 10.4319/lo.2001.46.5.1181
McClelland, J. W., Déry, S. J., Peterson, B. J., Holmes, R. M., and Wood, E. F. (2006). A pan-Arctic evaluation of changes in river discharge during the latter half of the 20th century. Geophys. Res. Lett. 33:GL025753. doi: 10.1029/2006GL025753
McClelland, J. W., Townsend-Small, A., Holmes, R. M., Pan, F., Stieglitz, M., Khosh, M., et al. (2014). River export of nutrients and organic matter from the North Slope of Alaska to the Beaufort Sea. Water Resour. Res. 50, 1823–1839. doi: 10.1002/2013WR014722
McGarigal, K. (2015). Biostats. Available online at: http://www.umass.edu/landeco/teaching/ecodata/labs/biostats (accessed February 25, 2015).
McTigue, N. D., and Dunton, K. H. (2014). Trophodynamics and organic matter assimilation pathways in the northeast Chukchi Sea, Alaska. Deep Sea Res. II Top. Stud. Oceanogr. 102, 84–96. doi: 10.1016/j.dsr2.2013.07.016
McTigue, N. D., Gardner, W. S., Dunton, K. H., and Hardison, A. K. (2016). Biotic and abiotic controls on co-occurring nitrogen cycling processes in shallow Arctic shelf sediments. Nat. Commun. 7, 1–11. doi: 10.1038/ncomms13145
Menzel, P., Ng, K. L., and Krogh, A. (2016). Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat. Commun. 7:11257. doi: 10.1038/ncomms11257
Morison, J., Kwok, R., Peralta-Ferriz, C., Alkire, M., Rigor, I., Andersen, R., et al. (2012). Changing Arctic Ocean freshwater pathways. Nature 481, 66–70. doi: 10.1038/nature10705
Müller, O., Wilson, B., Paulsen, M. L., Rumińska, A., Armo, H. R., Bratbak, G., et al. (2018). Spatiotemporal dynamics of ammonia-oxidizing thaumarchaeota in distinct arctic water masses. Front. Microbiol. 9:24. doi: 10.3389/fmicb.2018.00024
Nurk, S., Meleshko, D., Korobeynikov, A., and Pevzner, P. A. (2017). metaSPAdes: a new versatile metagenomic assembler. Genome Res. 27, 824–834. doi: 10.1101/gr.213959.116
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2019). vegan: Community Ecology Package. R package version 2.5-4. Available online at: https://CRAN.R-project.org/package=vegan (accessed April, 2018).
O’Leary, N. A., Wright, M. W., Brister, J. R., Ciufo, S., Haddad, D., McVeigh, R., et al. (2016). Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 44, D733–D745. doi: 10.1093/nar/gkv1189
Pandey, C. B., Kumar, U., Kaviraj, M., Minick, K. J., Mishra, A. K., and Singh, J. S. (2020). DNRA: a short-circuit in biological N-cycling to conserve nitrogen in terrestrial ecosystems. Sci. Total Environ. 738:139710. doi: 10.1016/j.scitotenv.2020.139710
Peterson, B. J. (2002). Increasing River Discharge to the Arctic Ocean. Science 298, 2171–2173. doi: 10.1126/science.1077445
Poretsky, R. S., Sun, S., Mou, X., and Moran, M. A. (2010). Transporter genes expressed by coastal bacterioplankton in response to dissolved organic carbon. Environ. Microbiol. 12, 616–627. doi: 10.1111/j.1462-2920.2009.02102.x
Rawlins, M. A., Steele, M., Holland, M. M., Adam, J. C., Cherry, J. E., Francis, J. A., et al. (2010). Analysis of the arctic system for freshwater cycle intensification: observations and expectations. J. Clim. 23, 5715–5737. doi: 10.1175/2010JCLI3421.1
Romanovsky, V. E., Smith, S. L., and Christiansen, H. H. (2010). Permafrost thermal state in the polar Northern Hemisphere during the international polar year 2007–2009: a synthesis. Permafrost Periglacial Process. 21, 106–116. doi: 10.1002/ppp.689
Sakshaug, E. (2004). “Primary and secondary production in the Arctic Seas,” in The organic carbon cycle in the Arctic Ocean, eds R. Stein and R. W. MacDonald (Berlin: Springer), 57–81.
Satinsky, B. M., Gifford, S. M., Crump, B. C., and Moran, M. A. (2013). “Use of internal standards for quantitative metatranscriptome and metagenome analysis,” in Methods in Enzymology, eds J. Abelson, M. Simon, G. Verdine, and A. Pyle (Amsterdam: Elsevier), 237–250. doi: 10.1016/b978-0-12-407863-5.00012-5
Serreze, M. C., and Barry, R. G. (2011). Processes and impacts of Arctic amplification: a research synthesis. Glob. Planet. Change 77, 85–96. doi: 10.1016/j.gloplacha.2011.03.004
Shen, M., Li, Q., Ren, M., Lin, Y., Wang, J., Chen, L., et al. (2019). Trophic status is associated with community structure and metabolic potential of planktonic microbiota in Plateau Lakes. Front. Microbiol. 10:2560. doi: 10.3389/fmicb.2019.02560
Sipler, R. E., Kellogg, C. T. E., Connelly, T. L., Roberts, Q. N., Yager, P. L., and Bronk, D. A. (2017). Microbial community response to terrestrially derived dissolved organic matter in the Coastal Arctic. Front. Microbiol. 8:1018. doi: 10.3389/fmicb.2017.01018
Sliekers, A. O., Haaijer, S., Schmid, M., Harhangi, H., Verwegen, K., Kuenen, J. G., et al. (2004). Nitrification and anammox with urea as the energy source. Syst. Appl. Microbiol. 27, 271–278. doi: 10.1078/0723-2020-00259
Smith, S. L., and Riseborough, D. W. (2010). Modelling the thermal response of permafrost terrain to right-of-way disturbance and climate warming. Cold Regions Sci. Technol. 60, 92–103. doi: 10.1016/j.coldregions.2009.08.009
Souza, A. C., Gardner, W. S., and Dunton, K. H. (2014). Rates of nitrification and ammonium dynamics in northeastern Chukchi Sea shelf waters. Deep Sea Res. II Top. Stud. Oceanogr. 102, 68–76. doi: 10.1016/j.dsr2.2013.12.017
Sperandio, V., Torres, A. G., and Kaper, J. B. (2002). Quorum sensing Escherichia coli regulators B and C (QseBC): a novel two-component regulatory system involved in the regulation of flagella and motility by quorum sensing in E. coli. Mol. Microbiol. 43, 809–821. doi: 10.1046/j.1365-2958.2002.02803.x
Stief, P., Kamp, A., Thamdrup, B., and Glud, R. N. (2016). Anaerobic nitrogen turnover by sinking diatom aggregates at varying ambient oxygen levels. Front. Microbiol. 7:98. doi: 10.3389/fmicb.2016.00098
Stocker, R., and Seymour, J. R. (2012). Ecology and physics of bacterial chemotaxis in the ocean. Microbiol. Mol. Biol. Rev. 76, 792–812. doi: 10.1128/MMBR.00029-12
Swensen, B., and Bakken, L. R. (1998). Nitrification potential and urease activity in a mineral subsoil. Soil Biol. Biochem. 30, 1333–1341. doi: 10.1016/S0038-0717(98)00015-7
Tang, Y., Horikoshi, M., and Li, W. (2016). ggfortify: unified interface to visualize statistical results of popular R packages. R J. 8, 474–485. doi: 10.32614/rj-2016-060
Teeling, H., Fuchs, B. M., Becher, D., Klockow, C., Gardebrecht, A., Bennke, C. M., et al. (2012). Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 336, 608–611. doi: 10.1126/science.1218344
Todd, J. D., Curson, A. R. J., Dupont, C. L., Nicholson, P., and Johnston, A. W. B. (2009). The dddP gene, encoding a novel enzyme that converts dimethylsulfoniopropionate into dimethyl sulfide, is widespread in ocean metagenomes and marine bacteria and also occurs in some Ascomycete fungi. Environ. Microbiol. 11, 1376–1385. doi: 10.1111/j.1462-2920.2009.01864.x
Tran, P., Ramachandran, A., Khawasik, O., Beisner, B. E., Rautio, M., Huot, Y., et al. (2018). Microbial life under ice: metagenome diversity and in situ activity of Verrucomicrobia in seasonally ice-covered Lakes. Environ. Microbiol. 20, 2568–2584. doi: 10.1111/1462-2920.14283
Traving, S. J., Rowe, O., Jakobsen, N. M., Sørensen, H., Dinasquet, J., Stedmon, C. A., et al. (2017). The effect of increased loads of dissolved organic matter on estuarine microbial community composition and function. Front. Microbiol. 8:351. doi: 10.3389/fmicb.2017.00351
Uhlig, C., Kirkpatrick, J. B., D’Hondt, S., and Loose, B. (2018). Methane-oxidizing seawater microbial communities from an Arctic shelf. Biogeosciences 15, 3311–3329. doi: 10.5194/bg-15-3311-2018
Underwood, G. J. C., Michel, C., Meisterhans, G., Niemi, A., Belzile, C., Witt, M., et al. (2019). Organic matter from Arctic sea-ice loss alters bacterial community structure and function. Nat. Clim. Change 9:170. doi: 10.1038/s41558-018-0391-7
Vallières, C., Retamal, L., Ramlal, P., Osburn, C. L., and Vincent, W. F. (2008). Bacterial production and microbial food web structure in a large Arctic river and the coastal Arctic Ocean. J. Mar. Syst. 74, 756–773. doi: 10.1016/j.jmarsys.2007.12.002
van Bloois, E., Torres Pazmiño, D. E., Winter, R. T., and Fraaije, M. W. (2010). A robust and extracellular heme-containing peroxidase from Thermobifida fusca as prototype of a bacterial peroxidase superfamily. Appl. Microbiol. Biotechnol. 86, 1419–1430. doi: 10.1007/s00253-009-2369-x
Vavourakis, C. D., Ghai, R., Rodriguez-Valera, F., Sorokin, D. Y., Tringe, S. G., Hugenholtz, P., et al. (2016). Metagenomic insights into the uncultured diversity and physiology of microbes in four hypersaline soda lake brines. Front. Microbiol. 7:211. doi: 10.3389/fmicb.2016.00211
Vila-Costa, M., Rinta-Kanto, J. M., Sun, S., Sharma, S., Poretsky, R., and Moran, M. A. (2010). Transcriptomic analysis of a marine bacterial community enriched with dimethylsulfoniopropionate. ISME J. 4, 1410–1420. doi: 10.1038/ismej.2010.62
von Biela, V. R., Zimmerman, C. E., and Moulton, L. L. (2011). Long-term increases in young-of-the-year growth of Arctic cisco Coregonus autumnalis and environmental influences. J. Fish Biol. 78, 39–56. doi: 10.1111/j.1095-8649.2010.02832.x
Wagner, G. P., Kin, K., and Lynch, V. J. (2012). Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 131, 281–285. doi: 10.1007/s12064-012-0162-3
Waleron, M., Waleron, K., Vincent, W. F., and Wilmotte, A. (2007). Allochthonous inputs of riverine picocyanobacteria to coastal waters in the Arctic Ocean. FEMS Microbiol. Ecol. 59, 356–365. doi: 10.1111/j.1574-6941.2006.00236.x
Wang, X., Tamiev, D., Alagurajan, J., DiSpirito, A. A., Phillips, G. J., and Hargrove, M. S. (2019). The role of the NADH-dependent nitrite reductase, Nir, from Escherichia coli in fermentative ammonification. Arch. Microbiol. 201, 519–530. doi: 10.1007/s00203-018-1590-3
Ward, B. B. (2005). Temporal variability in nitrification rates and related biogeochemical factors in Monterey Bay, California, USA. Mar. Ecol. Prog. Ser. 292, 97–109. doi: 10.3354/meps292097
Ward, B. B., Talbot, M. C., and Perry, M. J. (1984). Contributions of phytoplankton and nitrifying bacteria to ammonium and nitrite dynamics in coastal waters. Continent. Shelf Res. 3, 383–398. doi: 10.1016/0278-4343(84)90018-9
Wassmann, P., and Reigstad, M. (2011). Future arctic ocean seasonal ice zones and implications for pelagic-benthic coupling. Oceanography 24, 220–231. doi: 10.5670/oceanog.2011.74
Watteaux, R., Stocker, R., and Taylor, J. R. (2015). Sensitivity of the rate of nutrient uptake by chemotactic bacteria to physical and biological parameters in a turbulent environment. J. Theor. Biol. 387, 120–135. doi: 10.1016/j.jtbi.2015.08.006
Williams, T. J., Wilkins, D., Long, E., Evans, F., DeMaere, M. Z., Raftery, M. J., et al. (2013). The role of planktonic Flavobacteria in processing algal organic matter in coastal East Antarctica revealed using metagenomics and metaproteomics. Environ. Microbiol. 15, 1302–1317. doi: 10.1111/1462-2920.12017
Zhou, J., Lyu, Y., Richlen, M., Anderson, D. M., and Cai, Z. (2016). Quorum sensing is a language of chemical signals and plays an ecological role in algal-bacterial interactions. CRC Crit. Rev. Plant Sci. 35, 81–105. doi: 10.1080/07352689.2016.1172461
Zhu, Y. S., and Hearst, J. E. (1986). Regulation of expression of genes for light-harvesting antenna proteins LH-I and LH-II; reaction center polypeptides RC-L, RC-M, and RC-H; and enzymes of bacteriochlorophyll and carotenoid biosynthesis in Rhodobacter capsulatus by light and oxygen. PNAS 83, 7613–7617. doi: 10.1073/pnas.83.20.7613
Keywords: estuary, archaea, bacteria, omics, arctic national wildlife refuge
Citation: Baker KD, Kellogg CTE, McClelland JW, Dunton KH and Crump BC (2021) The Genomic Capabilities of Microbial Communities Track Seasonal Variation in Environmental Conditions of Arctic Lagoons. Front. Microbiol. 12:601901. doi: 10.3389/fmicb.2021.601901
Received: 01 September 2020; Accepted: 04 January 2021;
Published: 12 February 2021.
Edited by:
Lasse Riemann, University of Copenhagen, DenmarkReviewed by:
Linta Reji, Princeton University, United StatesBongkeun Song, College of William and Mary, United States
Copyright © 2021 Baker, Kellogg, McClelland, Dunton and Crump. 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(s) 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: Byron C. Crump, byron.crump@oregonstate.edu