Turbulence simultaneously stimulates small- and large-scale CO2 sequestration by chain-forming diatoms in the sea.

Chain-forming diatoms are key CO2-fixing organisms in the ocean. Under turbulent conditions they form fast-sinking aggregates that are exported from the upper sunlit ocean to the ocean interior. A decade-old paradigm states that primary production in chain-forming diatoms is stimulated by turbulence. Yet, direct measurements of cell-specific primary production in individual field populations of chain-forming diatoms are poorly documented. Here we measured cell-specific carbon, nitrate and ammonium assimilation in two field populations of chain-forming diatoms (Skeletonema and Chaetoceros) at low-nutrient concentrations under still conditions and turbulent shear using secondary ion mass spectrometry combined with stable isotopic tracers and compared our data with those predicted by mass transfer theory. Turbulent shear significantly increases cell-specific C assimilation compared to still conditions in the cells/chains that also form fast-sinking, aggregates rich in carbon and ammonium. Thus, turbulence simultaneously stimulates small-scale biological CO2 assimilation and large-scale biogeochemical C and N cycles in the ocean.

Introduction

Diatoms contribute ca. 20% to the total primary production on Earth and play a key role in class="Chemical">CO2 sclass="Chemical">pan class="Chemical">equestration in the ocean[1]. Chain-forming genera of diatoms, e.g. Chaetoceros, Skeletonema and Thalassiosira, are dominant primary producers in environments with high turbulent shear, silicate, and nitrate concentrations, e.g., spring blooms in the polar oceans and temperate regions, and in subtropical and tropical upwelling regions. During these blooms, the formation of millimeter-sized, fast-sinking diatom aggregates in the upper sunlit ocean are instrumental for connecting the surface and the deeper ocean with respect to aggregate-attached biota, carbon (C) and nutrient sources, the aerobic and anaerobic nitrogen (N) cycles, and for biological CO2 sequestration to the mesopelagic ocean[2-7].

Over the past half century, several theories have been developed to explain the global success of chain-forming diatoms, and how chain formation may provide adaptations to physical, chemical and biological constraints in the sea[8-17]. Lapan class="Chemical">boratory exclass="Chemical">periments with diatom cultures have shown that chain formation can facilitate nutrient uclass="Chemical">ptake under turbulent shear[10-13]. However, the role of turbulent shear for cell-sclass="Chemical">pecific C-, and class="Chemical">pan class="Chemical">nitrate assimilation by chain-forming diatoms within mixed field communities, and the relation between nutrient demand and supply in co-existing individual diatom cells and cell chains remain poorly understood due to technical and methodological limitations.

Mixing of class="Chemical">water is inefficient at a micrometer scale where viscous forces dominate inertial forces. Thus, transclass="Chemical">port of gases and nutrients between class="Chemical">phytoclass="Chemical">plankton cells and ambient class="Chemical">pan class="Chemical">water occurs by diffusion through the diffusive boundary layer (DBL) in the water adjacent to the cells and cell chains (Fig. 1a). Diffusion-limited primary production occurs when cellular demand and potential uptake rate of gases or nutrients are higher than the maximum diffusive supply from the ambient water to the cell, which is constrained by the concentration gradients within the DBL at the cell–water interface (Eq. (1)). Small phytoplankton cells or cell colonies can theoretically cover their C and nutrient demands faster at low concentrations in ambient water, primarily due to shorter diffusion distances, lower cellular C and nutrient demands, and larger surface area:volume ratios compared to larger phytoplankton cells or colonies. The active nutrient uptake across the cell membrane keeps concentration low at the cell surface relative to that in the ambient water, and a large surface area:volume ratio of small cells allows for more ion carriers across the cell membrane relative to their cell volume compared to large cells. Diffusion-limited primary production occurs when the concentration gradients across the DBL limit the supply of gases or nutrients from the ambient water, and therefore occurs more likely in large cells or colonies of cells than in small, free-living cells when ambient nutrient concentrations are low[18-20].

Fig. 1

Diffusive boundary layers (DBL) surrounding chain-forming diatoms. The radial diffusion-limited concentration gradient across the DBL from the ambient water to cell surface for the middle cell under low turbulent shear (a) and under high turbulent shear (b)

Diffusive class="Chemical">boundary layers (class="Chemical">pan class="Chemical">DBL) surrounding chain-forming diatoms. The radial diffusion-limited concentration gradient across the DBL from the ambient water to cell surface for the middle cell under low turbulent shear (a) and under high turbulent shear (b)

The Kolmogorov length scale describes the size of the smallest eddies due to turbulence and is usually on the order of 1 to 10 mm in the surface ocean[21]. Below this scale fluid flow is laminar and the turbulent energy is dissipated through viscous shear, which enhances diffusion-limited gas and nutrient exchange in large phytoplankton (cells >60 μm in radius)[8,9]. Turbulent shear is theoretically more significant for mass transfer to large cells compared to small cells (Fig. 1b; pan class="Chemical">Eq. (2)) and may alleviate diffusion-limited fluxes of gases and nutrients in large cells or colonies[9-13,18-22].

Secondary ion mass spectrometry (SIMS) combined with stable isotopic tracers allows insight into life at the single-cell level of micrometer-sized individuals in mixed microbial field populations[23-27]. SIMS combines the qualities of a microscope with those of a mass spectrometer and reveals elemental and isotopic compositions of single cells of known identity and their individual biological activity. Using SIMS, we tested three central hypotheses in biological oceanography and biogeochemistry in two field populations of chain-forming diatoms. The first hypothesis is that turbulence enhances C assimilation and nutrient fluxes in chain-forming diatoms. The second hypothesis is that small cells grow faster than large cells under still or low, turbulent shear at low ambient nutrient concentrations. The third hypothesis is that aggregates of chain-forming diatoms mediate significant class="Chemical">ammonium fluxes to class="Chemical">pan class="Chemical">nitrogen-depleted, ambient water. To test these hypotheses, we analysed cell-specific C-fluxes and N-fluxes in two co-existing field populations of the chain-forming diatoms Skeletonema and Chaetoceros using SIMS combined with stable isotopic tracers and compared our data with those predicted by mass transfer theory. By combining sensitive fluorometry and mass transfer theory, we quantified ammonium concentrations within aggregates and fluxes from aggregates to the ambient water.

Our study showed that turbulent shear indeed increases C assimilation in chain-forming diatoms as predicted by mass transfer theory. However, small cells do not always grow faster than large cells when comparing field populations of different diatom genera. Finally, aggregates of chain-forming diatoms do mediate a substantial class="Chemical">ammonium flux to class="Chemical">pan class="Chemical">nitrogen-depleted water.

Results

Phytoplankton community composition

Experiments were performed with Baltic Sea spring bloom populations during April 2013. The fraction of large cells (>2 µm) in the phytoplankton community was dominated by the diatom genera Skeletonema and Chaetoceros (Supplementary Fig. 1). Average chains consisted of seven cells for Skeletonema (maximum: 22 cells) and four cells for Chaetoceros (maximum: 12 cells chain−1) (Table 1). Less abundant phytoplankton were the diatoms class="Species">Coscinodiscus sp., class="Chemical">pan class="Species">Melosira arctica, Thalassiosira sp. and Navicula sp., and the dinoflagellate Peridiniella catenata. Picoplankton (<2 µm) was presumably also present, but was not quantified here.

Table 1

Characteristics of Skeletonema and Chaetoceros

Genera	Cell length(µm)	Cell width(µm)	Cell volume(µm3)	Cellular POC cont.a(pmol C cell−1)	Cellular PON cont.a(pmol N cell−1)	Chain length(cells chain−1)	Abundance(cells L−1)
Skeletonema	15	5	294	2.4	0.34	7.0	2.023 × 106
Chaetoceros	10	15	1766	10.3	1.6	4.3	1.486 × 106

Cell dimensions, cellular POC and PON content, chain length and abundance

POC: particulate organic carbon, PON: particulate organic nitrogen

aCalculated from the cellular dimensions (Menden-Deuer and Lessard[28])

POC: particulate organic pan class="Chemical">carbon, PON: class="Chemical">particulate organic class="Chemical">pan class="Chemical">nitrogen

Community primary production and nitrogen assimilation

During sampling, the spring bloom was well-developed as reflected by Chl. a concentrations, which had a mean value of 13 µg L−1. Suspended POC and PON concentrations were 57.8 µmol C L−1 and 5.9 µmol N L−1, respectively (Table 2), and the concentration of class="Chemical">silicate was high (5.6 µM) relative to inorganic class="Chemical">pan class="Chemical">nitrogen and phosphate concentrations (<0.3 µM). The POC:PON (mol:mol) ratio of total organic matter and the ratio of C to nitrate plus ammonium assimilation were in the range of 8.5 to 10.6 indicating N-limitation of primary production by the phytoplankton community. Nitrate was the primary source of N-uptake within the mixed community as the nitrate assimilation rate was 6.4-fold higher than the ammonium assimilation rate. C- and nitrate-assimilation rates by the total phytoplankton community were statistically similar under still conditions and under turbulent shear (two-way ANOVA, p > 0.13, DF:4).

Table 2

Chemical parameters and assimilation rates

Nitrite + Nitrate (µmol L−1)	0.28 ± 0.36
Ammonium (µmol L−1)	0.176 ± 0.107
Phosphate (µmol L−1)	0.067 ± 0.005
Silicate (µmol L−1)	5.58 ± 0.36
Chl. a (µg L−1)	13.0 ± 4.0
Chl. a: POC (w:w)	0.018
POC (µmol L−1)	57.8 ± 5.8
PON (µmol L−1)	5.9 ± 0.5
C:N ratio (mol:mol)	9.8
C-assimilation rate (µmol C L−1 h−1) under still conditionsa	1.56 ± 0.34
C-assimilation rate (µmol C L−1 h−1) under turbulent sheara	1.21 ± 0.15
Nitrate-assimilation rate (µmol N L−1 h−1) under still conditionsa	0.147 ± 0.016
Nitrate-assimilation rate (µmol N L−1 h−1) under turbulent sheara	0.137 ± 0.039
Ammonium-assimilation rate (µmol N L−1 h−1) under still conditions	0.023 ± 0.001
C:Nitrate assimilation ratio under still conditions	10.6
C:Nitrate assimilation ratio under turbulent shear	8.8
C: (nitrate + ammonium) assimilation ratio under still conditions	8.5
Ammonium production rate (µmol N L−1 h−1) under still conditions	0.019 ± 0.001

Nutrient concentrations, Chl.a, POC, and PON in ambient water. C-, nitrate and ammonium assimilation rates, and ammonium production rates of the phytoplankton community. The numbers are average values ± standard deviation of the mean value (n = 3)

POC: particulate organic carbon, PON: particulate organic nitrogen

aValues of community C- and nitrate-assimilation rates were statistically similar under still conditions and under turbulent shear (ANOVA two-way, p > 0.05)

Nutrient concentrations, Chl.a, POC, and PON in ambient class="Chemical">water. C-, class="Chemical">pan class="Chemical">nitrate and ammonium assimilation rates, and ammonium production rates of the phytoplankton community. The numbers are average values ± standard deviation of the mean value (n = 3)

POC: particulate organic pan class="Chemical">carbon, PON: class="Chemical">particulate organic class="Chemical">pan class="Chemical">nitrogen

aValues of community C- and pan class="Chemical">nitrate-assimilation rates were statistically similar under still conditions and under turbulent shear (ANOVA two-way, class="Chemical">p > 0.05)

Cell-specific C and N assimilation

C and N assimilation in single cells within individual cell chains were quantified from class="Chemical">13C:class="Chemical">pan class="Chemical">12C and 15N:14N SIMS images (Fig. 2). The stable C- and N-isotope composition in individual cells revealed that C- and nitrate-assimilation rates in individual cells of both diatom populations varied by more than one order of magnitude independent of still or turbulent conditions (Fig. 3). Under still conditions, the average C:N (mol:mol) assimilation ratio and the standard error of the average value measured in individual cells was 4.5 ± 0.3 and 5.6 ± 0.9 in Skeletonema and Chaetoceros, respectively. Under turbulent shear, it was 5.9 ± 0.7 and 4.5 ± 0.4 in Skeletonema and Chaetoceros, respectively. Thus, it was slightly below Redfield ratio of 6.6 independent of still conditions or turbulent shear. Thirty-one percent more C was assimilated in Skeletonema cells under turbulent shear relative to still conditions. We used two-way ANOVA to test if cell-specific C- and N-assimilation rates in Chaetoceros or Skeletonema, were significantly different under still and turbulent conditions as proposed by our first hypothesis. In the Skeletonema population, we found that C-assimilation was significantly higher under turbulent shear relative to still conditions (two-way ANOVA, p < 0.0001, DF: 598), but nitrate assimilation was not significantly increased under turbulent shear compared to still conditions. Turbulent shear caused no statistically significant changes in the cell-specific C- and nitrate-assimilation rates in the overall Chaetoceros population (two-way ANOVA, p > 0.23, DF: 598). However, C assimilation in individual cells of long Chaetoceros cell chains (>7 cells chain−1) differed significantly between still conditions and turbulent shear (Fig. 4). Under still conditions, the average C assimilation was 45% lower in cells positioned in the middle of a Chaetoceros chain compared to that in apical cells. Under turbulent shear C assimilation did not change in apical cells compared to still conditions, but increased by 59% in cells positioned in the middle of the chains. Consequently, under turbulent shear C assimilation was evenly distributed among Chaetoceros cells within the same chain. N assimilation showed a higher variability among the cells and was not statistically different between still and turbulent conditions. In contrast, C and nitrate assimilation rates were uniformly distributed within Skeletonema cell chains under still and turbulent conditions. The lower C assimilation in Chaetoceros cells positioned in the middle of chains relative to that of end cells under still conditions as well as the increase in C assimilation in the middle of chains relative to that of end cells under turbulent shear were statistical significant (two-way ANOVA, p < 0.05; DF:57). Thus, turbulent shear increased C assimilation both within Skeletonema and Chaetoceros cell chains relative to that under still conditions as proposed by the first hypothesis.

Fig. 2

SIMS images of Chaetoceros and Skeletonema cell chains. Chaetoceros (a, c, e) and Skeletonema (b, d, f) incubated with 13C-HCO3− and 15N-NO3−. The biomass distributions are demonstrated by the 12C14N- images (a, b). The unit represent counts per pixel (a, b). SIMS images of 15N:14N ratio (c, d) and 13C:12C ratio (e, f) reflect nitrate assimilation and C assimilation, respectively. The scale bars represent 20 μm

Fig. 3

C- and nitrate-assimilation rates in individual cells. Each symbol represents the rates measured in individual cells of Skeletonema (open symbols) and of Chaetoceros (closed symbols) under still conditions (a) and under turbulent shear (b). The solid lines represent the Redfield ratio

Fig. 4

C- and N-assimilation rates depending on position within cell chains. Cell-specific C- and nitrate (N) assimilation rate measured in cells positioned at the end or in the middle of cell chains (>7 cells chain−1). Skeletonema (a, c) and Chaetoceros (b, d) under still (black bars) and turbulent conditions (open bars). The bars represent the mean value and the error bars one standard deviation (n = 7–28)

SIMS images of Chaetoceros and Skeletonema cell chains. Chaetoceros (a, c, e) and Skeletonema (b, d, f) incubated with class="Chemical">13C-HCO3− and class="Chemical">pan class="Chemical">15N-NO3−. The biomass distributions are demonstrated by the 12C14N- images (a, b). The unit represent counts per pixel (a, b). SIMS images of 15N:14N ratio (c, d) and 13C:12C ratio (e, f) reflect nitrate assimilation and C assimilation, respectively. The scale bars represent 20 μm

C- and class="Chemical">nitrate-assimilation rates in individual cells. Each symclass="Chemical">pan class="Chemical">bol represents the rates measured in individual cells of Skeletonema (open symbols) and of Chaetoceros (closed symbols) under still conditions (a) and under turbulent shear (b). The solid lines represent the Redfield ratio

C- and N-assimilation rates depending on position within cell chains. Cell-specific C- and pan class="Chemical">nitrate (N) assimilation rate measured in cells class="Chemical">positioned at the end or in the middle of cell chains (>7 cells chain−1). Skeletonema (a, c) and Chaetoceros (b, d) under still (black bars) and turbulent conditions (oclass="Chemical">pen bars). The bars reclass="Chemical">present the mean value and the error bars one standard deviation (n = 7–28)

Mass transfer theory

The calculated, diffusion-limited DIC supply with respect to class="Chemical">bicarbonate to Chaetoceros was 167,000 fmol C cell−1 h−1 and more than 1000-times higher than the actual net C-assimilation rate of 67 to 120 fmol C cell−1 h−1, declass="Chemical">pending on the cell class="Chemical">position in the chain (class="Chemical">pan class="Chemical">Eq. (1); Fig. 4). Under still conditions, the diffusion-limited supply of nitrate from the ambient water to Chaetoceros cell chains was 38.4 fmol N cell−1 h−1, which was similar to the average measured nitrate assimilation rates of 37.6 fmol N cell−1 h−1 (Fig. 4).

The calculated diffusion-limited DIC supply to Skeletonema was 114,000 fmol C cell−1 h−1 (class="Chemical">Eq. (1)), but the measured net C assimilation was only 15 to 20 fmol C cell−1 h−1, i.e. more than 6000 times lower than the class="Chemical">potential DIC flux (Fig. 4). Measured class="Chemical">pan class="Chemical">nitrate assimilation rates (5.5 fmol C cell−1 h−1) were 20% of the diffusion-limited supply of nitrate from the ambient water under still conditions (26 fmol N cell−1 h−1), indicating that primary production and cellular growth in Skeletonema was N-limited, but possibly also limited by another element, e.g. phosphorus (P).

Assuming Redfield stoichiometry (N:P assimilation ratio of 16), the P assimilation rate was 2.4 fmol P cell−1 h−1 and 0.3 fmol P cell−1 h−1 for Chaetoceros and Skeletonema, respectively. By comparison, the diffusion-limited P-supply was 3.6 fmol P cell−1 h−1 and 2.5 fmol P cell−1 h−1 for Chaetoceros and Skeletonema, respectively. Hence, the estimated P assimilation rate was 67% of the diffusion-limited P-supply in Chaetoceros, while it was 14% of that in Skeletonema.

The Sherwood number (Sh) describes the relative increase in total flux of gases or nutrients due to turbulent shear compared to still conditions as predicted by mass transfer theory[9,22]. The predicted increase in total fluxes of class="Chemical">nitrate and class="Chemical">pan class="Chemical">ortho-phosphate due to turbulence relative to still conditions was 38% and 62%, respectively, in Chaetoceros, whereas they were 25% and 40%, respectively, in Skeletonema (cf. Eq. (2)). Our empirical observations of <59% and 31% in Chaetoceros and Skeletonema, respectively, were thus within the expected theoretical range for N and P coupled to the net C-assimilation rates in the diatoms (Discussion section).

Cell-specific contributions to community C and N assimilation

The average biovolume of single Chaetoceros cells was sixfold larger than that of Skeletonema cells (Table 1), and the average cell-specific C- and N-assimilation rates were sixfold to sevenfold higher in Chaetoceros compared to Skeletonema. The C-assimilation rate was similar during our class="Chemical">ammonium and class="Chemical">pan class="Chemical">nitrate assimilation experiments in both Chaetoceros and Skeletonema under still conditions. In both diatom genera, the cell-specific ammonium assimilation rate was ca. 16% of the nitrate assimilation rate (Fig. 5a). Under still conditions, the Chaetoceros population and the Skeletonema populations contributed 38% and 8% to the total nitrate assimilation by the phytoplankton community, respectively (Table 3). For ammonium and C assimilation, their relative contributions were 35% and 6%, and 10% and 2%, respectively. Under turbulent shear, their relative contribution to C assimilation was 15% and 4%, respectively, whereas their relative contribution to nitrate assimilation was 32% and 6%, respectively.

Fig. 5

C- and N-assimilation rates in the Skeletonema population and the Chaetoceros population. Cell-specific C-assimilation rates (a); Cell-specific N-assimilation rates (b); C-specific C-assimilation rates (c) and N-specific N-assimilation rates (d) in Chaetoceros (black bars) and Skeletonema (white bars) measured for nitrate as N sources and under still conditions or turbulent shear and at light intensities >200 μmol photons m−2 s−1. Ammonium assimilation rates were measured without turbulent shear and represent the average value in light and darkness. All values represent the average with one standard deviation (n = 70–262). Statistically significant differences are described in the text

Table 3

C and N assimilation associated to Skeletonema and Chaetoceros

Diatom	POC%	PON%	C assimilation%	NO3 assimilation%	NH4+ assimilation%
Skeletonema (still condition)	9	12	2a	8	6
Skeletonema (turbulent shear)			4a	6	nd
Chaetoceros (still conditions)	27	38	10	38	35
Chaetoceros (turbulent shear)			15	32	nd

POC: Particulate organic carbon, PON: particulate organic nitrogen

POC, PON, C-, nitrate assimilation and ammonium assimilation in % of the total phytoplankton community under still conditions and turbulent shear

aDifferences in cell-specific rates were statistically significant between still conditions and turbulent shear (two-way ANOVA, DF: 598; p < 0.0001)

C- and N-assimilation rates in the Skeletonema population and the Chaetoceros population. Cell-specific C-assimilation rates (a); Cell-specific N-assimilation rates (b); C-specific C-assimilation rates (c) and N-specific N-assimilation rates (d) in Chaetoceros (black bars) and Skeletonema (white bars) measured for pan class="Chemical">nitrate as N sources and under still conditions or turbulent shear and at light intensities >200 μmol class="Chemical">photons m−2 s−1. class="Chemical">pan class="Chemical">Ammonium assimilation rates were measured without turbulent shear and represent the average value in light and darkness. All values represent the average with one standard deviation (n = 70–262). Statistically significant differences are described in the text

POC: Particulate organic pan class="Chemical">carbon, PON: class="Chemical">particulate organic class="Chemical">pan class="Chemical">nitrogen

POC, PON, C-, pan class="Chemical">nitrate assimilation and class="Chemical">pan class="Chemical">ammonium assimilation in % of the total phytoplankton community under still conditions and turbulent shear

The sixfold difference in biovolume between Chaetoceros and Skeletonema can partly be explained by silicified cell walls as the C and N contents of Chaetoceros cells are approximately only 4-fold larger than those of Skeletonema cells[28]. C- and N-assimilation rates and doubling times can be expressed independent of cell size when C- and N-assimilation rates are normalised to cellular C and N contents as measured by SIMS. The average C-specific C-assimilation rate (C-growth rate) in Chaetoceros was 40% higher than in Skeletonema (cf. class="Chemical">Eq. (3); Fig. 5c, d), and the average N-sclass="Chemical">pecific N-assimilation rates (N-growth rate) in Chaetoceros, with resclass="Chemical">pect to class="Chemical">pan class="Chemical">nitrate and ammonium, were 67% and 97% higher, respectively, relative to those in Skeletonema. C-growth rates corresponded to average doubling times of 66 and 92 light hours under still conditions, and 62 and 71 light hours under turbulent shear in Chaetoceros and Skeletonema, respectively. N doubling times were 40 and 66 h under still conditions, and 37 and 56 h under turbulent shear in Chaetoceros and Skeletonema, respectively. We used two-way ANOVA to test if C and N doubling times were statistically significantly different in Chaetoceros vs. Skeletonema as proposed by our second hypothesis. Contrary to mass transfer theory, the C and N doubling times were significantly shorter in the larger Chaetoceros cells compared to the smaller Skeletonema cells both under still conditions and turbulent shear (p < 0.001; DF: 598).

NH4+ production in diatom aggregates and ambient water

Fast-sinking diatom aggregates are instrumental in class="Chemical">CO2 sclass="Chemical">pan class="Chemical">equestration in the ocean. Formation of mm-sized diatom aggregates from natural seawater in the laboratory also created ammonium-rich microenvironments that released ammonium to the ambient water (Fig. 6). NH4+ concentrations within aggregates were in the mid micromolar range, as modeled from the measured NH4+ release to the ambient water, and were <100-fold higher than in ambient water (176 nM) (Table 1). On average three mm-sized aggregates L−1 formed, and the measured average net flux of NH4+ from individual aggregate to the ambient water was 1.5 nmol NH4+ agg−1 h−1, or 4.5 nmol NH4+ L−1 h−1. The ammonium production rate measured in these aggregates (nmol N L−1 h−1) was equivalent to 24% of the ammonium production in the ambient water without aggregates (19 nmol N L−1 h−1), and equivalent to 6% of the nitrate assimilation rate (147 nmol N L−1 h−1) by the total phytoplankton community. Thus, diatom aggregates provided a substantial flux of ammonium to the ambient water as proposed in our third hypothesis. In ambient water, the ammonium production rate was similar to the ammonium assimilation rate showing that the plankton community efficiently assimilated the ammonium produced (Table 2).

Fig. 6

The radial distribution of ammonium concentrations in sinking diatom aggregates. The concentrations were modeled from the measured ammonium release to the ambient water. An example of such an aggregate is shown in the picture. The white size bar is 3 mm

The radial distribution of class="Chemical">ammonium concentrations in sinking diatom aggregates. The concentrations were modeled from the measured class="Chemical">pan class="Chemical">ammonium release to the ambient water. An example of such an aggregate is shown in the picture. The white size bar is 3 mm

Discussion

Our knowledge and understanding of physical-biological coupling and eco-physiology at the single-cell level within mixed field populations and communities have been largely limited due to technical constraints. The recent introduction of SIMS in biological oceanography has revealed a high variability of physiological rates at the single-cell level across field populations of various (micro)-organisms[23-27]. Our study demonstrated that C and class="Chemical">nitrate assimilation can vary by more than one order of magnitude among single cells in diatom class="Chemical">poclass="Chemical">pulations. However, average values of the assimilation rates stabilised at a reclass="Chemical">presentative mean value of each class="Chemical">poclass="Chemical">pulation when >50 cells were analysed (Fig. 7). class="Chemical">pan class="Chemical">Nitrate assimilation rates co-varied with C-assimilation rates in the majority of cells (Fig. 3). A similar large variability and co-variation of C-fixation and N2-fixation or NH4+ assimilation has been observed in field populations of cyanobacteria[24,27] and other phototrophic bacteria[23]. Intra-species variations in eco-physiology and nutrient demand at the cellular level provide plasticity for a population to adapt to regional heterogeneity and environmental constraints on a short time scale, as recently demonstrated at a single-cell level in cultures of heterotrophic N2-fixing bacterial populations[29]. Genetic diversity within a population composed of various strains is ultimately also a prerequisite for evolution on a long time scale.

Fig. 7

Average C- and N-assimilation rates as a function of cells analysed. Average cell-specific C-assimilation rate (closed symbols) or nitrate assimilation rate (open symbols) with the standard error of the mean as a function of cells analysed in Skeletonema (a) and in Chaetoceros (b) under turbulent shear

Average C- and N-assimilation rates as a function of cells analysed. Average cell-specific C-assimilation rate (closed symclass="Chemical">bols) or class="Chemical">pan class="Chemical">nitrate assimilation rate (open symbols) with the standard error of the mean as a function of cells analysed in Skeletonema (a) and in Chaetoceros (b) under turbulent shear

The combination of SIMS analysis with stable isotope tracers, and sensitive fluorometry enabled us to directly examine three basic hypotheses of cell-specific C and N assimilation by chain-forming diatoms in natural mixed communities at low-nutrient concentrations in the sea. Net C assimilation in class="Chemical">both diatom sclass="Chemical">pecies increased with turbulent shear as exclass="Chemical">pected from mass transfer theory, but surclass="Chemical">prisingly, no concurrent significant increase of class="Chemical">pan class="Chemical">nitrate assimilation was detected in the same cells as a response to turbulent shear. Contrary to predictions by mass transfer theory, large Chaetoceros cells assimilated N 70–94% faster than the smaller Skeletonema. Contributions of Skeletonema to C and N assimilation were modest, although its cell abundance was close to its maximum during this spring bloom[30]. Chl. a, suspended POC, and PON concentrations were typical of those occurring at the end of diatom spring blooms when nutrients are depleted in temperate coastal zones[5]. The larger, but less abundant Chaetoceros grew significantly faster than Skeletonema and contributed substantially to total C, nitrate, and ammonium assimilation by the entire phytoplankton community. Similar observations using nanoSIMS were recently reported in phototrophic microbial populations, where the largest, but less abundant phototrophic species contributed most to C and N assimilation in field communities[23,27]. Hence, smaller organisms do not always grow faster than larger ones do, even when nutrient concentrations are low in the field, and less abundant genera may dominate the total community production.

In our experiments, a relatively high shear rate was applied to enable significant C and N assimilation changes at the single-cell level where variability is large. The applied shear rate (20 s−1) corresponds to that close to the wave-breaking zone of the upper ocean[31]. Under these conditions, turbulent shear increased the net C-assimilation rate up to 59% within the longest diatom chains of Chaetoceros and by 31% in Skeletonema at low-nutrient concentrations compared to still conditions. The fluxes of nutrients with lowest molecular diffusion coefficients are generally most sensitive to turbulent shear, e.g. inorganic P fluxes are more sensitive to turbulent shear than those of inorganic N (pan class="Chemical">Eq. (2))[9,19]. At a shear rate of 1 s−1, often considered to reclass="Chemical">present the average of that occurring within the euclass="Chemical">photic zone, the theoretical enhancements in net C-assimilation rates due to turbulence at P-limitation are 14 and 9%, while they are 9 % and 6 % at N-limitation in Chaetoceros and Skeletonema, resclass="Chemical">pectively. Hence, the class="Chemical">predicted effect of turbulent shear may also be detectable below the wave-breaking zone, and it becomes substantial when integrated over a few days. Skeletonema cell chains are long, thin and stiff while those of Chaetoceros are shorter, but thicker and more flexible. Thus, the different resclass="Chemical">ponses of cell-sclass="Chemical">pecific C-assimilation rates to turbulent shear within chains of Chaetoceros and Skeletonema may be exclass="Chemical">plained by their different dimensions and morclass="Chemical">phology[32].

DIC fixation in Skeletonema and Chaetoceros is largely supplied by class="Chemical">bicarbonate uclass="Chemical">ptake (85% of total DIC fixation) during class="Chemical">photosynthesis[33,34]. In the ocean, class="Chemical">pan class="Chemical">bicarbonate concentrations are generally in the millimolar range whereas nutrient concentrations are in the nano- to micromolar range. Diffusion-limited fluxes of DIC and nutrients to cells are proportional to their respective concentrations in the ambient water (Eq. (1)), and the modeled DIC fluxes were three orders of magnitude higher than the measured net cellular C-assimilation rates. Hence, DIC fluxes do not limit photosynthesis in many diatoms[33,34], although gradients in pH, CO2 and bicarbonate concentrations occur at the cell–water interface[35]. Gross DIC fixation often exceeds C demands relative to nutrient uptake and cellular growth under nutrient limitation in chain-forming diatoms. A fraction of newly fixed carbon is excreted as dissolved organic carbon (DOC) and transparent exopolymer particles (TEP) while net C assimilation in biomass reflects nutrient-limited growth[36]. The average nitrate assimilation rate in Chaetoceros was similar to the upper limit determined by diffusion-limited supply through the DBL, but any short-term response in nitrate-assimilation rates to turbulent shear was masked by a high variability at the single-cell level. This may reflect storage and release of nitrate from the vacuole within the cells[37] as well as concurrent ammonium uptake. A (co)-limitation by P may also have occurred as indicated by the modeled P uptake. Nutrient limitation and uptake were thus more complex than described by diffusion limitation of one single nutrient.

Diffusion-limited nutrient uptake by individual phytoplankton, including chain-forming diatoms and other colony-forming phytoplankton, has previously been studied only using cultures and/or theoretical approaches[8-13]. In cultures, measured maximum cell-specific N-uptake rates were 53 fmol N cell−1 h−1 and 29 fmol N cell−1 h−1 in Chaetoceros calcitrans and class="Species">Skeletonema sp. resclass="Chemical">pectively, and a half-saturation constant of 0.4 µM for class="Chemical">pan class="Chemical">nitrate uptake has been reported in Skeletonema cultures[10]. In the present study, the average net nitrate assimilation rate within the field population of Chaetoceros cells was 37 fmol N cell−1 h−1 and close to the theoretical limit constrained by diffusion limitation at a concentration of 300 nM nitrate in the ambient water (Eq. (2)). By contrast, it was 5.5 fmol N cell−1 h−1 and ca. 5 times lower than that predicted under diffusion limitation in Skeletonema cells. Thus, the physiological state of Skeletonema may have been in the transition to stationary growth phase in the field population[30]. Carbon and nitrate, however, were assimilated close to or slightly below the Redfield ratio in both Skeletonema and Chaetoceros under still conditions, and some of the nitrate was presumably stored in the vacuole[37]. Cell-specific ammonium assimilation rates in non-aggregated Skeletonema and Chaetoceros in the ambient water were only 16% of the nitrate assimilation rates during daytime. In Chaetoceros, the rate increased by 39% during night compared to that measured during day (Students t-test: p < 0.01; n = 138, 129; Supplementary Table 1), while it increased by 86% during night in Skeletonema (Students t-test: p < 0.01; n = 93, 44; Supplementary Table 1). Thus, nitrate assimilation dominated over ammonium assimilation in these diatoms, which did not appear to be N-limited judged by their C:N-assimilation ratios. In contrast, the ratio of C to nitrate plus ammonium assimilation in the total phytoplankton community was 8 to 10 and thus well above the Redfield ratio. Hence, N-limitation of community primary production was apparently associated with other organisms that presumably lacked nitrate reductases and were limited by available ammonium[38].

Despite its high C and N growth rates, the cell abundance of class="Species">Chaetoceros sp. was lower than that of Skeletonema. Large cell size and long sclass="Chemical">pines facilitate the formation of fast-sinking aggregates in Chaetoceros as confirmed by our laclass="Chemical">pan class="Chemical">boratory experiments. In contrast, Skeletonema is well known to require relatively high cell concentrations for aggregate formation compared to other diatoms[5,36]. Turbulent shear stimulates aggregate formation in chain-forming diatoms, and small aggregates can form within a few hours while milimeter-sized aggregates may need 12 to 24 h to form depending on cell concentration, cell size and chain length, as was also the case in the present study[4,5,36,39]. Silicified frustules in diatoms act as biogenic ballast minerals, and sinking velocities of diatom aggregates observed in situ are on the order of 100 to 200 m d−1 or even faster in the laboratory[2,40]. Such aggregates are readily colonised by bacteria and protozoa that partly respire and remineralize the organic matter during their export out of the euphotic zone, linking processes in the upper ocean with those in the mesopelagic ocean[6,7,41,42]. Ammonium production in aggregates of the present study was likely driven by bacterial hydrolysis and remineralization of dead diatoms and other organic matter in the aggregates[3,6,7,41]. The half-life of organic C in aggregates during microbial hydrolysis and respiration is usually one to two weeks in the upper ocean at 15 °C but increases ca. 3 fold with decreasing temperature down to the 4 °C prevailing in the mesopelagic ocean[40-42]. Sinking aggregates are therefore important vehicles of organic C and nutrient transport in the ocean.

Our study demonstrates that chain-forming diatoms link small-scale class="Chemical">CO2 and class="Chemical">pan class="Chemical">nitrate assimilation as well as ammonium production in the euphotic zone to large-scale biogeochemical fluxes of C and N in the ocean interior through the formation of aggregates. Field populations of chain-forming diatoms had a high capacity to grow and assimilate nitrate at ambient concentrations of <0.3 µM. Turbulent shear increased C assimilation and growth in the longest and largest cell chains, which were relatively low in abundance. Larger cells and cell chains, however, are more efficiently exported out of the euphotic zone compared to smaller cells and cell chains because larger chains coagulate to form aggregates at lower cell abundance than smaller ones do. Thus, turbulence is instrumental for CO2 sequestration by large, chain-forming diatoms through simultaneous stimulation of cellular growth and formation of fast-sinking diatom aggregates, which drive biological CO2 sequestration in the ocean[42,43].

Methods

Study area

All class="Chemical">water samclass="Chemical">ples were collected at the sea surface at monitoring station B1 (N 58° 48’ 18, E 17° 37’ 52) of the Swedish National Monitoring Program (Swedish Meteorological and Hydrological Institute, SMHI). At the time of samclass="Chemical">pling salinity was 6 and the class="Chemical">pan class="Chemical">water temperature was 3 °C. Water analyses and experiments were performed at Askö laboratory.

Chlorophyll a analysis

Three replicates each of 500 mL of seaclass="Chemical">water were filtered onto GF/F filters. class="Chemical">pan class="Chemical">Chlorophyll a was extracted in 90% methanol at 4 °C during 6 h, and analysed on a spectrophotometer (Pharmacia KLB Ultrospec III)[44].

Nutrient analysis

Sea class="Chemical">water was filtered through 0.2 µm syringe filters into 5 mL class="Chemical">plastic vials and frozen until analysis. class="Chemical">pan class="Chemical">Silicic acid, nitrate, ammonium and phosphate concentrations were quantified with a continuous-flow automated system (TRAACS 800, Bran-Luebbe, Germany)[45]. Ammonium concentrations were analysed immediately after sampling (see below).

Microscopy

Lugol-preserved 150 mL seapan class="Chemical">water samclass="Chemical">ples for cell counting were examined under an inverted microscoclass="Chemical">pe (Zeiss, Axio) at ×200 magnification after the cells had settled in a 25 mL Utermöhl sedimentation chamber (Hydrobios, Germany). Number of cells, cell size and chain length of the dominant sclass="Chemical">pecies (Skeletonema and Chaetoceros) were quantified. Cell volumes were calculated as cylinders: V = πr2h.

The potential diffusive DIC and nutrient supply to cells within cell chains were calculated from the analytical solutions of diffusion to a cylinder under still conditions[46]:where Q is the quantity of substance diffusing to a cell within the chain with the length (L) and diameter (D) per unit time, t, D is the diffusion coefficient of the substance, C0 its concentration at the cell surface, r0 and C∞ the concentration in the ambient class="Chemical">water. The class="Chemical">pan class="Chemical">equation is considered accurate for L/D < 8. We assumed C0 to be zero, and used the molecular diffusion coefficients for HCO3− (6.0 × 10−6 cm2 s−1), NO3− (9.7 × 10−6 cm2 s−1) and ortho-phosphate (3.6 × 10−6 cm2 s−1) at 3 °C and a salinity of 6 as previously reported[47].

The Sherwood number, Sh, describes the relative increase in transport-limited fluxes due to fluid motion in the vicinity of a cell as compared to that at still conditions[22]. It depends on the Peclet number, Pe, and can be calculated as a function of shear rate due to turbulence for 0.1 < Pe < 100[9]:where E is the shear rate (s−1) and (Er02/D) is Pe.

In our experiments, we used couettes to create laminar shear[36] with a shear rate of 21 s−1 in order to measure statistically significant differences in nutrient uptake under still and under turbulent conditions (please see below). Pe numbers were 1 and 3 for Skeletonema and Chaetoceros, respectively, when calculations were based on the class="Chemical">equivalent sclass="Chemical">pherical radius of the colony volume. The average class="Chemical">pan class="Chemical">equivalent spherical radius of cell chains was 7.9 µm and 12.2 µm for Skeletonema and Chaetoceros, respectively.

Incubations with stable isotope tracers

The phytoplankton community was incubated with stable isotope tracers (class="Chemical">H13CO3−combined withclass="Chemical">pan class="Chemical">15NO3− or 15NH4+). The final excess 13C:12C isotope ratio of H13CO3− in the ambient water was 3.6% combined with 5.6% excess 15N:14N isotope ratio of NO3− or 80% excess 15N:14N isotope ratio of NH4+, respectively. Final concentrations of NO3− and NH4+ were similar (0.3 µM). Bottles were carefully shaken after addition of the stable isotope tracer. For the quantification of NH4+assimilation and NH4+ production, nine 1 L bottles were incubated by dual labelling with H13CO3− and 15NH4+ and one control without tracers in natural light in an outdoor mesocosm at in situ temperature at noon. The control bottle was used to measure the natural isotope ratios of 13C:12C and 15N:14N in the diatoms cells as well as in the phytoplankton community. An additional set of bottles were incubated around midnight. Triplicate incubation bottles were stopped by filtration at T0 and after 2 and 5 h. For the quantification of nitrate assimilation, the phytoplankton community was incubated with H13CO3− and 15NO3− in 1 L Duran bottles under still conditions or under laminar shear in couettes[36] in the laboratory. Three Duran bottles and three couettes with stable isotope tracers were incubated for 8.5 h in a thermo-stated room at in situ temperature and a light intensity of 220 µmol m−2s−1. Three 1 L Duran bottles and one couette without stable isotope tracers served as controls. All Duran bottles were placed horizontally next to the couettes, which were operated with a speed of 0.9 revolutions s−1.

Samples were fixed at 4°C for 24 h with 2% class="Chemical">paraformaldehyde (class="Chemical">pan class="Chemical">PFA) for SIMS analysis of C and N assimilation by Skeletonema sp. and Chaetoceros sp. and washed after filtration onto GTTP filters (pore size 0.22 µm; diameter 25 mm, Millipore, Eschborn, Germany). Furthermore, 500 mL samples were filtered onto pre-combusted GF/F filters for elemental analysis-isotope ratio mass spectrometry (EA-IRMS) analysis of C and N assimilation by the whole phytoplankton community (see below). Twelve millilitre samples were filled, headspace free into gas tight Exetainer©vials to which 100 µL 50% ZnCl2 solution was added for later analysis of 13C/12C-HCO3−, 15N/14N-NO3− and 15N/14N-NH4+. 13C-12C isotope ratios of bicarbonate were analysed at Stable Isotope Facilities, University of California Davis, USA, and the 15N/14N-NO3−, and 15N/14N-NH4+ isotope ratios were analysed at the Stable Isotope Laboratory Stockholm University (see below).

Elemental analysis-isotope ratio mass spectrometry

GF/F filters were freeze-dried and decalcified overnight in fuming 37% pan class="Chemical">HCl in a desiccator class="Chemical">prior to analysis by a Thermo Flash EA 1112 elemental analyzer couclass="Chemical">pled to an isotoclass="Chemical">pe ratio mass sclass="Chemical">pectrometer (Thermo Delta Plus XP, Thermo Fisher Scientific). class="Chemical">pan class="Chemical">Caffeine was used as standard for isotope correction and C/N quantification of bulk C and N assimilation and bulk C and N assimilations rates were calculated as previously described[48]. EA-IRMS analysis was completed before SIMS analysis was initiated.

Secondary ion mass spectrometry analysis

We used a large-geometry SIMS instrument (Cameca IMS 1280; Gennevilliers, France) at the NORDSIM facility, Swedish Museum of Natural History. The IMS1280 is a dynamic SIMS, with magnetic sector mass separation, and imaging capability similar to nanoSIMS. However, the spatial resolution of IMS1280 is ca. 1 µm whereas that of nanoSIMS is ca. 50 nm. The lower spatial resolution of IMS1280 allows for higher through-put compared to nanoSIMS when cells are large (>3 µm). Prior to analysis, the filters containing chemically fixed cells were coated with a nanometer thin gold layer and cut into ca. 4×4 mm2 pieces and mounted on a sample holder. Analysis was performed using a class="Chemical">Cs+ class="Chemical">primary beam with a sclass="Chemical">patial resolution of ca. 1 µm. We only analysed disclass="Chemical">persed cell chains on the filters. Diatom cell chains were class="Chemical">pre-sclass="Chemical">puttered with a beam of 10 nA for 5 min to remove the silicified cell wall, and then imaged using a 40–80 class="Chemical">pA class="Chemical">primary beam for 100 cycles. The class="Chemical">pre-sclass="Chemical">puttered area was larger than the imaged area in order to eliminate class="Chemical">possible slight offsets between the two beams (sclass="Chemical">putter and analytical) and edge effects. Analyses were automated into cells after a fixed class="Chemical">period of class="Chemical">pre-sclass="Chemical">puttering needed to reach the interior of cells with high CN had been determined. For each cell chain we recorded secondary ion (SIMS) images of class="Chemical">pan class="Chemical">13C14N− and 12C14N−, and 12C15N− using a peak-switching routine at a mass resolution of ca. 6000 (M/ΔM), which is less than the ~12,000 M/ΔM necessary to resolve BO− from 13C14N−47. BO− is a common surface contaminant and present as a trace element in the diatom frustules[49,50]. Also, ~7000 M/ΔM is necessary to fully resolve 12C14N- from 13C13C−, which can be significant for highly 13C-enriched samples. However, measurements on control samples suggest that these interferences did not significantly affect the isotope measurements inside the cells, but BO− slightly affected ratios in low yield regions, e.g. the frustules. The use of CN isotopomers have the advantage that they are bright and two isotope ratios can be measured from measurements of three isotopomers. The imaged area was 80×80 µm with a scanning ion imaging of 256×256 pixels. The dwell times were 1 s, 5 s and 2 s for 12C14N, 12C15N and 13C14N, respectively (wt time 0.8 s, 100 cycles). Isotope ratio data were extracted from individual cells based on cell morphology visualised in the 12C14N ion images. Image and data processing were performed as previously described for nanoSIMS analysis using the CAMECA software[23-26]. Each image of each cell was carefully examined by eye and the region of interest (ROI) was defined along the border of the 12C14N image and drawn by hand using the CAMECA software. The 13C/12C and 15N/14N ratios of 572 individual Skeletonema cells and of 306 individual Chaetoceros cells were determined for samples incubated with 13C and 15N tracers, whereas 13C/12C and 15N/14N ratios of 102 cells were determined from incubations without isotope tracers (controls). The small interference by 11B16O isotopes from the silicate frustule at a mass resolution of 6000 (M/ΔM) was corrected for when subtracting the natural apparent signal from 13C/12C in un-amended samples (controls) from that of amended samples. The average 13C/12C in un-amended samples (controls) was 0.0113 ± 0.0006 (stdev; n = 102), whereas that of amended samples was 0.0143 ± 0.0017 (stdev; n = 878).

C and N growth rates by cells

C-specific C-assimilation rates and N-specific N-assimilation rates, k (h−1), in the cells were calculated from the excess class="Chemical">13C/class="Chemical">pan class="Chemical">12C and 15N/14N isotopic ratios (IR) in individual cells relative to those measured in the un-amended controls, the excess %—isotopic ratios of bicarbonate or inorganics nitrogen in ambient water (Fbulk), and the incubation time (t1−t0) with stable isotopic tracers[24-26]:C and N content of diatom cells were calculated from their dimensions measured by light microscopy in parallel samples[28]. Doubling time was calculated as ln2/k. Cell-specific rates were calculated by multiplication of k with the cellular C or N content.

Ammonium analysis

Total class="Chemical">NH4+ concentration in the ambient class="Chemical">pan class="Chemical">water was measured in triplicate on a Trilogy® Laboratory Fluorometer (Turner Designs, USA)[51]. Triplicate subsamples for isotope analysis were transferred into 12 mL Exetainer©vials and biological activity stopped by adding 100 µL 50% ZnCl2 solution. The ratio of 15N/14N-NH4+ was analysed after helium de-gassing to remove 15N2 and 14N2 and subsequent chemical conversion of NH4+ to N2 with alkaline hypobromite (NaOBr)[52,53]. Isotopic ratios of 28N2, 29N2 and 30N2 were analysed by gas chromatography isotope ratio mass spectrometry on a Thermo Delta V isotope ratio mass spectrometer. NH4+ utilisation, assimilation and production rates were calculated from changes in 15NH4+ and 14NH4+ concentrations over time and the change in the isotope composition of particulate matter analysed by EA-IRMS[54]. The 15NH4+ concentration decreased exponentially over time, and 15NH4+ labelling percentage was averaged over the incubation period.

Ammonium production in aggregates

Seaclass="Chemical">water samclass="Chemical">ples were incubated in six 1.5 L roller tanks with a rotation sclass="Chemical">peed of 3 rclass="Chemical">pm within an hour after samclass="Chemical">pling[39]. After 24 h, 24 aggregates were collected from the roller tanks and incubated individually in class="Chemical">pan class="Disease">Falcon tubes with 0.2 µm filtered sea water (and six controls without aggregates)[6]. Ammonium concentrations in ambient water of the vials containing aggregates were significantly higher than those in the controls (Student’s t-test: p < 0.05; n = 6, 24). The ammonium distribution within and around aggregates were modeled as previously described[6].

Statistics

The statistical significance of differences in measured class="Chemical">carbon and class="Chemical">pan class="Chemical">nitrogen assimilation rates at the community level (EA-IRMS data) were analysed using two-way ANOVA test in SigmaPlot version 11. The replicates of cell-specific C and N-assimilation rates measured by SIMS were selected from the bottles with community nitrate assimilation rates closest to the mean value the nitrate assimilation measured by EAIRMS under still conditions and turbulent shear, respectively. Increasing numbers of the cellular C and N-assimilation rates measured by SIMS were calculated until the average values per cell were stable and the standard error was <5% of the average value to achieve representative average values for single cells of both genera of diatoms (Fig. 7). The C and N-assimilation rates calculated from the 13C/12C and 15N/14N ratios of the 572 individual Skeletonema cells and of the 306 individual Chaetoceros cells, and the 102 control cells were used in the statistical analysis.The C-assimilation rates followed a normal distribution as tested in detail using statistical analysis software (SAS) of their residuals. The C and N-assimilation data in the sub-groups of cells were also described by normal distributions and (non)-significant differences between these were analysed by two-way ANOVA test in SigmaPlot version 11 as well as by Student’s t-test in Microsoft Excell (2010). These two different statistical analyses gave similar results.

Data availability

The data that support the findings of this study are available and can be acquired from the corresponding author on rpan class="Chemical">equest.