- Open Access
Genomic and metatranscriptomic analyses of carbon remineralization in an Antarctic polynya
- So-Jeong Kim†1,
- Jong-Geol Kim†2,
- Sang-Hoon Lee3,
- Soo-Je Park4,
- Joo-Han Gwak2,
- Man-Young Jung5,
- Won-Hyung Chung6,
- Eun-Jin Yang3,
- Jisoo Park3,
- Jinyoung Jung3,
- Yoonsoo Hahn7,
- Jang-Cheon Cho8,
- Eugene L. Madsen9,
- Francisco Rodriguez-Valera10,
- Jung-Ho Hyun11 and
- Sung-Keun Rhee2Email author
© The Author(s). 2019
- Received: 22 August 2018
- Accepted: 1 February 2019
- Published: 20 February 2019
The Correction to this article has been published in Microbiome 2019 7:38
Polynyas in the Southern Ocean are regions of intense primary production, mainly by Phaeocystis antarctica. Carbon fixed by phytoplankton in the water column is transferred to higher trophic levels, and finally, to the deep ocean. However, in the Amundsen Sea, most of this organic carbon does not reach the sediment but is degraded in the water column due to high bacterial heterotrophic activity.
We reconstructed 12 key bacterial genomes from different phases of bloom and analyzed the expression of genes involved in organic carbon remineralization. A high correlation of gene expression between the peak and decline phases was observed in an individual genome bin-based pairwise comparison of gene expression. Polaribacter belonging to Bacteroidetes was found to be dominant in the peak phase, and its transcriptional activity was high (48.9% of the total mRNA reads). Two dominant Polaribacter bins had the potential to utilize major polymers in P. antarctica, chrysolaminarin and xylan, with a distinct set of glycosyl hydrolases. In the decline phase, Gammaproteobacteria (Ant4D3, SUP05, and SAR92), with the potential to utilize low molecular weight-dissolved organic matter (LMW-DOM) including compatible solutes, was increased. The versatility of Gammaproteobacteria may contribute to their abundance in organic carbon-rich polynya waters, while the SAR11 clade was found to be predominant in the sea ice-covered oligotrophic ocean. SAR92 clade showed transcriptional activity for utilization of both polysaccharides and LMW-DOM; this may account for their abundance both in the peak and decline phases. Ant4D3 clade was dominant in all phases of the polynya bloom, implicating the crucial roles of this clade in LMW-DOM remineralization in the Antarctic polynyas.
Genomic reconstruction and in situ gene expression analyses revealed the unique metabolic potential of dominant bacteria of the Antarctic polynya at a finer taxonomic level. The information can be used to predict temporal community succession linked to the availability of substrates derived from the P. antarctica bloom. Global warming has resulted in compositional changes in phytoplankton from P. antarctica to diatoms, and thus, repeated parallel studies in various polynyas are required to predict global warming-related changes in carbon remineralization.
- Carbon remineralization
Large polynyas developed in the Antarctic coast during the austral summer are ecological hotspots due to extremely productive phytoplankton bloom and intense biogeochemical cycling. In the spring, melting glaciers and sea ice stratify the water column and increase levels of light and nutrients for stimulating the growth of phytoplankton, especially Phaeocystis antarctica [1, 2]. Phaeocystis is one of the few microalgae that benefits strongly from eutrophication, resulting in almost uni-algal plankton blooms [3, 4]. This alga has a polymorphic life cycle (free-living cells and colonies composed of thousands of cells embedded in a polysaccharidic matrix). Organic matter is released from phytoplankton upon lysis by zooplankton grazing, virus infection, or senescence . In fact, the highest chromophoric dissolved organic matter (DOM) concentrations are frequently observed near the peak of the phytoplankton bloom within polynyas of the Amundsen Sea (the ASPIRE cruise in 2010–2011) .
Carbon fixed by phytoplankton in the water column is transferred to higher trophic levels, and finally, is exported from the system to the deep ocean . However, intriguingly, in the Amundsen Sea dominated by P. antarctica, most of the organic carbon from primary production does not reach the sediment; instead, it is actively recycled and mineralized in the water column via microbial loop [8–11]. Indeed, higher heterotrophic bacterial production, accounting for ca. 20% of the organic carbon produced by primary production, has been observed in the surface mixed layer of the polynya . These reports indicate that the Amundsen Sea polynya (ASP) is highly inefficient at transporting carbon to the deep sea via biological pump.
Increases in the pCO2 concentration and water column stratification  by global warming may have shifted the phytoplankton community structure from P. antarctica to diatoms in the Southern Ocean [12–14]. Diatoms possess relatively faster sinking rate [15, 16]. This makes export flux of organic carbon formed by diatom blooms two times faster than that by Phaeocystis blooms , which dramatically increases the contribution of diatoms to total export below the photic zone .
Bacterial assemblages may be responsible for the enhanced assimilation and remineralization of organic matter from phytoplankton blooms in polynyas. Associations of bacteria with phytoplankton blooms have been reported in polar oceans based on conventional clone library analyses . Using next-generation sequencing approaches with enhanced sequencing depth [19–21], distinct prokaryotic assemblages have been observed at the peak of the P. antarctica-dominated bloom in ASP [21, 22]. However, logistical constraints and safety concerns have limited the investigation of the temporal succession of microbial-bloom communities in Antarctic waters—particularly in the decline period (late February) following the peak in phytoplankton activity.
Furthermore, although bacterial community composition in various ocean provinces has been described, ecophysiological and biochemical properties of key bacteria are still hard to predict owing to their low culturability in artificial media. Since free-living bacterioplankton plays a pivotal role in the mineralization of algae-derived organic matter, free-living bacterioplankton was frequently used to explain bacterioplankton succession dynamics during the bloom . Genomic reconstruction techniques from metagenomic sequences can be used to characterize uncultivated lineages of microorganisms. Physiological properties and biochemical pathways can be predicted using genomes reconstructed by de novo assembly and binning from random shotgun genomic libraries that are obtained directly from environmental samples [24–26]. Further, environmental mRNA sequence reads can be mapped to metabolic pathway genes of reconstructed genome bins to evaluate in situ activities . In the present study, nucleic acids harvested from water samples gathered during three cruises were used to reconstruct genomes of key heterotrophic bacteria from a P. antarctica-dominant ASP. We also analyzed the expression of genes in the reconstructed genomes using metatranscriptomics. A comparative analysis of samples from different phases of the bloom (peak and decline) and from sea ice-covered ocean provides fundamental insights into the ecology of abundant bacteria associated with phytoplankton blooms in these Antarctic polynyas.
Sampling and geochemical analysis
During three cruises, four water samples from the Amundsen Sea (two each from the peak and decline phases) were gathered from the center of a polynya. The average bloom termination date is February 23 (± 5.38 days) . Satellite-based surface chlorophyll data was obtained from the GlobColour webpage (http://hermes.acri.fr/), which has been constructed from multi ocean color sensors. In addition, a control ocean water sample was gathered from below sea ice by ice core drilling (Additional file 1: Table S1). A rosette water sampler equipped with ten Niskin bottles (10 L) and a conductivity-temperature-depth unit (sensors for Chl-a fluorescence, pressure, temperature, and salinity) was deployed. The concentrations of Chl-a and inorganic nutrients (nitrate + nitrite, ammonium, phosphate, and silicate) were measured aboard, as described by Kim et al.  Briefly, the concentration of Chl-a extracted in 90% acetone was measured using a Trilogy fluorometer (Turner Designs, San Jose, CA, USA). The inorganic nutrients were analyzed using a QuAAtro Auto Analyzer (SEAL Analytical, Southampton, UK). Bacterial 16S rRNA gene copies were evaluated by real-time quantitative PCR (qPCR) as described by Kim et al.  For qPCR, the Bac518F and Bac786R primer set was used. To evaluate the phytoplankton composition, sampled water was immediately fixed using glutaraldehyde. Filter-harvested phytoplankton cells were stained with a staining solution, and the composition was manually checked using a microscope. The detailed methods are described in Additional file 2: Supplementary methods according to Lee et al. . The depth of the samples was set to the maximum chlorophyll layer.
Free-living bacterioplankton cells were harvested by sequential filtration of 10 L of water through a 1.2-μm pore-sized filter and 0.22-μm pore-sized filter (Supor polyethersulfone, Pall Life Sciences, Ann Arbor, MI, USA). Seawater samples were collected for approximately 1–2 h around midday during each sampling period. Within 30 min of collection, samples were filtered and maintained at under 10 °C. The filters were immediately stored at − 80 °C until further processing in the laboratory.
Nucleic acid extraction and pyrosequencing
The total nucleic acids were extracted from cells collected on filters with a 0.22-μm pore size using a previously described protocol . Briefly, harvested cells in the filter were grazed with quartz and liquid nitrogen. Grazed cells were lysed with CTAB buffer, and nucleic acids were separated using chloroform. Ethanol precipitation was used to purify and concentrate total nucleic acids. DNA and RNA were separated using the AllPrep DNA/RNA Mini Kit (Qiagen, Valencia, CA, USA) following the manufacturer’s protocols.
For the community analysis, the V1 to V3 regions of the 16S rRNA gene were amplified with the barcoded primer set V1-9F (5′-GAGTTTGATCMTGGCTCAG-3′) and V3-541R (5′-WTTACCGCGGCTGCTGG-3′) . PCR was conducted using the following conditions: 94 °C for 5 min, 30 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 45 s, elongation at 72 °C for 60 s, and 72 °C for 5 min for the final extension. The amplified PCR products were purified using the PCR Purification Kit (Cosmogentech, Seoul, Republic of Korea). PCR products for each sample were evenly mixed (0.4 μg per sample), quantified using the NanoDrop Spectrophotometer (Wilmington, DE, USA), and then sequenced by ChunLab Inc. (Seoul, Republic of Korea) using a Roche/454 GS FLX Titanium sequencer. Pyrosequencing data were processed as follows. Chimeric sequences were removed using UCHIME (v. 4.2) . OTUs at the 97% similarity level were obtained using uclust . A taxonomic classification was assigned to each read using Greengenes as a reference database . These steps were performed using the pick_de_novo_otus.py script in QIIME (ver. 1.9.1). Unassigned reads that showed no match to sequences in the Greengenes database were manually reassigned against NCBI Microbial 16S rRNA gene database using BLASTn (cutoff, e value 1 × E−5). Diversity and species richness indices were calculated using the alpha_diversity.py script with 5456 reads subsampled from each sample. Neighbor-joining trees were constructed using representative sequences of each OTU selected by QIIME. These sequences were compared against the NCBI nucleotide database and the EzBioCloud database to select the nearest neighbors. Selected sequences were aligned and edited using BioEdit, and a tree was reconstructed using MEGA 7.0  with the Kimura 2-parameter model with 1000 bootstrap replicates. Additionally, previous pyrosequencing data obtained during the 2010/2011 cruise  were reanalyzed using QIIME to compare bacterial communities.
Statistical analysis of environmental factors and the microbial community
To identify the relationships between major bacterial groups and the environmental parameters, a non-metric multidimensional scaling (NMDS) analysis was performed using R with the vegan package (v. 2.4).
Metagenome sequencing, assembly, and annotation
Three genomic DNAs obtained above (marked in Additional file 1: Table S1; where PK indicates the peak phase of the bloom, DC indicates the declining phase, and SI indicates the ocean sample under sea ice) were used to construct Illumina sequencing libraries with target insert sizes of 300 and 5000 bp using the TruSeq PCR Free Kit and Nextera Mate Pair Kit, respectively. Fragments of libraries were used for paired-end and mate-pair sequencing using the Illumina Hiseq2000 according to the manufacturer’s instructions at Macrogen (Seoul, Republic of Korea). Raw reads were trimmed using Sickle with default settings (https://github.com/najoshi/sickle/). Each paired-end read set was assembled using the IDBA-UD (v. 1.1.1) assembler since it provided the best results. The quality of assembly was determined by using N50. Default settings in IDBA_UD were used since those are normally suitable for most metagenome datasets (https://ggkbase-help.berkeley.edu/overview/data-preparation-metagenome/) .Scaffolding was performed by SSPACE (v. 2.1) using preassembled contigs from each sample with mate-pair reads . For taxonomic classification of metagenome and metatranscriptome reads, Centrifuge (v. 1.0.4) was used against NT database (supported by Centrifuge) .
For scaffolds of ≥ 1 kb, putative genes were predicted using MetaGeneAnnotator . Functions of the protein sequences were predicted using BLASTp similarity searches against the NCBI NR database, KEGG, COG, and TIGR (e value cutoff; 1 × E−5). The rRNAs and tRNAs were identified using RNAmmer  and tRNAscan-SE , respectively. Domain information was obtained using the Pfam database (pfam_scan.pl, default parameters). TonB-dependent transporters (TBDTs) were classified as described by Tang et al. . For the classification of ABC transporters and peptidases, putative genes were compared against the TCDB  and MEROPS databases  using BLASTp (e value cutoff; 1 × E−5), respectively. Signal peptides were predicted for the determination of extracellular proteins using SignalP version 4.1 .
Binning and phylogenetic analysis
The binning was performed with SI metagenomic scaffolds using the differential coverage plotting method in combination with the tetranucleotide frequency, as described by Albertsen et al.  based on R script (http://madsalbertsen.github.io/multi-metagenome//). Bowtie2 (v. 2.1.0) and SAMtools (v. 1.9) were used to calculate scaffold coverages [45, 46]. The coverage value was calculated according to Lander and Waterman  as LN/G, where L is the read length, N is the number of mapped reads, and G is the scaffold length (excluding Ns). After binning, unmatched scaffolds to the bin were manually removed according to a taxonomic analysis of scaffolds. The taxonomy of scaffold was determined based on an assignment of the majority of its genes to a specific taxon using BLASTp against the NR database. To estimate the genome completeness of reconstructed bins, CheckM (v. 1.0.5) was used . Using criteria (completeness and contamination) evaluated using CheckM, final bins were selected for further analysis. For the taxonomic classification of genomes, rRNA (rpoB; in case of GM5, recA was used because this bin lacked the rpoB gene) was used. Single marker genes were analyzed using BLASTp against the NR database. Hits for 16S rRNA genes were extracted from the total reads using BLASTn (alignment length ≥ 100 and e value ≤ 0.001) for comparisons of pyrosequencing data .
The TruSeq RNA Kit was used to generate the cDNA template library from extracted RNA, and transcript reads were obtained by Illumina paired-end Hiseq2000 sequencing according to the manufacturer’s instructions at Macrogen. Raw transcript reads were preprocessed using the FASTX-Toolkit (v. 0.0.14) for quality filtering (minimum quality score, 20; minimum percent of bases that must have a quality score exceeding 20, 90; http://hannonlab.cshl.edu/fastx_toolkit/). To remove rRNA reads, reads with matches in the SILVA database (SSU and LSU ver. 127)  and 5S rRNA database  were removed using Bowtie2 and SAMtools. Reads mapped to eukaryotic and viral genomes were removed using Centrifuge (v. 1.0.4). For the transcriptomic analysis, predicted bacterial scaffolds among assembled scaffolds exceeding 1 kb obtained from SI were also used as a custom nucleotide database. Non-rRNA cDNA reads from three samples were mapped against the predicted gene of the custom nucleotide database using Bowtie2 (under default parameters), and assigned reads were counted using SAMtools. Unmapped reads were aligned against UniRef90 using DIAMOND blastx search (query cover 90, e value 1, according to HUMAnN2) [52, 53]. Finally, mapped reads to metagenome and UniRef90 were used for mRNA normalization. Expression levels of individual genes were obtained using TPM value (transcripts per kilobase million).
Correlation between DNA-TPM and mRNA-TPM
The normality of the population distribution was evaluated using the Shapiro–Wilk test for samples of mRNA-TPM. As the data were not normally distributed, the Spearman rank correlation was determined using R (vegan package, v. 2.4) to determine correlations between PK-mRNA-TPM:DC-mRNA-TPM for each sample.
Description of the location and oceanographic properties of sampling sites
2010/2011 cruise; 07 Jan 2011
2011/2012 cruise; 14 Feb 2012
2011/2012 cruise; 29 Feb 2012
2013/2014 cruise; 03 Jan 2014
2010/2011 cruise; 05 Jan 2011
Sampling pointa (depth)
Polynya (30 m)
Polynya (20 m)
Polynya (25 m)
Polynya (20 m)
Under sea ice (2 m)
Abundance of bacteria (×106 16S rRNA gene copies/mL)
NO2- + NO3- (μmol/L)
Phytoplankton composition (μg C/L)
Total carbon (μg C/L)
Carbon to Chl-a
Bacterial community composition and dynamics
Based on an analysis of 16S rRNA gene sequences, alpha diversity indices in the decline phase were greater than those during the peak phase of the bloom and in the sea ice-covered ocean (Additional file 1: Table S1). Bacteroidetes and Proteobacteria were the most abundant phyla in all samples (Additional file 1: Table S1). Reads of 16S rRNA gene sequences of Polaribacter, the gammaproteobacterium Ant4D3 (Oceanospirillaceae), and Pelagibacter were abundant at ASP. Polaribacter sequences were particularly dominant (37.4–50.5% of total sequences) in the peak phase of the bloom. Despite differences between bloom phases, the major OTUs assigned to Polaribacter, SAR92, Ant4D3, and SAR11 were nearly identical (Additional file 3: Figure S2). The taxonomic distributions are described in detail in Additional file 1: Table S1.
A NMDS analysis of bacterial communities showed that the communities at the peak phase clustered with each other but were distinct from those at the decline phase or in the sea ice-covered ocean, although they were collected in different years (Additional file 3: Figure S3). The bacterial communities in the two different decline-phase samples clustered strongly together, despite differences in various factors, such as P. antarctica abundance, total algal carbon, and total bacterial abundance. The bacterial communities in the sea ice-covered oceans were variable but were distinct from those in polynya samples. The major environmental parameters influencing bacterial community structure were Chl-a and nitrite + nitrate concentration and depth, and the samples of peak phase were positively related to the Chl-a concentration.
Features of the twelve reconstructed genomes (genome bins) identified by an analysis of metagenomic sequences derived from DNA extracted from Antarctic waters
Total scaffold length (contig length) (Mb)
Number of scaffolds
Key genes for carbon assimilation (#)
Rank abundances of reconstructed genomes
mRNA transcript abundance
mRNA reads from PK, DC, and SI were classified as shown in Additional file 1: Table S5. Among non-rRNA (mRNA) reads, 66.2%, 42.1%, and 63.1% were matched to the total metagenomic scaffolds and uniref90 database for PK, DC, and SI, respectively (Additional file 1: Table S6). In the peak phase, approximately 25.1% and 6.4% of mRNA transcripts were assigned to genes of BC1_Pol and GM1_Ant, respectively, suggesting their high abundances. In SI, 3.4% of total bacterial mRNA reads matched to genes of AL1_Pel.
Pairwise comparisons of the expression levels of all genes (mRNA-TPM:mRNA-TPM) in PK and DC showed a low overall correlation (p = 0.16). However, in an individual genome bin-based pairwise comparison of gene expression between PK and DC, a high correlation of gene expression was observed, with regression slopes ranging from 0.44 to 0.90 (Additional file 3: Figure S5; Additional file 1: Table S7). The strong correlation was further supported by the low average fold changes in transcripts in most of the bins between PK and DC (Additional file 1: Table S8).
Traits overrepresented in each phase
To investigate the specialized metabolic potential for biogeochemical cycling, gene transcripts abundant during each phase of the bloom were analyzed using reconstructed genome-based approaches. Pairwise comparisons of the expression (mRNA-TPM: mRNA-TPM) of genes involved in biogeochemical processes indicated that gene transcripts overrepresented in PK and DC belong to the taxa with the highest abundance in PK and DC in the metagenome, respectively (Additional file 3: Figure S6a). The expression of distinct repertoires of genes for transporters, the glycoside hydrolase (GH) family, and the assimilation of compatible solutes (Additional file 3: Figure S6b–e) were overrepresented at each phase and were further investigated.
Although GM4_SAR92 belongs to Gammaproteobacteria, putative genes for TBDTs were enriched (Additional file 1: Table S9) and transcripts assigned to the putative genes were abundant. In particular, transcripts assigned to TBDT cluster 427 for the uptake of arabinose, a major component of the saccharide pools of Phaeocystis cells , were observed at the peak phase (Fig. 2b). In addition, putative genes for TBDTs of cluster 973 and 410 in GM4_SAR92 involved in the transport of vitamin B12 and various cofactors (especially iron), respectively, were similarly expressed in all three samples.
Putative genes encoding ABC transporters specific for the uptake of amino acids, branched amino acids, peptides, and monosaccharides were abundant in GM1_Ant and GM2_Ant (n = 73 and n = 97, respectively), and transcripts assigned to these putative genes were abundant in the peak phase (Table 2 and Fig. 2b). Transcripts of putative genes for the ABC transporter for glycerol uptake were also prominent in GM2_Ant (Fig. 2b). Further, putative TRAP genes (in particular, for mannitol uptake) were enriched (n = 12 and n = 16, respectively) in GM1_Ant and GM2_Ant, with high recruitment of transcripts in the peak phase (Fig. 2). Glycerol and mannitol are used as osmolytes  and might be widely available substrates in polynya blooms. Putative genes encoding mono- and disaccharide, peptide, and glycine betaine (GB) ABC transporters were abundant in GAM6_SUP05 and highly expressed in the decline phase (Fig. 2c). High abundances of transcripts for ABC and TRAP transporters in AL1_Pel for low molecular weight (LMW)-DOM, such as amino acids, taurine, mono- and disaccharides, were observed in the sea ice-covered ocean (Fig. 2d).
Utilization of polysaccharides
Putative genes in the GH family were abundant (n = 27) in BC1_Pol (Table 2), i.e., there were more than four copies of GH16, GH20, and GH92 for the degradation of laminarin, chitin, and mannose-containing polysaccharides, respectively. Many of these putative genes were co-localized with TBDTs and SusD genes in five putative polysaccharide utilization loci (PUL). In contrast to BC1_Pol, BC5_Pol contained six of ten putative sulfatase genes residing in PUL associated with the potential for sulfated polysaccharide assimilation.
Assimilation of dimethylpropionate
Putative genes for dimethylpropionate(DMSP) demethylase (dmdA) were widely present and expressed in the major bacteria in the PK, DC, and SI samples, i.e., GM1_Ant and GM4_SAR92 (Additional file 3: Figure S8a). While dmdA of GM1_Ant was highly expressed at PK, that of GM4_SAR92 was highly expressed in all three samples. Putative genes related to the assimilation of 3-methylmercaptopropionate (MMPA) were co-localized with DMSP demethylase in GM1_Ant, GM2_Ant, and GM3 (Additional file 3: Figure S8b). Assimilation of glycine betaine and other compatible solutes is described in the Additional file 2: Supplementary results.
Phytoplankton blooms in Antarctic polynyas are accompanied by various bacterial taxa, as revealed by 16S rRNA gene analyses . In order to investigate bacterial plankton succession dynamics during the bloom, we used samples for two phases, namely the peak and decline of bloom collected over different years. Due to the unique nature of blooms in ASP (i.e., a short period of opening in the sea ice, low temperature, the abundance of P. antarctica-derived organic matter), inter-annual variation in biodiversity may be limited, and thus, bacterioplankton dynamics likely exhibit recurrent patterns. Similarly, repeated patterns have been observed in the North Sea among specific bacterial clades of four consecutive blooms, despite some inter-annual variation in physicochemical conditions . Moreover, the dominance of Polaribacter at the peak of the bloom has been observed evident in previous studies of ASP [21, 57]. The obligate psychrophilic nature of Polaribacter spp. among members of Bacteroidetes isolated from Southern Ocean including ASP  may contribute to their adaptation to blooms of ASP . The high ANI value and conserved synteny between fosmid clone Ant4D3  and GM1_Ant, which was dominant in all phases of the polynya bloom (see Additional file 3: Figure S9), support the recurring pattern of the most well-adapted clades in the Southern Ocean.
For adaptation to the polynya blooms, each dominant bacterial clade may harbor distinct metabolic properties. Our genome reconstruction and metatranscriptomic analysis provide relevant information on the major taxa responsible for DOM utilization in the ASP. Polaribacter and SAR92 clades had the potential to utilize HMW-DOC, while the Ant4D3, SAR92, and SUP05 clades had the potential to utilize LMW-DOC. Further fine-scale taxonomic (ecotype) differences could also be observed.
In contrast to diatom-dominated spring bloom in the North Sea, where the expression of GHs most likely allowed diverse members of Bacteroidetes to decompose algal biopolymers , Polaribacter was the single predominant clade associated with the degradation of organic matter derived from P. antarctica at the peak bloom in the ASP (Table 1). The direct release of LMW-DOM by growing P. antarctica cells is low [62–64], and most DOM from intact P. antarctica cells is composed of polysaccharides. This corresponds with the predominance of novel Polaribacter ecotypes, BC1_Pol and BC5_Pol.
Chrysolaminarin, a β-1,3-glucan, laminarinase (GH16)-specific substrate, is considered one of the most abundant polysaccharides in Phaeocystis blooms . BC1_Pol might be specialized to degrade chrysolaminarin with high expression of a laminarinase (GH16) gene in PUL. Arabinose and xylose are predominant in the carbohydrate composition of Phaeocystis [15, 65], which is a crucial distinction from diatoms . Thus, pentose-containing polysaccharides might be widely available substrates during P. antarctica blooms, and BC5_Pol with the genetic potential for xylan degradation (GH10, endo-1,4-beta-xylanase) may be specialized for the utilization of these polysaccharides. Despite the high average nucleotide identity (ANI 91.8%) between BC1_Pol and BC5_Pol, our results indicate differentiation in the utilization of polysaccharides (e.g., chrysolaminarin and xylan with GH16 and GH10, respectively; Fig. 3 and Fig. 4), which are significant polymers in P. antarctica , using a distinct set of GHs. Similar to observations by Teeling et al. , frequencies of CAZymes of BC1_Pol and BC5_Pol in metagenomes and metatranscriptomes were significantly decreased as blooms subsided in DC. Detection of transcripts from putative genes encoding GHs in addition to the high representation of TBDT genes in GM4_SAR92 indicates a potential role of GM4_SAR92 in the degradation and assimilation of algae-derived polysaccharides, as previously suggested in a diatom bloom , which is a key distinction from other members of Gammaproteobacteria.
Phaeocystis accumulates large amounts of dimethylpropionate (DMSP) (~ 150 mM in cells) as a compatible solute, which can be released from decaying P. antarctica cells [67, 68]. One of the common traits of the abundant Ant4D3 and SAR92 clades in the Amundsen Sea was the metabolic potential of C1 compounds associated with the assimilation of these compatible solutes. DMSP can be converted to a highly reactive volatile sulfur species, methanethiol (MeSH), by DMSP demethylase (DmdA) (see Additional file 3: Figure S8c). MeSH is likely to remain in the ocean surface and is readily used by marine bacteria as a carbon and sulfur source . DMSP lyase (DddD) and CoA transferase (DddD and DddL, − W and − Y) catalyze key biochemical steps in the conversion of DMSP to dimethyl sulfide (DMS), an important climate change gas . Remarkably, representative genes encoding DMSP lyase and CoA transferase were not observed in the total assembled scaffolds. Since significant concentrations of DMS were observed at these study sites , DMSP lyase of Phaeocystis  or an unknown pathway might be involved in DMS production in the polynya.
Unique LMW-DOM-assimilating bacterial assemblages in the decline phase in the Amundsen Sea are dictated by gene expression in the SUP05, SAR92, and Ant4D3 clades of Gammaproteobacteria (Additional file 3: Figure S6d), while the RCA (Roseobacter clade-affiliated) cluster of Alphaproteobacteria assimilates LMW-DOM in spring algal blooms in other oceans [73–75]. A comparative genomic analysis indicated that the metabolic properties of the abundant clades of Gammaproteobacteria largely overlap with those of the RCA cluster, e.g., amino acid and C1 compound (e.g., DMSP) assimilation and cobalamin synthesis (Additional file 5: Table S11). Additionally, the differentiation of the clades of Gammaproteobacteria from the RCA cluster can be observed in the glyoxylate shunt, proteorhodopsin synthesis, CO oxidation, and urea utilization (Additional file 1: Table S10). Isocitrate lyase of GM6_SUP05 and GM1_Ant was highly expressed (ranked 229th and 137th in the expression in PK and DC, respectively (Additional file 1: Table S12). At polynyas with iron limitation, cells with the glyoxylate shunt can use fewer hemeproteins and make fewer electron carriers to decrease oxidative stress . Genes for CODH were only observed in RCA clades for which iron is required as an essential cofactor. Insufficient ATP or proton motive force can be supplemented by proteorhodopsin, which is an abundant gene transcript in the SUP05, SAR92, and Ant4D3 clades in both PK and DC (Additional file 3: Figure S6e). These genomic features may contribute to the dominance of gammaproteobacterial clades in Phaeocystis-dominant Antarctic blooms in which iron is depleted . The possible depletion of iron was supported by the high expression of genes encoding TBDT and ABC transporters involved in iron transport in most of the genomes (Fig. 2).
The polynya is a biogeochemical hotspot where high carbon biomass reached 10 mg C 1−1is sequestered via biological pump . Surprisingly, an exceptionally low amount of fixed carbon is exported to the ASP . The low recalcitrance of the DOM derived from the slowly sinking Phaeocystis is likely to be associated with rapid mineralization of the organic matter in the water column [79, 80]. Laboratory studies have revealed that carbohydrates derived from Phaeocystis were readily degraded by heterotrophic bacteria . Together, P. antarctica-dominated bloom may contribute relatively uniform and readily degradable organic matter, which can be efficiently degraded by a few specialized bacterial clades, contributing to high rates of remineralization of fixed carbon.
In this study, genomic reconstruction and in situ gene expression analyses were used to characterize the metabolic contributions of dominant bacteria native to the ASP to understand carbon remineralization processes. This analysis allowed us to evaluate the major findings of previous studies of the Southern Ocean at a finer taxonomic level. The metabolic potential of dominant bacterial clades of Bacteroidetes and Gammaproteobacteria can be used to predict temporal community succession linked to the availability of substrates derived from phases of the P. antarctica bloom. Global warming has resulted in compositional changes in phytoplankton from Phaeocystis to diatoms in Antarctic polynyas [12–14]. This change might be accompanied by bacterioplankton community changes, which may affect carbon remineralization and thus carbon sequestration. Repeated parallel studies of the carbon remineralization potential of bacterioplankton in diatom-dominated polynyas are required to predict global warming-related changes in biogeochemical cycles in polynyas.
We thank the captain and crew members of the Korean research icebreaker Araon, for their outstanding assistance during the cruises.
This work was supported by National Research Foundation of Korea grants (NRF-2015M3D3A1A01064881, NRF- 2015R1A4A1041869 and NRF-2018R1A2B6008861) funded by the Ministry of Science and ICT and a grant (PE19060) from the Korea Polar Research Institute. S.-J.K was supported by a grant from National Research Foundation of Korea grants (NRF-2016R1C1B1010946).
Availability of data and materials
Pyrosequencing data of 16S rRNA gene have been deposited in the NCBI Sequence Read Archive (SRA) under accession numbers SRX1887715-SRX1887719. Selected genome bins are available at IMG under accession number 3300009842. All of the cDNA reads are deposited in the NCBI SRA under accession numbers SRX1887523, SRX1887527, and SRX1887528.
S-JK, J-GK, and S-KR designed and performed the research. S-JK, J-GK, S-JP, J-HG, M-J, W-HC, YH, and S-KR analyzed the data. J-GK, S-HL, E-JY, JP, and JJ undertook the field work and processing of samples. J-CC, FR, and J-HH discussed the results and commented on the manuscript. S-JK, J-GK, ELM, and S-KR wrote the paper. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kim BK, Joo H, Song HJ, Yang EJ, Lee SH, Hahm D, Rhee TS, Lee SH. Large seasonal variation in phytoplankton production in the Amundsen Sea. Polar Biol. 2015;38(3):319–31.View ArticleGoogle Scholar
- Lee SH, Kim BK, Yun MS, Joo H, Yang EJ, Kim YN, Shin HC, Lee S. Spatial distribution of phytoplankton productivity in the Amundsen Sea, Antarctica. Polar Biol. 2012;35(11):1721–33.View ArticleGoogle Scholar
- Stefels J, Gieskes WW, Dijkhuizen L. Intriguing functionality of the production and conversion of DMSP in Phaeocystis sp, Biological and environmental chemistry of DMSP and related sulfonium compounds: Springer; 1996. p. 305–15.Google Scholar
- Riegman R, Noordeloos AA, Cadée GC. Phaeocystis blooms and eutrophication of the continental coastal zones of the North Sea. Mar Biol. 1992;112(3):479–84.View ArticleGoogle Scholar
- Schoemann V, Becquevort S, Stefels J, Rousseau V, Lancelot C. Phaeocystis blooms in the global ocean and their controlling mechanisms: a review. J Sea Res. 2005;53(1):43–66.View ArticleGoogle Scholar
- Yager PL, Sherrell L, Stammerjohn SE, Alderkamp A-C, Schofield O, Abrahamsen EP, Arrigo KR, Bertilsson S, Garay D, Guerrero R. ASPIRE: the Amundsen Sea Polynya international research expedition. Oceanography. 2012;25(3):40–53.View ArticleGoogle Scholar
- Henson SA, Sanders R, Madsen E. Global patterns in efficiency of particulate organic carbon export and transfer to the deep ocean. Glob Biogeochem Cy. 2012;26(1):GB1028.View ArticleGoogle Scholar
- Ducklow HW, Wilson SE, Post AF, Stammerjohn SE, Erickson M, Lee S, Lowry KE, Sherrell RM, Yager PL. Particle flux on the continental shelf in the Amundsen Sea Polynya and Western Antarctic Peninsula. Elem Sci Anthrepol. 2015;3(1):000046.View ArticleGoogle Scholar
- Hyun JH, Kim SH, Yang EJ, Choi A, Lee SH. Biomass, production, and control of heterotrophic bacterioplankton during a late phytoplankton bloom in the Amundsen Sea Polynya, Antarctica. Deep Sea Res. 2015;123:102–12.View ArticleGoogle Scholar
- Rinta-Kanto JM, Sun S, Sharma S, Kiene RP, Moran MA. Bacterial community transcription patterns during a marine phytoplankton bloom. Environ Microbiol. 2012;14(1):228–39.View ArticleGoogle Scholar
- Lee S, Hwang J, Ducklow HW, Hahm D, Lee SH, Kim D, Hyun JH, Park J, Ha HK, Kim TW. Evidence of minimal carbon sequestration in the productive Amundsen Sea polynya. Geophys Res Lett. 2017;44(15):7892–9.View ArticleGoogle Scholar
- Arrigo KR, Lowry KE, van Dijken GL. Annual changes in sea ice and phytoplankton in polynyas of the Amundsen Sea, Antarctica. Deep Sea Res. 2012;71:5–15.View ArticleGoogle Scholar
- Tortell PD, Payne CD, Li Y, Trimborn S, Rost B, Smith WO, Riesselman C, Dunbar RB, Sedwick P. DiTullio GR: CO2 sensitivity of Southern Ocean phytoplankton. Geophys Res Lett. 2008:35(4).Google Scholar
- Arrigo KR, Van Dijken GL. Phytoplankton dynamics within 37 Antarctic coastal polynya systems. J Geophys Res. 2003;108(C8):3271.View ArticleGoogle Scholar
- Alderkamp A-C, Buma AG, van Rijssel M. The carbohydrates of Phaeocystis and their degradation in the microbial food web. In: Phaeocystis, major link in the biogeochemical cycling of climate-relevant elements: Springer; 2007. p. 99–118.Google Scholar
- DeJong HB, Dunbar RB, Koweek DA, Mucciarone DA, Bercovici SK, Hansell DA. Net community production and carbon export during the late summer in the Ross Sea, Antarctica. Global Biogeochem Cycles. 2017;31(3):473–91.View ArticleGoogle Scholar
- Reigstad M, Wassmann P. Does Phaeocystis spp. contribute significantly to vertical export of organic carbon? Biogeochemistry. 2007;83(1–3):217–34.View ArticleGoogle Scholar
- Claire Horner-Devine M, Leibold MA, Smith VH, Bohannan BJ. Bacterial diversity patterns along a gradient of primary productivity. Ecol Lett. 2003;6(7):613–22.View ArticleGoogle Scholar
- Galand PE, Casamayor EO, Kirchman DL, Potvin M, Lovejoy C. Unique archaeal assemblages in the Arctic Ocean unveiled by massively parallel tag sequencing. ISME J. 2009;3(7):860–9.View ArticleGoogle Scholar
- Ghiglione JF, Murray AE. Pronounced summer to winter differences and higher wintertime richness in coastal Antarctic marine bacterioplankton. Environ Microbiol. 2012;14(3):617–29.View ArticleGoogle Scholar
- Kim JG, Park SJ, Quan ZX, Jung MY, Cha IT, Kim SJ, Kim KH, Yang EJ, Kim YN, Lee SH, et al. Unveiling abundance and distribution of planktonic bacteria and archaea in a polynya in Amundsen Sea, Antarctica. Environ Microbiol. 2014;16(6):1566–78.View ArticleGoogle Scholar
- Delmont TO, Eren AM, Vineis JH, Post AF. Genome reconstructions indicate the partitioning of ecological functions inside a phytoplankton bloom in the Amundsen Sea, Antarctica. Front Microbiol. 2015;6:1090.PubMedPubMed CentralGoogle Scholar
- Teeling H, Fuchs BM, Bennke CM, Kruger K, Chafee M, Kappelmann L, Reintjes G, Waldmann J, Quast C, Glockner FO, et al. Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms. Elife. 2016;5:e11888.View ArticleGoogle Scholar
- Castelle CJ, Hug LA, Wrighton KC, Thomas BC, Williams KH, Wu D, Tringe SG, Singer SW, Eisen JA, Banfield JF. Extraordinary phylogenetic diversity and metabolic versatility in aquifer sediment. Nat Commun. 2013;4:2120.View ArticleGoogle Scholar
- Seitz KW, Lazar CS, Hinrichs KU, Teske AP, Baker BJ. Genomic reconstruction of a novel, deeply branched sediment archaeal phylum with pathways for acetogenesis and sulfur reduction. ISME J. 2016;10:1696–705.View ArticleGoogle Scholar
- Wrighton KC, Thomas BC, Sharon I, Miller CS, Castelle CJ, VerBerkmoes NC, Wilkins MJ, Hettich RL, Lipton MS, Williams KH, et al. Fermentation, hydrogen, and sulfur metabolism in multiple uncultivated bacterial phyla. Science. 2012;337(6102):1661–5.View ArticleGoogle Scholar
- Klindworth A, Mann AJ, Huang S, Wichels A, Quast C, Waldmann J, Teeling H, Glockner FO. Diversity and activity of marine bacterioplankton during a diatom bloom in the North Sea assessed by total RNA and pyrotag sequencing. Mar Genomics. 2014;18:185–92.View ArticleGoogle Scholar
- Lee Y, Yang EJ, Park J, Jung J, Kim TW, Lee S. Physical-biological coupling in the Amundsen Sea, Antarctica: influence of physical factors on phytoplankton community structure and biomass. Deep Sea Res I. 2016;117:51–60.View ArticleGoogle Scholar
- Hurt RA, Qiu X, Wu L, Roh Y, Palumbo AV, Tiedje JM, Zhou J. Simultaneous recovery of RNA and DNA from soils and sediments. Appl Environ Microbiol. 2001;67(10):4495–503.View ArticleGoogle Scholar
- Hur M, Kim Y, Song HR, Kim JM, Choi YI, Yi H. Effect of genetically modified poplars on soil microbial communities during the phytoremediation of waste mine tailings. Appl Environ Microbiol. 2011;77(21):7611–9.View ArticleGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.View ArticleGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069–72.View ArticleGoogle Scholar
- Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.View ArticleGoogle Scholar
- Peng Y, Leung HC, Yiu SM, Chin FY. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28(11):1420–8.View ArticleGoogle Scholar
- Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27(4):578–9.View ArticleGoogle Scholar
- Kim D, Song L, Breitwieser FP, Salzberg SL. Centrifuge: rapid and sensitive classification of metagenomic sequences. Genome Res. 2016.Google Scholar
- Noguchi H, Taniguchi T, Itoh T. MetaGeneAnnotator: detecting species-specific patterns of ribosomal binding site for precise gene prediction in anonymous prokaryotic and phage genomes. DNA Res. 2008;15(6):387–96.View ArticleGoogle Scholar
- Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35(9):3100–8.View ArticleGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.View ArticleGoogle Scholar
- Tang K, Jiao N, Liu K, Zhang Y, Li S. Distribution and functions of TonB-dependent transporters in marine bacteria and environments: implications for dissolved organic matter utilization. PLoS One. 2012;7(7):e41204.View ArticleGoogle Scholar
- Saier MH Jr, Reddy VS, Tamang DG, Vastermark A. The transporter classification database. Nucleic Acids Res. 2014;42(Database issue):D251–8.View ArticleGoogle Scholar
- Rawlings ND, Barrett AJ, Finn R. Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2016;44(D1):D343–50.View ArticleGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.View ArticleGoogle Scholar
- Albertsen M, Hugenholtz P, Skarshewski A, Nielsen KL, Tyson GW, Nielsen PH. Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes. Nat Biotechnol. 2013;31(6):533–8.View ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.View ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome Project Data Processing S: the sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.View ArticleGoogle Scholar
- Lander ES, Waterman MS. Genomic mapping by fingerprinting random clones: a mathematical analysis. Genomics. 1988;2(3):231–9.View ArticleGoogle Scholar
- Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25(7):1043–55.View ArticleGoogle Scholar
- Kim SJ, Park SJ, Cha IT, Min D, Kim JS, Chung WH, Chae JC, Jeon CO, Rhee SK. Metabolic versatility of toluene-degrading, iron-reducing bacteria in tidal flat sediment, characterized by stable isotope probing-based metagenomic analysis. Environ Microbiol. 2014;16(1):189–204.View ArticleGoogle Scholar
- Yilmaz P, Parfrey LW, Yarza P, Gerken J, Pruesse E, Quast C, Schweer T, Peplies J, Ludwig W, Glockner FO. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014;42(Database issue):D643–8.View ArticleGoogle Scholar
- Szymanski M, Barciszewska MZ, Erdmann VA, Barciszewski J. 5S Ribosomal RNA Database. Nucleic Acids Res. 2002;30(1):176–8.View ArticleGoogle Scholar
- Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2014;12(1):59.View ArticleGoogle Scholar
- Franzosa EA, McIver LJ, Rahnavard G, Thompson LR, Schirmer M, Weingart G, Lipson KS, Knight R, Caporaso JG, Segata N. Species-level functional profiling of metagenomes and metatranscriptomes. Nat Methods. 2018;1.Google Scholar
- Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glockner FO. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41(1):e1.View ArticleGoogle Scholar
- Martens EC, Koropatkin NM, Smith TJ, Gordon JI. Complex glycan catabolism by the human gut microbiota: the Bacteroidetes Sus-like paradigm. J Biol Chem. 2009;284(37):24673–7.View ArticleGoogle Scholar
- Yancey PH. Organic osmolytes as compatible, metabolic and counteracting cytoprotectants in high osmolarity and other stresses. J Exp Biol. 2005;208(15):2819–30.View ArticleGoogle Scholar
- Delmont TO, Hammar KM, Ducklow HW, Yager PL, Post AF. Phaeocystis antarctica blooms strongly influence bacterial community structures in the Amundsen Sea polynya. Front Microbiol. 2014;5:646.View ArticleGoogle Scholar
- Gosink JJ, Woese CR, Staley JT. Polaribacter gen. nov., with three new species, P. irgensii sp. nov., P. franzmannii sp. nov. and P. filamentus sp. nov., gas vacuolate polar marine bacteria of the Cytophaga-Flavobacterium-Bacteroides group and reclassification of ‘Flectobacillus glomerat’ as Polaribacter glomeratus comb. nov. Int J Syst Bacteriol. 1998;48:223–35.View ArticleGoogle Scholar
- Choi S-B, Kim J-G, Jung M-Y, Kim S-J, Min U-G, Si O-J, Park S-J, Hwang CY, Park J, Lee S. Cultivation and biochemical characterization of heterotrophic bacteria associated with phytoplankton bloom in the Amundsen Sea polynya, Antarctica. Deep Sea Res II. 2016;123:126–34.View ArticleGoogle Scholar
- Grzymski JJ, Carter BJ, DeLong EF, Feldman RA, Ghadiri A, Murray AE. Comparative genomics of DNA fragments from six Antarctic marine planktonic bacteria. Appl Environ Microbiol. 2006;72(2):1532–41.View ArticleGoogle Scholar
- Teeling H, Fuchs BM, Becher D, Klockow C, Gardebrecht A, Bennke CM, Kassabgy M, Huang S, Mann AJ, Waldmann J, et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science. 2012;336(6081):608–11.View ArticleGoogle Scholar
- Van Boekel W, Hansen F, Riegman R, Bak R. Lysis-induced decline of a Phaeocystis spring bloom and coupling with the microbial foodweb. Mar Ecol Prog Ser. 1992:269–76.Google Scholar
- Brussaard C, Kuipers B, Veldhuis M. A mesocosm study of Phaeocystis globosa population dynamics: I. Regulatory role of viruses in bloom control. Harmful Algae. 2005;4(5):859–74.View ArticleGoogle Scholar
- Brussaard C, Riegman R, Noordeloos A, Cadée G, Witte H, Kop A, Nieuwland G, Van Duyl F, Bak R. Effects of grazing, sedimentation and phytoplankton cell lysis on the structure of a coastal pelagic food web. Mar Ecol Prog Ser. 1995:259–71.Google Scholar
- Mann AJ, Hahnke RL, Huang S, Werner J, Xing P, Barbeyron T, Huettel B, Stuber K, Reinhardt R, Harder J, et al. The genome of the alga-associated marine flavobacterium Formosa agariphila KMM 3901T reveals a broad potential for degradation of algal polysaccharides. Appl Environ Microbiol. 2013;79(21):6813–22.View ArticleGoogle Scholar
- Gugi B, Le Costaouec T, Burel C, Lerouge P, Helbert W, Bardor M. Diatom-specific oligosaccharide and polysaccharide structures help to unravel biosynthetic capabilities in Diatoms. Mar Drugs. 2015;13(9):5993–6018.View ArticleGoogle Scholar
- Keller MD, Bellows WK, Guillard RR. Dimethyl sulfide production in marine phytoplankton. In: Saltzman E, Copper W, editors. Biogenic sulfur in the environment, vol. 393. Washington, D.C: American Chemical society; 1989. p. 167–82.View ArticleGoogle Scholar
- Stefels J, van Leeuwe MA. Effects of iron and light stress on the biochemical composition of Antarctic Phaeocystis sp.(Prymnesiophyceae). I. Intracellular DMSP concentrations. J Phycol. 1998;34(3):486–95.View ArticleGoogle Scholar
- Buchan A, LeCleir GR, Gulvik CA, Gonzalez JM. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat Rev Microbiol. 2014;12(10):686–98.View ArticleGoogle Scholar
- Curson AR, Todd JD, Sullivan MJ, Johnston AW. Catabolism of dimethylsulphoniopropionate: microorganisms, enzymes and genes. Nat Rev Microbiol. 2011;9(12):849–59.View ArticleGoogle Scholar
- Tortell PD, Long MC, Payne CD, Alderkamp A-C, Dutrieux P, Arrigo KR. Spatial distribution of pCO2, ΔO2/Ar and dimethylsulfide (DMS) in polynya waters and the sea ice zone of the Amundsen Sea, Antarctica. Deep Sea Res. 2012;71:77–93.View ArticleGoogle Scholar
- Stefels J, Dijkhuizen L. Characteristics of DMSP-lyase in Phaeocystis sp. (Prymnesiophyceae). Mar Ecol Prog Ser. 1996;131(1):307–13.View ArticleGoogle Scholar
- Janse I, Zwart G, van der Maarel MJ, Gottschal JC. Composition of the bacterial community degrading Phaeocystis mucopolysaccharides in enrichment cultures. Aquat Microb Ecol. 2000;22(2):119–33.View ArticleGoogle Scholar
- Voget S, Wemheuer B, Brinkhoff T, Vollmers J, Dietrich S, Giebel HA, Beardsley C, Sardemann C, Bakenhus I, Billerbeck S, et al. Adaptation of an abundant Roseobacter RCA organism to pelagic systems revealed by genomic and transcriptomic analyses. ISME J. 2015;9(2):371–84.View ArticleGoogle Scholar
- Buchan A, Gonzalez JM, Moran MA. Overview of the marine Roseobacter lineage. Appl Environ Microbiol. 2005;71(10):5665–77.View ArticleGoogle Scholar
- Ahn S, Jung J, Jang IA, Madsen EL, Park W. Role of glyoxylate shunt in oxidative stress response. J Biol Chem. 2016;291(22):11928–38.View ArticleGoogle Scholar
- Bender SJ, Moran DM, McIlvin MR, Zheng H, McCrow JP, Badger J, DiTullio GR, Allen AE, Saito MA. Iron triggers colony formation in Phaeocystis antarctica: connecting molecular mechanisms with iron biogeochemistry. In: Biogeosciences Discuss; 2018. in review.Google Scholar
- Schofield O, Miles T, Alderkamp A-C, Lee S, Haskins C, Rogalsky E, Sipler R, Sherrell RM, Yager PL. In situ phytoplankton distributions in the Amundsen Sea Polynya measured by autonomous gliders. Elem Sci Anth. 2015;3:000073.View ArticleGoogle Scholar
- Wassmann P. Significance of sedimentation for the termination of Phaeocystis blooms. J Mar Syst. 1994;5(1):81–100.View ArticleGoogle Scholar
- Rousseau V, Becquevort S, Parent J-Y, Gasparini S, Daro M-H, Tackx M, Lancelot C. Trophic efficiency of the planktonic food web in a coastal ecosystem dominated by Phaeocystis colonies. J Sea Res. 2000;43(3–4):357–72.View ArticleGoogle Scholar
- Janse I, van Rijssel M, Ottema A, Gottschal JC. Microbial breakdown of Phaeocystis mucopolysaccharides. Limnol Oceanogr. 1999;44(6):1447–57.View ArticleGoogle Scholar
- Xing P, Hahnke RL, Unfried F, Markert S, Huang S, Barbeyron T, Harder J, Becher D, Schweder T, Glockner FO, et al. Niches of two polysaccharide-degrading Polaribacter isolates from the North Sea during a spring diatom bloom. ISME J. 2015;9(6):1410–22.View ArticleGoogle Scholar