Formation of Chromophoric Dissolved Organic Matter by Bacterial Degradation of Phytoplankton-Derived Aggregates

Kinsey, Joanna D.; Corradino, Gabrielle; Ziervogel, Kai; Schnetzer, Astrid; Osburn, Christopher L.

doi:10.3389/fmars.2017.00430

ORIGINAL RESEARCH article

Front. Mar. Sci., 04 January 2018

Sec. Aquatic Microbiology

Volume 4 - 2017 | https://doi.org/10.3389/fmars.2017.00430

This article is part of the Research Topic Extracellular enzymes in aquatic environments: exploring the link between genomic potential and biogeochemical consequences. View all 13 articles

Formation of Chromophoric Dissolved Organic Matter by Bacterial Degradation of Phytoplankton-Derived Aggregates

$\r\nJoanna D. Kinsey*$ Joanna D. Kinsey¹^*

Gabrielle Corradino¹

Kai Ziervogel²

Astrid Schnetzer¹

Christopher L. Osburn¹

¹Department of Marine, Earth, and Atmospheric Sciences, North Carolina State University, Raleigh, NC, United States
²Ocean Process Analysis Laboratory, Institute for the Study of Earth, Oceans, and Space, University of New Hampshire, Durham, NH, United States

Organic matter produced and released by phytoplankton during growth is processed by heterotrophic bacterial communities that transform dissolved organic matter into biomass and recycle inorganic nutrients, fueling microbial food web interactions. Bacterial transformation of phytoplankton-derived organic matter also plays a poorly known role in the formation of chromophoric dissolved organic matter (CDOM) which is ubiquitous in the ocean. Despite the importance of organic matter cycling, growth of phytoplankton and activities of heterotrophic bacterial communities are rarely measured in concert. To investigate CDOM formation mediated by microbial processing of phytoplankton-derived aggregates, we conducted growth experiments with non-axenic monocultures of three diatoms (Skeletonema grethae, Leptocylindrus hargravesii, Coscinodiscus sp.) and one haptophyte (Phaeocystis globosa). Phytoplankton biomass, carbon concentrations, CDOM and base-extracted particulate organic matter (BEPOM) fluorescence, along with bacterial abundance and hydrolytic enzyme activities (α-glucosidase, β-glucosidase, leucine-aminopeptidase) were measured during exponential growth and stationary phase (~3–6 weeks) and following 6 weeks of degradation. Incubations were performed in rotating glass bottles to keep cells suspended, promoting cell coagulation and, thus, formation of macroscopic aggregates (marine snow), more similar to surface ocean processes. Maximum carbon concentrations, enzyme activities, and BEPOM fluorescence occurred during stationary phase. Net DOC concentrations (0.19–0.46 mg C L⁻¹) increased on the same order as open ocean concentrations. CDOM fluorescence was dominated by protein-like signals that increased throughout growth and degradation becoming increasingly humic-like, implying the production of more complex molecules from planktonic-precursors mediated by microbial processing. Our experimental results suggest that at least a portion of open-ocean CDOM is produced by autochthonous processes and aggregation likely facilitates microbial reprocessing of organic matter into refractory DOM.

Introduction

Oceanic primary production (PP) contributes between 35 and 65 petagrams (Pg) of carbon annually to global net PP (Field et al., 1998; del Giorgio and Duarte, 2002; Carr et al., 2006; Chavez et al., 2011), mainly through autochthonous production by photosynthetic phytoplankton. Phytoplankton release 2–50% of photosynthetically fixed carbon as dissolved organic matter (DOM) through active exudation and passive leakage (Thornton, 2014), contributing to the largest reservoir of reduced carbon on earth (662 Pg C) (Hansell et al., 2009). Additional phytoplanktonic DOM is released through viral infection, sloppy feeding by predators such as zooplankton, and cell death, with zooplankton contributing their own DOM through excretion of organic matter and fecal pellets (Urban-Rich et al., 2006; Steinberg et al., 2008; Saba et al., 2011). A significant, yet poorly characterized, portion of primary production is exported out of the euphotic zone and into the twilight zone (Siegel et al., 2016). The vertical export has been estimated to range from 5 to more than 12 Pg C y⁻¹ (e.g., Boyd and Trull, 2007; Henson et al., 2011) spanning annual anthropogenic CO₂ emissions (~7.0 Pg C y⁻¹) (Siegenthaler and Sarmiento, 1993). Sinking aggregates (or marine snow) are an amalgam of intact phytoplankton cells, detritus, fecal pellets, silica and/or carbonate frustules (that provide ballast) and transparent exopolymeric material (TEP) (Alldredge and Silver, 1988). These aggregates exhibit a wide variability with respect to chemical composition (Wakeham and Lee, 1989, 1993; Minor et al., 2003) and may result in significant export of organic matter to depth (e.g., Hansell et al., 2009; Carlson et al., 2010) and contribution to refractory DOM (RDOM) formation (Lechtenfeld et al., 2015).

DOM is a complex mixture of organic compounds that play an important role in marine biogeochemical cycling. The optically active fraction, chromophoric DOM (CDOM), is essential in both physical and biological processes, controlling light attenuation and photochemical reactions in the surface ocean, impacting the depth of primary production and screening out harmful ultraviolet light (e.g., Arrigo and Brown, 1996; Blough and Del Vecchio, 2002; Mopper et al., 2015). A small fraction of DOM fluoresces, allowing for the identification of DOM sources using excitation-emission matrices (EEMs). Two types of fluorescence are typically described, amino acid-like compounds and humic-like compounds (Coble, 2007), and have been attributed to autochthonous (e.g., primary production) and allochthonous (e.g., terrestrial) production, respectively (Coble, 1996; Yamashita and Tanoue, 2003). Fluorescence of base-extracted particulate organic matter (BEPOM), has been established as a relatively new technique that allows for the optical properties and sources of intracellular fluorophores in POM to be determined (Brym et al., 2014). This work has suggested strong discrete fluorophores including amino acids in estuarine and coastal ocean POM, while studies investigating CDOM production from phytoplankton have confirmed strong amino acid-like fluorescence in axenic cultures with additional fluorescence in the region of humic-like fluorescence for some species (Romera-Castillo et al., 2010; Fukuzaki et al., 2014). Additionally, humic-like fluorescence has been found to increase in degradation studies (Rochelle-Newall and Fisher, 2002; Stedmon and Markager, 2005; Biers et al., 2007; Shimotori et al., 2009; Romera-Castillo et al., 2011) and with apparent oxygen utilization (AOU), especially in the Pacific and Southern oceans (Chen and Bada, 1992; Yamashita et al., 2007; Yamashita and Tanoue, 2008; Jørgensen et al., 2011; Catalá et al., 2015). This suggests a linkage to microbial oxidation and organic matter degradation (Hayase and Shinozuka, 1995; Yamashita et al., 2007; Jørgensen et al., 2011) and may be especially important as an indicator for biogeochemical cycling since DOM serves as a substrate for heterotrophic microbial remineralization of carbon and other trace elements. Microbes utilize extracellular enzymes to catalyze high-molecular weight DOM into smaller compounds that can be transported across cell membranes of bacteria (Arnosti, 2011) or through osmotic uptake by fungi and labyrinthulomycetes (Bochdansky et al., 2017). The carbon is then incorporated into biomass, respired to CO₂, or excreted as dissolved organic carbon (DOC) in the form of metabolically transformed products (Arnosti, 2011).

Given the potential importance of heterotrophic bacterial processing of organic matter for both biogeochemical cycling and contribution to deep-sea fluorescence, we aimed to investigate the role of microbes in the formation of open-ocean CDOM from phytoplankton-derived POM to better understand sources of open-ocean fluorescence. Our study builds upon previous culture studies that investigated CDOM production by incorporating fluorescence of base-extracted POM and hydrolytic enzyme activities in addition to traditional CDOM fluorescence measurements. Roller bottle incubations were used to form planktonic aggregates of four monocultures, sampling at initial, exponential, and stationary growth phases, and after 6 weeks of degradation in the dark. Bacterial abundance and two classes of hydrolytic extracellular enzymes indicative of carbohydrate and peptide hydrolysis (glucosidase and leucine-aminopeptidase, respectively) were measured in conjunction with carbon concentrations and CDOM and BEPOM fluorescence. We demonstrate that the production of CDOM fluorescence is directly related to microbial processing of phytoplankton-derived material.

Materials and Methods

Phytoplankton Growth Experiments and Sample Collection

Growth experiments were conducted with four phytoplankton taxa commonly found in coastal waters. Three strains were obtained from the NCMA Bigelow collection (https://ncma.bigelow.org/) including diatoms Skeletonema grethae (CCMP2801, isolated May 2004 off North Carolina) and Leptocylindrus hargravesii (CCMP1856, isolated Feb 1997 in the Gulf of Mexico), and one haptophyte Phaeocystis globosa (CCMP2754, isolated Aug 2003 in the Atlantic Bight). The fourth was a coastal non-axenic isolate of Coscinodiscus sp. (COS1, isolated April 2015 off North Carolina). Two growth experiments were conducted, one with Skeletonema, Leptocylindrus, and Phaeocystis spp. grown in parallel and the other with only Coscinodiscus sp. All cultures were grown in optically clear 0.2 μm-filtered artificial seawater (ASW) prepared with baked (550°C, 6 h) or ultrapure salts to minimize background CDOM fluorescence. A total of twelve 1.9 L borosilicate bottles per organism were inoculated with an initial concentration of ~21,000–66,000 cells mL⁻¹ with the exception of the much slower growing Coscinodiscus sp. with ~5 cells mL⁻¹. All inoculated bottles, along with control bottles of ASW containing no phytoplankton, were amended with f/20 growth media (one-tenth of f/2 media, Guillard and Ryther, 1962). Inoculated bottles and controls were grown on a Wheaton roller culture apparatus rotating at ~1 rpm (Shanks and Edmondson, 1989). Cultures were maintained at 18°C and a 12:12 L:D cycle at ~90 μEinstein m⁻² s⁻¹ under cool white fluorescent light in a culture incubator (Percival Scientific, Iowa).

Growth was monitored every other day by in vivo fluorescence and cell counts until particles formed. This information guided more comprehensive sampling for additional chemical and biological parameters (see further detail below) during each growth phase. Skeletonema sp., Leptocylindrus sp., and Phaeocystis sp. were grown for 3–4 weeks under the above conditions, sampling at the beginning of the experiment (T₀), during exponential growth (T₁), and early stationary phase (particle formation; T₂). Immediately after sampling for T₂, the remaining culture and control bottles were placed in the dark for 6 weeks to promote bacterial degradation of particles (T₃). The Coscinodiscus sp. experiment lasted 6 weeks with no degradation phase and was sampled at the beginning of the experiment (T₀), twice during exponential growth [early-exponential and late-exponential (T₁)], and early stationary phase (particle formation; T₂). Changes in cell densities were used to calculate specific growth rates of phytoplankton during exponential phases (Brand et al., 1981). Cell counts (mL⁻¹) were conducted using an Olympus BX53 compound microscope after preservation with acid Lugol's solution (5%) with settling volumes depending on growth phase (Utermöhl, 1958).

At each time point, three culture bottles of each organism and one control bottle (0.2 μm filtered ASW with f/20 nutrients) were sacrificed. Sample aliquots were gently vacuum filtered through pre-combusted (450°C, minimum 6 h) 0.7 μm (nominal mesh size) glass fiber filters. Filters were stored at −20°C until analyzed for particulate organic carbon (POC), POC/PON, and base-extracted POM (BEPOM) fluorescence followed by base-extracted POC (BEPOC). Filtrate was collected in polycarbonate bottles for absorbance and fluorescence and analyzed within 24 h or in pre-combusted borosilicate vials, acidified with 85% phosphoric acid (H₃PO₄) to pH 2, and stored at 4°C for dissolved organic carbon (DOC) analysis.

Optical Measurements

POM filters were base-extracted (BEPOM) in 10 mL of 0.1 M sodium hydroxide (NaOH) for 24 h at 4°C in the dark. After 24 h, extracts were neutralized with concentrated hydrochloric acid, then filtered using a 0.2 μm porosity Sterivex-GP polyethersulfone filter (EMD) (Brym et al., 2014). The filtrate was analyzed for absorbance and fluorescence in the same manner as CDOM.

Absorbance of CDOM and BEPOM was measured from 200 to 800 nm in a 1 or 10 cm quartz cell (Starna Cells, Inc.) on a Varian 300 UV spectrophotometer with ultrapure water as the reference blank for DOM samples and neutralized 0.1 M NaOH solution as the reference blank for BEPOM samples. Blank-corrected absorbance values were converted to Napierian absorption coefficients (α_λ) (Osburn et al., 2012).

CDOM and BEPOM fluorescence were measured using a Varian Eclipse spectrophotometer with excitation (Ex) from 240 to 450 nm at 5 nm intervals and emission (Em) from 300 to 600 nm every 2 nm. Samples were run with slit widths set at 5 nm for both Ex and Em modes, an integration time of 0.0125 s, scan rate of 2,400 nm m⁻¹, and at either 800 or 950 V, depending on sample response. Fluorescence measurements were corrected for excitation energy and emission detector response using correction factors supplied by the manufacturer and any inner-filter effects were corrected following Tucker et al. (1992). Final BEPOM fluorescence values were corrected for dilution and filtration volumes and both BEPOM and CDOM values were calibrated in quinine sulfate equivalents (QSE, where 1 QSE = 1 ppb quinine sulfate) (Lawaetz and Stedmon, 2009). Fluorescence results were concatenated into excitation-emission matrices (EEMs) for visualization as contour plots. The humification index (HIX) was determined by dividing the emission intensity in the 435–480 nm region by the intensity in the 300–345 nm region (Ohno, 2002).

Carbon Concentrations and POC/PON

DOC and BEPOC (from BEPOM samples) concentrations were measured on an OI Analytical 1030D TOC analyzer in combustion mode (Lalonde et al., 2014). Prior to analysis all acidified samples were sparged with ultrahigh purity argon. DOC concentrations were blank corrected using 18.2 MΩ ultrapure Milli-Q (Millipore) water and calibrated with caffeine standards (0–5 mg C L⁻¹) and Hansell deep sea reference (DSR) water. BEPOC concentrations were blank-corrected against ultrapure water and caffeine standards ranging from 0 to 20 mg C L⁻¹. Final BEPOC concentrations were corrected for extraction volume and initial filtration volume. Filters for POC concentrations and POC/PON ratios were dried overnight at 60°C followed by flash combustion to CO₂ and N₂ using a Thermo Flash 1112 elemental analyzer with acetanilide as a standard and corrected for volume filtered.

Bacteria Abundance

Bacterial cells were enumerated by flow cytometry (Gasol and Del Giorgio, 2000). At each sampling point 1 mL of experimental or control water was fixed with 0.1% glutaraldehyde (final concentration) for 10 min at room temperature in the dark, and stored at −80°C. Prior to analysis, thawed samples were pipetted through a cell strainer (Flowmi, 70 μm porosity) and stained with SYBR Green I for 15 min on ice in the dark. Counts were performed with a FACSCalibur flow cytometer (Becton-Dickinson) using fluorescent microspheres (Molecular Probes) of 1 μm in diameter as internal size standard. Cells were enumerated according to their right angle scatter and green fluorescence using the FloJo 7.6.1 software. This method quantifies free bacteria and does not account for bacteria attached to the aggregate particles.

Hydrolytic Enzyme Activity

Hydrolytic enzyme activities were determined using L-leucine-4-methylcoumarinyl-7-amide (MCA) hydrochloride, 4-methylumbelliferyl α-D-glucopyranoside, and 4-methylumbelliferone (MUF) β-D-glucopyranoside (Sigma-Aldrich) as substrate proxies for leucine-aminopeptidase, α-glucosidase, and β-glucosidase activities, respectively (Hoppe, 1983). For each bottle and substrate proxy, 196 μL of unfiltered experimental or control water was added in duplicate to a pure-grade black 96-well plate (Brand Life Sciences) containing a single substrate proxy at saturation levels (final concentration 200 μM). Fluorescence (excitation 370 nm, emission 440 nm) was measured in a Tecan Infinite 200 Pro microplate reader immediately following the addition of the substrate and several more times over 7–20 h. The well plates were incubated in the dark at ~20°C. MCA and MUF standard solutions prepared in ASW were used to determine hydrolysis rates. Killed controls (boiled sample water) and ultrapure water samples showed little change over the incubations.

Statistical Analysis

All statistical analyses were performed using SigmaPlot Version 12.5 (Systat Software Inc.). Results were not normally distributed thus non-parametric Spearman correlation coefficients (r) were computed to examine relationships between all variables. An α-value of ≤0.05 was used for all statistical analyses.

Results

Phytoplankton Growth

Exponential growth (T₁) was reached within 15–23 days in all experiments. Specific growth rates (μ) for each of the strains were 0.25 d⁻¹ for Skeletonema sp., 0.21 d⁻¹ for Leptocylindrus sp., 0.17 d⁻¹ for Phaeocystis sp., and 0.11 d⁻¹ for Coscinodiscus sp. Macroscopic particles formed within 3–8 days of the onset of exponential growth. Particle characteristics varied for each organism, from lightly aggregated particles for Skeletonema sp., to small well suspended particles for Leptocylindrus and Phaeocystis spp. Coscinodiscus sp. aggregates started out stringy and suspended becoming rounded, dense aggregates over time.

Fluorescence Characteristics

BEPOM samples for genera Skeletonema, Leptocylindrus, and Phaeocystis during exponential growth (T₁) were dominated by low-intensity fluorescence (<0.15 QSE) in the region of 275/340 (excitation/emission, peak T) with additional diffuse secondary emission occurring in the regions of 260/400–460 (peak A) and 320–360/420–460 (peak C) (Figures 1, 2, Table 1). BEPOM peaks A and C became more distinct and shifted to longer wavelengths over the duration of the experiment (Figure 3). Peak A red-shifted in emission wavelengths to 260/440–500 (peak A*) and peak C red-shifted in both excitation and emission to 350–400/450–500 (peak C*), resulting in a three-peak pattern with peak T over growth and degradation. For example, BEPOM peak A emission fluorescence in Skeletonema sp. occurred at 400–410 nm for exponential and stationary growth phases but shifted to longer wavelengths (~490 nm) after 6 weeks of degradation (Figure 3). Coscinodiscus sp. showed only faint secondary emission during both exponential growth samplings, but still followed the three-peak pattern observed in the other cultures. By stationary phase (T₂), Coscinodiscus sp. had minimal secondary emission in the region of peaks C and C*. Maximum fluorescence for all fluorescent peaks occurred during stationary phase (T₂) with a sharp decrease (ca. 70%) over the 6 weeks of degradation (T₃, Figures 1, 2).

FIGURE 1

Figure 1. Representative excitation-emission matrices (EEMs) of BEPOM fluorescence for genera Skeletonema, Leptocylindrus, and Phaeocystis at initial, exponential, stationary, and degradation phases and for Coscinodiscus sp. at initial, early exponential, late exponential, and stationary growth phases. Color bars are presented in quinine sulfate equivalence (QSE) where 1 QSE = 1 ppb quinine sulfate equivalents. Note different scales for color bars for each EEM. Major peak regions of EEMs are indicated in the initial time point for Skeletonema sp.

FIGURE 2

Figure 2. BEPOM fluorescence intensity of Coble EEM peaks over growth and degradation for genera Skeletonema (n = 3), Leptocylindrus (n = 2–4), Phaeocystis (n = 3), and Coscinodiscus (n = 3–4). Note the different growth stages for Coscinodiscus sp. All data are presented as mean ± standard error.

TABLE 1

Table 1. Central regions of EEM fluorescence attributed to different sources of organic matter, modified from Coble (2007) by Stedmon and Nelson (2015).

FIGURE 3

Figure 3. BEPOM (A–C) and CDOM (D–F) F_max-normalized emission spectra at 250 nm excitation for genera Skeletonema, Leptocylindrus, and Phaeocystis at exponential (dotted line), stationary (solid line), and degradation (dashed line) growth phases.

In contrast to the distinct fluorophores in the BEPOM samples, CDOM samples from exponential growth (T₁) had low (< 0.3 QSU), unstructured fluorescence with slightly enhanced fluorescence in the region of peak T and B (275/305), most commonly associated with protein-like fluorescence (Figures 4, 5, Table 1). Higher overall fluorescence intensities occurred during stationary (T₂) and degradation (T₃) phases and demonstrated a three-peak pattern similar to the BEPOM samples, though CDOM maximum fluorescence intensities were about half of BEPOM maximum fluorescence intensities. While CDOM shared the three-peak features observed in the BEPOM, the three peaks were much broader and less distinct in CDOM than in BEPOM and increased in fluorescence throughout growth and degradation for all areas of fluorescence (Figures 4, 5). Overall, CDOM fluorescence patterns were similar between all the phytoplankton species, with peaks A and T being the most prominent. CDOM peaks A and C were significantly correlated with BEPOM peaks A and C for many of the cultures, and all were significantly correlated when degradation (T₃) time points were excluded (Table 2).

FIGURE 4

Figure 4. Representative excitation-emission matrices (EEMs) of CDOM fluorescence for genera Skeletonema, Leptocylindrus, and Phaeocystis at initial, exponential, stationary, and degradation phases and for Coscinodiscus sp. at initial, early exponential, late exponential, and stationary growth phases. Color bars are presented in quinine sulfate equivalence (QSE) where 1 QSE = 1 ppb quinine sulfate equivalents. Note different scales for color bars for each EEM. Major peak regions of EEMs are indicated in the initial time point for Skeletonema sp.

FIGURE 5

Figure 5. CDOM fluorescence intensity of Coble EEM peaks over growth and degradation for genera Skeletonema (n = 3), Leptocylindrus (n = 2–4), Phaeocystis (n = 3), and Coscinodiscus (n = 3–4). Note the different growth stages for Coscinodiscus sp. All data are presented as mean ± standard error.

TABLE 2

Table 2. Spearman correlation coefficient (r) between select CDOM fluorescence peaks and select BEPOM fluorescence peaks (T, M, A, and C) for all data (n = 12) and excluding degradation phase data in parenthesis (n = 9).

The humification index (HIX) is a florescence index of the degree of organic matter degradation, with higher values characteristic of higher molecular weight, aromatic compounds (Huguet et al., 2009). CDOM HIX values increased over the duration of the experiment for all of the cultures, whereas BEPOM HIX values were more variable with no clear trend among the cultures (Figure 6).

FIGURE 6

Figure 6. Humification index (HIX) for BEPOM and CDOM for genera Skeletonema (Skel, n = 3), Leptocylindrus (Lepto, n = 2–4), Phaeocystis (Phaeo, n = 3), and Coscinodiscus (Cos, n = 3–4). Note the different growth stages for Coscinodiscus. All data are presented as mean ± standard error. Ratios from the initial time point were excluded due to low fluorescence at this stage.

Carbon Concentrations and POC/PON

Carbon concentrations were similar between all three carbon pools (particulate, base-extracted, and dissolved) for all four cultures (Figure 7, Table 3). In general all carbon concentrations increased through stationary phase and decreased over the 6 weeks of degradation. The exception to this was Phaeocystis sp. DOC concentrations that leveled off between stationary and degradation phases, and for Coscinodiscus sp. which did not include a degradation stage. The largest increase occurred in the particulate pool with POC increasing between initial (T₀) and stationary phase (T₂) by 1.56 mg C L⁻¹ for genera Skeletonema and 5.08 mg C L⁻¹ for Leptocylindrus. For base-extracted POC, Phaeocystis sp. had the smallest increase (0.33 mg C L⁻¹) compared to Leptocylindrus sp. (1.33 mg C L⁻¹) over the experiment. DOC concentrations had a similar range, with Coscinodiscus sp. increasing by 0.29 mg C L⁻¹ and Leptocylindrus sp. increasing by 1.14 mg C L⁻¹ between initial (T₀) and stationary phase (T₂).

FIGURE 7

Figure 7. Dissolved (DOC), base-extracted particulate (BEPOC), and particulate (POC) organic carbon concentrations for genera Skeletonema (Skel, n = 3), Leptocylindrus (Lepto, n = 2–4), Phaeocystis (Phaeo, n = 3), and Coscinodiscus (Cos, n = 3–4). Samples were collected at initial, exponential, stationary, and after 6 weeks of degradation, except for Coscinodiscus sp. which had time points collected at initial, early exponential, late exponential, and stationary phases of growth. All data are presented as mean ± standard error.

TABLE 3

Table 3. Bacteria abundance, hydrolytic enzyme activities (α-glucosidase, β-glucosidase, and aminopeptidase), dissolved (DOC), base-extracted (BEPOC), and particulate (POC) organic carbon concentrations, particulate organic carbon to nitrogen ratios (POC/PON), and BEPOM extraction efficiencies for the four phytoplankton cultures and associated controls.

Base-extraction efficiencies (BEPOC concentration divided by POC concentration times 100) varied between 4 and 88%, with Leptocylindrus sp. having a slightly narrower range (4–36%) (Table 3). These extraction efficiencies were consistent with those reported in Brym et al. (2014). Greater extraction efficiencies occurred earlier in growth (e.g., exponential phase 21–78%) and decreased with the duration of the experiment with the lowest extraction efficiencies occurring after 6 weeks of degradation (4–24%, T₃). Overall BEPOC concentration was correlated to POC concentration (r = 0.75, p < 0.001).

POC/PON ratios ranged from 5.2 to 10.9 for time points between exponential (T₁) and degradation growth phases (T₃) (Table 3). Higher POC/PON ratios occurred for Leptocylindrus sp. (6.9–9.6) and Coscinodiscus sp. (6.9–10.9), while Skeletonema sp. (6.7–7.2) and Phaeocystis sp. (5.2–6.8) had slightly lower POC/PON ratios; however, these differences were not significant.

Bacterial Abundance and Activities

Bacterial Cell Counts

The greatest increase in bacterial cell numbers occurred for genera Phaeocystis and Leptocylindrus, which had maximum bacterial cell numbers during stationary phase of growth (T₂), reaching 5.5 ± 2.4 × 10⁶ cells mL⁻¹ and 4.3 ± 1.0 × 10⁶ cells mL⁻¹, respectively (Figure 8, Table 3). By the end of the 6 weeks of degradation (T₃), bacterial cell numbers decreased by nearly half to 3.6 ± 0.7 × 10⁶ cells mL⁻¹ and 2.4 ± 0.5 × 10⁶ cells mL⁻¹, respectively. Bacterial cell numbers in the Skeletonema culture increased throughout the growth experiment, but only increased from 0.02 × 10⁶ at the initial sampling (T₀) to 1.0 × 10⁶ during degradation (T₃). The Coscinodiscus culture had the lowest overall bacterial cell numbers, only reaching 0.2 × 10⁶ cells mL⁻¹ during stationary growth (T₂).

FIGURE 8

Figure 8. Hydrolytic enzyme activities (α-glucosidase, β-glucosidase, and aminopeptidase) and bacterial abundance for genera Skeletonema (n = 3), Leptocylindrus (n = 2–4), Phaeocystis (n = 3), and Coscinodiscus (n = 3–4). Samples were collected at initial, exponential (Expo), stationary (Stat), and after 6 weeks of degradation (Deg), except for Coscinodiscus sp. which had time points collected at initial and early exponential (Early); no enzyme activity data (ND) were available for late exponential (Late) and stationary. The number of days following inoculation are given in parentheses. All data are presented as mean ± standard error.

Hydrolytic Enzyme Activities

Aminopeptidase enzyme activities were typically three to five orders of magnitude higher than either glucosidase activities, with the Phaeocystis culture having the highest activities during stationary phase (T₂, Figure 8, Table 3). Activities of β-glucosidase were typically a factor of two greater than α-glucosidase activities, except for Phaeocystis sp. which had β-glucosidase activities a factor of six to seven higher than α-glucosidase. Due to sampling constraints, no data was available for late-exponential (second T₁) or stationary (T₂) phases for Coscinodiscus spp.

Enzyme activities were highly correlated to fluorescent components, especially for glucosidase activities with genera Leptocylindrus and Phaeocystis and for aminopeptidase with Skeletonema sp. (Table 4). α- and β-glucosidase activities were also correlated to all carbon concentrations for genera Leptocylindrus and Phaeocystis.

TABLE 4

Table 4. Spearman correlation coefficient (r) between bacterial abundance or enzyme activity and DOC, BEPOC, POC and select BEPOM and CDOM fluorescence peaks (T, M, A, and C).

Discussion

Formation of CDOM by Marine Plankton in Culture

Previous laboratory experiments have relied on high nutrients (e.g., Rochelle-Newall and Fisher, 2002; Romera-Castillo et al., 2010) or additions of carbon substrates (e.g., Gruber et al., 2006; Goto et al., 2017) to observe CDOM formation by phytoplankton and/or bacteria. In our experiments, we produced CDOM fluorescence and phytoplankton densities on the same magnitude as open ocean waters using moderate nutrient additions and rotating bottles to keep cells in suspension. Using this approach, our maximum CDOM fluorescence for all cultures was low and ranged between 0.5 and 2 QSE (Figure 4), comparable to fluorescence intensities in several oceanic regions (0.5–6 QSE) (Nelson and Gauglitz, 2016; Netburn et al., in press). Additionally, trans-oceanic Pacific and Atlantic CDOM fluorescence ranged between 0.1 and 1 QSE for humic-like components and 0.1 to 11 QSE for protein-like components when data was modeled by parallel factor analysis (PARAFAC) (Murphy et al., 2008). Together, our results and previously published data, demonstrate that fluorescence intensities produced in culture experiments are similar to oceanic fluorescence intensities, even allowing for instrument variability. The low fluorescence intensities produced further suggest future culture experiments should use media with minimal background fluorescence to observe the low CDOM production by plankton.

Base-extraction of chromophoric and fluorophoric material from marine particles is a relatively new way to analyze POM and allowed for direct comparison between the two organic matter pools. The available field studies have been conducted in estuaries and coastal waters, but show remarkable consistency in BEPOM fluorescence. Overall fluorescence was low, typically less than 0.25 QSE, with maximum fluorescence corresponding to chlorophyll maxima and decreasing with depth (Brym et al., 2014; Ziervogel et al., 2016; Netburn et al., in press). For our cultures, maximum BEPOM fluorescence ranged between 0.4 and 3 QSE during stationary phase and decreased to 0.1–0.9 QSE during degradation (Figure 1). The elevated BEPOM fluorescence in cultures relative to natural samples was likely the result of measuring the cultures at the peak of phytoplankton biomass in a controlled environment compared to natural phytoplankton blooms that have more complex growth dynamics and influenced by other processes such as aggregation, grazing, and hydrodynamics.

Spectral Properties of Phytoplankton-Derived CDOM and BEPOM

All four monocultures produced a three-peak pattern that was dominated by discrete BEPOM fluorophores, suggestive of distinct compounds such as aromatic amino acids (Figure 1). This three-peak pattern has been described previously in natural BEPOM samples of estuarine (Brym et al., 2014) and oceanic (Ziervogel et al., 2016; Netburn et al., in press) origin. The prominence of the protein-like peak strongly resembles tryptophan, which has been shown to be dominant during phytoplankton growth (Stedmon and Markager, 2005; Murphy et al., 2008; Jørgensen et al., 2011). Distinct fluorescence in the region of peaks A and C have been attributed to fluorophores present in derivatives of benzoic acids (including phenols) and flavins, respectively (Wolfbeis, 1985; Wünsch et al., 2015). Other possibilities include siderophores (Fukuzaki et al., 2014). More generally, peaks A and C have previously been attributed to, and are abundant in, terrestrially-derived humic substances. Given our algal growth medium lacked material of terrestrial origin, as it was prepared with artificial seawater with f/20 nutrients, the observed formation and subsequent red-shift of peaks A and C were likely the result of microbially-transformed planktonic material rather than terrestrially-derived fluorescence.

CDOM fluorescence showed similar peak regions to BEPOM, but were generally broader and unstructured in emission and increased through growth and degradation phases (Figure 4). We interpret these findings to suggest that degradation processes occurred on the algal-derived material once in the dissolved phase. In previous studies, CDOM generated in phytoplankton and bacterial cultures showed similar fluorescence patterns to those presented here. Culture work on diatoms, dinoflagellates, and prasinophytes illustrated CDOM production by phytoplankton at intensities observed in the ocean (Romera-Castillo et al., 2010). Fukuzaki et al. (2014) presented EEM spectra measured on DOM from axenic diatoms, dinoflagellates, chlorophytes, cryptophytes, haptophytes, and raphidophytes monocultures. Though the dominant fluorophores differed between individual species, the ubiquitous fluorescence patterns produced overall by a range of genera further supports planktonic-sources of open ocean CDOM. However, these studies were only able to suggest that CDOM was formed contemporaneously with plankton growth as no fluorescence measurements on POM material were performed (Romera-Castillo et al., 2010; Fukuzaki et al., 2014). Additionally, these previous studies as well as our own still show elevated fluorescence in the regions of peaks A, C, and T, unlike open-ocean fluorescence which is low overall with maximum fluorescence in the area of peak A (Jørgensen et al., 2011). This suggests longer degradation experiments than performed in the present study and others (e.g., Gruber et al., 2006; Fukuzaki et al., 2014) are required to fully transform the protein and humic-like fluorescence into the low CDOM fluorescence observed in the ocean.

Results from our phytoplankton growth experiments, viewed through BEPOM to represent the particle phase, in addition to CDOM, strongly indicate that the unstructured CDOM fluorescence signal was a direct product of the microbial loop operating on planktonic biomass. Correlations between CDOM and BEPOM fluorescence peaks A and C (Table 2) were strongest when samples from degradation time points were excluded, due to the continual increase in CDOM fluorescence and concurrent decrease in BEPOM fluorescence (Figures 4, 5, Tables 1, 2). The overall increase in CDOM fluorescence and humification index (Figure 6), along with general decreases in all three carbon pools (Figure 7) and high microbial activities (Figure 8) between stationary and degradation time points, suggest microbial alteration of the DOM pool resulted in its compositional alteration to more fluorescent components. CDOM fluorescence increased in all peak regions throughout growth and degradation, including peaks A and C that have previously been described as humic-like substances arising from terrestrial material (Coble, 2007; Stedmon and Nelson, 2015). However, with no terrestrial material in our culture experiments, the humic-like fluorescence has to be a direct product of autochthonous production and microbial transformation. Similar processes occurring in the open-ocean water column would explain the global observations of elevated fluorescence and correlations with AOU in the deep ocean (Chen and Bada, 1992; Yamashita et al., 2007; Yamashita and Tanoue, 2008; Jørgensen et al., 2011; Catalá et al., 2015) that has been linked to the remineralization of organic matter and formation of RDOM (Yamashita et al., 2007; Jiao et al., 2010; Hansell, 2013; Lechtenfeld et al., 2015).

One advantage of incorporating BEPOM measurements into field sampling would be to observe subtle fluorescence shifts that are more distinct than in the dissolve phase and indicate microbial transformation of autochthonous material (Figure 3). A few previous studies have observed red-shifts in peaks A and C in CDOM, but these signals could easily be misinterpreted as terrestrial humic-like material or overshadowed by strong terrestrial influence, especially in coastal and freshwater systems. Murphy et al. (2008) hypothesized degradation of their coastal and open-ocean CDOM samples resulted in a PARAFAC component (C3) with excitation maxima at 260 and 370 nm and an emission maxima at 490 nm, similar to the locations of our red-shifted peaks A and C. Additionally, a similar PARAFAC component (C2, Ex 240 and 370 nm, Em 480) was determined from a meridional transect in the Atlantic Ocean by Kowalczuk et al. (2013), with its intensities doubling below the mixed layer. Burdige et al. (2004) also observed red-shifts in humic-like fluorescence of peaks A (239/429) and M (328/422) to longer wavelengths, resulting in peaks A' (248/461) and C (360/460) in sediment pore waters, and attributed the red-shift to diagenetic transformations of the original fluorophores. The red-shift in natural samples were all explained as a result of microbial degradation processes and most were at locations in the open ocean, far removed from terrigenous inputs of DOM. Thus, these studies support our interpretation that the appearance of red-shifted peaks in our plankton experiments were due to microbial degradation and that structurally complex substances were formed from planktonic precursors. The structural complexity may arise as a direct result of marine snow formation and microbial processing. Aggregation during particle formation could bring aromatic aldehydes and ketones in close proximity with hydroxylated benzoic acid derivatives through the microbial reprocessing into more complex molecules such as RDOM (Lechtenfeld et al., 2015). The presence of these molecules in the same moiety could facilitate charge transfer interactions between electron acceptors and donors, which has been suggested for lignin derivatives (Del Vecchio and Blough, 2004; Baluha et al., 2013; Sharpless and Blough, 2014). While the biochemical pathways and ecological ramifications of the planktonic CDOM phenomenon remain elusive, we can be certain that not all oceanic CDOM is terrestrially-derived and plankton are important sources of RDOM to the deep ocean.

Organic Matter Remineralization

Positive correlations between bacterial cell numbers and the different carbon pools (Table 4) suggest a close link between bacterial growth and activities in our experiments, with the microbes actively transforming organic matter compounds, such as carbohydrates, to shape the CDOM and BEPOM pools. Furthermore, decreases in BEPOM fluorescence (ca. 70%), BEPOM extraction efficiencies, and BEPOC concentrations during the 6 weeks of degradation suggest that substantial amounts of base-extractable POM components were degraded or solubilized to the dissolved phase by bacteria. This solubilized material was then rapidly removed from the dissolved phase, likely though organic matter remineralization to CO₂, resulting in net increases in DOC concentrations ranging between 0.19 and 0.46 mg C L⁻¹, regardless of the maximum carbon concentrations produced in each culture (Figure 7). This net DOC production was on the same order of magnitude as RDOC concentrations observed in the deep ocean (0.41–0.96 mg C L⁻¹) (Hansell et al., 2009).

Strong correlations between bulk glucosidase activities and organic matter fluorescence indicate that either carbohydrate hydrolysis products were a major source of CDOM likely resulting from POM sources, or that bacterial formation of CDOM is an energetically demanding process. Goto et al. (2017) demonstrated that a strain of Alteromonas grown solely on glucose produced humic-like CDOM fluorescence with emission properties similar to what we measured in our cultures. In contrast, amino acids and other nitrogen compounds apparently played a minor role in the formation of the pool of fluorescing organic matter as indicated by fewer correlations between aminopeptidase activities and organic matter parameters (Table 4). Rather, intermediate reactions, such as the reaction of aldehydes with organic amines may incorporate N into structures that can produce charge transfer, thus creating “humic-like” CDOM fluorescence on these purely planktonic precursors (Kieber et al., 1997; Brandes et al., 2004; Del Vecchio and Blough, 2004). Furthermore, the overall high rates of aminopeptidase (Figure 8, Table 3) indicate rapid turnover of nitrogen compounds and recycling of nutrients throughout the phytoplankton growth experiments. The coupling of BEPOM analysis of fluorescence along with measurements of enzyme activities has clearly opened a new dimension of observation of processes by which CDOM in the ocean may be formed by plankton.

Conclusions

Overall, our study demonstrates that phytoplankton-derived organic matter directly results in the formation of unstructured “humic-like” CDOM fluorescence, facilitated by aggregation and microbial processing of precursor POM biomolecules containing discrete, yet unidentified, fluorophores. We have demonstrated the ability to measure fluorescence from the cellular biomass of POM, extracted into dilute base, but without the need to use resin extraction techniques to concentrate the signal, which may alter the composition (and optical properties) of CDOM. The overall trends in growth and production of fluorescent CDOM follows global ocean observations in fluorescence patterns and intensity, as well as, net production of DOC. Furthermore, the significant correlations between BEPOM and CDOM fluorescence, especially peaks A and C, challenges the paradigm that humic-like fluorescence in the open ocean is the result of terrestrially-derived material (Andrew et al., 2013; Jørgensen et al., 2014). Finally, hydrolytic enzymatic rates suggest that the bacterial community actively transforms organic matter with significant turnover of carbon and nitrogen, making bacterially-mediated production of RDOM from planktonic biomass a key pathway in global element cycles (e.g., Jiao et al., 2010; Lechtenfeld et al., 2015; Moran et al., 2016).

Author Contributions

All authors contributed to the design of the study and data interpretation. JK and GC performed the incubation experiment. JK, GC, and KZ ran sample analyses. JK wrote the initial draft of the manuscript and all authors contributed to its revision.

Funding

Financial support for this work was provided by the National Science Foundation Chemical Oceanography award 1459406.

Conflict of Interest Statement

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.

Acknowledgments

The authors thank Daniel Wiltsie, Mackenzie Fiss, and Alexandra Gizzi for their assistance in sample collection and analyses.

References

Alldredge, A. L., and Silver, M. W. (1988). Characteristics, dynamics and significance of marine snow. Prog. Oceanogr. 20, 41–82. doi: 10.1016/0079-6611(88)90053-5