Uptake of Dimethylsulfoniopropionate (DMSP) by Natural Microbial Communities of the Great Barrier Reef (GBR), Australia

Dimethylsulfoniopropionate (DMSP) is a key organic sulfur compound that is produced by many phytoplankton and macrophytes and is ubiquitous in marine environments. Following its release into the water column, DMSP is primarily metabolised by heterotrophic bacterioplankton, but recent evidence indicates that non-DMSP producing phytoplankton can also assimilate DMSP from the surrounding environment. In this study, we examined the uptake of DMSP by communities of bacteria and phytoplankton within the waters of the Great Barrier Reef (GBR), Australia. We incubated natural GBR seawater with DMSP and quantified the uptake of DMSP by different fractions of the microbial community (>8 µm, 3–8 µm, <3 µm). We also evaluated how microbial community composition and the abundances of DMSP degrading genes are influenced by elevated dissolved DMSP levels. Our results showed uptake and accumulation of DMSP in all size fractions of the microbial community, with the largest fraction (>8 µm) forming the dominant sink, increasing in particulate DMSP by 44–115% upon DMSP enrichment. Longer-term incubations showed however, that DMSP retention was short lived (<24 h) and microbial responses to DMSP enrichment differed depending on the community carbon and sulfur demand. The response of the microbial communities from inside the reef indicated a preference towards cleaving DMSP into the climatically active aerosol dimethyl sulfide (DMS), whereas communities from the outer reef were sulfur and carbon limited, resulting in more DMSP being utilised by the cells. Our results show that DMSP uptake is shared across members of the microbial community, highlighting larger phytoplankton taxa as potentially relevant DMSP reservoirs and provide new information on sulfur cycling as a function of community metabolism in deeper, oligotrophic GBR waters.


Introduction
Dimethylsulfoniopropionate (DMSP) is an important and ubiquitous organic sulfur compound in the marine environment [1]. Phytoplankton are the dominant producers of DMSP, although intracellular concentrations vary from species to species, ranging from extremely low in cyanobacteria and most diatoms to very high,~300 pg DMSP cell −1 , in haptophytes and dinoflagellates [2][3][4][5]. In DMSP-producing phytoplankton, DMSP is central to cell metabolism, satisfying up to 71% of the sulfur and 15% of the cells carbon demand [3,6]. Likewise, for DMSP-associated bacterioplankton, DMSP can account for up to 95% of the sulfur and 15% of the carbon demand [7]. Several physiological roles have been attributed to DMSP, including osmoregulation, cryoprotection and antioxidation [8][9][10], while ecologically, DMSP has been shown to function as a grazing deterrent and have bacteriocidal properties [11][12][13].
Intracellular or particulate DMSP (DMSPp) is transported from the synthesising cell into the external seawater environment via exudation, viral lysis or cell senescence [14][15][16]. Once of DMSP in deep water tropical ecosystems. To better understand microbial DMSP sinks and sources within the oligrotrophic waters of the GBR, we investigated DMSP uptake among natural microbial communities from two distinct reef locations, with the aim to quantify and compare the uptake of DMSP by the bacterial and phytoplankton communities using a size fractionation approach and to evaluate DMSP-induced changes in microbial composition and DMSP-degradation gene abundance over time.

Study Site and Experimental Design
This study was conducted during an oceanographic voyage on the RV Investigator in October 2016. Two sites within the Great Barrier Reef were chosen: a southern open ocean (depth ~400 m) site (outer reef site; OR) (−20°8′ S, 153°1′ E) and a shallower (~45 m) northern inner-reef site (−18°6′ S, 146°9′ E), located between the coast and the reefs (inner reef site; IR). The inner reef site was located ~13 km from the closest reefs of John Brewer Reef, Rib Reef and Fore and Aft Reefs; the deeper outer reef was located ~25 km from Bills Reef (Figure 1). At both sites, surface seawater samples (5 m) were collected with Niskin bottles from CTD casts and filtered through a 210 µm mesh to remove large grazers before being sub-sampled into experimental units. Samples for nutrients, DMS and microbial communities (16S and 18S) were taken from initial samples prior to incubations. Seawater temperature and salinity were taken from the CTD data. At each site, two separate incubation experiments were conducted: (1) a short-term study (7 h) to quantify DMSP uptake by different size fractions of the microbial community, and (2) a longer (120 h) experiment to investigate the influence of DMSP on community structure and DMSP degradation gene regulation. At each site, two separate incubation experiments were conducted: (1) a short-term study (7 h) to quantify DMSP uptake by different size fractions of the microbial community, and (2) a longer (120 h) experiment to investigate the influence of DMSP on community structure and DMSP degradation gene regulation.

Experiment 1: Quantification of DMSP Uptake by Size Fractions of the Microbial Community
To study the uptake of DMSP by the microbial community of the Great Barrier Reef, 28 L of seawater was collected and transferred to 4 L-polycarbonate bottles with no headspace. Bottles were separated into three treatments, in triplicate: controls (unamended seawater), DMSP addition (+10 nM DMSP) and fixed (dead). The fixed sample was incubated with glutaraldehyde (final concentration 1%) for 1 h, before DMSP (10 nM Microorganisms 2021, 9,1891 4 of 23 final concentration) was added, to control for any passive uptake of DMSP by the cells. Following DMSP addition, the bottles were closed with screw caps and inverted gently to dissolve the added compound. All flasks were then incubated for up to 7 h in an ondeck flow-through incubator. At four time points (T0, T2, T4 and T7 h) each bottle was subsampled in triplicate for quantification of DMSP total and dissolved (DMSPt, DMSPd) concentrations, either as a whole water sample or a gravity filtered (GF/F, nominal pore size 0.7 µm) sample (3 mL), respectively. Sub-samples were acidified with H 2 SO 4 (0.5%, pH~1.1) and kept at −80 • C until processing [50]. After 7 h, two litres from each replicate was size-fractionated via serial filtration onto 8, 3 and 0.22 µm polycarbonate filters for analysis of particulate DMSP (DMSPp) and the filters were snap frozen and kept at −80 • C until analysis. Filters were defrosted and hydrolysed with NaOH prior to GC-FPD analysis, described below. During filtration, aliquots (800 µL) were taken from the whole water fraction, fixed with glutaraldehyde (final concentration 1%) and snap frozen for determination of microbial abundance via flow cytometry (see below).

Experiment 2: Effect of DMSP on Community Structure and Gene Regulation
To study microbial community changes caused by DMSP enrichment, 30 L of water were divided across 15 × 2 L polycarbonate bottles (no headspace). Bottles were separated into two treatments: six controls (unamended) and nine amended with DMSP (10 nM final concentration). Bottles were then incubated in an on-deck flow-through incubator until sampling at three time points (T24, T72 and T120 h). At each time point, subsamples (3 mL) for DMSP total and dissolved (gravity filtration, GF/F) were taken and acidified as described above.

Quantification of Biogenic Sulfur Compounds
Biogenic sulfur compounds (DMSPt, DMSPd) were quantified as total DMS after conversion with 100 mg of NaOH and measured using a gas chromatograph (GC-2010 Plus, Shimadzu, Kyoto, Japan) coupled with a flame photometric detector (FPD). Samples were purged with He (60 mL/min for 4 min) while cryo-trapped in liquid N 2 and subsequently eluted onto a capillary column (DB-1, Agilent; injector: 120 • C, column: 110 • C, FPD: 130 • C, column flow: 2.1 mL min −1 ). For experiment 1, particulate DMSP (DMSPp) was measured directly on the filtered size fractions, where filters were placed into vials with NaOH and quantified as total DMS as described above. Particulate DMSP for experiment 2 was calculated by subtracting DMSPd from DMSPt. Integrated peak areas for each sample were quantified against a calibration curve obtained from solutions prepared with known amounts of DMSP hydrolysed with NaOH and analysed using the same method as the samples. For determination of DMSP lyase activity (DLA) rates of the starting communities, subsamples (600 mL) were filtered first through a GF/C filter (nominal pore size 1.2 µm) for the phytoplankton fraction rates (DLAp), and then 300 mL of the phytoplankton-free sample filtrate was filtered onto a 0.22 µm polycarbonate filter for bacterial DLA (DLAb). Both filters were placed into cryotubes, flash frozen with in liquid nitrogen and stored in −80 • C until analysis. DMSP lyase activity was determined via direct injection of 500 µL of headspace (column flow: 3.66 mL min −1 , FPD: 160 • C) as previously described [51].

Quantification of Microbial Composition
Quantification of microbial populations via flow cytometry was performed using a Beckman Coulter Inc flow cytometer (Indianapolis, IN, USA). Populations of Prochlorococcus, Synechococcus and picoeukaryotes were discriminated using side scatter (SSC) and red and orange fluorescence [52]. Samples for quantification of heterotrophic bacteria were stained with SYBR Green I nucleic acid stain (1:10,000 final dilution; Invitrogen, Waltham, MA, USA) and populations were discriminated according to green fluorescence and side scatter properties [53,54]. Data were analysed using CytExpert software (Beckman Coulter, Brea, CA, USA).

Nutrients
Nutrients, measured in the initial sample (i.e., at time of collection) and at T24, T72 and T120 h from the incubations were carried out on board as previously described [55]. Data were subsequently obtained from the Australian Ocean Data Network (AODN).

Microbial Community Composition and DMSP-Related Genes
In both incubation experiments, the remaining sample (~1.9 L) at the final time-point was filtered through a 0.22 µm polycarbonate filter, which was subsequently snapped frozen and kept at −80 • C. DNA was extracted using the Powerwater DNA isolation kit (QIAGEN, Hilden, Germany) following manufacturer's instructions, and the DNA concentration was quantified using a Nanodrop spectrophotometer (Nanodrop 2000, Thermo Fisher Scientific, Waltham, MA, USA). Extracted DNA was stored at −20 • C until analysis. Prior to analysis the DNA was split into two fractions, with one used for the analysis of prokaryotic and eukaryotic microbial diversity using (16S and 18S rRNA) amplicon sequencing, and the second fraction used for the quantification of DMSP degrading genes using real time quantitative Polymerase Chain Reaction (qPCR) described below.
Gene sequencing reads were processed to an abundance table of Amplicon Sequence Variants (ASVs) from demultiplexed paired fastq files using DaDa2 (version 1.14.1 [56]). Briefly, after visual inspection of the read qualities the primers were removed from the forward and reverse reads using Cutadapt ( [57]). Sequences were then filtered for N and trimmed with the following specific parameters Bacterial 16S: 'truncLen = c(255,250), maxN = 0, maxEE = c(2,6), truncQ = 6, rm.phix = TRUE'. These parameters were identical for Eukaryote 18S with the exception of truncLen = c (200,195). Error rates were estimated from 1e8 bases with the option 'max consist = 20'. Reads were then dereplicated and merged using the pool = "pseudo" option and potential chimeras were then removed using 'removeBimeraDenovo' with 'minFoldParentOverAbundance' = 4. The ASVs in the resulting abundance table were assigned to taxonomic lineages using 'assignTaxonomy' using Silva SSU taxonomic training data formatted for DADA2 (DOI: 10.5281/zenodo.3986799; Silva version 138) for 16S and PR2 version 4.13.0 for 18S (https://zenodo.org/record/3362765#.YL2O9y0RrUI, accessed on 19 April 2021). The dataset was further cleaned by removing sequence variants (ASVs) identified as chloroplasts, mitochondria and ASVs (abundance cut-off of 0.03%). Sequences were then rarefied to the same depth to remove the effect of sampling effort upon analysis. Raw FASTQ data files for the 16S rRNA and 18S rRNA assay were deposited in NCBI Sequence Read Archive (SRA) under Bioproject number PRJNA755150.

Data Analysis
Differences in microbial abundances from flow cytometry and particulate DMSP for each size fraction were analysed by two-factor ANOVA, with site and experiment (microbial abundance) or site and treatment (DMPSp), as fixed factors. A Tukey's post hoc test was performed to locate differences and account for multiple comparisons. Using a resemblance matrix based on Euclidean distance, a two factor Permutational Multivariate Analysis of Variance (PERMANOVA) with pair-wise comparisons was used to determine significant changes or differences in DMSPt and DMSPd concentrations over time and between treatments and to analyse the differences between treatments over time for the relative abundance of dmdA/A1, dmdA/Dall and DddP genes data. Differences in bacterioplankton (16S) and phytoplankton (18S) community composition between treatments over time (exp. 2) were analysed using PRIMER v6 statistical package [61]. Data were square root transformed and a resemblance matrix generated using Bray-Curtis similarity. Data were analysed using non-parametric Analysis of Similarities (ANOSIM) and discriminatory microbes identified using a two-way nested (time within treatment) Analysis of Similarity Percentages (SIMPER). Relationships between relative abundance of DMSP degrading genes and bacterial abundance (including all time points-T0, T24, T72 and T120 h and treatments (controls and DMSP-enriched)) were analysed using linear regression. All data analyses were performed in R 4.0.4 [62], using packages tidyverse [63], vegan [64] and ggpubr [65].

Characteristics of Initial Water Masses
The physicochemical characteristics of the initial GBR water masses at the time of sampling were similar across both reef sites and experiments (Table 1). Initial water temperatures ranged between 26.3 and 27.1 • C for the two sites, while salinity did not vary (~35.4). Dimethyl sulfide (DMS) concentrations were~0.9 nM at both sites, except for the inner reef site for experiment 1, where it was below the detection limit. Nutrient concentrations were generally very low, with nitrates (NO X ) consistently low (≤0.02 µM) at both sites and for both incubation experiments. Phosphate concentrations were similar at both sites, ranging from 0.038 to 0.050 µM in the inner reef, and 0.049 to 0.066 µM in the outer reef, while silicate concentrations were 0.798 and 0.882 µM in the outer reef compared with 0.500 and 0.623 µM at the inner reef site for the two incubation experiments (Table 1). Picoplankton abundances varied between sites and experiments, with generally higher abundances at the inner reef site ( Table 2). There were higher abundances (F = 13.6, p = 0.006) of Synechococcus inside the reef compared with outside the reef (differing by an order of magnitude), while Prochlorococcus abundances were highest for experiment 1 (F = 23.54, p = 0.001) at the inner reef site (6.56 ± 2.7 × 10 5 cell mL −1 ) and experiment 2 (5.59 ± 0.53 × 10 5 cell mL −1 ) at the outer reef site. Heterotrophic bacterial abundances were higher at the inner reef site (0.79 ± 0.16 × 10 5 cell mL −1 ) compared to the outer reef site (0.30 ± 0.05 × 10 5 cell mL −1 ), with differences in abundance between experiments for the inner reef site (F = 25.8, p < 0.0001) and Picoeukaryote abundance was more than double inside (1490 ± 250 cell mL −1 ) than outside the reef (F = 37.8, p = 0.0002). Table 2. Microbial abundance (cells mL −1 ) measured using flow cytometry. Data represent mean ± SD (n = 3). Superscript letters denote statistical differences at p < 0.05 between sites and experiments. Unique lettering indicates significantly different from all other treatments, whereas shared superscript letters denote non-significance between treatments.

Inner Reef
Outer Reef   As expected, DMSP enrichment was immediately evident in the dissolved fraction, being significantly higher in the +DMSP treatment compared to the controls for both the inner (F = 37.4, p = 0.001) and outer (F = 217, p = 0.001) reef sites over time ( Figure 2B). For the inner reef site, concentrations of DMSPd were 4.7 and 15.6 nM in the control and +DMSP samples, respectively, and for the control, DMSPd concentrations at the final time point were only marginally higher (6 nM) than at the outset (T0 h)-the abnormally high value of 10.88 nM at 4 h was likely due to filtration artifacts. For the outer reef site, DMSP remained constant in both the control and +DMSP treatments over time. There was no As expected, DMSP enrichment was immediately evident in the dissolved fraction, being significantly higher in the +DMSP treatment compared to the controls for both the inner (F = 37.4, p = 0.001) and outer (F = 217, p = 0.001) reef sites over time ( Figure 2B). For the inner reef site, concentrations of DMSPd were 4.7 and 15.6 nM in the control and +DMSP samples, respectively, and for the control, DMSPd concentrations at the final time point were only marginally higher (6 nM) than at the outset (T0 h)-the abnormally high value of 10.88 nM at 4 h was likely due to filtration artifacts. For the outer reef site, DMSP remained constant in both the control and +DMSP treatments over time. There was no change in DMSPt and DMSPd concentrations over time for the fixed samples (data not shown).
Particulate DMSP levels varied according to size fraction and treatment ( Figure 2C). Within Inner reef samples, the largest size fraction (>8 µm) contained the highest DMSPp concentrations (2.1 ± 0.5 nM), followed by the medium fraction (>3 µm; 0.99 ± 0.26 nM), and then the picoplankton (<3 µm; 0.28 ± 0.18 nM) fraction ( Figure 2C). In contrast, DMSPp levels in the outer reef samples did not differ significantly across all three fractions in the absence of enrichment, but addition of DMSP led to significant increases (F = 80.43, p = 0.001) in DMSPp in the largest (>8 µm) fraction ( Figure 2C). For the smallest fraction (0.22-3.0 µm), DMSPp levels were significantly higher in the outer reef community compared to the inner reef (F = 16.9, p = 0.0026), but with no effect of DMSP enrichment. In the 3-8 µm fraction, DMSPp levels were higher in the inner reef (F = 20.2, p = 0.0015), but enrichment did not result in a further increase in DMSPp. For the largest fraction (>8 µm), the inner reef had higher levels of DMSPp, and enrichment increased DMSPp at both sites (F = 8.47, p = 0.008). There was no change in DMSPp in fixed cells with DMSP enrichment, which remained low in all fractions ( Figure 2C).
Rates of DMSP lyase activity (DLA) differed between sites ( Figure 3). The inner reef was dominated by phytoplankton DLA (DLAp), with rates three-times higher (0.06 µM DMS min −1 ) than those of the bacterial (DLAb) fraction (0.017 µM DMS min −1 ; Figure 3A). In contrast, a reverse pattern was found for the outer reef site, where DLAb rates (0.23 µM DMS min −1 ) were three times higher than DLAp rates (0.07 µM DMS min −1 ; Figure 3B). Interestingly, DLAb rates in the outer reef were more than an order of magnitude higher than those measured for the inner reef, while DLAp rates were similar across both sites.
Of the DMSP catabolising genes tested, DmdA/Dall occurred in the highest abundance at both sites, with 7-10 × 10 −2 copies per 16S copies for the inner reef site and 6.5-7 × 10 −2 copies per 16S copies for the outer reef site ( Figure 3C). The common lyase gene DddP was approximately twice as abundant (6.7 × 10 −3 copies/16S) in the bacterial population of the inner reef compared with the outer reef site (3.6 × 10 −3 copies/16S; Figure 3D). As with DddP, the abundance of the gene DmdA/A1 was more abundant in the bacterial populations of the inner reef (9.4 × 10 −3 copies/16S) than the outer reef site (2.5 × 10 −3 copies/16S; Figure 3E). change in DMSPt and DMSPd concentrations over time for the fixed samples (data not shown).
Particulate DMSP levels varied according to size fraction and treatment ( Figure 2C). Within Inner reef samples, the largest size fraction (>8 µm) contained the highest DMSPp concentrations (2.1 ± 0.5 nM), followed by the medium fraction (>3 µm; 0.99 ± 0.26 nM), and then the picoplankton (<3 µm; 0.28 ± 0.18 nM) fraction ( Figure 2C). In contrast, DMSPp levels in the outer reef samples did not differ significantly across all three fractions in the absence of enrichment, but addition of DMSP led to significant increases (F = 80.43, p = 0.001) in DMSPp in the largest (>8 µm) fraction ( Figure 2C). For the smallest fraction (0.22-3.0 µm), DMSPp levels were significantly higher in the outer reef community compared to the inner reef (F = 16.9, p = 0.0026), but with no effect of DMSP enrichment. In the 3-8 µm fraction, DMSPp levels were higher in the inner reef (F = 20.2, p = 0.0015), but enrichment did not result in a further increase in DMSPp. For the largest fraction (>8 µm), the inner reef had higher levels of DMSPp, and enrichment increased DMSPp at both sites (F = 8.47, p = 0.008). There was no change in DMSPp in fixed cells with DMSP enrichment, which remained low in all fractions ( Figure 2C).
Rates of DMSP lyase activity (DLA) differed between sites ( Figure 3). The inner reef was dominated by phytoplankton DLA (DLAp), with rates three-times higher (0.06 µM DMS min −1 ) than those of the bacterial (DLAb) fraction (0.017 µM DMS min −1 ; Figure 3A). In contrast, a reverse pattern was found for the outer reef site, where DLAb rates (0.23 µM DMS min −1 ) were three times higher than DLAp rates (0.07 µM DMS min −1 ; Figure 3B). Interestingly, DLAb rates in the outer reef were more than an order of magnitude higher than those measured for the inner reef, while DLAp rates were similar across both sites.
Of the DMSP catabolising genes tested, DmdA/Dall occurred in the highest abundance at both sites, with 7-10 × 10 −2 copies per 16S copies for the inner reef site and 6.5-7 × 10 −2 copies per 16S copies for the outer reef site ( Figure 3C). The common lyase gene DddP was approximately twice as abundant (6.7 × 10 −3 copies/16S) in the bacterial population of the inner reef compared with the outer reef site (3.6 × 10 −3 copies/16S; Figure 3D). As with DddP, the abundance of the gene DmdA/A1 was more abundant in the bacterial populations of the inner reef (9.4 × 10 −3 copies/16S) than the outer reef site (2.5 × 10 −3 copies/16S; Figure 3E). Bacterioplankton community composition was similar at both sites, with up to eleven genera of >1% of relative abundance making up more than 80% of the community ( Figure  4A). Reef waters were dominated by Synechococcales (44-47%) and SAR11 clade (16-22%). The remaining genera present in >1% relative included abundance consisted of genera from the orders of SAR86, Puniceispirillales, Flavobacteriales, Rhodospirillales, Bacterioplankton community composition was similar at both sites, with up to eleven genera of >1% of relative abundance making up more than 80% of the community ( Figure 4A). Reef waters were dominated by Synechococcales (44-47%) and SAR11 clade (16-22%). The remaining genera present in >1% relative included abundance consisted of genera from the orders of SAR86, Puniceispirillales, Flavobacteriales, Rhodospirillales, Cellvibrionales, Rhodobacterales, Rickettsiales, Alteromonadales and Actinomarinales (1.03%). At both sites, the dominant genus within the Synechococcales order was the photosynthetic cyanobacterium from the genus Prochlorococcus, while the other important photosynthetic cyanobacterium Synechococcus was present at 3.6% in the inner reef site, yet <1% at the outer reef. Phytoplankton community composition was similar between sites ( Figure 4B) with dinoflagellates (Dino-Group I-V, 48-56%; Gymnodiniales, 11-17%; Dinophyceae, 3.7-6.8%; Prorocentrales, 3-4.6%; Peridiniales 2-5.6%) as the dominant algal class (65% IR, 85% OR) at both sites. The only notable difference in community composition was the presence of diatoms at the inner reef site, which had a relative abundance of 13%, and were completely absent from the outer reef site.
tosynthetic cyanobacterium from the genus Prochlorococcus, while the other important photosynthetic cyanobacterium Synechococcus was present at 3.6% in the inner reef site, yet <1% at the outer reef. Phytoplankton community composition was similar between sites ( Figure 4B) with dinoflagellates (Dino-Group I-V, 48-56%; Gymnodiniales, 11-17%; Dinophyceae, 3.7-6.8%; Prorocentrales, 3-4.6%; Peridiniales 2-5.6%) as the dominant algal class (65% IR, 85% OR) at both sites. The only notable difference in community composition was the presence of diatoms at the inner reef site, which had a relative abundance of 13%, and were completely absent from the outer reef site. Relative abundance of (A) bacteria and (B) phytoplankton for inner (IR) and outer (OR) reef sites. Graphs display orders with >1% of relative abundance for clarity purposes (n = 1).

Experiment 2: Effect of DMSP on Community Structure and Gene Regulation
To evaluate the effect of DMSP-enrichment on microbial community composition and DMSP catabolising gene abundance in natural seawaters of the GBR, we incubated water samples with trace levels of DMSP (final concentration 10 nM) and monitored changes in microbial communities and sulfur compounds over five days. Macronutrient concentrations over time differed between the inner and outer reef experiments, where the inner reef, phosphate dropped below detection limit by 120 h, and silicate concentrations increased from initial concentrations (Supplementary Table S2). For the outer reef silicate also increased by 120 h, but there was no change in phosphate (Supplementary Table S3). Concentrations of NOX and ammonium remained low throughout.

Experiment 2: Effect of DMSP on Community Structure and Gene Regulation
To evaluate the effect of DMSP-enrichment on microbial community composition and DMSP catabolising gene abundance in natural seawaters of the GBR, we incubated water samples with trace levels of DMSP (final concentration 10 nM) and monitored changes in microbial communities and sulfur compounds over five days. Macronutrient concentrations over time differed between the inner and outer reef experiments, where the inner reef, phosphate dropped below detection limit by 120 h, and silicate concentrations increased from initial concentrations (Supplementary Table S2). For the outer reef silicate also increased by 120 h, but there was no change in phosphate (Supplementary Table S3 Contrasting trends in DMSP degradation gene abundance were evident between the two sites, with significant increases over time for the inner reef community, compared with a stable or declining gene abundance over time for the outer reef community ( Figure  5D-F Figure 5B). Like the inner reef, the highest rate of loss occurred during the first 24 h with 7.74 nM day −1 (controls) and 18.56 nM day −1 (+DMSP). This loss of DMSPd over 120 h equates to a loss of the initial DMSPd of approximately 73% in the controls and~91% in DMSP-enriched bottles. As with DMSPt, the rapid decline meant that initial differences in DMSPd between treatments for the inner (F = 71.8, p = 0.001) and outer (F = 21.3, p = 0.006) reefs disappeared after only 24 h.
Particulate DMSP for the inner reef samples, calculated by subtracting DMSPd from DMSPt ( Figure 5C) remained relatively constant (~12 nM) for the DMSP-enriched samples, while declining significantly from 14 nM to 9 nM in the controls (F = 7.19, p = 0.015) with a loss rate of 1.07 nM day −1 . Conversely, DMSPp remained constant at~6 nM for the controls yet significantly diminished from 10.46 nM to 4.94 nM for the DMSPamended samples (F = 5.01, p = 0.020) at a rate of 1.10 nM day −1 in the outer reef samples ( Figure 5C).
Contrasting trends in DMSP degradation gene abundance were evident between the two sites, with significant increases over time for the inner reef community, compared with a stable or declining gene abundance over time for the outer reef community ( Figure 5D-F Figure 5D-F). The inverse pattern was observed outside the reef, where the higher relative abundances of all three genes were in the +DMSP samples and DddP (F = 143.4, p = 0.001) and DmdA/A1 (F = 18.4, p = 0.001) decreased significantly over 120 h ( Figure 5E,F). Despite an apparent spike in DmdA/Dall at 24 h for the DMSP-enriched population, gene abundance increased over the 120 h (F = 10.5, p = 0.003) and maintained a higher relative abundance (F = 49.4, p = 0.001) than the controls ( Figure 5D). DMSP lyase activity for the phytoplankton (DLAp) and bacterial (DLAb) fractions were similar (between 0.47-0.79 µM DMS min −1 ) for the inner reef site. In contrast, DLAp rates were approximately four times higher than rates of DLAb for the outer reef site (0.04 and 0.01 µM DMS min −1 , respectively), and were an order of magnitude lower than those measured for the inner reef community (data not shown).
The key changes across treatments from the starting community were the dramatic increase in the relative abundance of Rhodobacterales (predominantly from the genus Sulfitobacter), which reached 50% in the controls, compared with only 8% in the DMSP-enriched samples (F = 6.8, p = 0.025; Figure 6B). While controls saw a decline in the relative abundance of SAR11 clade, DMSP-enriched communities exhibited an increase in the relative abundances, making up~28% of the final community (F = 5.9, p = 0.038). In both treatments we saw a decline in relative abundance of Synechococcales after 120 h (F = 33.9, p = 0.001), where final abundances were lower in the controls (<10%) compared with the +DMSP (~28%) treatment (F = 8.3, p = 0.018; Figure 6B).
For the outer reef site, the prokaryotic community resembled the initial inner reef waters ( Figure 7A), with strong changes in community structure over 120 h that differed between treatments (ANOSIM; global R = 0.617, p = 0.001). For the outer reef, the major groups contributing to 53.61% of the cumulative difference were Rhodobacterales, Synechococcales, Actinomarinales and Sphingomonadales (Supplementary Table S4). As with the inner reef, both treatments saw an increase in diversity, with several new orders (Alteromonadales, Caulobacterales, Chitinophagales, Microtrichales, Pseudomonadales and Sphingomonadales) appearing at >1%. Of the initial community, we saw a strong initial (24 h) increase in the relative abundance of Rhodobacterales in both treatments (F = 8.03, p = 0.008). In the +DMSP treatment this shift towards Rhodobacterales, (predominantly Sulfitobacter) was retained for 120 h, becoming the dominant order (~42%) within the community at 120 h, whereas in the control, the relative abundance returned to low numbers (~6%) by the end of the experiment ( Figure 7B). In both treatments, SAR11 clade showed no change in abundance through time, while the relative abundance of Synechococcales declined in both treatments (F = 39.1, p = 0.001), but remained higher (F = 35.3, p = 0.001) in the control (~27%) than the +DMSP (~7%) treatment ( Figure 7B). fitobacter), which reached 50% in the controls, compared with only 8% in the DMSP-enriched samples (F = 6.8, p = 0.025; Figure 6B). While controls saw a decline in the relative abundance of SAR11 clade, DMSP-enriched communities exhibited an increase in the relative abundances, making up ~28% of the final community (F = 5.9, p = 0.038). In both treatments we saw a decline in relative abundance of Synechococcales after 120 h (F = 33.9, p = 0.001), where final abundances were lower in the controls (<10%) compared with the +DMSP (~28%) treatment (F = 8.3, p = 0.018; Figure 6B). For the outer reef site, the prokaryotic community resembled the initial inner reef waters ( Figure 7A), with strong changes in community structure over 120 h that differed between treatments (ANOSIM; global R = 0.617, p = 0.001). For the outer reef, the major groups contributing to 53.61% of the cumulative difference were Rhodobacterales, Synechococcales, Actinomarinales and Sphingomonadales (Supplementary Table S4). As with the inner reef, both treatments saw an increase in diversity, with several new orders (Alteromonadales, Caulobacterales, Chitinophagales, Microtrichales, Pseudomonadales and Sphingomonadales) appearing at >1%. Of the initial community, we saw a strong initial (24 h) increase in the relative abundance of Rhodobacterales in both treatments (F = 8.03, p = 0.008). In the +DMSP treatment this shift towards Rhodobacterales, (predominantly Sulfitobacter) was retained for 120 h, becoming the dominant order (~42%) within the community at 120 h, whereas in the control, the relative abundance returned to low numbers (~6%) by the end of the experiment ( Figure 7B). In both treatments, SAR11 clade showed no change in abundance through time, while the relative abundance of Synechococcales declined in both treatments (F = 39.1, p = 0.001), but remained higher (F = 35.3, p = 0.001) in the control (~27%) than the +DMSP (~7%) treatment ( Figure 7B). Phytoplankton community composition of the inner reef was dominated by dinoflagellates belonging to seven different taxonomic groupings (Figure 8). Dino-Group (I, II, III and V) made up the highest relative abundance (38%), followed by diatoms from the order Centrales (18%). Dinoflagellates from the order Gymnodinales (13%) had the next highest relative abundance within the initial inner reef phytoplankton community. The remaining Phytoplankton community composition of the inner reef was dominated by dinoflagellates belonging to seven different taxonomic groupings (Figure 8). Dino-Group (I, II, III and V) made up the highest relative abundance (38%), followed by diatoms from the order Centrales (18%). Dinoflagellates from the order Gymnodinales (13%) had the next highest relative abundance within the initial inner reef phytoplankton community. The remaining phytoplankton taxa consisted of chlorophytes from the order Mamiellales (5.6%), dinoflagellates from the orders Prorocentrales (4.8%) and Peridiniales (2.3%), uncultivated marine Stramenopiles-MAST-group (3.4%), pennate diatoms (2.9%) and Pelagomonodales (1.5%). Changes in community structure were observed for both the control and DMSP-enriched communities over the 120 h (ANOSIM; global R = 0.556, p = 0.002) with five orders-Diatoms (centrales), Dino-Group (I, II, III, V), Gymnodiniales, Suessinales and Diatom (pennales)-contributing to 51.41% of the dissimilarity between treatments (Supplementary Table S4).  Table S4). Relative abundance of phytoplankton taxa that make up >1% of the total community at T0, 24 (T1), 72 (T2) and 120 h (T3) for control and +DMSP incubations from the inner reef site. Data are based on eukaryotic (18S) sequences. * In cases where taxonomic identification could not be made to lower level than class, class is shown.
To test for relationships between DMSP catabolising taxa and the relative gene abundance, linear regressions were performed on selected prokaryotic orders known to harbour relevant taxa ( Figure 10). Specifically, qPCR primer target groups SAR11 (DmdA/Dall; [66]) and Rhodobaterales, which harbour Roseobacters, the target group for A1 and DddP primers [58,66] were tested. We found significant relationships between Rhodobacterales and both 16S normalised DmdA/A1 (adjR 2 = 0.67; p = 0.0003) and DddP (adjR 2 = 0.52; p = 0.002), but only for the inner reef site ( Figure 10A,B). This relationship appeared to be driven by the abundance of the Roseobacter genus Sulfitobacter (Supplementary Figure S1), where the significant relationships explained 67% (DmdA/A1) and 52% (DddP) of the variation. We found no relationships between the relative abundance of SAR11 clade and DmdA/Dall at either reef site ( Figure 10C). Relative abundance of phytoplankton taxa that make up >1% of the total community at T0, 24 (T1), 72 (T2) and 120 h (T3) for control and +DMSP incubations from the inner reef site. Data are based on eukaryotic (18S) sequences. * In cases where taxonomic identification could not be made to lower level than class, class is shown.
For the outer reef, the phytoplankton community structure strongly resembled that of the inner reef, with a dominance of dinoflagellates, but a notable absence (<1%) of diatoms ( Figure 9). As with the inner reef, Dino-Group (I, II, III and V) made up the highest relative abundance (51%), followed by dinoflagellates from the Gymnodinales (23%), Dinophyceae (7%), Prorocentrales (5%), Peridinales (4.6%) and Suessiales (1.3%). The only non-dinoflagellate taxa present at >1% in the initial community belong to the MAST-Group (2.7%). As with the inner reef community, there were shifts in community structure over 120 h (ANOSIM; global R = 0.521, p = 0.001), with Dino-Group (I, II, III, V), Diatoms (centrales), MAST-Group, Gymnodiniales and Stramenopiles as the major contributors to 53.49% of the dissimilarity between treatments (Supplementary Table S4). Relative abundance of phytoplankton taxa that make up >1% of the total community at T0, 24 (T1), 72 (T2) and 120 h (T3) for control and +DMSP incubations from the outer reef site. Data are based on eukaryotic (18S) sequences. * In cases where taxonomic identification could not be made to lower level than class, class is shown.

Discussion
Corals and reef lagoons of the GBR are well characterised as DMS/P hot spots, due to high local production by corals [47]; however, limited information exists on DMSP concentration and cycling outside shallow reef waters. Dimethylsulfoniopropionate (DMSP), which in the open ocean is predominantly produced by marine phytoplankton and taken up by several groups of marine heterotrophic bacteria [6,7] plays an important role in marine ecosystems. While there is considerable evidence for bacterial uptake and processing of DMSP [29,66], little research exists on the uptake of DMSP by non-DMSP producing phytoplankton [33,34]. Using a combination of short-and longer-term DMSP enrichment experiments, we investigated which members of the marine microbial community take up available DMSP from the environment and the effect that increasing DMSP availability has on the abundance of DMSP degradation genes and microbial community composition in natural sea waters of the Great Barrier Reef, Australia. Relative abundance of phytoplankton taxa that make up >1% of the total community at T0, 24 (T1), 72 (T2) and 120 h (T3) for control and +DMSP incubations from the outer reef site. Data are based on eukaryotic (18S) sequences. * In cases where taxonomic identification could not be made to lower level than class, class is shown.
To test for relationships between DMSP catabolising taxa and the relative gene abundance, linear regressions were performed on selected prokaryotic orders known to harbour relevant taxa ( Figure 10). Specifically, qPCR primer target groups SAR11 (DmdA/Dall; [66]) and Rhodobaterales, which harbour Roseobacters, the target group for A1 and DddP primers [58,66] were tested. We found significant relationships between Rhodobacterales and both 16S normalised DmdA/A1 (adjR 2 = 0.67; p = 0.0003) and DddP (adjR 2 = 0.52; p = 0.002), but only for the inner reef site ( Figure 10A,B). This relationship appeared to be driven by the abundance of the Roseobacter genus Sulfitobacter (Supplementary Figure S1), where the significant relationships explained 67% (DmdA/A1) and 52% (DddP) of the variation. We found no relationships between the relative abundance of SAR11 clade and DmdA/Dall at either reef site ( Figure 10C).

Discussion
Corals and reef lagoons of the GBR are well characterised as DMS/P hot spots, due to high local production by corals [47]; however, limited information exists on DMSP concentration and cycling outside shallow reef waters. Dimethylsulfoniopropionate (DMSP), which in the open ocean is predominantly produced by marine phytoplankton and taken up by several groups of marine heterotrophic bacteria [6,7] plays an important role in

Discussion
Corals and reef lagoons of the GBR are well characterised as DMS/P hot spots, due to high local production by corals [47]; however, limited information exists on DMSP concentration and cycling outside shallow reef waters. Dimethylsulfoniopropionate (DMSP), which in the open ocean is predominantly produced by marine phytoplankton and taken up by several groups of marine heterotrophic bacteria [6,7] plays an important role in marine ecosystems. While there is considerable evidence for bacterial uptake and processing of DMSP [29,66], little research exists on the uptake of DMSP by non-DMSP producing phytoplankton [33,34]. Using a combination of short-and longer-term DMSP enrichment experiments, we investigated which members of the marine microbial community take up available DMSP from the environment and the effect that increasing DMSP availability has on the abundance of DMSP degradation genes and microbial community composition in natural sea waters of the Great Barrier Reef, Australia.

Quantification of DMSP Uptake by Size Fractions of the Microbial Community-Experiment 1
The seawater characteristics of both sites were typical of coral reef waters, with nutrients and temperature values within ranges of previous studies made on the GBR [67][68][69]. Measured DMS/P concentrations from both reef sites fell within the range previously found in the GBR, albeit at the lower end. Previous studies have measured DMS concentrations of~2.0 (range; BDL-54) nM and~1.3 (BDL-3.9) nM in the reefs of Orpheus and Magnetic islands, respectively [48], and DMSP values of~3.8 (0.  nM at One Tree Reef [46]. The relatively low values of DMS/P measured in this study are likely due to the location of the sampling points and time of collection; as peak DMS concentrations have been measured in high summer above the reef flat [48], whereas our data were sampled in October and collected from open water~13 km from the closest reef. The patterns of DMSP uptake by the different fractions of the marine microbial communities of the GBR determined over 7 h (experiment 1) showed a loss of DMSPt, for both sites and treatments, indicating that DMSP was escaping from the system. This could occur via bacterial conversion to DMS, DMSO and MeSH or by photochemical oxidation to DMSO [70][71][72]. As DMSPp accumulated with enrichment, we suggest that a significant fraction of available DMSPd was taken up and transformed, either to satisfy the sulfur demand in form of MeSH, to cleave into DMS and be lost to the atmosphere, or to oxidise to DMSO to protect the cells from ROS [10,21,37]. At both reef sites, DMSP enrichment led to greater DMSP uptake in the largest (8 µm) fraction, supporting the findings that phytoplankton, as well as bacteria take up DMSP [18,34,35,73], providing new evidence that the larger members of the microbial community can act as important DMSP sinks in coral reef waters. The differences in DMSP uptake between sites may be due to differences in the eukaryotic composition, where the presence of diatoms at the inner reef site may account for the greater uptake with DMSP enrichment [34].
At both reef sites, prokaryotic communities were dominated by the photosynthetic cyanobacteria Prochlorococcus, typical of warm oligotrophic areas, due to its adaptability to thrive in nutrient poor waters [74,75]. Because of its high abundance and known ability to take up DMSP [33,76], it is possible that this taxon plays an important role in DMSP cycling on the reef and may be an important DMSP sink within the smallest size fraction (>0.22 µm fraction). While previous work has shown that cyanobacteria take up DMSP, little is known whether they catabolise it. To date, only one sequenced cyanobacterial strain, Synechococcus sp. KORDI-100, has been found to have a DMSP lyase gene [77]. Given that Prochlorococcus and SAR11 were present in high relative abundances at both sites, the fact that DMSP showed higher accumulated in the small fraction of the outer reef, suggests that rather than differences in microbial community composition, it was the metabolic state of the community driving this difference, whereby the bacterial community of the outer reef site were preferentially taking up DMSPd to meet their sulfur and carbon demands for growth [21,27]. Given the difference in depth of the two sites and therefore potentially different mixing regimes, the influence of mechanistic controls, such as UV-A dose [58], on the bacterial switch cannot be ruled out as a potential driver of the differences between sulfur demand for the two sites.
We found dinoflagellates dominated the reef waters from both sites (Figure 4). Given that dinoflagellates are known to be some of the highest-DMSP producers [2,5], it is likely that species from this class constituted the primary source of DMSP in these waters. The higher concentrations of DMSP inside the reef, suggests a greater absolute abundance of these cells, with the possibility of additional contributions from dinoflagellate symbionts associated with corals [46][47][48]. In support of higher biomass, the fluorescence trace from the CTD revealed the chlorophyll a signature was twice as high at the inner reef site compared with that of the outer reef, indicating higher phytoplankton biomass in the inner reef surface waters. As such, we propose that the dinoflagellates accounted for the high proportion of DMSPp in the >8 µm fraction, but that the main taxa responsible for the uptake of additional DMSP in the +DMSP treatment be attributed to the diatoms, which have been shown to take up and accumulated DMSP in high concentrations [34]. The absence of diatoms from the outer reef site, suggests that other large non-DMSP producing phytoplankton may act as DMSP sinks, but further work is needed to confirm uptake and identify these groups.
At both sites, the most abundant DMSP degrading pathway was demethylation, represented by the gene DmdA (Dall and A1). Similar relative abundances of these DMSP degradation genes using qPCR have been reported for the North Pacific subtropical gyre [78]. The DmdA gene is responsible for the first step of the demethylation pathway and so, it is likely that the disappearance of DMSP from the incubations over 7 h for both reef sites was through MeSH production [29]. Of the bacterioplankton, the groups most likely to demethylate DMSP in theses water are SAR11 and members of the Rhodobacterales [29]. However, as SAR11 were more abundant than Rhodobacterales at both reef sites, we speculate that they may form the dominant DMSP consumer in the surface waters of the GBR, a finding commensurate with those for the Sargasso Sea and North Atlantic Ocean [79], Tropical and Subtropical Pacific Ocean [80], and East China sea [81]. With respect to DMSP cleavage via DddP genes, these processes can be attributed to several genera from within the Rhodobacterales, many of which have been shown to have strong correlations with DMS production [81]. It is important to note however, that the highly conserved nucleotide sequences of the primers used in qPCR means that copy numbers recorded in this study likely underestimate total gene abundance in the GBR reef waters, as many sequences are left out of the analysis. For example, in the case of the demethylation gene DmdA, Dall primers target SAR11-like sequences, while A1 targets Roseobacter-like sequences [29], and while both clades D and A have been shown to make up the majority of sequences when a universal DmdA primer has been used [78], the primer sets used in this study do not cover all DmdA subclades. Similarly, it is impossible to rule out that some bacterioplankton lacking DMSP degradation genes may take up DMSP for other reasons, such as for osmolytic benefit, as shown in a previous study, in which coastal seawater filtrates containing mostly bacteria, diminished their MeSH production and retained up to 54% of their intracellular DMSP when under osmotic stress [82].

Dynamic Changes in Microbial Community Composition and DMSP Degradation Gene Abundance in Reef Waters-Experiment 2
Longer-term incubations designed to follow DMSP-induced community changes and shifts in catabolising gene abundance revealed strong differences between sites and treatments. For both sites, concentrations of DMSPt and DMSPd decreased over time in both control and DMSP-enriched treatments, with DMSP consumption rates (IR: 5-12, and OR 9-20 nM d −1 ) within previously reported values (7-80 nM d −1 ; [83,84]). These data indicate that DMSP was rapidly (within hours) taken up from the solution and then lost within 24 h via production of DMS or MeSH [21], in accordance with experiment 1. The lack of intracellular accumulation (increase in DMSPp) over several days suggests that the DMSP that was taken up within the first few hours and was being utilised or converted, regardless of DMSP availability, a finding that is in accordance with previous results where added 35 S-DMSPd was transformed into several products over time scales of minutes to one day [66].
To date, DMSP lyase activity (DLA) in phytoplankton has only been shown in Emiliania huxleyi [85], Phaeocystis spp. [86] and different Symbiodiniaceae genera [87,88], all high DMSP producers. As such, we hypothesise that the DLAp measure in these GBR waters is likely attributable to the high DMSP-producing dinoflagellates that dominated the phytoplankton communities at both sites. For the bacterial fractions, the high rates of DLAb of the inner reef suggest that DMSP concentrations may have exceeded the necessary levels to cover the sulfur demands of the bacterioplankton community, and therefore much of the excess was cleaved to DMS [21,26,27]. In contrast, outside the reef, it is likely that most of the added DMSP was utilised to meet the sulfur demands of the bacterioplankton community. Taken together, our data suggest that the inner reef site had higher inputs of organic carbon and sulfur and more efficient nutrient recycling systems, while the outer reef site was more representative of the oligotrophic open ocean, whereby DMSP enrichment provided cells some relief from carbon and sulfur limitation.
As with the short-term experiments, DmdA/Dall was the dominant DMSP degradation gene throughout the experiment (120 h) with similar relative abundances at both sites at the outset when most of the DMSP was lost from the system, suggesting that most of the DMSP was degraded through demethylation [22,23]. These data are congruent with previous studies, where DmdA subclade D is present at an order of magnitude higher than DmdA/A1 and DddP [29,30,58,78,81]. The relative abundances of the different DMSP degradation genes over time showed variable behaviours for the two locations. In agreement with the DLAb data, all degradation genes increased in relative abundance with time at the inner reef site, suggesting that sulfur demands were met quickly [27], after which cleavage to DMS dominated. For the outer reef site, the relative abundance of genes was higher in the DMSP-amended samples compared to the controls, indicating that the microbial community from the outer reef was partially advantaged by DMSP enrichment, but overall abundances declined over time. These data are supported by the low DLA levels measured for this site, strengthening the idea that the DMSP concentrations in the oligotrophic waters outside the reef were insufficient to completely satisfy the sulfur and carbon demands of the outer reef microbial community.
In favouring DMSP utilisers, enrichment with DMSP can disrupt competitive hierarchies in microbial communities, allowing some taxa to proliferate over others. We saw DMSP-related changes in the phytoplankton community over the 120 h at both sites, predominantly between the diatom and dinoflagellate groups, indicating an interplay between DMSP sources and sinks within the phytoplankton community. For the prokaryotic community we measured key changes over time, some of which were in response to DMSP enrichment and could be linked with changes in gene abundance ( Figure 10). In particular, the proliferation of SAR11 clade with DMSP addition at the inner reef site suggest that these taxa were able to monopolise on increased DMSP availability, despite finding no relationship with demethylation gene abundance. The absence of a significant relationship between SAR11 clade and DmdA/Dall, indicates that SAR11 clade may not be the major contributors to MeSH production via demethylation in GBR waters. It does not preclude the possibility, however, that SAR11 clade may have contributed to the high DLAb activity measured at the inner reef site, via the DddK gene [32,89], but this was not measured directly here. Contrastingly, Rhodobacterales, known to harbour several genera with the capacity to utilise DMSP [81], showed no benefit from DMSP enrichment at the inner reef site, yet exhibited strong correlations with DMSP catabolising gene abundances (both DmdA and DddP), suggesting that taxa within this order are the dominant DMSP degrading bacteria in these inner coral reef waters. Rhodobacterales have been previously found to positively correlate with relative abundances of both DddP and DmdA [80,81], and are one of the major groups responsible for DMSP demethylation in the marine environment [29]. In particular, Sulfitobacter, a genus within the Rhodobacter order and the dominant genus in our study, have been previously linked with these genes and shown to use DMSP as their sole carbon source for MESH production [81]. These observed differences in community responses to DMSP enrichment between the two sites are compatible with the theory of the bacterial switch, whereby enrichment of the inner reef community resulted in the sulfur demand being satisfied and the subsequent loss of DMSP from the system via cleavage to DMS, whereas for the more oligotrophic waters of the outer reef, even with DMSP enrichment, the DMSP availability was insufficient to completely satisfy the sulfur and carbon requirements of the community.

Conclusions
In this study, we evaluated the uptake of DMSP by different fractions of marine microbial communities of the GBR and revealed that both bacteria and phytoplankton from natural reef waters can take up DMSP over short time scales. Specifically, DMSP enrichment revealed the main sink for DMSP comprised of taxa from the largest microbial fraction. By investigating the main DMSP catabolising genes over 120 h, we were able to show that despite similarities in microbial community composition, strong differences in DMSP catabolism and thus influence on sulfur cycling are possible. We found that taxa from the Rhodobacterales order best explained the increases in DMSP catabolising genes, but this relationship was only found for the inner reef site. The differences in DMSP lyase activity, catabolism and community structure observed between the two reef sites suggests that at the inner reef site sulfur and carbon demands were largely satisfied by existing DMSP availability, meaning that DMSP addition had minimal impact on the microbial community from these waters, with lyase activity dominating the conversion of DMSP to DMS. On the other hand, the absence of a strong response to DMSP enrichment for the outer reef bacterioplankton community (minimal increases in abundance and no positive relationships with DMSP catabolising genes), when taken in the context of the low DLA rates measured, suggests minimal DMSP cleavage to DMS occurring, and therefore we may presume that any DMSP in the system was being converted to MeSH or lost from the system via oxidation. Overall, our findings contribute to an improved understanding of the fate and function of DMSP in tropical marine ecosystems and identify microbial DMSP sinks and sources within the deeper, oligotrophic waters of the GBR.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/microorganisms9091891/s1, Figure S1: Regression analysis of the relative abundance of DMSPdegrading genes and Sulfitobacter abundance for the inner (red) and outer (yellow) reef sites, Table S1: Forward and reverse primer pairs for DMSP catabolizing genes and 16S normalisation gene, Table S2: Macronutrient concentrations for control and DMSP-enriched samples during experiment 2 (IR), Table S3: Macronutrient concentrations for control and DMSP-enriched samples during experiment 2 (OR), Table S4: SIMPER dissimilarity output for control and +DMSP treatment over time for inner reef (IR) and outer reef (OR) bacterioplankton (16S) and phytoplankton (18S) communities.
Author Contributions: Conceptualization, K.P. and J.R.S.; formal analysis, E.F. and M.O.; investigation, E.F. and N.S.; data curation, E.F., M.O. and K.P.; writing-original draft preparation, E.F. and K.P.; writing-review and editing, E.F., J.R.S. and K.P.; visualization, E.F. and K.P.; supervision, K.P.; funding acquisition, K.P. and J.R.S. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: Sequence data files for the 16S rRNA and 18S rRNA are available in the NCBI Sequence Read Archive (SRA) under Bioproject number PRJNA755150. All remaining data presented in this study are available on reasonable request from the corresponding author.