ORIGINAL RESEARCH article

Front. Mar. Sci., 10 January 2020

Sec. Aquatic Microbiology

Volume 6 - 2019 | https://doi.org/10.3389/fmars.2019.00797

A Unique Benthic Microbial Community Underlying the Phaeocystis antarctica-Dominated Amundsen Sea Polynya, Antarctica: A Proxy for Assessing the Impact of Global Changes

  • 1. Department of Marine Sciences and Convergent Technology, Hanyang University, Seoul, South Korea

  • 2. School of Earth and Environmental Sciences and Research Institute of Oceanography, Seoul National University, Seoul, South Korea

  • 3. Department of Microbiology, Chungbuk National University, Cheongju, South Korea

  • 4. Department of Biological Sciences, Eastern Illinois University, Charleston, IL, United States

  • 5. Max Planck Institute for Marine Microbiology, Bremen, Germany

  • 6. Korea Institute of Ocean Science & Technology, Busan, South Korea

  • 7. Korea Polar Research Institute, Incheon, South Korea

Abstract

Polynyas in the polar seas are regarded as windows through which ecosystem responses associated with global climate changes are to be noticed. However, little information is available on benthic microbial communities in the Amundsen Sea polynya (ASP), where environmental changes due to global warming are occurring rapidly, from which future climate change-induced ecosystem responses could be assessed. We performed high-throughput sequencing of 16S rRNA genes and quantitative PCR in combination with biogeochemical analyses and metabolic rate measurements to determine the composition, diversity and controls of major microbial communities in sediments of the ASP. A large fraction of the sequenced benthic microbial community (40% on average) in the polynya was uniquely affiliated with the phylum Planctomycetes, whereas Thaumarchaeota (51% on average) predominated in non-polynya areas. The relative abundance of Planctomycetes correlated significantly with organic carbon (Corg) content in the polynya sediment underlying the Phaeocystis-dominated water column. These results suggest that Planctomycetes comprise a major bacterial group utilizing relatively recalcitrant Corg produced primarily by Phaeocystis blooms. In contrast, the predominance of chemolithoautotrohic Thaumarchaeota in the sea-ice zone was attributed to low Corg supply due to low primary productivity in the ice-covered water column. The Planctomycetes-dominated microbial communities in the ASP is in stark contrast to that Proteobacteria (Delta- and Gamma-proteobacteria) occupy ecological niches as primary mineralizers of organic materials in most benthic systems in the Southern Ocean, where organic materials in the sediments mostly originate from diatom blooms. Given that microbial communities respond quickly to environmental changes, and that global climate change is proceeding rapidly in the ASP, our results suggest that any modifications in the Planctomycetes-dominated microbial communities will provide valuable insight into changes in organic matter transport to the seafloor that may result from shifts in planktonic communities from Phaeocystis to diatoms associated with climate change.

Introduction

In marine surface sediments, microorganisms occur in enormous numbers (approximately 1.7 × 1028 cells worldwide) (), are extremely diverse (), and are notably involved in global carbon, nitrogen, and sulfur cycles (; ). Because microbial community composition, diversity, and metabolic activity are significantly influenced by environmental change (; ; ; ; ), characterization of microbial distribution provides relevant information on spatial and temporal variations in environmental conditions (; ; ). However, little is known about benthic microbial communities in the Southern Ocean (SO), where environmental changes due to global warming are occurring rapidly (; ; ).

The SO plays a profound role in regulating global carbon cycles, accounting for approximately 20–30% of global ocean CO2 uptake (; ). As atmospheric pCO2 increases (), the role of the SO as an atmospheric carbon sink has received more attention. During the austral summer, the coastal zone of the SO is typically characterized by the occurrence of polynyas, areas of seasonally recurring open water surrounded by sea ice (; ). Because of the combined effects of enhanced light conditions and iron supply resulting from melting sea ice, polynyas are one of the most productive marine ecosystems (; ; ), and thus are regarded as significant sinks for atmospheric CO2 (; ). In addition, due to the significant gas exchange (CO2, DMS, etc.) and heat flux between the sea surface and atmosphere, polynyas are regarded as windows through which we can anticipate ecosystem responses associated with climate changes ().

Among the 37 known coastal polynyas around Antarctica, the Amundsen Sea polynya (ASP) is reported to be the most productive (, ), with a primary productivity per unit area of ∼220 g C m–2 y–1 (; ). The prymnesiophyte Phaeocystis antarctica is the main contributor to the high primary productivity in the central polynya (Figure 1B) (; ; ; ; ). Due to the inflow of warm circumpolar deep water (CDW), glaciers near the Amundsen Sea (AS) are undergoing the highest rates of melting and thinning on the Antarctic continent (; ; ). Consequently, as global warming continues, the heat flux intensity of the CDW may stimulate ice melting, which results in changes in phytoplankton productivity and community composition (), thereby regulating the function of the ASP in carbon sequestration (; ). The composition and metabolic activities of benthic microbial communities are ultimately determined by the quality and quantity of the organic matter supplied from the overlying water column (; ; ; ). Therefore, given that microbial communities respond quickly to environmental change (; ), quantitative and qualitative information on benthic microbial communities and their metabolic activities is essential for assessing the response of these ecosystems to climate change-induced variations in primary productivity in the water column of the SO. The benthic microbial community has been studied in several regions of the Southern Ocean, including the Mertz Glacier polynyas, Ross Sea, Bellingshausen Sea, and Australian-Antarctic ridge (; ; ; ). In those studies, Proteobacteria were reported as a predominant bacterial group in the sediment and are regarded as major organic carbon (Corg) oxidizers. However, there is no information on the composition and diversity of entire microbial communities in the sediment of the ASP.

FIGURE 1

The objectives of this paper were: (1) to identify microbial communities in the sediment underlying the Phaeocystis-dominated water column, (2) to elucidate the major factors controlling microbial communities, with a special emphasis on the Corg content in sediment across the marginal ice zone (MIZ), polynya, and ice shelf zone, and (3) to suggest an appropriate microbial proxy for assessing the response of benthic ecosystem to ongoing climate changes in the Amundsen Sea.

Materials and Methods

Study Area

The Amundsen Sea is located in western Antarctica between the Ross Sea and Bellingshausen Sea (69°S–74°S; 100°W–135°W, Figure 1), and is characterized by a large polynya from November to February (). The Korean Amundsen Sea Expedition was conducted during the austral summer, from February 9 to March 10 in 2012, aboard the Korean icebreaker research vessel RV Araon. Water depth ranged from 530 to 1,064 m, and bottom water temperature ranged from −1.8 to −1.1°C (Table 1).

TABLE 1

Polynya
Ice shelf zone
Sea-ice zone
StationStn 10Stn 17Stn 19Stn 83
Latitude73.250° S73.496° S74.202° S71.699° S
Longitude114.997° W114.008° W112.51° W114.037° W
Water depth (m)8257301064530
Temperature (°C)a−1.1−1.2−1.5−1.8
Salinity (psu)a33.533.433.633.6
TOC (%, dry wt.) (0–1 cm)1.030.710.620.41
TN (%, dry wt.) (0–1 cm)0.150.100.080.07
OPD (cm)1.8 ± 0.12.0 ± 0.23.6 ± 0.13.5 ± 0.3
Total oxygen uptake (mmol O2 m–2 d–1)2.443.111.581.57
Denitrification (nmol N cm–3 sed. d–1)b1.44–4.320.96–7.20n.d.n.d.
Sulfate reduction rate (mmol S m–2 d–1)c0.070.060.050.04
Anammox (nmol N cm–3 sed. d–1)b3.12–3.843.84–6.24n.d.n.d.

Oceanographic parameters and sediment properties, February 18–March 07, 2012.

OPD, oxygen penetration depth; n.d., not detected. aBottom water; bPotential N2 production rate within 0–6 cm depth; cDepth-integrated inventories of sulfate reduction down to 10 cm.

Sediment samples were collected using a box corer at four stations at three contrasting sites: the MIZ (Stn 83), inside the polynya (Stns 10 and 17), and in direct proximity to the ice shelf site (Stn 19) (Table 1). P. antarctica was the major planktonic algae in this highly productive polynya area, whereas diatoms were more abundant in the relatively less-productive MIZ (Stn 83) (Figure 1B) (). Subsamples for DNA extraction were taken from the center portion of the box corer using acryl sub-core liners (6 cm in diameter). Cores were sliced at 1-cm intervals to a depth of 11 or 18 cm, and immediately frozen at −80°C.

DNA Extraction, Quantitative PCR, and Pyrosequencing of 16S rRNA Genes

Total genomic DNA was extracted from the different sediment layers using a PowerMax DNA Isolation kit (Mo Bio Laboratories, Carlsbad, CA, United States), following the manufacturer’s instructions. Quantitative real-time polymerase chain reaction (q-PCR) was used to determine the copy number of 16S rRNA genes and archaeal ammonia monooxygenase A (amoA) genes using a TaqMan assay and a SYBR Green I assay, respectively (Supplementary Table S1). Details are described in Supplementary Material.

For each of the 48 sampled sediment layers at the polynya and non-polynya sites, we generated a 16S rRNA gene amplicon library, with one primer set covering the V5–V8 region of both bacterial and archaeal taxa (). PCR amplification of the 16S rRNA genes was performed in triplicate using a primer set of Uni787F () and Uni1391R () according to . Resulting amplicons were sequenced by Macrogen Corporation (Korea) using the 454 GS FLX + system (Roche). The raw data have been deposited to the NCBI SRA database under accession number SRX3405376. Raw flowgrams of pyrosequencing reads were filtered and de-noised by a PyroNoise algorithm () implemented in MOTHUR (ver. 1.36.1) (). Singleton OTUs were removed prior to analysis. Chimeric sequences were identified and removed by ChimeraSlayer. A total of 132,914 pyrosequencing reads from 48 samples were qualified for further processing. The sequences were then clustered into operational taxonomic units (OTUs) that met the criteria of a 97% similarity threshold and a minimum cluster size of 2 using a QIIME pipeline (ver. 1.9.1) (). Taxonomy for each OTU0.97 was assigned using the RDP classifier method () with the Greengenes database (ver. 13_8) (). To avoid the effects of different sample sizes for estimating diversity, comparison sequences were randomly subsampled to the smallest library size (), which were 1,628 sequences in the present study. Chao1 estimates were created using QIIME software to assess diversity ().

Geochemical Characterization

Geochemical constituents (, NOX [ + ], , and Fe2+) in the pore water, oxygen penetration depth (OPD), total oxygen uptake (TOU), sulfate reduction rates, and N2 removal rates by denitrification and anaerobic ammonia oxidation (anammox) were adopted from the results reported by and . Total organic carbon (TOC) and total nitrogen (TN) were analyzed with Elemental Analyzer (Carlo Erba, NA-1500) () (see Supplementary Material for detailed methods).

Statistical Analyses

Spatial difference of geochemical constituents was assessed using a Mann–Whitney U-test, with a statistical significance level of 0.05. The relative abundance of microbial populations was determined using a Kruskal–Wallis test. Tukey’s honest significance difference (HSD) test was used for multiple comparison of means at a 95% confidence interval. The OTU table from QIIME and geochemical measurements were analyzed using R (v. 3.3.2) () with custom scripts and several packages including vegan (v. 2.4-2) and lmtest (v. 0.9-35). Exploratory data analysis was carried out for both microbial communities and geochemical measurements data using non-metric multidimensional scaling (NMDS), diversity measures, and hierarchical clustering. Microbial community was analyzed at the OTU0.97 level as well as at the phylum and order levels. Compositional difference of microbial communities among stations was tested by multivariate analysis of variance–like non-parametric tests (analysis of similarities [ANOSIM] and permutational multivariate analysis of variance [PERMANOVA]). Ordination of microbial communities was fitted with geochemical measurements by vector fitting. Constrained ordination models were constructed by redundancy analysis (RDA) in an iterative fashion considering collinearity among constraining geochemical variables. Community structures were compared among stations, and between bacteria and archaea using a Procrustes test on RDA ordination configuration and a Mantel test.

Results

Environmental Parameters and Microbial Abundance

The distribution of geochemical constituents and community structures showed distinct spatial variation between the polynya sites and non-polynya sites of the Amundsen Sea. Contents of TOC and TN were approximately 1.25 to 2.1 times higher at the polynya than at the ice shelf (Stn 19) and sea-ice zone (Stn 83) (Figure 2 and Table 1) (P < 0.001). Accordingly, OPD in the polynya sites (1.8–2.0 cm) was shallower than in non-polynya sites (3.5–3.6 cm) (P < 0.001). Pore-water analysis revealed that the concentration of NOX at the polynya sites decreased with depth from approximately 30 μM at the top to 8 μM at 3–5 cm depth, and then remained constant down to a depth of 10–20 cm (Figure 2). In contrast, NOX concentration at the non-polynya sites (ice shelf and sea-ice zone) was high (>20 μM) at all depths. concentration was higher at the polynya sites than at the non-polynya sites (Figure 2) (P < 0.001). Concentration of Fe2+ in the pore-water was low at all sites (<10 μM), but the average concentration of dissolved Fe2+ was higher at the polynya sites (4.9 μM) than at the non-polynya sites (1.98 μM) (P = 0.003). Metabolic activities such as TOU rate and anaerobic respiration by sulfate reduction were consistently higher in the polynya sites with relatively higher Corg content compared to those measured in non-polynya sites (Figure 2 and Table 1) (P = 0.002).

FIGURE 2

Total prokaryotic abundance determined by q-PCR of 16S rRNA genes ranged from 2.3 × 105 to 3.3 × 108 copies cm–3 per each sample (Table 2). Although the 16S rRNA gene copy numbers appeared higher at Stn 17 than those at Stn 10, the variation between polynya sites was not significant (P = 0.06). However, 16S rRNA gene copies in the sediment of Stn 17 were clearly higher than those of Stn 83 (P < 0.05). Total prokaryotic abundances were higher in the surface sediments and decreased with depth at all sites (Figure 2). The relative abundance of bacterial 16S rRNA gene copies to total prokaryotic cells was highest at Stn 10 (96.7 ± 3.3%), and then decreased to 84.8 ± 3.7% (Stn 17), 77.6 ± 7.0% (Stn 19) and 70.3 ± 9.7% on average (Stn 83). In contrast, archaeal 16S rRNA gene proportion comprised 3.3 ± 0.7% at Stn 10, and then gradually increased with distance from Stn 10 to 15.2 ± 3.7% (Stn 17), 22.4 ± 7.0% (Stn 19), and 29.7 ± 9.7% on average (Stn 83) (Table 2). The percentage of bacterial and archaeal 16S rRNA gene copies of the total 16S rRNA genes determined by q-PCR (Table 2) corresponded with the relative abundance of Bacteria versus Archaea as estimated from pyrosequencing data (Figure 2). On the other hand, copy numbers of archaeal amoA gene ranged from 2.8 × 104 to 5.7 × 106 and from 1.7 × 104 to 8.3 × 106 cm–3 for polynya sites and non-polynya sites, respectively. The depth profiles of amoA gene copy numbers in the ASP sediments showed a similar trend to the archaeal 16S rRNA gene copy numbers (Supplementary Figure S1).

TABLE 2

Depth16S rRNA genesNo. of totalNo. ofGood’s
Station(cm)(copies cm–3)§readsOTUscoverageChao1
PolynyaStn 100–1239780566.41239
1–2320595870.11317
3–42580876661359
4–528941014651640
5–6231673568.31200
6–73041118960.91816
7–830491159621801
8–92973107663.81561
10–12185780956.41557
12–14171370958.61392
14–16241495960.31622
16–18315495969.61348
Stn 170–1315686172.71322
1–2258976770.41314
2–3274371673.91341
3–4264498762.71666
4–5162870556.71415
5–63526118666.41793
6–72844107662.21784
7–82706108060.11909
8–9237879266.71418
9–10239887263.61418
10–123103106765.61601
Ice-shelfStn 190–1253084966.41361
1–2300283172.31424
2–35309115378.31457
3–4326977376.41326
4–53154977691617
5–64887114276.61560
6–7305097468.11609
7–8238760874.5991
8–9230681864.51462
9–10251671771.51239
10–12296174374.91126
12–14208956073.2931
14–16334657782.8793
Sea-iceStn 830–1340280376.41182
1–2348068380.41125
3–4255375070.61243
4–5239773169.51229
5–6218160272.41194
6–7233373468.51179
7–8185962066.61277
8–9219550377.1779
9–10212150576.2735
10–12258764375.1963
12–14261850880.6781
14–16307446784.8704

Prokaryotic abundance estimated by q-PCR and estimates of the richness and coverage.

§Numbers in parenthesis indicates the percentage of bacterial and archaeal of the 16S rRNA genes copies. Based on 97% similarity clustering. Good’s coverage (%) = [1−(n/N)] × 100 (n, the number of OTUs; N, the total number of reads).

Microbial Community Structure

In total of 132,914 pyrosequencing reads with an average length of 601 bp, 9,852 reads were unique. Chao1 indices are as high as an average of 1,909 in the 7–8 cm depth of Stn 17, and were lowest at 704 in the 14–16 cm depth of Stn 83 (Table 2). Chao1 and observed OTU0.97 counts of the polynya and ice-shelf sites were higher than those of the sea-ice zone. These two indices were estimated to be highest within the sub-oxic layers (3–8 cm depth) of all sites.

Non-metric multidimensional scaling of OTUs resulted in a separate clustering between polynya and non-polynya samples, and the separation was mainly driven by differences in geochemical constituents (, Fe2+, and ), TOC, TN, and sulfate reduction rate (Figure 3). Both archaeal and bacterial communities were quite distinctive between polynya sites and non-polynya sites (P < 0.001 from both PERMANOVA and ANOSIM) (Supplementary Figure S2). Similar patterns in community structure were observed at order-level resolution (Supplementary Figure S2D). From RDA analysis (Supplementary Figure S3), NOX concentrations correlated with surface microbial communities at all stations. Corg contents correlated with microbial communities in surface sediments of Stn 10 and at intermediate depths of Stns 17 and 19, respectively. Archaeal communities were more distinctive among stations, while bacterial communities were more similar overall (Supplementary Table S2). Compositional similarity at the phylum level better reflected ecological or geographic settings.

FIGURE 3

The relative contributions of each phylum were dramatically different between sites (Figures 2, 4). At the polynya site (Stn 10), the most dominant phyla were Planctomycetes (35–71%), Proteobacteria (4.8–22%), Thaumarchaeota (1.2–12.1%), and Chloroflexi (1.9–10.8%). Candidate Division GN02 (0.5–4.4%), Candidate Division SBR1093 (0.2–8.3%), Acidobacteria (1.2–5.0%), and Bacteroidetes (0.4–3.2%) were of minor abundance. At Stn17, the proportion of Planctomycetes decreased slightly (14.8–43%), whereas the percentage of Thaumarchaeota (13–51%) appeared to be more abundant than at Stn 10. Other phyla included Proteobacteria (10.6–31%), Chloroflexi (2.1–8.7%), Candidate Division SBR1093 (0.5–8.3%), and Acidobacteria (2.6–4.1%). The proportion of archaeal sequences prominently increased at Stns 19 and 83, which is consistent with the archaeal cell abundance estimated from 16S rRNA gene quantification (Figure 2 and Table 2). Relative abundance of Thaumarchaeota accounted for 41–59% and 41–67% at Stn 19 and Stn 83, respectively. Proteobacteria (4.1–27.8%), Planctomycetes (5.8–24%), Chloroflexi (1.6–8.1%), Acidobacteria (2–8%), and Candidate Division SBR1093 (0.3–8.1%) appeared in similar proportions at both non-polynya sites (Figure 2). The relative abundance of Planctomycetes in the polynya sediments showed a high positive correlation with (r = 0.64, P < 0.001), Fe2+ (r = 0.46, P < 0.001), TOC (r = 0.75, P < 0.001), TN (r = 0.73, P < 0.001), and SRR (r = 0.48, P < 0.001) (Figure 5). In contrast, the relative abundant of Thaumarchaeota in the non-polynya sediments showed a negative correlation with TOC (r = −0.75, P < 0.001) (Supplementary Figure S4).

FIGURE 4

FIGURE 5

). The scale bar indicates the number of substitutions per nucleotide position. The color of the square bar represents the relative abundance of the sequences for each site (Stn 10, Stn 17, Stn 19, and Stn 83). Node support estimated using 1,000 bootstrap replicates. Bootstrap values above 50% are shown.

Bacterial Community Composition

The most abundant Planctomycetes in the microbial communities were divided into three clades: a Pirellula-like group, candidate order MSBL-9 (Mediterranean Sea Brine Lake-9) (), and Candidatus Brocadiae (Supplementary Figure S5). The Pirellula-like group and the candidate order MSBL-9 were the two most dominant bacterial groups at Stn 10, representing 34% and 32% of the total 16S rRNA gene sequences, respectively (Figure 4 and Supplementary Figure S5). In contrast, at Stn 83, the relative abundance of the Pirellula-like group decreased to < 17%, and the presence of the candidate order MSBL-9 group was not discernible. The members in Candidatus Brocadiae, which are capable of anaerobic ammonium oxidation (anammox) using nitrite as an electron acceptor (), were most abundant at Stn 17, and comprised a maximum of 9.2% of the total prokaryotic sequences at the 9–10 cm depth (Figure 4 and Supplementary Figure S5).

The relative abundance of Proteobacteria did not exhibit significant spatial variation among sites [Kruskal–Wallis, χ2(3) = 1.375, P = 0.711]. Similarly, the spatial distribution of Gammaproteobacteira was similar between sites. Most gammaproteobacterial sequences in all sites were affiliated with the orders Thiotrichales and Chromatiales (Supplementary Figure S6). On the other hand, based on the order level, delta- and alphaproteobacterial compositions showed spatial variations (Supplementary Figure S6). Desulfobacterales and Desulfuromonadales, which are well-known as sulfate- and sulfur (S0)-reducing bacterial groups in marine sediments, appeared to be most abundant at Stn 10, but were rarely detected at Stn 83. In contrast, the deltaproteobacterial sequences in the candidate order NB1-j appeared higher at Stn 19 and Stn 83, but showed low relative abundance at Stn 10 (Figure 4 and Supplementary Figure S6). Members of the Alphaproteobacteria made up 0.4–6.9% of total reads and were divided into three major orders (Supplementary Figure S6): Rhodobacterales, Rhodospirillales, and Rhizobiales. The sequences falling into the Rhodobacteraceae appeared more frequently in Stn 10, while members of Rhodospirillacaea were detected more frequently at non-polynya sites (Stns 19 and 83).

The OTUs in Chloroflexi and Bacteroidetes appeared to vary between polynya sites and non-polynya sites (Supplementary Figure S7). The Anaerolineae were a major member of Chloroflexi at Stn 10, whereas sequences clustered in class SAR202 and TK17 were more abundant at Stns 17, 19, and 83. Members of Bacteroidetes that have been reported as major organic matter decomposers in SO sediment (; ; ) were substantially low with relative abundance less than 5% of total 16S rRNA gene sequences in the polynya (Figure 2).

Archaeal Community Composition

The Thaumarchaeota were a predominant microbial group in sediments from the Amundsen Sea, except at the polynya center (Stn 10) (Figure 2). Most Thaumarchaeota sequences could be classified into three subgroups termed Alpha, Theta, and Upsilon (). The Thaumarchaeota subgroup Alpha was dominant in the top layers (0–5 cm) at all sites. The community composition of the Thaumarchaeota subgroup changed gradually from the top to the bottom at the non-polynya sites (Supplementary Figure S8). The relative abundance of archaeal sequences associated with Theta and Upsilon subgroups in Thaumarchaeota increased with increasing depth, especially in the marginal sea ice zone (Stn 83). Marine Benthic Group B (MBGB) and Miscellaneous Crenarchaeotic Group (MCG), which have been proposed to be Thorarchaeota () and Bathyarchaeota (), respectively, were only detected at Stn 10 (∼9.5% and ∼5% of total sequences, respectively). They have been reported as putative heterotrophic microorganisms (; ; ). Likewise, archaeal sequences related to the methane cycle, including Methanobacteria, Methanomicrobia, and Thermoplasmata, were only detected at Stn 10 (<2% of the total sequences).

Discussion

Planctomycetes-Dominated Microbial Communities in ASP Sediments

One of the prominent features revealed by 16S rRNA genes pyrosequencing was that the members of Planctomycetes were the most abundant microbial members detected in the highly productive polynya sites, especially at Stn 10 (Figure 2). Many Planctomycetes have been found attached to sinking marine aggregates in the water column (; ). At the ASP, however, none were detected in the water column (; ). Thus, the highly abundant Planctomycetes sequences in the ASP sediment could not have originated from the water column, confirming the results reported by Probandt and coworkers who showed that Planctomycetes in subtidal, sandy sediments differed from those in the overlaying water column (). Pirellula members in Planctomycetes have been reported often as heterotrophic bacteria that degrade organic matter produced by algae (; ; ), while little information is available to speculate about the ecological role of uncultured Planctomycetes (Figure 5). Most cultivated Planctomycetes are known as aerobes or facultative aerobes (). However, environmental sequences from Planctomycetes are often retrieved from anoxic zones, such as methane hydrate-bearing sediment () and subsurface seafloor sediment (). Therefore, detection of these sequences in oxygen-depleted layers is unsurprising (Figure 4).

In fine-grained marine sediment receiving high organic material input, Delta- and Gamma-proteobacteria have been reported as the predominant bacterial groups (; ; ; ). Similarly, in most benthic systems of the SO, where organic materials primarily originate from diatom blooms in the water column, proteobacterial groups such as Delta- and Gamma-proteobacteria occupy ecological niches as primary mineralizers of organic materials in the sediments (; ; ; ). Therefore, our observation of high Planctomycetes abundance, an average of 40% of total sequences (Figure 2) in ASP sediments, is intriguing. Further discussion on the dominance of Planctomycetes is presented in next section.

Organic Substrates Control Spatial Distribution of Planctomycetes

To the best of our knowledge, this is the first report suggesting the predominance of Planctomycetes in the ASP, which contrasts with the microbial communities in other benthic systems underlying diatom-dominated water columns. Major controls regulating the distribution of Planctomycetes in the sediments remain unknown, largely because no culture-based studies are available (; ). However, our statistical analysis (Figure 6 and Supplementary Figure S9) revealed that relative abundance of Planctomycetes in the ASP showed a significant positive correlation with TOC contents and inorganic constituents (, , and Fe2+), presumably resulting from the benthic Corg mineralization. Most labile Corg produced by P. antarctica blooms are rapidly decomposed by heterotrophic bacteria in the water column before reaching the sediment of the ASP (). The export flux of Corg formed by a Phaeocystis bloom is twofold slower than that formed by a diatom bloom (), and thus the contribution of P. antarctica cells to total export flux below the photic zone declines dramatically (). Indeed, in the ASP, most photosynthetically-produced particulate organic carbon (POC) (>95%) exported from the surface layer is converted to suspended POC and/or dissolved carbon within the top 400 m of the water column (; ). In addition, direct rate measurement conducted in this study also demonstrated that, despite high primary production (110 mmol C m–2 d–1) in the Phaeocystis-dominated water column of the ASP, benthic Corg mineralization (average, 2.1 mmol C m–2 d–1) accounted for only 1.9% of primary production, which was strikingly lower than that measured in other less productive polar regions (). Therefore, organic materials accumulated in surface sediments of the deep ASP (>700 m) would be intrinsically recalcitrant. Meanwhile, it has been well-established that Phaeocystis colonies excrete a mucous matrix containing both carboxylated and sulfated heteropolysaccharides as main constituents (; ). Based on metagenomic information from Namibia and Oregon coastal upwelling systems, revealed that all marine Planctomycete genomes, except for Candidatus Kuenenia stuttgartiensis, possess a high number of sulfatase genes, which suggests that marine Planctomycetes may be able to break down recalcitrant sulfated heteropolysaccharides (; ). Consequently, the high relative abundance of Planctomycetes at the center of the polynya (Stns 10 and 17) suggests that the members of Planctomycetes constitute a significant heterotrophic bacterial group utilizing recalcitrant organic materials produced by Phaeocystis blooms in the water column of the ASP. Recently, also suggested that Planctomycetes play a key role in the degradation of high molecular weight compounds and recalcitrant materials entering surface sediment from the water column of the Wadden Sea.

FIGURE 6

High Archaeal Abundance in Marginal Sea Ice Zone

In addition to the Planctomycetes-dominated bacterial communities in the sediments of the polynya sites (Stns 10 and 17), another interesting finding revealed from the quantification of 16S rRNA gene was that archaeal abundance occupied more than half (30–71%) of the total prokaryotic abundance in the MIZ (Stn 83), in which most archaeal 16S rRNA gene sequences were assigned to Thaumarchaeota (Figure 2 and Table 2). Members of Thaumarchaeota are major contributors to aerobic ammonia and nitrite oxidation in aquatic environments (). To examine the potential of sedimentary Thaumarchaeota to oxidize ammonia, we quantified the archaeal amoA gene, which is known for a genetic marker for the ammonia oxidation. The depth profiles of archaeal amoA and 16S rRNA gene copy numbers were very similar to each other in the samples except for some layers (Supplementary Figure S1), which implies that most archaeal members at all depths have a gene encoding ammonia monooxidase.

Both the proportion of archaeal abundance in total prokaryotic abundance (Supplementary Figure S10) and the relative abundance of Thaumarchaeota in total 16S rRNA gene sequences (Supplementary Figure S4B) showed a negative correlation with TOC contents. Previously, the environmental members of Thaumarchaeota have been reported in oligotrophic marine sediments as an important chemolithotrophic microbial assemblage (; ). Similarly, cultivated thaumarchaeal ammonia oxidizers have been shown to be adapted to oligotrophic conditions (). Recently, based on continuous succession of an active microbial group belonging to Thaumarchaeota in the absence of an external Corg supply, suggested that Thaumarchaeota play a significant role in sustaining oligotrophic bathypelagic ecosystems where the Corg supply is limited. Likewise, metaproteomic data indicate that MG-I are abundant and metabolically active at the surface water of the West Antarctic Ocean during the ice-covered winter (), and chemoautotrophic carbon fixation by Thaumarchaeota contributes up to 9% of bacterioplankton production during the Antarctic winter (). Indeed, in the present study, primary production was low in the water column of the ice-covered sea-ice zone (Stn 83) (; ), and the Corg content in the sediment was consequently lowest (<0.5%, Table 1). Our results thus suggest that the Thaumarchaeota predominating in the sediment of non-polynya sites are a significant chemolithoautotrophic group, sustaining the benthic ecosystem where input of Corg from the water column is limited.

Planctomycetes-Dominated Communities as a Proxy for Assessing Climate Change-Induced Ecosystem Response in the ASP

The structure and function of benthic heterotrophic microbial communities are ultimately controlled by the quantity and quality of organic matter from the water column (; ; ; ; ). Thus, the predominance of Planctomycetes in ASP sediment provides ecological and environmental baseline information for assessing the responses of benthic ecosystems to transitions in phytoplankton communities that may result from ongoing and future climate change in the Amundsen Sea. Due to the inflow of warm CDW, the western Antarctic Ocean including the Amundsen Sea has experienced significant surface warming and loss of sea ice (; ; ; ). In particular, glaciers near the Amundsen Sea are undergoing the highest rates of melting and thinning on the Antarctic continent (). Sea ice melting increases the supply of nutrients and intensifies water column stratification (, ; ), which consequently results in increased primary productivity (; ) and shifts in phytoplankton community structure from P. antarctica to diatoms (; , ). As discussed in previous sections, while Corg originating from Phaeocystis blooms settles slowly, the particulate organic material produced by diatoms possesses a relatively faster sinking rate and a lower C:N ratio (i.e., they are relatively more labile) than those produced by Phaeocystis (; ). Therefore, the transition of dominant phytoplankton from Phaeocystis to diatoms would deposit more labile organic matter into the benthic system of the ASP, which ultimately results in a shift in the benthic microbial community structure from Planctomycetes to Proteobacteria. Likewise, in the marginal sea-ice zone, the loss of sea ice due to global warming may cause an increase in primary production in the water column (; ). As a result, benthic microbial communities in the MIZ can undergo a shift from Thaumarchaeota to other heterotrophic bacteria. Consequently, in the Amundsen Sea, a spatial distribution of the benthic microbial community reflecting the composition and production of dominant phytoplankton groups in the water column can act as a relevant proxy for detecting possible variation in primary productivity and phytoplankton community that is associated with sea ice melt due to the global warming.

Conclusion

Characterization of benthic microbial communities is underrepresented in the Southern Ocean, especially in the Amundsen Sea where environmental changes due to global warming are occurring rapidly. We conducted high-throughput sequencing of 16S rRNA genes and quantitative PCR in combination with biogeochemical analyses and metabolic rate measurements to determine the composition, diversity and controls of major microbial communities in sediments of the ASP. Planctomycetes predominated benthic bacterial communities in the polynya sediment underlying the Phaeocystis-dominated water column, which suggests that Planctomycetes comprise a major bacterial group utilizing relatively recalcitrant Corg produced primarily by Phaeocystis blooms. Planctomycetes-dominated microbial communities in the ASP sediment are in stark contrast to that Proteobacteria (Delta- and Gamma-proteobacteria) occupy ecological niches as primary mineralizers of organic materials in most benthic systems in the Southern Ocean, where organic materials in the sediments mostly originate from diatom blooms. Sea ice melting by global warming increases the supply of nutrients and intensifies water column stratification, which consequently results in increased primary productivity and shifts in phytoplankton community structure from P. antarctica to diatoms. Given that microbial communities respond quickly to environmental changes, and that global climate change is proceeding rapidly in the ASP, our results suggest that any modifications in the Planctomycetes-dominated microbial communities will provide valuable insight into changes in organic matter transport to the seafloor that may result from shifts in planktonic communities from Phaeocystis to diatoms associated with climate change.

Statements

Data availability statement

The datasets generated for this study can be found in the SRX3405376, https://www.ncbi.nlm.nih.gov/sra/SRX3405376.

Author contributions

HC and J-HH as the first and corresponding author, designed the original experiments and conducted most writing. HC, CH, J-GK, and SK conducted molecular and statistical analysis. S-HK, AC, and J-HH participated in onboard research activities and performed biogeochemical analysis. J-GK, SK, KK, S-KR, EY, and SL participated in writing and discussion of the results. All authors read and approved the final manuscript.

Funding

Analysis for microbial communities was supported by the Korea Polar Research Institute (PE19060) and the Mid-career Researcher Program funded by the Korean Ministry of Science and ICT (No. 2018R1A2B2006340). Biogeochemical analysis and metabolic rate measurements were supported by the Korea-Long-Term Marine Ecological Research (K-LTMER) funded by the Korean Ministry of Oceans and Fisheries.

Acknowledgments

We thank the captain and crew of the Korean Research Icebreaker Araon, for their outstanding assistance during the cruises.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2019.00797/full#supplementary-material

References

Summary

Keywords

Amundsen Sea polynya, Antarctica, benthic microbial community, global climate change, organic carbon mineralization, Phaesocystis antarctica, Planctomycetes, Thaumarchaeota

Citation

Cho H, Hwang CY, Kim J-G, Kang S, Knittel K, Choi A, Kim S-H, Rhee S-K, Yang EJ, Lee S and Hyun J-H (2020) A Unique Benthic Microbial Community Underlying the Phaeocystis antarctica-Dominated Amundsen Sea Polynya, Antarctica: A Proxy for Assessing the Impact of Global Changes. Front. Mar. Sci. 6:797. doi: 10.3389/fmars.2019.00797

Received

27 September 2019

Accepted

11 December 2019

Published

10 January 2020

Volume

6 - 2019

Edited by

David A. Walsh, Concordia University, Canada

Reviewed by

Jiwen Liu, Ocean University of China, China; Yong Jiang, Ocean University of China, China; André M. Comeau, Dalhousie University, Canada

Updates

Copyright

*Correspondence: Jung-Ho Hyun,

This article was submitted to Aquatic Microbiology, a section of the journal Frontiers in Marine Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics