Woeseiales transcriptional response to shallow burial in Arctic fjord surface sediment.

Distinct lineages of Gammaproteobacteria clade Woeseiales are globally distributed in marine sediments, based on metagenomic and 16S rRNA gene analysis. Yet little is known about why they are dominant or their ecological role in Arctic fjord sediments, where glacial retreat is rapidly imposing change. This study combined 16S rRNA gene analysis, metagenome-assembled genomes (MAGs), and genome-resolved metatranscriptomics uncovered the in situ abundance and transcriptional activity of Woeseiales with burial in four shallow sediment sites of Kongsfjorden and Van Keulenfjorden of Svalbard (79°N). We present five novel Woeseiales MAGs and show transcriptional evidence for metabolic plasticity during burial, including sulfur oxidation with reverse dissimilatory sulfite reductase (dsrAB) down to 4 cm depth and nitrite reduction down to 6 cm depth. A single stress protein, spore protein SP21 (hspA), had a tenfold higher mRNA abundance than any other transcript, and was a hundredfold higher on average than other transcripts. At three out of the four sites, SP21 transcript abundance increased with depth, while total mRNA abundance and richness decreased, indicating a shift in investment from metabolism and other cellular processes to build-up of spore protein SP21. The SP21 gene in MAGs was often flanked by genes involved in membrane-associated stress response. The ability of Woeseiales to shift from sulfur oxidation to nitrite reduction with burial into marine sediments with decreasing access to overlying oxic bottom waters, as well as enter into a dormant state dominated by SP21, may account for its ubiquity and high abundance in marine sediments worldwide, including those of the rapidly shifting Arctic.

Introduction

Seafloor sediments host a significant portion of Earth’s biomass [1]. Microorganisms must have physiological adaptations to be active in marine sediments. Steep gradients in oxygen, redox active compounds, organic matter compositions, n class="Chemical">metal distributions, and metabolite pools are some of the many characteristics to which a microorganism must be able to respond [2-4]. Alternatively, increasing energy limitation with depth in the subseafloor environment may promote a strategy of metabolic standby until conditions again become optimal [5, 6], perhaps through exhumation through bioturbation, sediment slumping, or storm events in shallow waters. The particular strategy used by sedimentary microbes will likely depend on the sedimentological regime, degree of bioturbation, and depositional history.

One group that may have adaptive strategies for dealing with diverse sediment environments is the order Woeseiales (formerly ‘JTB255’) [7]. The Woeseiales are a globally-distributed group of Gammaproteobacteria that have been shown through 16S rRNA gene amplicon sequence surveys, as well as catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH), to be highly abundant in both coastal and deep-sea surface sediments [7-11]. Metagenomic studies [7] and 16S rRNA gene oligotyping studies [12] have revealed that this clade is composed of distinct lineages that inhabit either coastal or deep-sea sediments, suggesting that subclades within Woeseiales may interact with their sediment envn class="Chemical">ironment differently. Single-cell amplified genomes (SAGs), metagenome-assembled genomes (MAGs), transcriptomics, and incorporation of 13C-labeled inorganic carbon showed the potential for chemolithoautotrophy via the oxidation of reduced sulfur compounds and the Calvin Benson cycle in coastal subclades [13]. A stronger potential for heterotrophy via oxidation of carbohydrates, peptides, aromatics [7] and molecular hydrogen [10] has been suggested for deep-sea subclades [12]. Interestingly, the only cultured representative, Woeseia oceani XK5, from a coastal sediment sample, appears incapable of autotrophic growth [14]. This supports the idea that different populations of Woeseiales may have fundamentally different physiologies. Even single representatives of Woesiales have been hypothesized to be metabolic generalists, with evidence for a high degree of metabolic flexibility [15].

The physiological responses of individual populations to burial in marine sediments, such as transcriptional regulation and dormancy strategies, remain largely unexamined and may be key to understanding how Woeseiales taxa modulate their behavior in different sediment envn class="Chemical">ironments. We hypothesize that a high depth-resolution metatranscriptomic analysis will demonstrate metabolic plasticity as well as strategies for long-term survival and subsistence in the sediment environment. Altogether, this work will shed light on the ability to extrapolate the functions made previously from genomic inferences of Woeseiales.

Woeseiales have also been found in Arctic fjord sediments in Svalbard, Norway [11, 13]. The Arctic is currently experiencing climatic changes from increased temperatures [16], frequency and severity of storm events [17], and glacial retreat [18]. The ability to predict whether members of a globally important microbial taxon such as Woeseiales will be resilient to this change requires understanding the mechanisms by which it responds to the stresses of burial in marine sediments. We therefore examined changes in gene transcript recruitment with depth to MAGs of uncultured Woeseiales over the upper 5 to 6 cm of sediments from two Svalbard fjords: Kongsfjorden and Van Keulenfjorden. Here, glaciers terminate in marine-influenced fjord systems, where high rates of sediment loading (5–10 cm/year [19, 20]) derive directly from glacial run-off. However, as glaciers retreat, this sediment loading may decrease [20, 21]. We examined two sites from Kongsfjorden (sites F and P) and one site from Van Keulenfjorden (site AC) with heavy loading of glacially-derived sediment, as well as another site from Van Keulenfjorden (site AB) with minimal input of glacially-derived sediments. Five novel medium to high quality Woeseiales n class="Chemical">MAGs were reconstructed from metagenomic sequences from these four sites. Mapping of transcripts to these MAGs demonstrated broad decreases in transcription with depth for the majority of genes, while one gene, hspA (encoding spore protein SP21), comparatively increased transcription with depth by multiple orders of magnitude.

Results

Sites within Van Keulenfjorden (sites AB and AC) and Kongsfjorden (sites P and F), Svalbard, were investigated for their geochemical properties and relative abundance of Woeseiales with depth (Fig 1a–1c). Previous geochemical analysis on porewaters at these sites showed higher relative concentrations of Fe (up to ~300 nM) at site AB compared to <50 nM at site AC in the top 5 cm of sediment [22]. Within shallow sediments at site AB, Mn concentrations never exceed 65 nM, while at site AC, Mn concentrations exceed 250 nM [22]. Similarly high n class="Chemical">Fe concentrations (~300 nM) have been observed in shallow sediments at site P in Kongsfjorden, with Mn concentrations not exceeding 20 nM [23].

Fig 1

Physical and geochemical setting.

Overview map of Spitzbergen (a), containing red box indicating Kongsfjorden (b) and green box indicating Van Keulenfjorden (c). Landsat imagery courtesy of NASA Goddard Space Flight Center and U.S. Geological Survey. The relative abundance of Woeseiales by 16S rRNA gene amplicon analysis for sites AB (d) and AC (e) in Van Keulenfjorden are shown along with measurements of hydrogen (f), total organic carbon (g), carbon to nitrogen ratios (h), and stable isotopes of organic carbon (i). Outliers were determined with Cook’s distance.

Overview map of Spitzbergen (a), containing red box indicating Kongsfjorden (b) and green box indicating Van Keulenfjorden (c). Landsat imagery courtesy of NASA Goddard Space Flight Center and U.S. Geological Survey. The relative abundance of Woeseiales by 16S rRNA gene amplicon analysis for sites AB (d) and AC (e) in Van Keulenfjorden are shown along with measurements of n class="Chemical">hydrogen (f), total organic carbon (g), carbon to nitrogen ratios (h), and stable isotopes of organic carbon (i). Outliers were determined with Cook’s distance.

The relative abundance of Woeseiales 16S rRNA gene amplicons remained steady with depth in Van Keulenfjorden sediment, between 1–3% of 16S rRNA gene amplicon libraries down to ~20 cm at sites AB and AC in Van Keulenfjorden (Fig 1d and 1e). 16S rRNA genes did not amplify for sites P and F in Kongsfjorden, likely due to the high concentrations of Fe, which can intern class="Chemical">fere with DNA extraction and amplification [24].

In Kongsfjorden, hydrogen, total organic carbon (TOC), C/N ratios, and stable isotopes of organic carbon were variable in the upper five cm (Fig 1f–1i). This is similar to sites AB and AC in Van Keulenfjorden [22, 25]. High C/N ratios and low TOC values at sites AC, F, and P indicate the deposition of refractory terrestrial organic material relative to site AB where corresponding values were lower [22, 25].

Recruitment of metagenomic and metatranscriptomic reads to Woeseiales MAGs

Metagenomic sequence libraries were generated from the depth-integrated upper five cm at three of the four sites (AB, AC, and F), with read numbers ranging from 30,615,734 to 83,705,875 per library. The metagenomic sequencing failed for site P, perhaps due to the same coextraction of inhibitors that probably prevented amplification of 16S rRNA genes. Five Woeseiales MAGs were reconstructed from the metagenomic data (S1 Fig) and found to be closely related to the Woeseiales Lineage 11 representative JSS Woes1 [12]. Amino acid identity analysis shows that the majority of proteins in our n class="Chemical">MAGs match unclassified Gammaproteobacteria and/or Chromatiales (S2 Fig). Less than half of the proteins in the Woeseiales MAGs have average amino acid identity to the proteins from the isolate Woeseia oceani, pointing to proteomic novelty in our environmental MAGs compared to the isolate.

MAG completeness ranged from 80 to 94%, with 2 to 9% contamination (S1 Table). Woeseiales collectively recruited from 24% of 6,145,182 reads mapped to the metagenome at site F to 33% of 15,328,660 reads mapped to the n class="Species">metagenome at site AC, highlighting the importance of this clade to the sediment microbial community. Across all transcriptomes (which ranged from 272,522 reads to 1,307,568 reads after removal of rRNA), Woeseiales persistently recruited from 0.5 to 3% of mapped reads. A high similarity in predicted gene functions with KEGG-decoder [26] between the five Woeseiales MAGs permitted a combined analysis for transcripts (S3 Fig).

Most genes that recruited transcripts (54% of total transcribed genes) were recruited only from one of the four sites (S4 and S5 Figs). However, the majority of these transcripts that only appeared at one site (93%) occurred only in a single depth. We therefore restricted analyses to transcripts occurring at more than one depth, to filter out potential noise. This showed transcripts unique to a single site included electron transport proteins, proteins involved in cellular signaling, and chaperones (S5 Fig). Specifically at site AB, we identified transcripts for the oxidoreductase complex protein subunits RnfE (Em = -500 mV [27]) and RsxB, which are involved in electron transn class="Chemical">fer to regulatory proteins that respond to oxidative stress [28]. At site AC, we identified exclusive transcripts for bifunctional (p)ppGpp synthase/hydrolase SpoT. This enzyme catalyzes the synthesis and degradation of guanosine pentaphosphate ((p)ppGpp), an alarmone involved in signaling and bacterial stringent response to nutritional starvation and the transition to a state of dormancy under stress [29]. At site F, we detected unique transcripts for chaperone ClpB. This chaperone is involved in recovery after heat-induced damage and works in concert with the DnaK chaperone system (transcribed at F, AB, and P) to refold denatured protein aggregates [30]. Finally, at site P, we identified unique transcripts for RNA polymerase sigma factor SigV, an extracytoplasmic function sigma factor involved in transcription under cell envelope stress [31].

Among the genes that were transcribed at all sites (8% of total transcribed genes; S4 Fig), about 20% encoded proteins involved in energy production and conversion (S2 Table). This included transcripts for the central metabolism genes encoding citrate synthase, acetyl-coenzyme A synthase, and isocitrate lyase. In addition, several genes involved in sulfide oxidation metabolism were found at all sites. This included the cytochrome subunit of n class="Chemical">sulfide dehydrogenase (fccA, present in all MAGs) and the beta subunit of dissimilatory sulfite reductase (dsrB). Sequence analysis of the alpha and beta subunits (dsrAB) within two Woeseiales MAGs (Woeseia sp. stnF and Woeseia sp. stnAC-2) confirms that members of this clade encode a reverse dissimilatory sulfite reductase involved in oxidation of sulfur (Fig 2) [7]. In Woeseiales, this may facilitate thiosulfate oxidation through a truncated SOX pathway involving SoxXA, SoxYZ, SoxCD, SoxB and a c-type cytochrome (c-555) (Fig 3a). The oxidation of reduced sulfur and sulfur intermediates is potentially linked to carbon fixation through the Calvin-Benson Cycle using an encoded Type II RubisCO (Fig 3a, S6 Fig; S1 File); however, we were not able to detect transcripts for sox genes. Transcription of dsr subunits decreased with depth (Fig 3b).

Fig 2

Phylogeny of DsrA encoded on Woeseiales MAGs.

The CIPRES Science Gateway v. 3.3 online server was used for RAxML and the tree was visualized with the Interactive Tree of Life (iToL) visualization tool. Nodes with bootstrap values equal to or greater than 70% are indicated with a black circle. Reductive-type archaeal and bacterial sequences are represented in pink and orange fields, respectively, with solid tree branches. Oxidative-type bacterial sequences are represented in the blue field with dashed tree branches. Woeseiales DsrA sequences are highlighted in red text.

Fig 3

Sulfur oxidation and denitrification in Woeseiales MAGs.

Schematic representation of metabolic pathways reconstructed in Woeseiales MAGs (a). The gray/black bar at the top corresponds to the redox potential for the key cytochromes depicted in the figure (sulfide quinone reductase, cbb3-type c oxidase, cyt bc1 complex, NADPH oxidoreductase, and the Rnf complex). Colors in circles next to each protein represent the presence in the MAGs. Depth trends in transcriptional coverage (reported in TPM) for nirS and the alpha and/or beta subunits of dsr are reported (b).

Schematic representation of metabolic pathways reconstructed in Woeseiales MAGs (a). The gray/black bar at the top corresponds to the redox potential for the key cytochromes depicted in the figure (n class="Chemical">sulfide quinone reductase, cbb3-type c oxidase, cyt bc1 complex, NADPH oxidoreductase, and the Rnf complex). Colors in circles next to each protein represent the presence in the MAGs. Depth trends in transcriptional coverage (reported in TPM) for nirS and the alpha and/or beta subunits of dsr are reported (b).

Another gene that was transcribed at all sites was nirS. This gene encodes nitrite reductase, a key enzyme in denitrification. Transcripts for nirS were detected at every sediment depth (Fig 3b). Similar to rdsrAB, transcription of nirS decreased with depth. Transcripts of nirS were detected in deeper sediment depths compared to transcripts of dsrAB, which did not penetrate past 4–5 cmbsf (Fig 3b). While transcription levels of both genes decreased overall, transcription of nirS remained elevated below 2 cm, after which dsr transcripts were nearly zero at most sites. This may indicate a transcriptional response to decreased availability of respective substrates, with n class="Chemical">oxygen and intermediates of sulfur becoming depleted at shallower depths relative to nitrite [32-34].

Transcriptional trends of hspA (SP21)

The gene hspA, encoding SP21, also called heat-shock protein A, recruited at least ten times more transcripts in our Woeseiales MAGs than the next most highly transcribed gene all sites (chaperone protein DnaK, 1,916 TPM; S7 Fig). SP21 is a small, membrane-associated stress protein related to heat shock proteins [35, 36] (S8 Fig). This gene had very high transcriptional coverage, with values an order of magnitude greater on average (TPM = 2,770) compared to all other genes (average TPM = 531). This gene also had the largest transcription observed for any gene at any site, with 309,186 TPM at Site P, 5–6 cmbsf (S7 Fig).

Depth analysis of SP21 coverage with post-hoc Tukey testing showed statistically-significant increases with depth at sites AB and P (ANOVA p < 0.05; Fig 4). At sites AC and F, transcriptional coverage did not show statistically-supported trends with depth (ANOVA p > 0.05). Instead, at these sites, transcription of spore protein SP21 either remains high at all depths (site F) or in the low-to-middle-range of coverage at all depths (site AC).

Fig 4

Transcriptional coverage of SP21 with depth at each site.

Each point represents a copy of the gene encoding SP21 (hspA) in a Woeseiales MAG. Some MAGs have more than one copy of the gene (see S9 Fig). The upper and lower hinges correspond to the 75th and 25th percentiles, respectively, and the median is represented by a horizontal bar. Statistical differences in mean values across depths for a single site was calculated with Tukey means testing. Depths that do not share a letter/color have statistically significant different mean TPM values across the different copies of the gene.

Each point represents a copy of the gene encoding SP21 (hspA) in a Woeseiales MAG. Some MAGs have more than one copy of the gene (see S9 Fig). The upper and lower hinges correspond to the 75th and 25th percentiles, respectively, and the median is represented by a horizontal bar. Statistical difn class="Chemical">ferences in mean values across depths for a single site was calculated with Tukey means testing. Depths that do not share a letter/color have statistically significant different mean TPM values across the different copies of the gene.

In addition to SP21, transcripts for other stress proteins, transcriptional regulators, and starvation genes were detected (S9 Fig). However, transcript recruitment to these genes was much lower than that observed for SP21, at 89 TPM on average. In addition, none of these genes showed increases with depth, even those co-localized on the same contig as the spore protein SP21 gene hspA (S10 Fig). In fact, transcripts for the co-localized chaperonins encoded by groS and groL decreased with depth (S9 Fig), while other chaperones either showed no trend (co-localized chaperon DnaK) or were only detected at the surface (chaperone SurA and chaperone HtpG; S9 Fig).

Other genes flanking the SP21 gene hspA included those encoding the protein-methionine-sulfoxide reductase MsrPQ system (msrP and msrQ; S10 Fig). This membrane-associated system works with the periplasmic chaperone SurA to repair oxidatively damaged periplasmic proteins containing methionine sulfoxide residues, thereby maintaining cell envelope integrity [37]. Transcripts for chaperone SurA were identified in surface sediments of sites AB and AC (S9 Fig). We identified several additional nearby genes involved in membrane-associated processes. These included genes that participate in toxin removal/efflux (tolC transcribed at AB, P, and F; xcpT and ephA) [38-40], protein repair (degQ transcribed at AC, F, and P; pcm) [41, 42], and cell wall reorganization (uppP and gluP) [43].

Transcriptional coverage and richness with depth

To get a broader understanding of transcriptional responses to burial, we classified Woeseiales transcripts into Clusters of Orthologous Groups (COGs) (Fig 5). We took this approach because low overall transcriptional coverage and inconsistent detection of transcripts for an individual gene across several depths precluded analysis at the gene-level. Spearman’s rho was calculated to assess the degree of association between depth and TPM for each COG category. For this analysis, SP21 was removed to allow equitable comparison.

Fig 5

Coverage trends and transcriptional richness.

Changes in transcription for each of the three main COGs (information storage and processing, cellular processing and signaling, and metabolism) at each site (a). For equitable comparison, SP21 transcripts were removed (because of its anomalously-high coverage values) as well as the two libraries representing the bottom 10% of library sizes (2–3 cm and 4–5 cm at site F). Metagenome coverage is not reported for site P, as we did not sequence a metagenome for this site. Changes in transcriptional richness, or the number of transcribed genes, with depth (b). Lines correspond to a subCOG category within a main category, each rendered to share the same color palette: khaki lines represent information storage and processing transcripts, pink lines represent cellular processes and signaling transcripts, and green lines represent metabolism transcripts. Spearman’s rho was calculated for both COG-specific and site-wide trends in transcriptional richness. Only site-wide trends had significant p-values at alpha = 0.05.

Changes in transcription for each of the three main COGs (information storage and processing, cellular processing and signaling, and metabolism) at each site (a). For equitable comparison, SP21 transcripts were removed (because of its anomalously-high coverage values) as well as the two libraries representing the bottom 10% of library sizes (2–3 cm and 4–5 cm at site F). Metagenome coverage is not reported for site P, as we did not sequence a n class="Species">metagenome for this site. Changes in transcriptional richness, or the number of transcribed genes, with depth (b). Lines correspond to a subCOG category within a main category, each rendered to share the same color palette: khaki lines represent information storage and processing transcripts, pink lines represent cellular processes and signaling transcripts, and green lines represent metabolism transcripts. Spearman’s rho was calculated for both COG-specific and site-wide trends in transcriptional richness. Only site-wide trends had significant p-values at alpha = 0.05.

Information storage and processing genes (mostly involved in transcription/translation) had the highest average transcription amongst the three main COG categories (average TPM = 91) compared to cellular processes and signaling (average TPM = 65) and metabolism (average TPM = 60). With the exception of site P, all sites had significant decreases (p < 0.05) in transcripts for all three of these main COG categories (Fig 5a). Metabolism genes had the greatest decrease (Spearman’s rho from -0.5 at site AC to -0.6 at site F), followed by cellular processes and signaling (Spearman’s rho from -0.24 at site AB to -0.58 at site F) and information storage and processing (Spearman’s rho from -0.25 at site AB to -0.47 at site F).

The COG with the largest number of genes that were transcribed was energy production and conversion, followed at most sites by translation, ribosomal structure, and biogenesis (Fig 5b). We tested how richness within these COGs and others changed along a depth gradient through Spearman analysis. This analysis revealed that the number of transcribed genes at site AB decreased overall (Spearman’s rho = -0.26, p value = 0.018), while at sites AC and P, there was a slight increase (Spearman’s rho = 0.30 and 0.41, respectively, p ≤ 0.005). By contrast, site F had no statistically-supported depth trend in transcriptional richness, either when evaluating the entire data set or sub-setting by COG category (Spearman’s rho = 0.11, p value = 0.31).

Combined, these analyses suggest that overall transcription generally decreased, while trends in transcript richness were not uniform across sites. Instead, trends in transcriptional richness broadly intersected with increased SP21 transcription (Figs 4 and 5). For example, at site AB, overall decreases in transcriptional richness coincided with increased SP21 transcription from 0–1 cm to 2–3 cm. At site P, increased richness from 0–1 cm to 2–3 cm was followed by a dramatic decrease from 2–3 cm to 3–4 cm (Fig 5b). This decrease coincided with statistically elevated amounts of SP21 transcription (Fig 4). At site AC, SP21 transcripts remained low at all depths and were accompanied by an increase in total transcript richness. Interestingly, at site F, we observed no change in either transcriptional richness or in SP21 transcription. Altogether, this suggests that mRNA generation may be attenuated once a threshold of SP21 transcription is reached.

Discussion

Woeseiales populations made up 1–3% of the microbial community by 16S rRNA gene amplicon analysis, and had high metagenomic read recruitment and transcriptional activity in the sediments of Kongsfjord and Van Keulenfjord, Svalbard. This is in line with previous studies that showed Woeseiales is an abundant lineage in marine sediments worldwide [7, 12, 13], although the properties that determine such ecological success are poorly understood. Previous work has suggested metabolic plasticity to be at the root of the global distribution and high abundance continuously observed for Woeseiales worldwide [7, 13, 15]. Our high-resolution transcriptional evidence supports this idea and suggests that the Woeseiales in these Arctic fjords are able to sense and adapt to their environment upon burial, most likely in response to changes in redox chemistry. In addition to metabolic plasticity, transcriptional evidence also supports an additional mechanism that is potentially driving their ecological success: the transcription of a small, membrane-associated stress protein, SP21 (hspA) [35, 36].

Metabolic response to burial

In these surface sediments, where redox cycling of carbon, n class="Chemical">nitrogen, and sulfur compounds provide key feedbacks to the overlying water [22, 23, 25], we found metabolic potential for chemolithoautotrophy with sulfur oxidation paired to carbon fixation (Fig 3a). This agrees with previous work that hypothesized Woeseiales as members of the sulfur-oxidizing community in marine sediments [7, 12, 13]. We also found strong transcriptional evidence that Woeseiales contribute to denitrification by performing the NO2 reduction step in the examined sediment surface layers down to 5–6 cm (Fig 3b; S1 File). Transcriptional evidence suggests that these metabolisms are occurring in situ and suggests sulfur oxidation is down-regulated quickly with depth in response to changing redox conditions (Fig 3b). The electron sources driving nitrite reduction need to be examined in the future. Metabolic plasticity within the Woeseiales was first proposed by Dyksma et al. [13]. Using combined molecular and isotopic approaches, they found support for Woeseiales to switch between sulfur and hydrogen oxidation. They hypothesized that the global distribution and ubiquity of the Woeseiales clade may in part be attributable to their metabolic flexibility. Evidence for the potential for metabolic plasticity was supported by an analysis of metabolic features within all available Woeseiales genomes [7]. Likewise, analysis of MAGs from permeable marine sediments in Australia suggest Woeseiales encode genes for aerobic heterotrophy, denitrification, sulfur oxidation, and hydrogenotrophic sulfur reduction [15]. Further, meta-analysis of deep-sea versus coastal Woeseiales taxa found chemoheterotrophic metabolism is a deep-sea trait [12]. Our analysis extends these findings of metabolic flexibility to include Arctic fjord sediments, suggesting that metabolic flexibility likely drives the ecology and distribution of Woeseiales in the marine sediment realm.

Dormancy and quiescence strategies

An additional mechanism that may facilitate the worldwide distribution of Woeseiales is its transcriptional response to burial, suggesting cellular quiescence (Fig 4). The extremely high number of transcripts of the single gene encoding SP21, or hspA (6–161 times higher than the next most abundant transcript; S7 Fig) that increases with depth (1.5 to 9 times increase over the upper 5 cm at three sites) introduces a potentially important and overlooked strategy available to Woeseiales. SP21 is transcribed alongside other chaperones and genes for repair, perhaps as a general response to stressful conditions that occur as Woeseiales is buried. While the transcription of chaperones, repair proteins, and other genes related to stress was observed, transcription did not increase with depth and remained several orders of magnitude lower than transcription for spore protein SP21 on average (S9 Fig). Therefore, the transcription of spore protein SP21 may indicate a quiescence strategy related to stress-induced dormancy [44]. The conditions that trigger this response may not be as pronounced or had as much time to develop at site AC, where SP21 transcript abundance remained steady and diversity of transcripts did not decrease with depth. Swath-bathymetry of Van Keulenfjorden indicates that this part of the fjord contains submarine landform assemblages, such as dewatering structures and pockmarks, that result from rapid sediment deposition and debris flows [21]. These structures are absent at site AB. This suggests that that difn class="Chemical">ferences in the physical sediment environment, such as rapid deposition, at site AC may attenuate the stressors that elicit SP21 transcription in shallow sediment layers.

In Stigmatella aurantiaca, the expression of membrane-localized clusters of SP21 is inducible under stress conditions like heat shock, n class="Disease">anoxia, and nutrient deprivation [35, 45]. When expressed under heat shock, antigens for this stress protein have area densities ranging from ~400 to >800 conjugates/μm2 within the cell wall [35], suggesting that the protein is in extremely high concentrations in the cell. Unlike the sulfate-reducing Firmicutes whose spores are inducible under high temperatures in Svalbard fjord sediment [46], our transcriptional evidence for putative spore-like dormancy in Woeseiales is limited to just one protein located in the cell wall. Thus, it is very unlikely that Woeseiales can form spores sensu stricto as is observed in the Firmicutes, such as Bacillus subtilis. This is supported by the inability of the isolate Woeseia oceani XK5 to form spores [14]. However, the extremely high transcript abundance of SP21 relative to all other genes, and its sharp increase with depth suggest that it might also have a high protein abundance, as it does in S. aurantiaca, where it provides a spore-like quiescence.

Indeed, decreases in abundance and diversity of transcripts for metabolism, information storage and processing, and cellular processes may represent a stress response pattern indicating entry into a stage of quiescence. The transcriptional behavior observed for the Woeseiales MAGs is similar to quiescence strategies noted for other organisms, such as n class="Species">B. subtilis and Mycobacterium tuberculosis [44]. For example, the overall high transcription of a single protein coupled to much lower overall mRNA diversity and abundance agrees with a recent RNA-seq study on B. subtilis spores conducted by Korza et al. [47]. B. subtilis spores had overall low transcriptional diversity and abundance, with the exception of a few high-abundance transcripts that were ~4,000-fold more abundant than the least abundant mRNAs [47]. These high-abundance transcripts encoded proteins involved in stringent response and spore coat formation are thought to remain abundant post-sporulation because of mRNA stability conferred by dehydration upon spore maturation [47, 48]. Stabilization of mRNA is also observed for sporulating Actinobacteria [49, 50] and Mycobacteria entering into quiescence [51]. For most sites, mRNA diversity and abundance were highest in uppermost sediment layers prior to build up of SP21 transcripts. This is consistent with increased macromolecular synthesis observed to occur during entry and exit from quiescence in E. coli and M. tuberculosis [52] in order to facilitate cellular remodeling between growth states [53]. After quiescence is established, transcription has been shown to decrease 3- to 5- fold in Saccharomyces cerevisiae [54]. And so, while Woeseiales may not have the ability to form true spores, it may have a strategy that is similar to cellular quiescence that has been observed for both pathogenic bacteria [55] and environmental bacteria [56].

Altogether, high transcription of SP21, coupled to the transcription of chaperones, transcriptional regulators, membrane-associated stress genes, spore coat protein A (S9 Fig), and peptidases related to sporulation factors (M50B), suggests that Woeseiales has a strategy for cellular quiescence for surviving periods of stress. Such metabolic standby of microbial populations, wherein growth and metabolism are extremely depressed or put on hold entirely until conditions again become optimal, has been proposed as a key strategy of lin class="Chemical">fe in marine sediments [6, 57]. These dormant populations act as seed banks that can potentially be reactivated under different conditions [58], such as those resulting from resuspension events or bioturbation [59]. Such a subsistence strategy may confer resilience to Woeseiales populations, allowing their persistence in the environment and contributing, in part, to their global distribution and abundance.

Conclusions

Our results suggest that the Woeseiales MAGs reconstructed from Arctic fjord sediments responded physiologically to burial over the upper 5 cm of sediments. Transcriptional evidence for metabolic plasticity was found with shallow termination of transcription of rDSR likely involved in n class="Chemical">sulfur oxidation tied to carbon fixation and deeper transcription of the nirS gene for denitrification. A single membrane-associated stress protein, SP21, was highly transcribed at all depths and at all sites. Increases in transcripts for SP21 at three sites coincided with decreased investment in mRNA for metabolism, cellular processes and signaling, and information storage and processing genes with depth. The site where SP21 did not increase with depth did not show decreases in the abundance and diversity of cellular transcripts. Altogether, this suggests that in addition to metabolic flexibility, the success of Woeseiales as a ubiquitous and abundant marine sediment clade may be attributable at least in part to a novel subsistence strategy that involves SP21. With multiple subsistence strategies (metabolic flexibility and cellular quiescence) available to Woeseiales both within and outside the reach of large glacial deposits, they may be well-prepared to withstand shifts in the sub-benthic conditions of Arctic fjords due to climate change.

Materials and methods

Sediment collection

Sediment was collected in the summer of 2016 from four different sites within Svalbard fjords (79°N) with a HAPs corer [60] (Fig 1a). Sediments included in this study are from sites AB (77°35.249’ N, 15°05.121’E) and AC (77°32.260’ N, 15°39.434’ E) in Van Keuelenfjorden (outer and middle sites, respectively; Fig 1b) and sites F (78°55.075’ N, 12°15.929’ E) and P (78°57.915’N, 12°15.600’E) in Kongsfjorden (both located at areas of glacier outflow) (Fig 1c).

Geochemical measurements

Acrylic liners were used for collecting cores from each site only cm apart within each larger HAPs corer. These cores were kept at ~5°C until they were brought back to land and processed 1–3 days later. Cores were subsectioned every 3 cm down to a depth of ~25 cm. Samples for hydrogen analysis followed a previously published headspace equilibration method [61]. Briefly,1 mL of sediment was placed into a dark glass serum vial which was then crimp sealed and gassed with n class="Chemical">N2 for 15 min prior to storage at 4°C. Headspace was measured with glass syringes on a Peak Performer GC with a mercuric chloride detector (San Jose, California) at the University of Tennessee, Knoxville (UTK) after 2 days. Sediment for analysis of organic matter was freeze-dried after thawing from -80°C and subjected to acid fumigation overnight before analysis [62]. Total organic carbon (TOC) and isotope composition of carbon (δ13Corg) from bulk organic matter was measured using a Thermo-Finnigan Delta XL mass spectrometer coupled to an elemental analyzer at UTK. Carbon to nitrogen (C/N) ratios were calculated by dividing percent C by percent N. Isotopic values were calibrated against the USGS40 and USGS41 international standards. In-house standard sets were run every 12 samples. Across multiple runs, one standard deviation was 0.1–0.2‰ for δ13Corg, 1.1–1.8% for mgN, and 1.0–2.2% for mgC. Values for Van Keulenfjorden are published elsewhere [22].

DNA extraction and 16S rRNA gene amplicon libraries

Cores were collected with acrylic liners a few cm apart from those taken for geochemical analysis and were subsectioned at 1 cm depth intervals in the Ny-Ålesund Marine Lab down to ~20 cm depth. Sediment was frozen immediately on dry ice and remained frozen during transport using a dry shipper. Sediment was stored at -80°C until processing. Nucleic acids were extracted by the Lloyd Lab (sites AB, AC, P, F) using the Qiagen RNeasy Powersoil® Kit (Hilden, Germany) for RNA with the DNA extraction accessory according to manufacturer’s instructions. DNA from site AC at 18 cm depth was extracted by the Loy lab using the DNeasy PowerSoil Kit (Qiagen), following manufacturer’s instructions. 16S rRNA gene amplification from sites P and F in Kongsfjorden failed. The Phusion Master Mix (Thermo Fisher) was used with the primer set 515F/806R [63] at the Center for n class="Species">Environmental Biology at UTK for amplification. Reads were sequenced with Illumina MiSeq and trimmed for quality with Trimmomatic [64] using a window 10 base pairs wide and a minimum phred score of 28. Trimmed reads were then processed in mothur 1.35.1 [65] using the computational cluster at the Bioinformatics Resource Facility (BRC) at UTK. Operational taxonomic units (OTUs) were clustered with mothur as previously described [22]. All 16S rRNA gene amplicon library sequences can be accessed at NCBI under Bioproject PRJNA493859.

Metatranscriptomic and metagenomic sequencing

RNA extracts were treated with DNAse I at UTK (Qiagen) and further DNase treatment was performed at MRDNA (Shallowater, TX), followed by sequencing of metatranscriptomic libraries with Illumina HiSeq, PE 2x250 bp. Individual 1 cm-depth resolved metatranscriptomic libraries were generated with RNA extracts from the first 5 to 6 cm of sediment from sites AB, AC, F, and P, for a total of 20 metatranscriptomes. Metagenomic libraries were generated from the combined extracts from the first 5 cm (spanning 0 to 5 cm downcore) in sites AB and F in the Center for n class="Species">Environmental Biotechnology, Knoxville, using Illumina HiSeq, PE 2x250 bp. The decision to make high depth resolution metatranscriptomes, and only a single metagenome for each site was based on financial constraints. The Loy Lab sequenced DNA extracts from Site AC at 18 cm depth using Illumina HiSeq 3000, PE 2x150 by the Biomedical Sequencing Facility, at The University of Vienna.

Prior to metagenomic assembly, raw reads were trimmed for quality and adapters were removed using in-house scripts in the Loy lab (removing the leading eight 5’ bases, bases with QC < 15 and reads below 50 bp in length) and Illumina adapters were removed using AdapterRemoval using default settings [66]. The Lloyd Lab used Trimmomatic [64] for trimming both metagenomes and metatranscriptomes, with a sliding window of 10 and a phred cut off score of 28 for all reads above 90 bp. The quality of trimmed reads was assessed with Quast 4.5 [67]. Metagenomic sequences can be accessed in the NCBI SRA under BioSamples SAMN10372305, SAMN10372307, and SAMN10372306.

Metagenomic assembly: MetaSPAdes

Metagenomic assembly was accomplished using metaSPAdes [68] both in command line for Linux and in KBase, a browser interface with bioinformatics modules and applications [69]. Site F was the only metagenome that was assembled via metaSPAdes version 3.11 in the command line, with kmer size set to 21, 33, 55, 77, 99, and 127 to find the best assembly. All other metaSPAdes assemblies were completed on KBase with the default parameters (1000 bp length and kmer sizes of 21, 33, and 55). Assembled contigs were then filtered to contain only contigs with more than 5x coverage and 1000 bp length using in-house scripts.

Metagenomic assembly: IDBA and Megahit

To reduce RAM utilization and wall clock time, larger sequence datasets were normalized to 100X with BBN prior to assembly with IDBA version 1.1.3 [70] and Megahit version 1.0.5 [71]. Assemblies were completed in either the command line or on KBase [69]. All assemblies were generated with 1000 minimum contig length, except IDBA assembly in the Loy lab, which used a 500 bp cutoff.

Taxonomic binning of contigs into metagenome assembled genome (MAGs)

Contig binning was carried out with MaxBin2 version 2.2.3 [72], CONCOCT (version 0.4.1) [73] and MetaBAT (version 2.12.1) [74] in command line. Each of these binning tools utilizes genomic signatures within contigs, such as coverage and n class="Chemical">tetranucleotide frequency, to identify discrete clusters of contigs that likely represent a population’s genome. MaxBin2 and CONCOCT binning were performed with default parameters and 1000 bp minimum contig length. MetaBat binning was conducted with an interactive pipeline that applies decreasing levels of stringency within each successive iteration, collecting the best bins and their contigs each time. The remaining contigs are then placed into the next step with a lower stringency.

DASTool (version 1.1.0) [75] was applied to each set of binned contigs to create medium to high quality MAGs. This tool takes several difn class="Chemical">ferent MAGs as input and identifies a consensus, non-redundant genome for each MAG, leading to higher quality genomes. Quality was determined for our MAGs according to Bowers et al. [76]. Medium quality drafts have completeness and contamination values of ≥ 50% and <10%, respectively whereas high quality MAGs have completeness and contamination values >90% and contamination <5%. Completeness and contamination were determined with CheckM [77]. MAG taxonomy was determined through phylogenetic analysis using a concatenated alignment of proteins from single copy marker genes included in the CheckM suite and a tree was built with FastTree [78].

Transcript mapping and annotation

Ribosomal sequences were identified and removed using SortmeRNA version 2.1 [79]. Filtered transcripts were mapped to MAGs using Bowtie2 version 2.3.4 [80] with sensitive local mapping, which requires neither end-to-end alignment nor mismatches when evaluating seed substrings 20bp in length. Genomes were annotated with Prokka version 1.12 [81]. Mapping files, Prokka gene calls, and fasta files for each MAG were then used in the metagenomics pipeline of Anvi’o version 5.1 [82] for COG identification, visualization, and transcript abundance estimates. Transcript abundance is the coverage of the gene normalized by gene length. Transcripts per million was calculated by dividing transcript abundance by the number of mapped reads for each library and multiplying by 1,000,000.

COG annotations were computed using eggNOG-mapper [83] based on eggNOG 4.5 ontology data [84]. Spearman correlation coefficient analysis of changes in transcriptional richness associated with depth was performed in R [85]. Intersection analysis was performed with the UpsetR package in R [86]. TPM values for each hspA gene copy were log-transformed and normality was tested with a shapiro-wilks test in R [86]. ANVOA (alpha = 0.05) was performed in R to facilitate post-hoc Tukey testing with the agricolae package (alpha = 0.05). Co-localization analysis was performed in CLC Genomics Workbench.

Amino acid identity analysis for each Woeseiales MAG was conducted with AAI-profiler [87]. This tool uses SANSparallel [88] to find homology matches within the Uniprot database (update Apr2020).

RubisCO, rDSR, and hspA sequence analyses

Publicly available full length (> 480 aa) sequences for large and small chains of ribulose bisphosphate carboxylase (RubisCO) were compared with MAG RubisCO sequences. Mega version 7 [89] was used for Clustal alignments and neighbor-joining tree construction.

The alpha subunits of Woeseiales rDSR sequences were aligned in the dsrAB ARB database [90]. RAxML was used within the CIPRES Science Gateway v. 3.3 server [91] and the Interactive Tree of Life (iToL) server was used for tree visualization [92]. For spore protein SP21, publicly available amino acid sequences for hspA were acquired from n class="Species">environmental samples and isolates and compared with Woeseiales sequences. RAxML was used within the CIPRES Science Gateway v. 3.3 server [91] and the Interactive Tree of Life (iToL) server was used for tree visualization [92].

(XLSX)

Click here for additional data file.

Genome statistics for the Woeseiales MAGs in this study.

Completeness and contamination were determined for each genome with CheckM.

(DOCX)

Click here for additional data file.

Genes with transcripts detected at all sites.

C = Energy production and conversion; D = Cell cycle control, cell division, chromosome partitioning; E = Amino acid transport and metabolism; F = Nucleotide transport and metabolism; G = Carbohydrate transport and metabolism; H = Coenzyme transport and metabolism; I = Lipid transport and metabolism; J = Translation, ribosomal structure and biogenesis; K = Transcription; L = Replication, recombination and repair; N = Cell motility; O = Post-translational modification, protein turnover and chaperones; P = Inorganic ion transport and metabolism; Q = Secondary metabolites biosynthesis, transport and catabolism; S = Function unknown; T = Signal transduction mechanisms; U = Intracellular trafficking, secretion and vesicular transport; V = Defense mechanisms.

(DOCX)

Click here for additional data file.

Phylogenetic tree for Woesiales MAGs placement.

Maximum likelihood was calculated in Mega v. 7 with 1000 bootstraps on a concatenated alignment of 49 ribosomal proteins. Only nodes with >75% support are shown. MAGs from this study are indicated with blue circles. Tree was rooted using Roseobacter species as the outgroup.

(TIF)

Click here for additional data file.

Comparison of genomic contents across the 5 Woeseiales MAGs to show similarity.

Figure generated with KEGG-decoder [26].

(TIF)

Click here for additional data file.

Amino acid identity analysis.

Scatterplots, Krona diagrams, and histograms relating amino acid identity to the cultured isolate, Woeseia oceani, are presented for the n class="Chemical">MAGs Woeseia_stnAB (a, b, c), Woeseia_stnAC (d, e, f), Woeseia2_stnAC (g, h, i), Woeseia_stnF (j, k, l), and Woeseia2_stnF (m, n, o).

(TIF)

Click here for additional data file.

Intersection analysis of all transcripts.

Set intersections of transcripts visualized with UpsetR [86]. The total number transcribed genes (“set size”) is indicated as horizontal bars. Vertical bars represent the intersection size, or number of transcribed genes that fall into the 15 different sets of site combinations. Depth analysis of transcript abundance for unique sets of transcripts are depicted in S5 Fig.

(TIF)

Click here for additional data file.

Recruitment of unique transcripts.

Transcriptional coverage (transcript abundance normalized to gene length) is reported for the unique transcripts detected at each site for the COGs Cellular Processing and Signaling, Information Storage and Processing, Metabolism, and Function Unknown. Those with predicted function only were removed from this analysis. Dot size scales with metagenome coverage and color saturation scales from lighter to darker blue for increased transcriptional coverage. Site P did not have a n class="Species">metagenome.

(TIF)

Click here for additional data file.

Neighbor-joining tree of RubisCO sequences.

Only nodes with >80% support after 1000 bootstraps are shown.

(TIF)

Click here for additional data file.

Coverage of all Woeseiales MAG transcripts.

The transcripts per million (TPM) value is reported for each transcribed gene (along x-axis). Names of genes are not included to ease illustration and interpretation. COG assignment of genes is represented by color. The gene encoding for spore protein SP21 is indicated with arrows.

(TIF)

Click here for additional data file.

Maximum likelihood tree of hspA genes encoding for spore protein SP21.

Bootstrap support after 1000 bootstraps is reported for nodes with > 70% support. Woeseiales sequences are indicated in red, cultured microorganisms in purple, and experimentally-verified spore protein SP21 indicated in green [36, 45]. Archaeal sequences of spore protein SP21 were used as the outgroup. The scale bar represents the number of substitutions per site.

(TIF)

Click here for additional data file.

Coverage with depth for transcripts of stress proteins, chaperones, and transcriptional regulators.

TPM values reported for genes associated with starvation response, stress mitigation, protein repair/folding, and transcriptional regulation. Yellow border labels indicate genes that co-occur on the same contigs as hspA (encoding spore protein SP21; S9 Fig).

(TIF)

Click here for additional data file.

Co-localization of genes on the same contig as the gene for spore protein SP21 (hspA), highlighted with orange rectangle.

Visualizations from CLC Genomics Workbench.

(TIF)

Click here for additional data file.

9 Jul 2020

PONE-D-20-16591

Woeseiales transcriptional response in Arctic fjord surface sediment

PLOS ONE

Dear Dr. Buongiorno,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.

Please submit your revised manuscript by Aug 23 2020 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at plosone@plos.org. When you're ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

Please include the following items when submitting your revised manuscript:

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols

We look forward to receiving your revised manuscript.

Kind regards,

Peter R. Girguis, Ph.D.

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and

https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. We note that Figure 1 in your submission contain map images which may be copyrighted. All PLOS content is published under the Creative Commons Attribution License (CC BY 4.0), which means that the manuscript, images, and Supporting Information files will be freely available online, and any third party is permitted to access, download, copy, distribute, and use these materials in any way, even commercially, with proper attribution. For these reasons, we cannot publish previously copyrighted maps or satellite images created using proprietary data, such as Google software (Google Maps, Street View, and Earth). For more information, see our copyright guidelines: http://journals.plos.org/plosone/s/licenses-and-copyright.

We require you to either (1) present written permission from the copyright holder to publish these figures specifically under the CC BY 4.0 license, or (2) remove the figures from your submission:

2.1. You may seek permission from the original copyright holder of Figure 1 to publish the content specifically under the CC BY 4.0 license.

We recommend that you contact the original copyright holder with the Content Permission Form (http://journals.plos.org/plosone/s/file?id=7c09/content-permission-form.pdf) and the following text:

“I request permission for the open-access journal PLOS ONE to publish XXX under the Creative Commons Attribution License (CCAL) CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). Please be aware that this license allows unrestricted use and distribution, even commercially, by third parties. Please reply and provide explicit written permission to publish XXX under a CC BY license and complete the attached form.”

Please upload the completed Content Permission Form or other proof of granted permissions as an "Other" file with your submission.

In the figure caption of the copyrighted figure, please include the following text: “Reprinted from [ref] under a CC BY license, with permission from [name of publisher], original copyright [original copyright year].”

2.2. If you are unable to obtain permission from the original copyright holder to publish these figures under the CC BY 4.0 license or if the copyright holder’s requirements are incompatible with the CC BY 4.0 license, please either i) remove the figure or ii) supply a replacement figure that complies with the CC BY 4.0 license. Please check copyright information on all replacement figures and update the figure caption with source information. If applicable, please specify in the figure caption text when a figure is similar but not identical to the original image and is therefore for illustrative purposes only.

The following resources for replacing copyrighted map figures may be helpful:

USGS National Map Viewer (public domain): http://viewer.nationalmap.gov/viewer/

The Gateway to Astronaut Photography of Earth (public domain): http://eol.jsc.nasa.gov/sseop/clickmap/

Maps at the CIA (public domain): https://www.cia.gov/library/publications/the-world-factbook/index.html and https://www.cia.gov/library/publications/cia-maps-publications/index.html

NASA Earth Observatory (public domain): http://earthobservatory.nasa.gov/

Landsat: http://landsat.visibleearth.nasa.gov/

USGS EROS (Earth Resources Observatory and Science (EROS) Center) (public domain): http://eros.usgs.gov/#

Natural Earth (public domain): http://www.naturalearthdata.com/

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. n class="Species">participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: This a well designed, performed, and written study tracking the genomic diversity and metabolic activity of bacteria through fjord sediments. I have no issues with it in its current form and think it should be published.

Reviewer #2: This manuscript examines five novel metagenome assembled genomes (n class="Chemical">MAGs) from the Order Woeseiales, a bacterial clade which is abundant in marine sediments, coupled with depth-resolved metatranscriptomic data, from four sampling sites in Arctic fjord surface sediment. The authors report that these MAGs had the capability for chemoautotrophic carbon fixation via the Calvin-Benson-Bassham pathway, coupled to the oxidation of sulfite or thiosulfate, as well as the ability to reduce nitrite to nitric oxide. At three of the four sampling sites, transcription of most of the central metabolism genes dropped with sediment depth, while the transcription of a single gene, spore protein SP21, rose sharply. The authors propose that the Woeseiales in this environment may enter a quiescent state with burial, allowing these microbes to ride out stressful conditions by entering dormancy, only to be revived when improved conditions present themselves (if, for example, they are brought to the surface of the sediment through bioturbation).

While it is risky to extrapolate data from five MAGs to describe the metabolic capabilities of an entire clade, the power of the data presented in this manuscript lies in the coupling of the data to the metatranscriptomes. In this way the authors demonstrate how the Woeseiales found in their samples might survive changing conditions with sediment burial.

The science presented here is well-executed and the manuscript well-written. The methods are logical, well-explained, and justified. I found very little to criticize, but I do have a few minor points:

The authors state towards the end of the discussion that SP21 increased with depth at three of the four sampling sites where they were able to generate metagenomes/metatranscriptomes, and that at the one site where they did not observe this, they also did not observe a concurrent decrease in the abundance and diversity of cellular transcripts. (I assume from the figures that the outlier site is site AC, although this is not stated in the text.) Do the authors have any insight into why this site may be difn class="Chemical">ferent from the others? Could there be some geochemical parameter at play here? There may of course be no obvious explanation, but a comment or two about what may be driving this difference would be useful.

Line 454 (methods): "extracted by both the Lloyd" is a phrase missing here? Do you mean "both the Lloyd and the Loy labs"?

Figure 5: Some of the depths appear to be missing from sites AB and AC (unless I misinterpreted something in the text). The figure legend mentions that certain libraries from site F were not included- at first it was difficult to tell that points were missing from the line graph at the bottom of the figure, but I can see that those are not there. What happened with the data from sites AB and AC?

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: Yes: Lauren M Seyler

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

15 Jul 2020

We have re-made Figure 1 to include only images from NASA Earth Observatory. We have changed the figure caption to reflect the proper attribution of image data source.

Reviewer comments:

Reviewer #1: This a well designed, performed, and written study tracking the genomic diversity and metabolic activity of bacteria through fjord sediments. I have no issues with it in its current form and think it should be published.

The authors thank this reviewer for taking the time to reviewer our work and for their positive review.

Reviewer #2: This manuscript examines five novel metagenome assembled genomes (n class="Chemical">MAGs) from the Order Woeseiales, a bacterial clade which is abundant in marine sediments, coupled with depth-resolved metatranscriptomic data, from four sampling sites in Arctic fjord surface sediment. The authors report that these MAGs had the capability for chemoautotrophic carbon fixation via the Calvin-Benson-Bassham pathway, coupled to the oxidation of sulfite or thiosulfate, as well as the ability to reduce nitrite to nitric oxide. At three of the four sampling sites, transcription of most of the central metabolism genes dropped with sediment depth, while the transcription of a single gene, spore protein SP21, rose sharply. The authors propose that the Woeseiales in this environment may enter a quiescent state with burial, allowing these microbes to ride out stressful conditions by entering dormancy, only to be revived when improved conditions present themselves (if, for example, they are brought to the surface of the sediment through bioturbation).

While it is risky to extrapolate data from five MAGs to describe the metabolic capabilities of an entire clade, the power of the data presented in this manuscript lies in the coupling of the data to the metatranscriptomes. In this way the authors demonstrate how the Woeseiales found in their samples might survive changing conditions with sediment burial.

The science presented here is well-executed and the manuscript well-written. The methods are logical, well-explained, and justified. I found very little to criticize, but I do have a few minor points:

The authors state towards the end of the discussion that SP21 increased with depth at three of the four sampling sites where they were able to generate metagenomes/metatranscriptomes, and that at the one site where they did not observe this, they also did not observe a concurrent decrease in the abundance and diversity of cellular transcripts. (I assume from the figures that the outlier site is site AC, although this is not stated in the text.) Do the authors have any insight into why this site may be difn class="Chemical">ferent from the others? Could there be some geochemical parameter at play here? There may of course be no obvious explanation, but a comment or two about what may be driving this difference would be useful.

Thank you for inviting us to speculate about this. We have added text (L376-383) that addresses a potential cause for the differences observed at site AC.

Line 454 (methods): "extracted by both the Lloyd" is a phrase missing here? Do you mean "both the Lloyd and the Loy labs"?

We have fixed this in the text. (The Lloyd Lab extracted DNA and RNA using the RNeasy kit with DNA accessory kit and the Loy lab extracted DNA using the DNeasy kit).

Figure 5: Some of the depths appear to be missing from sites AB and AC (unless I misinterpreted something in the text). The figure legend mentions that certain libraries from site F were not included- at first it was difficult to tell that points were missing from the line graph at the bottom of the figure, but I can see that those are not there. What happened with the data from sites AB and AC?

Thanks for checking this. All libraries are reported in the original figure, it just looks as though some are missing as you suggest because other sites had successful metatranscriptomes at different depths.

Site AB had four metatranscriptomes (0.5, 1.5, 2.5, 3.5 cmbsf)

Site AC had four metatranscriptomes (1.5, 2.5, 3.5 4.5 cmbsf)

Submitted filename: Response to Reviewers.docx

Click here for additional data file.

5 Aug 2020

Woeseiales transcriptional response to shallow burial in Arctic fjord surface sediment

PONE-D-20-16591R1

Dear Dr. Buongiorno,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

Peter R. Girguis, Ph.D.

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

Reviewers' comments:

11 Aug 2020

PONE-D-20-16591R1

Woeseiales transcriptional response to shallow burial in Arctic fjord surface sediment

Dear Dr. Buongiorno:

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

If we can help with anything else, please email us at plosone@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Peter R. Girguis

Academic Editor

PLOS ONE