Abstract
The pelagic ocean receives terrigenous inputs of a range of organic compounds; however, the role that this terrigenous material plays in the ocean carbon cycle and biological pump is not entirely understood, and questions remain as to how oceanic cycles of terrigenous and autochthonous carbon interact. A significant portion of organic carbon that cannot be utilized by marine microbes in the epipelagic ocean escapes microbial remineralization to be sequestered in the deep ocean as refractory dissolved organic matter (DOM). Lignin, a “model” terrigenous compound, is thought to be refractory in the open ocean unless chemically altered. However, in this study, incubation experiments performed using lignin-amended oligotrophic seawater from the Sargasso Sea exhibited bacteria and archaea growth that doubled compared to unamended control treatments. The increase in bacteria and archaea cell abundance in lignin-amended treatments coincided with a 21–25% decrease in absorbance (250–400 nm) of chromophoric dissolved organic matter (CDOM), suggesting that certain microbes may be capable of altering fractions of this ostensibly recalcitrant organic matter. Furthermore, the microbial response to the lignin-amended treatments appears to be taxon-specific. Two phyla of Archaea, Euryarchaeota and Thaumarchaeota, exhibited an increase in abundance of 7-fold and 28-fold (from 2.42 × 106 cells L–1 to 1.72 × 107 cells L–1, and from 1.60 × 106 cells L–1 to 4.54 × 107 cells L–1, respectively), over 4 days of incubation in lignin-amended treatments. Additionally, an increase of 11-fold and 13-fold (from 2.93 × 106 cells L–1 to 3.30 × 107 cells L–1, and from 3.26 × 106 cells L–1 to 4.28 × 107 cells L–1, respectively), was observed in the abundance of these phyla in treatments containing lignin with added nitrogen and phosphorus, thus raising questions regarding primary and/or secondary responses to lignin degradation. Our findings indicate that marine bacteria and archaea play a role in the transformation of the optical properties of lignin in the open ocean and that they may serve as a potential sink for a portion of the lignin macromolecule.
Introduction
Dissolved organic matter (DOM) is an important component of many ocean biogeochemical cycles (; ; ) and can serve as a substrate for heterotrophic archaea and bacteria (; ). The majority of DOM in the open ocean is of marine origin (; ); however, terrigenous inputs of compounds such as lignin may be important to the oceanic carbon cycle. Lignin, a cross-linked macromolecule, unique to woody vascular plants, has been used as a model compound and tracer of terrigenous carbon (; ). On the basis of measurements of lignin phenol concentrations, it has been shown that less than 2% of oceanic DOC is terrigenous (); yet, if estimates of annual riverine export to the ocean and atmospheric deposition of dissolved or water soluble organic carbon to the surface ocean are considered (∼0.4 and ∼0.65 Pg organic C, respectively), the concentrations of lignin phenols should be higher (; ; ; ). This raises a number of questions about the transport and fate of terrigenous DOM in the open ocean.
Furthermore, a portion of DOM has been shown to be light-absorbing, or chromophoric (CDOM), and is found ubiquitously throughout the global ocean. CDOM accounts for up to 50% of blue light absorption in the open ocean, thereby effectively regulating photosynthesis and influencing remote sensing algorithms used to predict primary production. Despite its impact on biogeochemical cycles, the composition and sources of CDOM in the open ocean remain largely unresolved. Dark microbial incubation studies have demonstrated that production of CDOM is positively correlated with an increase in bacterial cell abundance, providing evidence that CDOM is a byproduct of heterotrophic bacterial production (; ). The same studies also demonstrated that some fraction of CDOM is subsequently removed, presumably due to microbial remineralization. It is thought that the chemical composition and degradation rate of marine CDOM varies depending on the quality of the organic nutrient substrate available to bacterioplankton (). Lignin is also chromophoric, but its role in the pelagic CDOM pool is poorly understood. A number of studies have investigated the relationship between lignin and microbes in soil or riverine water, and suggest that some fraction of lignin is biologically available and that distinct microbial lineages under specific nutrient conditions are capable of degrading lignin (; ; ; ). Additionally, recent genomics studies have shown a wide range of organisms with novel enzymes that have the potential to degrade lignin (; ). However, the ability of microbial communities to utilize lignin in nutrient-depleted oceanic regions and the relationship to CDOM variability has not been fully investigated.
Here, we evaluate the response of the bacteria and archaea communities to lignin amended oligotrophic seawater from the Sargasso Sea to determine if specific lineages respond to lignin amended treatments and to assess whether terrigenous organic matter can serve as an intermediary in the production of autochthonous CDOM.
Materials and Methods
Sample Collection
Samples were collected at the Bermuda Atlantic Time Series Study (BATS) site in the North Atlantic (31° 40′ N, 64° 10′ W) on board the R/V Atlantic Explorer in July, 2014. Water was collected from depths of 1 and 60 m, using Niskin bottles on a rosette containing a conductivity, temperature and depth (CTD) sensor package. An inoculum of unfiltered “whole water” containing natural bacteria and archaea assemblages was added to grazer-diluted incubation media. Water samples for the incubation media were collected from 1 m depth, to ensure low background concentrations of CDOM as a result of photobleaching, low lignin concentrations due to photooxidation and low nutrients (, ; ). Seawater collected from 1 m depth was filtered through a 0.2-μm filter (Whatman Polycap AS) and 14-L aliquots were transferred into 20-L polycarbonate carboys (Nalgene, United States) that were pre-cleaned with 10% hydrochloric acid. Water for the inoculum was collected from a depth of 60 m (where maximum bacterioplankton and CDOM concentrations can typically be observed at BATS) (; ; ), and 7 L of inoculum was added, unfiltered, to each carboy, thereby diluting the filtered surface seawater to 70% of the total volume (Supplementary Table S1). This diluted seawater culture approach has been shown to alleviate grazing pressure by protistan grazers (, ; ). In situ concentrations of NO3 and PO4 were <0.02 and <0.03 μmol L–1, respectively, at both 1 m and 60 m depth (); CDOM absorbance (325 nm) at 1 m was 0.075 m–1 and 0.155 m–1 at 60 m. Although in situ DOC concentrations were not measured on the cruise, TOC concentrations ranged from 67.1 μmol L–1 at 60 m to 67.4 μmol L–1 at the surface. Previous studies have indicated that DOC comprises the majority of TOC in the epipelagic zone at BATS and have evidenced DOC concentrations of ∼68 μmol L–1 from June to October ().
Nutrient Amendments
Replicate incubations (n = 2) were performed in the dark and at the in situ temperature of the inoculum (23°C). Seawater dilution cultures included unamended controls and treatments that included a combination of organic and inorganic amendments (Supplementary Table S1). The DOM amendments (shown in Supplementary Table S1) included: (1) lignin (96% purity, alkali-extracted with low sulfonate content; Sigma-Aldrich) (Lignin); (2) lignin with the addition of 1 μmol L–1 ammonium chloride (99.998% trace metal basis; Sigma-Aldrich) and 0.1 μmol L–1 dipotassium phosphate (BioUltra anhydrous ≥ 99.0%; Sigma-Aldrich) (LNP); or (3) glucose (D-(+ -Glucose ≥ 99.5% (Sigma-Aldrich) with the addition of 1 μmol L–1 ammonium chloride and 0.1 μmol L–1 dipotassium phosphate (GNP). The inclusion of both GNP and LNP treatments in this study allowed for the direct comparison of the microbial response to a model labile substrate, compared with a complex terrigenous organic substrate. The commercial lignin was chosen according to previous work by . Concentrations of N and P were chosen to be similar to those used in previous work at the BATS site (). D-glucose was added to the carboys to obtain a final concentration of 10 μmol L–1 C. The solution was prepared such that less than 20 mL of the aqueous glucose solution was added to the 20 L sample, and thus did not affect the ionic strength of the sample. A solution of lignin was prepared by the addition of solid lignin material to de-ionized water (MilliQ, Millipore, United States), with subsequent sonication and filtration (GF/F filters, Whatman, and a pre-rinsed 0.2-μm pore size track-etched polycarbonate filter, Whatman). The total organic carbon (TOC) content of the solution was determined according to prior to inoculation and the solution was pipetted directly into the appropriate carboys to target a concentration of 10 μmol L–1 C. Samples for bacteria and archaea cell abundance, DOC, CDOM, and fluorescence in situ hybridization (FISH) were drawn by spigot directly into sampling containers at discrete time intervals over a period of 35 days.
Determination of CDOM Absorption Coefficient
Changes in CDOM concentrations were determined by UV-Visible spectroscopy. Subsamples for each amendment were gravity filtered through a 0.2-μm filter (polycarbonate track-etched membrane filter; Whatman, Chicago, IL, United States) that had been pre-flushed with 500 mL of de-ionized water and collected into muffled glass vials (I-Chem, Thermo Fisher Scientific, Waltham, MA, United States) fitted with PTFE-lined polypropylene caps. Samples were stored in the dark at 4°C until analysis. Absorbance measurements were made using a dual-beam spectrophotometer (Perkin Elmer Lambda-18), equipped with a photo-multiplier tube. Samples were analyzed in quartz cuvettes with a 10-cm pathlength and were blank corrected against de-ionized water. Blank-corrected spectra (250–700 nm) were reported as absorption coefficients (m–1) using the equation:
where a is the absorption coefficient (m–1), A is the absorbance and l is the pathlength (m). Absorbance values at discrete wavelengths in the near-UV spectrum (325 and 340 nm) were selected from individual spectra to measure changes in concentration and changes in CDOM components over time, particularly lignin (; ; ).
Enumeration of Bacteria and Archaea Cell Abundance
Bacteria and archaea abundance samples (40 mL) were fixed with formalin (10% final concentration), and stored at −80°C until slide preparation. Slides were prepared according to , , and . Slides were then enumerated using an AX70 epifluorescence microscope (Olympus, Tokyo, Japan) under ultraviolet excitation at 1000× magnification. At least 500 cells per slide (10 fields) were counted for bacteria and archaea abundance.
Quantification of Dissolved Organic Carbon
Dissolved organic carbon was sampled in duplicate and filtered through 0.7-μm pre-combusted glass fiber filters (Whatman, GF/F) into pre-combusted borosilicate glass vials (40 mL; I-Chem) fitted with acid-washed Teflon-lined caps and stored frozen at –20°C until analysis. All DOC samples were sent to the University of California for analysis. Samples were analyzed using high-temperature combustion on Shimadzu TOC-V analyzers, with slight modifications from the manufacturer model, as described in . In order to minimize the machine blanks, expansive conditioning of the combustion tube with repeated injections of low carbon water (LCW) as well as deep seawater was performed. A daily standardization of the system response was performed using a four-point standardization calibration curve of glucose solution in LCW. In-house surface and deep references were included in every analytical run in order to ensure comparability between sample sets. The in-house surface and deep references used in this study were collected from the Santa Barbara Channel in April 2014 and calibrated against consensus reference material () and historical house reference materials.
The average surface and deep reference values were 66.6 ± 0.9 and 38.0 ± 0.7 μmol L–1 C, respectively. Further details of this method are described in , .
Enumeration of Specific Lineages by Fluorescence in situ Hybridization (FISH) and Catalyzed Reporter Deposition-FISH (CARD-FISH)
Fluorescent in situ hybridisation was utilized to quantify the abundance of the major bacteria and archaea phylotypes present in the seawater. We followed the protocols described in and for all FISH and CARD-FISH. The FISH probes used for this study (Supplementary Table S3) have previously been validated in silico for specificity using Probe Match on the Ribosomal Database Project () and TestProbe and Probebase on the SILVA website (). The bacterial and archaeal groups quantified included the SAR11 clade (152R, 441R, 542R, 732R probes), Alteromonas spp. (AC137R), Flavobacteria II (CFB563R), Rhodobacteraceae (536R), SAR202 (103R, 311R), Euryarchaeota (Eury806), and Thaumarchaeota (Thaum537) ().
Results
Bacteria and Archaea Cell Abundance
The cell abundance data show that in all treatments where lignin was added (Lignin or LNP), the microbial response was greater than that observed in the unamended Controls (Paired t-test for Control and Lignin (p < 0.05), and for Control and LNP (p < 0.05); Figure 1). The unamended Controls exhibited a 2.5-fold increase in bacteria and archaea cell abundance (from 1.875 × 108 cells L–1 to 4.42 × 108 cells L–1) between days 0 and 12. The treatments that were amended solely with commercial lignin showed an average increase in cell abundance of ∼3-fold between days 0 and 5 (from 1.67 × 108 cells L–1 to 4.95 × 108 cells L–1), after which the population remained relatively unchanged until day 12, when it increased by another 1.2-fold (to 6.59 × 108 cells L–1). The treatments amended with lignin in addition to nitrogen and phosphorus (LNP) showed an average increase in cell abundance of 3-fold between days 0 and 5 (from 1.99 × 108 cells L 1 to 6.20 × 108 cells L–1) and then remained in stationary growth for the duration of the incubation. The glucose amendments with added nitrogen and phosphorus (GNP) demonstrated an increase in cell abundance (from 1.81 × 108 cells L–1 to 1.10 × 109 cells L–1) between days 0 and 2, then approximated stationary growth through day 6 (at ∼ 1.40 × 109 cells L–1), after which cell densities declined through day 12.
FIGURE 1
Changes in CDOM Absorbance
The absorbance of the unamended Control (absorption coefficient at 325 nm; m–1) ranged from 0.098 m–1 at day 0 (T0) to 0.134 m–1 at day 12 (T15; Figure 2, top panel). In the GNP amendments, CDOM absorbance at 325 nm increased continuously over the first 3 days from 0.129 m–1 to 0.182 m–1, then decreased to 0.123 m–1 by day 5 (T9) (Figure 2, top panel) and remained relatively unchanged (paired t-test for T9 and T14, p > 0.05) for the remainder of the experiment (Figure 2). CDOM 325 nm absorption coefficients at the absorbance maximum in the GNP treatments were almost 70% higher than for Controls at the same time point. Both the Lignin and LNP amended samples exhibited a higher initial absorbance (between 250 and 450 nm) than the Controls at day 0 (T0), due to the chromophoric nature of the amendments themselves (T0 = 0.664 m–1 and 0.650 m–1, respectively). Spectra (250–450 nm) of both Lignin and LNP amended samples (Figure 3, top and center panels, respectively), exhibited a discrete peak at approximately 340 nm, attributed to the amendment of lignin, as well as a less pronounced shoulder from 260–280 nm. Both of these peaks decreased substantially between days 2 and 10, indicating a loss of absorbance at these wavelengths of approximately 30%, and resulting in spectra that resembled the shape of oceanic CDOM spectra (Figure 3, bottom panel). In both treatments, absorbance at 325 nm increased until day 2 (Figure 2, bottom panel), indicating the production of CDOM in situ and then exhibited an almost 25% decrease in absorbance, falling below initial absorbance values (Figure 2, bottom panel). Both LNP and Lignin absorbance at 325 nm remained relatively unchanged from day 10 to 35 (paired t-tests, p = 0.06 and p = 0.25, respectively; Supplementary Figure S1).
FIGURE 2
FIGURE 3
Dissolved Organic Carbon Utilization
Dissolved organic carbon analysis confirmed an addition of ∼10 μmol L–1 C to each of the Lignin, LNP, and GNP treatments, compared to the control at the time of inoculation (Control = 66.6 μmol L–1 C ± 1.2; Lignin = 78.2 μmol L–1 C ± 1.0; LNP = 78.2 μmol L–1 C ± 1.1; GNP = 77.6 μmol L–1 C ± 1.1.). There was no resolvable change in DOC throughout the incubation for either the unamended control or the LNP treatments. We observed a loss of ∼ 4 μmol L–1 and 11 μmol L–1 in the Lignin and GNP treatments, respectively (Figure 4).
FIGURE 4
Specific Lineages Measured Using FISH and CARD-FISH
The initial microbial community (T0) in all treatments was dominated by members of the SAR11 lineage, with cell abundances between 4.18 and 5.82 × 107 cells L–1, while cell abundances of all other targeted groups were at least an order of magnitude lower (Supplementary Figures S2A–D). In the unamended Control samples, SAR11 continued to increase throughout the incubation, reaching 1.61 × 108 cells L–1 by day 6 (T10), with cell abundances still above 1 × 108 cells L–1 on day 10 (T13; Figure 5 and Supplementary Figure S2A). In contrast, bacteria and archaea populations in the amended treatments shifted by day 4, and SAR11 abundance reached stasis, while other lineages increased (Supplementary Figures S2B–D).
FIGURE 5
In the GNP treatments, Alteromonas spp. increased from >2 × 106 cells L–1 (T0) to almost 9 × 108 cells L–1 by day 4, comprising 56.3% of the total microbial community (Figure 5B and Supplementary Figure S2B). The increase in Alteromonas was less pronounced in both lignin-amended samples (a net increase of 2.60 and 4.76 × 107 cells L–1, respectively, compared to an increase of 8.87 × 108 cells L–1 in the GNP treatments), whereas Archaea became more pronounced in these treatments (Figure 5A and Supplementary Figures S2C,D). During days 0 through 4, the cell population of Thaumarchaeota increased from 1.60 × 106 cells L–1 to 4.54 × 107 cells L–1 (a factor of ∼28) in the Lignin treatments and from 3.26 × 106 cells L–1 to 4.28 × 107 cells L–1 (a factor of ∼13) in the LNP treatments. Over the same time period, Euryarchaeota increased from 2.42 × 106 cells L–1 to 1.72 × 107 cells L–1 (a factor of ∼7) in the Lignin treatments and from 2.93 × 106 cells L–1 to 3.30 × 107 cells L–1 (a factor of ∼11) in the LNP treatments (Figure 5A and Supplementary Figures S2C,D). In fact, Thaumarchaeota increased from 3.26 × 106 cells L–1 to 4.54 × 107 cells L–1 between days 0 and 4 (Supplementary Figure S2C). However, by day 6, Thaumarchaeota abunances in all lignin-amended treatments had decreased to approximately 2 × 107 cells L–1. Euryarchaeota cell abundance decreased from 1.71 to 1.16 × 107 cells L–1 between days 4 and 6 in the Lignin treatments, and decreased from 3.30 × 107 cells L–1 to only 5.70 × 106 cells L–1 between days 4 and 6 in the LNP treatments (Supplementary Figures S2C,D). In the unamended Control treatments, the abundance of Thaumarchaeota and Euryarchaeota decreased by day 4 (a loss of 5.02 × 106 cells L–1 and 4.47 × 107 cells L–1, respectively. Cell densities of Thaumarchaeota also decreased in the GNP treatments throughout the first 4 days (from 9.87 × 105 cells L–1 to 7.29 × 105 cells L–1), while Euryarcheota cell abundances were low but exhibited an increase (from 8.65 × 105 cells L–1 to 1.41 × 106 cells L–1) (Supplementary Figures S2A,B).
In addition to increases in Archaea, cell growth was observed for SAR202, which increased by between day 0 and day 4 for all amended samples (Figure 5 and Supplementary Figures S2A–D). Cell densities of Rhodobacteracea in the unamended controls did not exceed an average of 1.55 × 107 cells L–1 (Figure 5A), but demonstrated a significant increase in cell abundance in all of the amended treatments (paired t-tests for growth between day 4 and 6 for each amendment, p < 0.05 for all tests; Supplementary Figures S2B–D). Cell growth was highest in the treatments with added inorganic nutrients (GNP = 1.30 × 108 cells L–1 and LNP = 1.08 × 108cells L–1, respectively, on day 6), whereas it only reached a maximum of 6.33 × 107 cells L–1 (on day 6) in the Lignin treatments. Notably, in all of the amended samples and the controls, the maximum growth for Rhodobacteracea did not occur until day 6. Additionally, growth of Flavobacteriales was investigated, but although it was detected at all time points for all treatments, in the lignin amended treatments, it reached a minimum on day 4, while like Rhodobacteracea, it reached a maximum on day 6.
Discussion
Microbial Growth and DOC Utilization
The GNP amended treatments yielded the greatest growth response and coincident drawdown of DOC (Figures 4, 5B) of all treatments, which is consistent with previous studies that have demonstrated that copiotrophic populations grow rapidly in the presence of a labile carbon source and available macronutrients (, , ; ; ; ). Microbial cell growth was also evidenced in the Lignin and LNP treatments (Figure 1), which suggests that some marine heterotrophs were able to utilize a portion of the added lignin; however, the question of which components of the lignin molecule were used remains unresolved. Unlike the GNP amendments, where DOC drawdown accounted for all of the added organic carbon, no significant DOC drawdown was observed in the Lignin or LNP amendments (Figure 4). However, DOC removal in the Lignin treatment was greater than that observed in the unamended control (Student’s t-test, p = 0.174), where DOC removal was also not resolved. As the DOC measurements were not resolvable below the micromolar level, we cannot say definitively whether DOC was utilized in the Lignin and LNP treatments. However, the lignin amended dark incubations revealed greater prokaryote growth compared to the unamended controls, suggesting that marine heterotrophs were able to utilize a portion of the DOM in the lignin amended treatments to support anabolism. A such, we utilized CDOM measurements to further investigate changes in DOM.
Changes in CDOM Absorbance
The microbial growth over the first few days, coincident with an increase in CDOM absorption at 325 nm in all amended samples (Figures 1, 2), is consistent with the idea that as organic matter is remineralized by microbes, chromophoric byproducts are produced, either through direct release or by transformation of organic matter to a molecular structure with enhanced chromophoric properties (). While CDOM production in the GNP amendments was followed by a decrease in the 325 nm absorption coefficient to levels similar to those at T0 (Figure 2, top panel; Supplementary Table S2), the CDOM absorbance at 325 nm in the Lignin and LNP amendments decreased to values 20% lower than the T0 values (Figure 2, bottom panel), further suggesting an alteration of chromophoric organic matter in the lignin amended treatments. Nevertheless, the final CDOM values between 250 and 450 nm in the Lignin and LNP treatments remained well above those of the unamended control (Figure 3), suggesting that portions of the lignin remained unaltered in these treatments. Absorbance measurements alone do not conclusively prove that microbial remineralization of the amended carbon took place in the incubated samples; however, a decrease in CDOM absorbance coincident with an increase in some lineages of bacteria and archaea in the lignin amended treatments does suggest that these microbes altered the chromophoric components of lignin throughout the incubation. Thaumarchaeota and Euryarchaeota, in particular, reached maximum cell abundances in the Lignin and LNP treatments in conjunction with minimum CDOM 325 nm values (Figures 2, 5). Specifically, there was a coinciding decrease in spectral absorbance at the lignin excitation maxima (∼340 nm; Figure 2, bottom panel; Figure 3, top panel and middle panel) as Thaumarchaeota cell abundance increased (Supplementary Figures S2C,D). The more pronounced growth of Euryarchaeota in the LNP treatment than in the Lignin treatment (2.91 × 107 cells L–1 vs. 1.27 × 107 cells L–1, respectively; Supplementary Figures S2C,D) may indicate the additional requirements of nitrogen and phosphorus for cell growth of this lineage. Furthermore, metagenomic analysis has revealed that Marine Group II Euryarcheota have protein degradation pathways (; ) and may have a role in DON transformation.
The featureless absorption spectrum of CDOM in natural waters, which exhibits an exponential increase in absorbance with decreasing wavelength, has been the subject of much discussion, as it has been shown that this spectrum cannot be the result of the superposition of spectra of individual chromophores (; , ), and it has been suggested that chemical alteration of chromophores must occur to produce the CDOM spectrum. , proposed that charge-transfer complexes formed from the partial oxidation of lignin precursors were responsible for the absorbance spectra in laboratory experiments using dissolved humic substances and that these complexes would be found in natural waters, where donor-acceptor complexes arise. Upon initial amendment of the incubations with lignin, a distinct absorbance peak was observed at ∼340 nm (which is not observed in monthly in situ CDOM measurements at BATS), along with a broad “shoulder” between 265–280 nm (Figure 3, top and middle panels). However, as absorbance was measured over time, the enhanced absorbance at these wavelengths gradually decreased, resulting in a spectrum which more closely resembled the shape of the featureless “textbook” CDOM spectrum, as described by and (Figure 3, bottom panel). The timing of this CDOM transformation coincided with an increase in cell abundance and could indicate that microbes may play a previously unknown role in creating the charge-transfer complexes proposed to govern CDOM absorbance.
Taxon-Specific Microbial Response
The taxa targeted by the probes used in this study have been found within the North Atlantic Subtropical Gyre, which is dominated by the oligotrophic alphaproteobacterium, SAR11, representing ∼40% of the prokaryotic community in the euphotic zone (, ). In the amended incubations, the absolute cell abundance of SAR11 was relatively constant over time (Supplementary Figures S2C,D), while several other lineages increased. Members of the SAR202 clade of marine bacteria exhibited a maximum increase in cell abundance (Figure 5 and Supplementary Figures S2C,D) in the amended treatments at the same time when CDOM decreased. SAR202 have been estimated to account for ∼30% of all bathypelagic plankton (; ), and studies have shown that multiple SAR202 clade genomes encode families of enzymes such as flavin-dependent monooxygenases, enolases and ring-hydroxylating dioxygenases, that may allow for the degradation of recalcitrant organic matter (; ). As such, the ability of SAR202 to degrade lignin warrants further study. While maximum cell abundance was reached on day 4 for most of the targeted lineages, the copiotrophic alphaproteobacterium Rhodobacteracea exhibited maximum growth on day 6 in all treatments (Supplementary Figures S2A–D). The maxima for Rhodobacteracea occurred after CDOM had already decreased and therefore, it remains unclear if members of this lineage were consuming lignin directly or if they were utilizing byproducts of its remineralization. Members of Rhodobacteracea have been shown to comprise approximately 5% of the total bacterioplankton community within the surface 300 m of the Sargasso Sea (). However, in this study, Rhodobacteracea comprised up to 13% of the population in the Lignin amendment and 17% of the total population in the LNP amendment. The maximum cell abundance of Rhodobacteracea in the LNP treatments was almost double that observed in the Lignin treatments (Figure 5), and was comparable to the maximum observed for Rhodobacteracea in the GNP treatments (Figure 5). This presents the possibility that Rhodobacteracea growth was influenced by added N and P.
It is interesting to note that Alteromonas exhibited higher cell abundances in the LNP than the Lignin treatment (Figure 5 and Supplementary Figures S2C,D), as Alteromonas is a copiotrophic organism, and has been shown to outcompete other species in the presence of labile carbon substrates, especially when nitrogen and phosphorus are not limited. Additionally, the added N and P in the LNP amended samples could account for the fact that Alteromonas cell abundance was almost double that in the Lignin treatment. Several possible mechanisms could account for its growth in the lignin amended samples. It is possible that Alteromonas was able to utilize the DOM that was originally present in the seawater before the carbon substrate was added. Conversely, labile breakdown products could have been generated by the degradation of lignin by other lineages, which Alteromonas was then able to utilize. Alternately, Alteromonas may be able to directly utilize a portion of the lignin macromolecule.
The growth of Thaumarchaeota in both lignin amended treatments is of particular interest, as this lineage is known for its demonstrated ability to oxidize ammonia and as such, is largely known as a chemoautotroph (; ). However, it has also been shown to take up amino acids, and genes linked to the utilization of organic carbon have been identified in the same lineage of marine Archaea as that identified in our study, thus revealing the potential mixotrophic nature of this lineage (; ; ). It has been reported that approximately 60% of Archaea in samples collected from the Mediterranean and Pacific, including Marine GI Crenarchaeota-1 and Crenarchaeota-2, could uptake dissolved amino acids (). This study concluded that a portion of Crenarchaeota are heterotrophic and coexist with bacteria. Furthermore, a 2006 study of Crenarchaeota in the North Atlantic also demonstrated their ability to uptake amino acids, with certain lineages actively incorporating D- and L- aspartic acid (). Notably, while the Crenarchaeota phylum has since been divided, creating the new Thaumarchaeota phylum, the oligonucleotide probe used by both and to identify Crenarchaeota (Cren537 [5-TGACCACTTGAGGTGCTG-3′]) is the same as was used in the present study to detect Thaumarchaeota.
Proposed Mechanisms of Lignin Degradation by Bacteria and Archaea
As a component of cell walls in vascular plants, lignin serves to protect plants from pathogens and is therefore naturally resistant to degradation from many types of microbes. The complex nature of the cross-linked lignin macromolecule, comprising of three phenolic monomers, creates a natural hindrance to radical decoupling and therefore can impede microbial degradation of the compound. While enzymatic breakdown of lignin has not previously been demonstrated for oceanic microbes, particularly archaea, two species of bacteria have been found to degrade lignin in marine pulp mill effluent (), and there are groups of soil-associated bacteria, as well as many fungi (; ; ; ) that use peroxidases and laccases () to degrade lignin. Although fungi typically utilize heme-based lignin peroxidase (LiP) or manganese peroxidase (MnP) to depolymerize lignin, the mechanisms for microbial degradation of lignin are less clear (; ). It is interesting to note that archaeal ammonia monooxygenases belong to a family of copper-containing membrane-bound monooxygenases (CuMMOs) () and the ability of these enzymes to oxidize lignin into phenolic breakdown products has not yet been investigated.
Conclusion
This study provides evidence that transformation of terrigenous DOM in dark incubations by marine heterotrophs is taxon-specific, with some marine lineages of bacteria and archaea altering what were previously thought to be biologically recalcitrant carbon sources, on a timescale of days. Furthermore, these results indicate that in addition to the photochemical interactions speculated as being responsible for the distinct CDOM absorbance spectrum, that biologically driven alterations of chromophores may occur.
Additionally, this study demonstrates a clear growth response to added lignin by multiple bacteria and archaea lineages, including Thaumarchaeota, warranting further investigation as to whether this is a primary response consistent with previous evidence that Thaumarcheota is mixotrophic, or a secondary response to reduced remineralization byproducts of other organisms. These results also provide further insight into the mechanisms by which cycles of terrigenous and autochthonous DOM in the pelagic ocean are inter-connected. The subsequent in situ production of CDOM (at 325 nm) in lignin-amended incubations of marine heterotrophs suggests a cycle in which terrigenous and autochthonous CDOM are linked, with specific lineages of bacteria and archaea capable of transforming portions of lignin and then producing autochthonous CDOM, either through the production of DOM or the molecular rearrangement of existing DOM. Because up to 25% of lignin remains undegraded by photolysis, and lignin has been detected in the bathypelagic region of the North Atlantic, questions remain about the role of microbial degradation in the deep ocean. It is increasingly evident that the chemical composition and optical properties of DOM are linked to microbial community composition, and that these factors play a large role in determining the lability and fate of organic carbon in the ocean. It is by better understanding these links that a more comprehensive view of DOM cycling in the ocean can be achieved.
Statements
Data availability statement
All datasets generated for this study are included in the article/Supplementary Material.
Author contributions
NM, CC, RP, and JM-B designed the study. NM, SL, JM-B, and RP carried out the experimentation and analysis. EA, MG, and NN provided expertise and advice in analyzing data and interpreting the results. NM, EA, CC, MG, SL, JM-B, NN, and RP contributed to the writing of the manuscript.
Funding
This study was supported in part by the NSF Oceanic Microbial Observatory grant (award OCE-0801991) as well as by the Canadian Associates of BIOS (CABIOS) and Simons Foundation International’s BIOS-SCOPE Program. Additionally, travel support was granted by Future Ocean.
Acknowledgments
We would like to thank Ellie Halewood and Keri Opalk for their valuable contributions and assistance with DOC analysis, as well as Dr. Rod Johnson and BATS technicians for facilitating at-sea sample collection, and the captain and crew of the R/V Atlantic Explorer.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2019.00743/full#supplementary-material
References
1
BatesN. R.JohnsonR. J. (2016). Bermuda Atlantic Time Series (Bats) Bottle Data.St. George’s: Institute of Ocean Sciences.
2
BloughN. V.Del VecchioR. (2002). “Chromophoric DOM in the coastal environment,” in Biogeochemistry of Marine Dissolved Organic Matter, edsHansellD.CarlsonC. (San Diego, CA: Academic Press).
3
BoyleE. S.GuerrieroN.ThialletA.Del VecchioR.BloughN. V. (2009). Optical properties of humic substances and Cdom: relation to structure.Environ. Sci. Technol.432262–2268. 10.1021/es803264g
4
BrownM. E.ChangM. C. Y. (2014). Exploring bacterial lignin degradation.Curr. Opin. Chem. Biol.191–7. 10.1016/j.cbpa.2013.11.015
5
BrownM. E.WalkerM. C.NakashigeT. G.ChangM. C. (2011). Discovery and characterization of heme enzymes from unsequenced bacteria: application to microbial lignin degradation.J. Am. Chem. Soc.13318006–18009. 10.1021/ja203972q
6
BuggT. D. H.AhmadM.HardimanE. M.RahmanpourR. (2011). Pathways for degradation of lignin in bacteria and fungi.Nat. Prod. Rep.281883–1896. 10.1039/c1np00042j
7
CarlsonC.HansellD. (2015). “Doc sources, sinks, reactivity, and budgets,” in Biogeochemistry of Marine Dissolved Organic Matter, 2nd Edn, edsHansellD.CarlsonC. (San Diego, CA: Academic Press).
8
CarlsonC. A.DucklowH. W. (1996). Growth of bacterioplankton and consumption of dissolved organic carbon in the Sargasso Sea.Aquat. Microb. Ecol.1069–85. 10.3354/ame010069
9
CarlsonC. A.DucklowH. W.HansellD. A.SmithW. O. (1998). Organic carbon partitioning during spring phytoplankton blooms in the Ross Sea polynya and the Sargasso Sea.Limnol. Oceanogr.43375–386. 10.4319/lo.1998.43.3.0375
10
CarlsonC. A.GiovannoniS. J.HansellD. A.GoldbergS. J.ParsonsR.OteroM. P.et al (2002). Effect of nutrient amendments on bacterioplankton production, community structure, and Doc utilization in the northwestern Sargasso Sea.Aquat. Microb. Ecol.3019–36. 10.3354/ame030019
11
CarlsonC. A.GiovannoniS. J.HansellD. A.GoldbergS. J.ParsonsR.VerginK. (2004). Interactions among dissolved organic carbon, microbial processes, and community structure in the mesopelagic zone of the northwestern Sargasso Sea.Limnol. Oceanogr.491073–1083. 10.4319/lo.2004.49.4.1073
12
CarlsonC. A.HansellD. A.NelsonN. B.SiegelD. A.SmethieW. M.KhatiwalaS.et al (2010). Dissolved organic carbon export and subsequent remineralization in the mesopelagic and bathypelagic realms of the North Atlantic basin.Deep Sea Res. Part Ii-Top. Stud. Oceanogr.571433–1445.10.1016/j.dsr2.2010.02.013
13
CarlsonC. A.MorrisR.ParsonsR.GiovannoniS. J.VerginK. (2009). Seasonal dynamics of Sar11 populations in the euphotic and mesopelagic zones of the northwestern Sargasso Sea.Isme J.3283–295. 10.1038/ismej.2008.117
14
ColeJ. R.WangQ.CardenasE.FishJ.ChaiB.FarrisR. J.et al (2009). The ribosomal database project: improved alignments and new tools for rrna analysis.Nucleic Acids Res.37D141–D145. 10.1093/nar/gkn879
15
CottrellM. T.KirchmanD. L. (2000). Natural assemblages of marine proteobacteria and members of the Cytophaga-Flavobacter cluster consuming low and high molecular weight dissolved organic matter.Appl. Environ. Microbiol.661692–1697. 10.1128/aem.66.4.1692-1697.2000
16
CraggS. M.BeckhamG. T.BruceN. C.BuggT. D. H.DistelD. L.DupreeP.et al (2015). Lignocellulose degradation mechanisms across the tree of life.Curr. Opin. Chem. Biol.29108–119. 10.1016/j.cbpa.2015.10.018
17
Del VecchioR.BloughN. V. (2002). Photobleaching of chromophoric dissolved organic matter in natural waters: kinetics and modeling.Mar. Chem.78231–253. 10.1016/s0304-4203(02)00036-1
18
Del VecchioR.BloughN. V. (2004). On the origin of the optical properties of humic substances.Environ. Sci. Technol.383885–3891. 10.1021/es049912h
19
EwartC. S.MeyersM. K.WallnerE. R.McgillicuddyD. J.CarlsonC. A. (2008). Microbial dynamics in cyclonic and anticyclonic mode-water eddies in the northwestern Sargasso Sea.Deep Sea Res. Part Ii-Top. Stud. Oceanogr.551334–1347. 10.1016/j.dsr2.2008.02.013
20
FashamM. J. R.BalinoB. M.BowlesM. C.AndersonR.ArcherD.BathmannU.et al (2001). A new vision of ocean biogeochemistry after a decade of the joint global ocean flux study (Jgofs).AMBIO104–31.
21
GlennJ. K.GoldM. H. (1985). Purification and chracterization of an extracellular Mn(Ii)-dependent peroxidase from the lignin-degrading basidomycete, phanerochaete chrysosporium.Arch. Biochem. Biophys.242329–341. 10.1016/0003-9861(85)90217-6
22
GlennJ. K.MorganM. A.MayfieldM. B.KuwaharaM.GoldM. H. (1983). An extracellular H2O2-requiring enzyme preparation involved in lignin biodegradation by the white rot basidiomycete phanerochaete-chrysosporium.Biochem. Biophys. Res. Commun.1141077–1083. 10.1016/0006-291x(83)90672-1
23
GoldbergS. J.CarlsonC. A.BrzezinskiM.NelsonN. B.SiegelD. A. (2011). Systematic removal of neutral sugars within dissolved organic matter across ocean basins.Geophys. Res. Lett.38:L17606.
24
GongG.WooH. M.UmY.ParkT. H. (2014). Lignin utilization by Bacillus sp. associated with the growth enhancement and the molecular weight distribution change of lignin.New Biotechnol.31:S102.
25
GonzalesJ. M.MayerF.MoranM. A.HodsonR. E.WhitmanW. B. (1997). Microbulbifer hydrolyticus gen. nov., sp. nov., and Marinobacterium georgiense gen. nov., sp. nov., two marine bacteria from a lignin-rich pulp mill waste enrichment community.Int. J. Syst. Bacteriol.47369–376. 10.1099/00207713-47-2-369
26
HallamS. J.MincerT. J.SchleperC.PrestonC. M.RobertsK.RichardsonP. M.et al (2006). Pathways of carbon assimilation and ammonia oxidation suggested by environmental genomic analyses of marine crenarchaeota.PLoS Biol.4:e95. 10.1371/journal.pbio.0040437
27
HansellD. (2005). Dissolved organic carbon reference material program.Eos Trans. Am. Geophys. Union86:318. 10.1016/j.aca.2013.06.015
28
HedgesJ. I.KeilR. G.BennerR. (1997). What happens to terrestrial organic matter in the ocean?Organ. Geochem.27195–212. 10.1016/s0146-6380(97)00066-1
29
HernesP. J.BennerR. (2003). Photochemical and microbial degradation of dissolved lignin phenols: implications for the fate of terrigenous dissolved organic matter in marine environments.J. Geophys. Res. Oceans108:3291.
30
HernesP. J.BennerR. (2006). Terrigenous organic matter sources and reactivity in the North Atlantic Ocean and a comparison to the Arctic and Pacific oceans.Mar. Chem.10066–79. 10.1016/j.marchem.2005.11.003
31
HuangX. F.SanthanamN.BadriD. V.HunterW. J.ManterD. K.DeckerS. R.et al (2013). Isolation and characterization of lignin-degrading bacteria from rainforest soils.Biotechnol. Bioeng.1101616–1626. 10.1002/bit.24833
32
IversonV.MorrisR. M.FrazarC. D.BerthiaumeC. T.MoralesR. L.ArmbrustE. V. (2012). Untangling genomes from metagenomes: revealing an uncultured class of marine Euryarchaeota.Science335587–590. 10.1126/science.1212665
33
JanuszG.PawlikA.SulejJ.Świderska-BurekU.Jarosz-WilkołazkaA.PaszczyńskiA. (2017). Lignin degradation: microorganisms, enzymes involved, genomes analysis and evolution.FEMS Microbiol. Rev.41941–962. 10.1093/femsre/fux049
34
KinseyJ.CorradinoG.ZiervogelK.SchnetzerA.OsburnC. (2017). Formation of chromophoric dissolved organic matter by bacterial degradation of phytoplankton-derived aggregates.Front. Mar. Sci.4:430. 10.3389/fmars.2017.00430
35
KochB. P.WittM. R.EngbrodtR.DittmarT.KattnerG. (2005). Molecular formulae of marine and terrigenous dissolved organic matter detected by electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry.Geochim. Cosmochim. Acta693299–3308. 10.1016/j.gca.2005.02.027
36
KonnekeM.BernhardA. E.De La TorreJ. R.WalkerC. B.WaterburyJ. B.StahlD. A. (2005). Isolation of an autotrophic ammonia-oxidizing marine archaeon.Nature437543–546. 10.1038/nature03911
37
LandryZ.SwanB. K.HerndlG. J.StepanauskasR.GiovannoniS. J. (2017). SAR202 genomes from the dark ocean predict pathways for the oxidation of recalcitrant dissolved organic matter.mBio8:e00413–17. 10.1128/mBio.00413-17
38
LavorivskaL.BoyerE. W.DewalleD. R. (2016). Atmospheric deposition of organic carbon via precipitation.Atmos. Environ.146153–163. 10.1016/j.atmosenv.2016.06.006
39
LiuS.ParsonsR.OpalkK.BaetgeN.GiovannoniS.BolanosL.et al (2019). Different carboxyl-rich alicyclic molecules (CRAM) proxy compounds select distinct bacterioplankton for oxidation of DOM in the mesopelagic Sargasso Sea.Limnol. Oceanogr.(in press).
40
LoyA.MaixnerF.WagnerM.HornM. (2007). probeBase - an online resource for rrna-targeted oligonucleotide probes: new features 2007.Nucleic Acids Res.35D800–D804.
41
McNallyS. P.ParsonsR. J.SantoroA. E.ApprillA. (2017). Multifaceted impacts of the stony coral Porites astreoides on picoplankton abundance and community composition.Limnol. Oceanogr.62217–234. 10.1002/lno.10389
42
MopperK.KieberD. J.StubbinsA. (2015). “Marine photochemistry of organic matter,” in Biogeochemistry of Dissolved Organic Matter, 2nd Edn, edsHansellD.CarlsonC. (San Diego, CA: Academic Press).
43
MorrisR. M.ConnonS. A.RappeM.VerginK. L.SieboldW. A.CarlsonC. A.et al (2002a). High cellular abundance of the Sar11 picoplankton clade in seawater.Nature420806–809.
44
MorrisR. M.RappeM. S.ConnonS. A.VerginK. L.SieboldW. A.CarlsonC. A.et al (2002b). Sar11 clade dominates ocean surface bacterioplankton communities.Nature420806–810. 10.1038/nature01240
45
MorrisR. M.RappeM.UrbachE.ConnonS. A.GiovannoniS. J. (2004). Prevalence of the Chloroflexi-related Sar202 bacterioplankton cluster throughout the mesopelagic zone and deep ocean.Appl. Environ. Microbiol.702836–2842. 10.1128/aem.70.5.2836-2842.2004
46
NelsonC. E.CarlsonC. A. (2012). Tracking differential incorporation of dissolved organic carbon types among diverse lineages of Sargasso Sea bacterioplankton.Environ. Microbiol.141500–1516. 10.1111/j.1462-2920.2012.02738.x
47
NelsonN. B.CarlsonC. A.SteinbergD. K. (2004). Production of chromophoric dissolved organic matter by Sargasso Sea microbes.Mar. Chem.89273–287. 10.1016/j.marchem.2004.02.017
48
NelsonN. B.SiegelD. A. (2002). “Chromophoric Dom in the open ocean,” in Biogeochemistry of Marine Dissolved Organic Matter, ed.HansellD. A. C. C. (San Diego, CA: Academic Press).
49
NelsonN. B.SiegelD. A.CarlsonC. A.SwanC.SmethieW. M.KhatiwalaS. (2007). Hydrography of chromophoric dissolved organic matter in the North Atlantic.Deep Sea Res. Part I Oceanogr. Res. Pap.54710–731. 10.1016/j.dsr.2007.02.006
50
OrsiW. D.SmithJ. M.LiuS.LiuZ.SakamotoC. M.WilkenS.et al (2016). Diverse, uncultivated bacteria and archaea underlying the cycling of dissolved protein in the ocean.ISME J.102158–2173. 10.1038/ismej.2016.20
51
OuverneyC.FuhrmanJ. (2000). Marine planktonic archaea take up amino acids.Appl. Environ. Microbiol.664829–4833. 10.1128/aem.66.11.4829-4833.2000
52
ParsonsR. J.BreitbartM.LomasM. W.CarlsonC. A. (2011). Ocean time-series reveals recurring seasonal patterns of virioplankton dynamics in the northwestern Sargasso Sea.ISME J.20111–12. 10.1038/ismej.2011.101
53
ParsonsR. J.NelsonC. E.CarlsonC. A.DenmanC. C.AnderssonA. J.KledzikA. L.et al (2015). Marine bacterioplankton community turnover within seasonally hypoxic waters of a subtropical sound: Devil’s Hole, Bermuda.Environ. Microbiol.173481–3499. 10.1111/1462-2920.12445
54
PengR. H.XiongA. S.XueY.FuX. Y.GaoF.ZhaoW.et al (2008). Microbial biodegradation of polyaromatic hydrocarbons.FEMS Microbiol. Rev.32927–955. 10.1111/j.1574-6976.2008.00127.x
55
PesterM.SchleperC.WagnerM. (2011). The Thaumarchaeota: an emrging view of their phylogeny and ecophysiology.Curr. Opin. Microiol.14300–306. 10.1016/j.mib.2011.04.007
56
PorterK. G.FeigY. S. (1980). The use of Dapi for identifying and counting aquatic microflora.Limnol. Oceanogr.25943–948. 10.4319/lo.1980.25.5.0943
57
StedmonC. A.NelsonN. B. (2015). “The optical properties of Dom in the ocean,” in Bioegeochemistry of Marine Dissolved Organic Matter, ed.HansellD. A. C. C. (San Diego, CA: Academic Press).
58
SteinbergD. K.CarlsonC. A.BatesN.JohnsonR. J.MichaelsF. A. (2001). Overview of the Us Jgofs Bermuda atlantic time-series study (Bats): a decade-scale look at ocean biology and biogeochemistry.Deep Sea Res. Ii481405–1447. 10.1016/s0967-0645(00)00148-x
59
TeiraE.Van AkenH.VethC.HerndlG. J. (2006). Archaeal uptake of enantiomeric amino acids in the meso- and bathypelagic waters of the North Atlantic.Limnol. Oceanogr.5160–69. 10.4319/lo.2006.51.1.0060
60
TienM.KirkT. (1984). Lignin-degrading enzyme from phanerochaete chrysosporium: purification, characterization, and catalytic properties of a unique H2O2-requiring oxygenase.Proc. Natl. Acad. Sci. U.S.A.812280–2284. 10.1073/pnas.81.8.2280
61
VarelaM. M.Van AkenH. M.HerndlG. J. (2008). Abundance and activity of Chloroflexi-type Sar202 bacterioplankton in the meso- and bathypleagic waters of the (sub)tropical Atlantic.Environ. Microbiol.101903–1911. 10.1111/j.1462-2920.2008.01627.x
62
VarelaM. M.Van AkenH. M.SintesE.ReinthalerT.HerndlG. J. (2011). Contribution of Crenarchaeota and Bacteria to autotrophy in the North Atlantic interior.Environ. Microbiol.131524–1533. 10.1111/j.1462-2920.2011.02457.x
63
WalkerC. B.de la TorreJ. R.KlotzM. G.UrakawaH.PinelN.ArpD. J.et al (2010). Nitrosopumilus maritimus genome reveals unique mechanisms for nitrification and autotrophy in globally distributed marine crenarchaea.Proc. Natl. Acad. Sci. U.S.A.1078818–8823. 10.1073/pnas.0913533107
64
WilliamsP. M.DruffelE. R. M. (1987). Radiocarbon in dissolved organic-matter in the Central North Pacific-Ocean.Nature330246–248. 10.1038/330246a0
Summary
Keywords
thaumarchaeota, lignin, bacterioplankton abundance, CDOM, fluorescent in situ hybridisation (FISH)
Citation
McDonald N, Achterberg EP, Carlson CA, Gledhill M, Liu S, Matheson-Barker JR, Nelson NB and Parsons RJ (2019) The Role of Heterotrophic Bacteria and Archaea in the Transformation of Lignin in the Open Ocean. Front. Mar. Sci. 6:743. doi: 10.3389/fmars.2019.00743
Received
15 August 2019
Accepted
14 November 2019
Published
10 December 2019
Volume
6 - 2019
Edited by
Carol Robinson, University of East Anglia, United Kingdom
Reviewed by
Xosé Anxelu G. Morán, King Abdullah University of Science and Technology, Saudi Arabia; David L. Kirchman, University of Delaware, United States
Updates
Copyright
© 2019 McDonald, Achterberg, Carlson, Gledhill, Liu, Matheson-Barker, Nelson and Parsons.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Natasha McDonald, natashamc1@gmail.com
This article was submitted to Marine Biogeochemistry, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.