- 1Department of Microbiology, Oregon State University, Corvallis, OR, United States
- 2Department of Statistics, Oregon State University, Corvallis, OR, United States
- 3Lawrence Livermore National Laboratory, Livermore, CA, United States
Introduction: Volatile organic compounds (VOCs) are small, low-vapor-pressure molecules emitted from the surface ocean into the atmosphere. In the atmosphere, VOCs can change OH reactivity and condense onto particles to become cloud condensation nuclei. VOCs are produced by phytoplankton, but the conditions leading to VOC accumulation in the surface ocean are poorly understood.
Methods: In this study, VOC accumulation was measured in real time over a 12 h day−12 h night cycle in the model diatom Phaeodactylum tricornutum during exponential growth.
Results: Sixty-three m/z signals were produced in higher concentrations than in cell-free controls. All VOCs, except methanol, were continuously produced over 24 h. All VOCs accumulated to higher concentrations during the day compared to the night, and 11 VOCs exhibited distinct accumulation patterns during the morning hours. Twenty-seven VOCs were associated with known metabolic pathways in P. tricornutum, with most VOCs involved in amino acid and fatty acid metabolism.
Discussion: Patterns of VOC production were strongly associated with diel shifts in cell physiology and the cell cycle. Diel VOC production patterns give a fundamental understanding of the first steps in VOC accumulation in the surface ocean.
Introduction
Phytoplankton are the primary source of the diverse array of labile dissolved organic compounds (DOC) in the surface ocean (Hansell, 2013; Moran et al., 2022), where DOC concentrations are maintained at 10–30 μM via uptake and metabolism by heterotrophic bacterioplankton (Hansell, 2013). Volatile organic compounds (VOCs) are chemicals with relatively low molecular weight and high vapor pressures that can comprise up to 30% of DOC (Hansell, 2013; Ruiz-Halpern et al., 2014). VOC physicochemical properties allow them to readily diffuse out of phytoplankton cells into seawater (Achyuthan et al., 2017; Zuo, 2019), where they can be taken up by bacterioplankton (Dawson et al., 2021; Moore et al., 2022; Simó et al., 2022), photochemically oxidized (Kieber et al., 1990; De Bruyn et al., 2020), or emitted into the atmosphere where they impact climate (Lelieveld et al., 2009). VOC air-sea emissions alter atmospheric chemistry and can influence the formation of secondary organic aerosols (SOA) (Quinn and Bates, 2011; Liu and Matsui, 2022) and cloud condensation nuclei (Charlson et al., 1987; Jackson et al., 2020; Mayer et al., 2020; Meskhidze and Nenes, 2006). The roles of VOCs in the marine microbial carbon cycle and climate highlight the need to understand the metabolic processes underlying VOC accumulation in the surface ocean.
VOC production is often associated with phytoplankton stress, which can occur when cells are exposed to light intensities or nutrient concentrations that are much higher or lower than the cells' acclimation status (Zuo et al., 2012; Meskhidze et al., 2015; Wang et al., 2021). For example, the well-studied VOCs, DMS and isoprene, are thought to provide antioxidant properties under high light stress (Zuo, 2019; Hackenberg et al., 2017; Dani et al., 2020; Fisher et al., 2020), iron limitation and UV exposure (Sunda et al., 2002), and nitrogen limitation (Sunda et al., 2007). Production of benzene and xylene was stimulated by ozone exposure in seawater communities, and those VOCs were suggested to protect cells from oxidative damage (Rocco et al., 2021). DMS is also produced by phytoplankton under grazing pressure (Fink, 2007; Fisher et al., 2020; Pozzer et al., 2022) and has recently been shown to serve as an info-chemical that attracts grazers to their phytoplankton prey (Shemi et al., 2023), potentially causing a positive feedback loop. The microalga, Microchloropsis salina, produced several volatile hydrocarbons that could serve as biomarkers of its grazing by a marine rotifer (Fisher et al., 2020). During harmful cyanobacterial algal blooms, oxidative stress in a freshwater lake caused saturated fatty aldehydes to increase in concentrations (Collart et al., 2023). Methanol production was highest at the onset of the stationary phase in cultures of phytoplankton from a range of taxonomic groups, suggesting nutrient limitation and cell envelope restructuring trigger methanol release (Mincer and Aicher, 2016), similar to methanol production associated with pectin demethylation in higher plants (Fall and Benson, 1996).
VOC concentrations are about ten-fold higher during optimal, non-stressed phytoplankton growth than during stationary phase (Shaw et al., 2003; Padaki et al., 2024) and may be a significant component of gross carbon production (Moore et al., 2020). Ethane, propane, and hexane increased in parallel with chlorophyll concentrations in diatom and dinoflagellate cultures (McKay et al., 1996) and in seawater (Ratte et al., 1993; Broadgate et al., 1997; Kameyama et al., 2009), suggesting their production by healthy, actively growing phytoplankton. Isoprene production rates were variable in a wide range of actively growing phytoplankton taxa (Colomb et al., 2008; Bonsang et al., 2010; Shaw et al., 2010) and correlated with photosynthetic activity in diatoms (Dani et al., 2020). In terrestrial plants, VOCs help maintain plant homeostasis during daily light fluctuations (Vivaldo et al., 2017; Liu et al., 2023) and are differentially produced depending on environmental stimuli (Broz et al., 2010; Mccormick et al., 2014). Some VOCs, such as acetaldehyde, acetone, and formaldehyde, are metabolic intermediates in biochemical pathways (Dudareva et al., 2013; Moore et al., 2020) active during normal growth. Acetaldehyde is both a product and a substrate in multiple pathways involving fatty acids, amino acids, carbohydrates, and precursors of some nucleic acids (Lees and Jago, 1978; Mazeed et al., 2021), which may explain its abundance in some VOC assays (Halsey et al., 2017). VOCs that are secondary metabolites, such as eugenol, and the isomers estragole and anethole, are antioxidants produced in healthy, growing cells (Wang et al., 2021; Rico et al., 2024).
Phytoplankton photosynthesis and cell division exhibit periodicity across the light-dark diel cycle (Bruyant et al., 2005; Mella-Flores et al., 2012; Su et al., 2024) and are tightly aligned with changes in the cell cycle (Su et al., 2024; Halsey et al., 2013). Phytoplankton cell cycles include three main phases: G1 encompasses most of the cell's biosynthetic activities and therefore spans most of the daytime (Brzezinski, 1992; Huysman et al., 2014). The S phase, or DNA synthesis, typically commences before dusk, while the G2-M phase, in which the cell replicates its organelles and divides into daughter cells, occurs during the night (Kim et al., 2017; Brzezinski, 1992). Photosynthesis increases intracellular metabolite concentrations during G1 to facilitate carbon precursors available for cell growth and division (Prézelin, 1992; Uchimiya et al., 2022). DOC in the surface waters of Suruga Bay, Japan, was often higher in concentration during the day compared to pre-dawn (Shinomura et al., 2005). These differences may be caused by phytoplankton primary production surpassing heterotrophic bacterial uptake of DOC (Shinomura et al., 2005). Daytime increases in DOC suggest that VOC production may also vary in concert with the diel cycle. During summertime in the North Atlantic Ocean, net VOC production (i.e., the rate of VOC production minus the rate of biological VOC losses) of methanol, DMS, and methanethiol were highest at midday (Davie-Martin et al., 2020). Isoprene and methanethiol exhibited clear diel variation in the diatom, Thalassiosira pseudonana (Kameyama et al., 2011), and higher day vs night concentrations of DMS, isoprene, and halogenated VOCs in a coral reef were attributed to biogenic and photochemical processes (Masdeu-Navarro et al., 2024).
Daytime accumulation of VOCs in the surface ocean would increase the likelihood of air-sea emissions if the accumulating VOCs are not consumed by other biological and chemical sinks (Halsey and Giovannoni, 2023). Therefore, the association of individual VOCs to the metabolic pathways from which they are produced is important to understanding the relationships between VOC production and air-sea emissions.
Diatoms contribute up to 40% of marine primary production (Field et al., 1998) and are known to produce a wide array of VOCs (Moore et al., 2020; Rocco et al., 2021). Here, we investigated diel variation in VOC production in the model diatom, P. tricornutum CCMP 2561. Cells were fully acclimated to sinusoidal light that mimicked the radiation in the surface ocean. VOCs were collected in real-time using proton transfer reaction time of flight mass spectrometry (PTR-TOF MS). A wide array of VOCs was detected, with VOC concentrations per cell peaking during the light phase. Many of the detected VOCs were associated with metabolic pathways to help explain their production patterns in the context of diel and cell cycle physiology. Knowledge of light-driven VOC production by phytoplankton is important to inform mechanisms of VOC accumulation in the surface ocean and their potential for air-sea VOC emissions.
Materials and methods
Culture growth conditions
Axenic P. tricornutum CCMP 2561 was cultivated under a sinusoidal 12 h:12 h light: dark cycle (0 to 400 μmole photons m−2 s−1) in f/2+Si artificial seawater medium (ASW) (Guillard and Ryther, 1962). Axenicity was confirmed prior to and during diel experiments by staining culture samples daily using SyBr Green, and fluorescence distributions analyzed by GUAVA flow cytometer (Millipore; Billerica, MA, USA) for any indication of bacterial contaminants <3 μm. Second, culture samples were analyzed prior to each diel experiment via fluorescent microscopic examination of DAPI-stained cells. The culture growth light was delivered by a programmable phyto-panel composed of white and blue light-emitting diodes (LEDs) (Photo Systems Instruments). Irradiance was measured with a 4π spherical quantum meter (Biospherical Instruments QSL-100). Cultures were maintained at 19 +/- 0.5 °C in an acrylic incubator with submerged copper tubing connected to a recirculating water bath (Figure 1).

Figure 1. The diel VOC experimental setup with four 1-L cultures of axenic P. tricornutum. One culture was designated for real-time VOC measurements over 24 h (marked with a red arrow), and the remaining three were used for supporting measurements of cell counts, Chl-a, and photosynthetic efficiency. Temperature was maintained at 19 +/- 0.5 °C using copper coils circulating cold water around the acrylic incubator. Cultures were bubbled with breathing-grade air at 100 ml min−1. This setup was repeated four times with axenic P. tricornutum and four times with ASW-only controls.
Diel experimental setup
For each of four replicated experiments, the incubator held four P. tricornutum cultures (0.6 L) in 1-L polycarbonate Nalgene® Bottles. One culture was dedicated to VOC measurements, and the remaining three were used for supporting measurements (Figure 1). Abiotic control experiments (n = 4), in which the bottles contained only sterile ASW, were conducted after each of the four replicated culture experiments. Abiotic controls were conducted in bottles and with tubing identical to the experimental cultures. Experimental and abiotic control cultures were aerated with Ultra Zero Grade Air, and the headspace was vented through 0.2-μm Polytetrafluoroethylene (PTFE) filters (Figure 1). Cultures were kept in the exponential growth phase, using semi-continuous culturing, under these conditions for 2 weeks and then used as inocula for the diel VOC experiment. ASW prepared less than 2 days prior to initiating the diel VOC experiment was aerated for at least 6 h to remove dissolved VOCs originating from lab air and then inoculated with the axenic culture or used for the abiotic control. VOC measurements were initiated when cell densities reached 2.25–4.5 x105 cells mL−1.
Supporting measurements
Supporting measurements were collected every 2 h from culture samples pulled from ports in the bottle caps. Cell densities were measured by GUAVA flow cytometer (Millipore; Billerica, MA, USA). Chlorophyll concentration (Chla) was determined by filtering 2–5 mL of culture (GF/F, Whatman, 25 mm), extracting in 5 mL 90% acetone at −20 °C for 24 h, and then measuring by spectrophotometer (Shimadzu; Kyoto, Japan) as described in Ritchie (2006). Photosynthetic efficiency (Fv/Fm) was measured by a fast repetition-rate fluorometer (FRR) following 10 min dark acclimation (Kolber et al., 1998).
Proton transfer reaction time of flight mass spectrometry (PTR-TOF MS)
Real-time VOC measurements were initiated by removing the filter from the vent line and attaching the line directly to the PTR-TOF MS inlet. VOCs were measured via soft ionization with H3O+. The mass spectra (30–240 atomic mass units) were acquired in real-time for 24 h and then binned into 2-h intervals to match up with supporting measurements. The first 15 min in the first hour of the 24 h acquired data were removed from the analysis because of the lag time between the culture headspace and detector. PTR-TOF MS data were analyzed using PTR viewer 3.4.3 (Ionicon Analytik). Raw data in .h5 format were first calibrated against ions of known mass and internal standards (nitrosonium, NO+, at m/z 29.998 and 1,3-diiodobenzene at m/z 203.943 and m/z 330.848) using the following criteria: 3-point calibration mode, cycle 350, 0.2 m/z search range, and three spectra averaging. The accuracy of the mass calibration was confirmed against the primary ions (hydronium ions, O+ at m/z 21.022, and H2O+ at m/z 39.035). In the PTR Viewer 3.4.3 software, each .h5 file shows an RMS error ranging from 0.006 to 0.008.
The mass spectra were visualized in the PTR-Viewer 3.4.3 (Ionicon Analytik) using the Gaussian approach, ion-mass percentage, and correctness percentage. A peak table previously collected from P. tricornutum (Padaki et al., 2024) was used as a reference to select m/z signals. Each peak (m/z signal) was tentatively identified using GLOVOCS (Yáñez-Serrano et al., 2021) and the ChemSpider database (integrated into PTR Viewer 3.4.3). Integrated peak signals were normalized to H3O+ concentration. Isotopic signatures of m/z signals (M) with known chemical formulas were assessed for M+1 (13C- and 15N-containing VOCs) and M+2 (18O- and 34S-containing VOCs) isotope peaks. M+1 and M+2 peak concentrations were added to M concentrations. M+1 peaks were detected for the hydrocarbons, benzene, toluene, and ethylbenzene/xylenes, C3H6, C9H10. Both 13C and 15N M+1 isotope peaks were detected for acetonitrile. The M+2 peak was detected for acetaldehyde. Detected analyte concentrations (ppbv) were calculated using the simple-reaction-kinetics approach with 30% accuracy for E/N values >100Td (E/N for this study = 126Td) (Holzinger et al., 2019).
Paired t-tests were used to identify m/z signals in each culture that were present in higher concentrations than the abiotic controls (n = 4). Benjamini-Hochberg procedure was used to control the false detection rate of m/z signals, with q-values >0.1 removed from further analysis. m/z signals from each replicate were adjusted by subtracting the corresponding m/z signal measured in abiotic controls. The adjusted m/z signals were averaged and then converted from ppbv to molar concentrations using the formula:
Where C* is the molarity of the analyte in the culture, C-B is the adjusted analyte mixing ratio (ppbv), Q is the flow rate of air bubbled through the culture (3 L h−1), p is the atmospheric pressure (1 atm), t is the time period of collection (1.75 h for the first time period and 2 h thereafter), V is the volume of the culture (0.6 L), R is the gas constant (0.0821 L atm K−1 mol−1), and T is the incubated culture temperature (292 K). Daily production rates for each m/z signal were determined by summing the concentrations (fmol per cell) over 12 time points, i.e., the 24h period.
Statistics
Statistically significant differences in peaks were analyzed using t-tests and q-values. Figures were created in R Studio 2024.12.1 + 563, R v.4.1.1, utilizing the fda, stats, ggplot2, and ComplexHeatmap packages (Gu et al., 2016; Gu, 2022), and the VOLCALC R-package (Meredith et al., 2023). The arrangement of figure panels was completed in Adobe Illustrator.
B-spline curve fits
A smooth curve was fit for the temporal change in the mean concentrations of each peak via a regression spline function. To this end, a set of B-Spline basis functions were generated using the timepoints of the 24-h cycle. A linear combination of the B-spline basis functions gave a smooth curve to fit the mean concentrations of each peak. The coefficients of the above linear combination were estimated by balancing a goodness-of-fit criterion and a smoothness penalty to avoid underfitting or overfitting. Spline regression allows the incorporation of local and non-linear effects with a continuous covariate, which was observed in the temporal change in the peak concentrations. In this application, a B-spline basis function of the 3rd order with 3 internal knots was used, resulting in 6 estimated B-spline coefficients per peak. Subsequent cluster analysis was performed using K-means on the six estimated B-spline coefficients to form groups of peaks with similar temporal changes in their mean concentrations (Supplementary Figure 1).
Integrating m/z putative IDs to CHEBI IDs, KEGG, and RVIs
All m/z identifications extracted from the GLOVOCs database were associated with their ChEBI identifiers from the Chemical Identities of Biological Interest Database (ChEBI) (Hastings et al., 2012). The ChEBI identities were used to locate metabolites in the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases for P. tricornutum or plants (Kanehisa, 2000; Kanehisa et al., 2016). Finally, ChEBI identities were input into the VOLCALC R-package (Meredith et al., 2023) to determine their relative volatility indices (RVIs).
Results and discussion
Cell physiology
P. tricornutum grown on a 12 h:12 h light:dark sinusoidal regime exhibited a specific growth rate of 0.84 ± 0.08 d−1, consistent with exponential growth. The cell density stayed constant from dawn to 2 h prior to dark and then increased 2.1-fold over the remaining 14 h (Figure 2). Cell Chla content doubled from dawn to midday, reaching 0.41 ± 0.04 pg cell−1 and then slowly decreased to the initial amount by dawn. Fv/Fm was maximal (0.67 ± 0.01) at dawn and during the night but decreased at higher light intensities, consistent with photoinhibition during the midday (Domingues et al., 2012; Jallet et al., 2016) (Figure 2).

Figure 2. Cell physiology of P. tricornutum during the 24 h diel experiment. Light intensity is shown at the top, and the shaded area from 12-24 h corresponds to the night.
Temporal changes in VOCs in P. tricornutum cultures
A total of 63 m/z signals were detected in the P. tricornutum cultures. Sixty-two of these signals were present at concentrations higher than media controls throughout the diel cycle, and one, m/z 33.03, corresponding to methanol, was only detected at concentrations higher than the controls from 0-8 h (Figure 3, FDRB − H <0.1, n = 4).

Figure 3. Temporal changes in relative VOC concentrations in P. tricornutum cultures. (a) Heatmaps showing concentrations of m/z signals normalized to their maximum values (rows) or (b), concentrations of m/z signals normalized to cell density and then to their maximum values (rows). Relative concentrations ranged from low (light purple) to high (dark purple). Temporal patterns of m/z signal concentrations clustered into 6 different groups are shown in the leftmost column of each heatmap. n = 4 independent replicated cultures of P. tricornutum.
To explore the data for temporal variations in VOC accumulation, we normalized the concentration of each m/z signal collected in real time during each 2 h interval in each of the four culture replicates to its maximum concentration over the 24 h experiment, and then took the average of the four replicates for each 2 h time period (Figure 3a). To decipher differences in the temporal patterns of each VOC, we fit a B-spline function with three internal knots positioned at 6, 12, and 18 h to the concentration data. These positions were selected based on the major changes in cell physiology associated with Fv/Fm and cell division (Figure 1). Spline regression was used to fit curves, using third-order B-spline basis functions. The resulting curves were grouped into six clusters using K-means on their respective estimated B-spline coefficients (Figure 3a, Supplementary Figure 1).
VOCs in three clusters (1, 3, 5) were maintained at relatively constant levels throughout the diel cycle, albeit with slight variations in the early morning (Cluster 3), midday (Cluster 1), or late night (Cluster 5). Forty-six of the 63 m/z signals were in Clusters 3 and 5 (Table 1). VOCs in the remaining three clusters (2, 4, 6) exhibited a wider range of variation over the 24 h cycle. Cluster 2, consisting of only methanol (m/z 33.03), was highest at dawn and early morning and decreased to zero during the early night (Figure 3a). Cluster 4 increased from dawn until 6 h into the night. VOCs in Cluster 6 increased from dawn to midday and then remained at relatively constant levels through the night (Figure 3a).

Table 1. Distribution of m/z signals in each cluster identified using B-spline coefficients and K-means clustering.
To evaluate cell-specific patterns in VOC concentrations over the diel cycle, the concentration of each m/z signal in each of the four culture replicates was divided by the cell density at each time point and then normalized to the maximum value over the 24 h period. The same 3-knot B-spline regression framework was used to fit curves for each m/z signal, and the m/z signal curves were again grouped into 6 clusters using K-means on their respective estimated B-spline coefficients, designated as Clusters 1* - 6* (Figure 3b). Fifty-one m/z signals were grouped into the same clusters as the un-normalized clusters (Figure 3a), and 12 m/z signals were grouped into different clusters (Figure 3, Table 1). Cluster 3* included 32 m/z signals, and 52 of the 63 m/z signals were categorized into Clusters 1*, 3*, or 5*.
Relative VOC concentrations in Clusters 1*, 3*, and 5* decreased by approximately half over the diel cycle, reflecting increasing P. tricornutum cell density over the course of the experiment. Cluster 1* showed two slight peaks at 6–8 h and 18–20 h. Cluster 3* was maximal from 0–6 h, after which relative concentrations gradually decreased over the remaining 24 h. Cluster 5* decreased steadily throughout the 24 h period (Figure 3b). The patterns exhibited by Clusters 2 and 2*, which consisted only of m/z 33.03, methanol, were independent of cell density. Methanol temporal dynamics were dominated by the changes in its concentrations, which showed strong signals between 0 and 8 h and ranged from 0.9 to 1.5 nM in the cultures (Figure 4, Supplementary Tables 1, 2). Between 10–22 h, the concentrations of methanol were lower than the detection limit [0.24 nM or 7.87 ppbv at 19 °C, (Davie-Martin et al., 2020)], which may indicate that methanol production and its magnitude of diel variation were underestimated in this experiment. The three m/z signals in Cluster 4* exhibited some intra-cluster variation, but these m/z signals all increased during the daytime and nearly returned to the dawn concentrations by 24 h. Cluster 6*, consisting of 7 m/z signals, peaked at midday and then returned to the dawn concentrations by 24 h.

Figure 4. Temporal changes in absolute concentrations of m/z signals (rows) over the diel cycle in P. tricornutum. (left) Heatmap showing concentrations of m/z signals in the cultures ranging from low (light blue) to high (dark blue), and (right) heatmap showing cell-specific concentrations of m/z signals in the cultures ranging from low (light red) to high (dark red). m/z values are shown in the column between the two heatmaps. A total of 27 m/z values were found P. tricornutum KEGG pathways (magenta font) or plant KEGG (yellow font) pathways. The cell-specific clusters (1* - 6*) into which the m/z values were grouped are shown in the rightmost column.
Isostatic VOC production and light-driven VOC accumulation
The clustering of VOCs determined from their concentrations over the diel cycle was influenced by cell division and the light-dark phases (Figure 3). The general decreases in relative VOC concentrations from 10 to 24 h in Figure 3b were caused by accounting for cell densities that increased during that period. In contrast, absolute m/z signals in Figure 3a (Clusters 1, 3, 5) were maintained at approximately the same concentrations over the entire 24 h period even as the P. tricornutum population increased during the nighttime. We call this behavior “isostatic VOC production.” In this case, passive diffusion of VOCs allows them to move from inside the cell, through the cell envelope, and into the bulk medium, where they achieve equilibrium with intracellular VOC concentrations. In the absence of a strong biological sink (i.e., bacterial VOC consumption) (Sun et al., 2011; Beale et al., 2013; Sargeant et al., 2018; Moore et al., 2020) or physical sink (more rapid bubbling that strips VOCs out of the culture or aeration by shaking), VOC concentrations remained stable. As the cell density increased over the nighttime, VOC production continued, albeit at slower cell-specific rates, to maintain the same VOC concentrations in the culture medium. VOCs in Clusters 4 and 6 accumulated rapidly during the light phase until they peaked in concentrations. We call this behavior “light-driven VOC accumulation.” VOCs in Clusters 4 and 4* reached maximum concentrations after 12 h, and VOCs in Clusters 6 and 6* reached maximum concentrations after 6 h. After reaching their maxima, VOC concentrations in Clusters 4* and 6* gradually decreased as the cell density increased.
In addition to relative VOC concentrations, we also examined the absolute concentrations of the 63 m/z signals over the diel cycle (Figure 4). Individual VOCs accumulated in the cultures to concentrations ranging from 0.03 to 13.68 fmol per cell over the 24 h period, corresponding to rates of VOC production over 24h ranging from 0.5 to 125 fmol (cell·d)−1 (determined by summing the 2 h average concentrations over the 24 h period) (Figure 4b). Seven m/z signals accumulated to concentrations at the high end of the range of 7-13 fmol cell−1. Twenty-three m/z signals reached concentrations exceeding 1 fmol cell−1, including m/z 43.05, 45.03, 47.04, 59.05, 75.04, and 91.03 (Figure 4, Supplementary Table 2). Methanol (m/z 33.033) concentration varied up to 5-fold over the diel cycle, with cell-specific concentrations of 3–5 fmol cell−1 at 0 to 6 h and 22 h but below the limit of detection between 6 and 22 h (Figure 4). Other m/z signals, including toluene (m/z 93.07) and styrene (m/z 105.06), varied up to 2-fold over the 24 h cycle.
To further explore VOC temporal dynamics, we focused on clusters determined from cell density-normalized m/z signal concentrations (Clusters 1*-6*; Figure 3b) because they represent VOC dynamics related to diel-dependent physiology (Sartori et al., 2021; Ledford and Meredith, 2024), including the cell cycle. Each cluster contained at least one m/z signal produced at >50 fmol (cell·day) −1. Forty-eight of the 63 m/z signals were putatively identified using the GLOVOCs database, and twenty-two of those 48 m/z signals were found in P. tricornutum KEGG pathways (Table 2; Figure 5, Supplementary Table 3). Five of the m/z signals not found in P. tricornutum KEGG pathways were present in KEGG pathways in plants. The remaining 21 putatively identified m/z signals represent metabolites not previously known from P. tricornutum or plant KEGG pathways. An additional 15 m/z signals were unknown (Figure 4, Supplementary Table 3).

Figure 5. Putatively identified VOCs in KEGG and their metabolic pathways in P. tricornutum. In the first column, m/z signals described in P. tricornutum KEGG pathways are colored magenta and plant KEGG pathways are colored yellow. Putative VOC identifications for m/z signals with two or three isomers are listed with slashes. Each black box indicates the presence of the VOC in the pathway indicated at the top of each column. RVI, relative volatility index.
Eight P. tricornutum metabolic pathways contained the 27 putatively identified VOCs found in KEGG (Figure 5). Some VOCs are metabolites in multiple pathways. For example, acetaldehyde (m/z 45.03) is a product or substrate in six of the eight VOC pathways (Figures 5, 6), highlighting its pivotal role in the interchange of carbon atoms amongst different cellular processes. The vast majority of VOCs were found in either amino acid metabolism (16 VOCs) or fatty acid synthesis (seven VOCs) (Figure 5). Twenty of the VOCs identified in KEGG were grouped into the isostatic VOC production Clusters 1*, 3*, and 5* and were distributed across all eight P. tricornutum pathways (Figures 3b, 5). Seven VOCs, categorized into the light-driven accumulation Clusters 4* and 6*, were associated with amino acid metabolism and secondary metabolite biosynthesis.

Figure 6. Proposed metabolic network of isostatic VOCs in central carbon metabolism. VOCs in Clusters 1*, 3*, and 5* are involved in Glycolysis/Gluconeogenesis, TCA cycle, and Fatty acid metabolism. Acetaldehyde is involved in nine different reactions and is a central node of intersecting pathways. The association of these isostatic VOCs with methanol in Cluster 2* is shown for clarity. Clusters numbered with an asterisk are from Figure 3b and Figure 4. Solid arrows with an asterisk represent acetaldehyde as an intermediate product; solid arrows represent single-step reactions; dashed arrows represent multistep reactions; and dotted arrows represent metabolites transferred between respective pathways.
We computed the relative volatility index (RVI) for each VOC identified in KEGG. The higher the RVI, the higher the vapor pressure of the compound. RVI may indicate how readily a VOC diffuses across the membrane. RVI will also impact the movement of VOCs into the gas phase and out of the culture in an open system like those used in this study. All of the VOCs identified in KEGG had RVI values >2 and are thus classified as VOCs with ‘high volatility' (Meredith et al., 2023). “Isostatic” VOCs exhibited RVIs in the range of 8–10, whereas VOCs that accumulated in the light had slightly lower RVIs of 5–7 (Figure 5). This difference suggests that diffusive loss of VOCs with high RVIs requires cells to continually replenish their intracellular VOC pools, resulting in the “isostatic production” observed in Clusters 1*, 3*, 5*. VOCs in clusters 4* and 6*, with lower RVIs, are less readily diffused out of the culture. These VOCs were upregulated during the light phase, indicating their increased demand during light-driven cell metabolism (Figure 5).
VOCs in central carbon metabolism and fatty acid metabolism exhibited isostatic production
Isostatic VOCs present in the highest concentrations were part of central carbon metabolism. Most of these highly produced VOCs are metabolic intermediates connecting glycolysis and the TCA cycle (Figure 6). Ethanol (m/z 47.05) and acetaldehyde (m/z 45.03), products from glycolysis/gluconeogenesis, can be metabolized into lactic acid (m/z 91.03) and butanal (m/z 73.06), which feed into fatty acid metabolism. Similarly, the isomers acetone/propanal (m/z 59.05) produced by glycolysis feed into fatty acid metabolism. PTR-TOF MS cannot distinguish between isomers of the same mass. Acetone/propanal is generated as a by-product of acetyl-CoA or acetoacetate decarboxylation in plants (Suganuma et al., 1993) and subsequently in P. tricornutum can enter either fatty acid metabolism or be converted to methyl acetate and hydrolyzed to methanol (Figure 6). Methanol oxidation to formaldehyde (m/z 31.01) sources C1 units in the tetrahydrofolate (THF) pathway (Figure 6, Supplementary Table 2). Additionally, m/z 89.06 corresponding to either acetoin or butanoic acid, depending on the molecular structure, are end products of glycolysis or converted to 2,3 butanediol. Acetoin buffers against pH fluctuation in plants, and its conversion to 2,3-butanediol creates a signal for plant growth (Forlani et al., 1999; Xiao and Lu, 2014). Acetaldehyde, ethanol, and acetone are commonly described as products of stress in plants (Kelsey and Westlind, 2017; Jardine and McDowell, 2023) and algae (Catalanotti et al., 2013; Yu and Li, 2021). For instance, in response to zinc-sulfate deficiency (salt stress), acetaldehyde and ethanol production increased ~3-fold in Chlamydomonas reinhardtii compared to replete controls (Van Lis et al., 2017). In our experiments, cells were in balanced, exponential growth under light conditions typical of the summer surface ocean mixed layer. Even at the peak light intensity (400 μE), when photoinhibition was observed, concentrations of acetaldehyde, ethanol, and acetone remained nearly constant (Figure 2). Thus, although VOC production by algae is observed during stress, many VOCs are produced, even in relatively high concentrations, in healthy growing cells (Halsey et al., 2017; Moore et al., 2020).
A goal of this research was to understand how VOC production is related to the diel cycle. The cell-specific production of isostatic VOCs in Clusters 1*, 3*, and 5* decreased by half over the course of the night. Even at night, some production of these VOCs occurred to maintain their intra-extracellular equilibrium as they were consumed either through metabolism or removal from the system by gentle bubbling. Carbon metabolism in P. tricornutum is closely coordinated with the diel cycle (Chauton et al., 2013). Carbon fixation, gluconeogenesis, and fatty acid biosynthesis are upregulated during the early morning when VOCs in clusters 1*, 3*, and 5* were at their highest concentrations. The TCA cycle and fatty acid beta-oxidation are more highly expressed during the late afternoon (Smith et al., 2016) when VOC production begins to slow. Glycolysis and the TCA cycle are active during the day and night, but the expression of these pathways is partly controlled by their spatial localization in the cytosol, chloroplast, or mitochondria (Smith et al., 2016; Broddrick et al., 2019) where glycolytic isozymes are differentially expressed during the day and night (Chauton et al., 2013; Smith et al., 2016). The amphibolic isozyme, glucose-6-phosphate isomerase, GPI_3, functions in glycolysis/gluconeogenesis in the chloroplast and is upregulated during the day (Smith et al., 2016). In contrast, glucose-6-phosphate isomerase, GPI_1, functions in the cytosol and is most highly expressed at night in tandem with genes in the TCA cycle (Kroth et al., 2008; Chauton et al., 2013). The cell's constant requirement for carbon precursors requires non-stop activity of central carbon metabolism. Consequently, isostatic VOCs in Clusters 1*, 3*, and 5* are produced continuously, at levels that depend on the activities of their respective pathways at different times of the day.
VOCs in amino acid metabolism and secondary metabolite biosynthesis accumulated during the light phase
VOCs in Clusters 4* and 6* accumulated during the light phase and are produced from aromatic amino acid metabolism and the biosynthesis of secondary metabolites (Figures 5, 7). In Cluster 4*, styrene (m/z 105.06) is a by-product of p-coumaric acid and caffeic acid metabolism in secondary metabolite biosynthesis (Figure 7). Secondary metabolite production and amino acid biosynthesis co-occur with protein and carbohydrate accumulation in plants (Liebelt et al., 2019). Benzenoid compounds from Cluster 6*, including toluene (m/z 93.07) and the isomers ethylbenzene/xylenes (m/z 107.08), are produced by aromatic amino acid metabolism. Their production during the daytime is consistent with the daytime upregulation of biosynthesis genes for tyrosine and phenylalanine (Smith et al., 2016). Also, in Cluster 6*, the isomers anethole/estragole (m/z 149.06, Figure 7), metabolized from p-coumaric acid, protect against high light stress and have antimicrobial activity in plants (Koeduka et al., 2009).

Figure 7. Metabolic network of VOCs exhibiting light-driven accumulation. VOCs from clusters 4* and 6* are linked to aromatic amino acid metabolism and secondary metabolite biosynthesis. Clusters numbered with an asterisk are from Figure 3b and Figure 4. Solid arrows represent a single-step reaction; dashed arrows represent a multistep reaction; and dotted arrows represent predicted metabolites from respective pathways.
The roles of mono-aromatic benzenoids, including benzene, toluene, ethylbenzene/xylenes, and styrene, in algae remain largely unknown. They have been proposed to serve as precursors for aromatic compound synthesis in plants (Misztal et al., 2015). In the cytosol of plants, approximately 30% of the carbon derived from photosynthesis is directed toward phosphoenolpyruvate and erythrose-4-phosphate synthesis, which feed the shikimate pathway (Maeda and Dudareva, 2012) used for biosynthesis of aromatic amino acids phenylpropanoid and benzenoid production (benzyl alcohol, benzaldehyde, benzyl benzoate, methyl benzoate, phenylacetaldehyde, and phenyl ethanol) (Orlova et al., 2007). Alternatively, some benzenoids may be products of protein degradation (Dudareva et al., 2013). Suppression of benzenoid compound production in plants affected their morphology and decreased transport of auxin, a hormone that facilitates cell division and growth (Orlova et al., 2007). It is possible that benzenoids have roles associated with the cell cycle in diatoms. Unexplained is why benzene (m/z 79.05) from Cluster 1*, likely produced from aromatic amino acid metabolism, appears to be isostatically produced rather than showing accumulation during the light phase as the other VOCs functioning in aromatic amino acid metabolism pathways.
Methanol (m/z 33.03) temporal behaviors may reflect P. tricornutum cell cycle dynamics
The temporal pattern exhibited by methanol (m/z 33.033) (Cluster 2*; Figures 3, 4) suggests its production is strongly associated with the diatom's mitotic cell cycle. In P. tricornutum, the S-phase, DNA synthesis, begins 2 h before dark (Figure 8) (Brzezinski, 1992; Kim et al., 2017). In our experiment, methanol was high in concentration during the daytime and dropped below the limit of detection 2 h before dusk, aligning with S-phase (Figure 8). This temporal correspondence suggests the intracellular methanol pool was rapidly consumed to generate C1 equivalents for nucleic acid biosynthesis during DNA synthesis. Methanol can be oxidized to formaldehyde (Figure 6) and spontaneously react with tetrahydrofolate (THF) to form 5,10-methylene-tetrahydrofolate, which donates C1 units for nucleotide biosynthesis (Nuccio et al., 1995; He et al., 2020; Bou-Nader et al., 2021). Methanol concentrations increased as the cells entered another growth phase. Formaldehyde was grouped in Cluster 1*, which showed little change over the diel cycle, possibly reflecting the constant need for C1 units in RNA biosynthesis and energy production throughout the cell cycle.

Figure 8. Relationships between VOC metabolism and the cell cycle in P. tricornutum. (a) Activities of pathways involved in VOC production are indicated by the shading of the purple bars (dark purple shows periods of greatest activity and most VOC accumulation, and light purple is lower activity and VOC production). Purple bars for each metabolic pathway are aligned with the diel cycle at the top, and cell cycle phases are shown in pink. (b) Metabolism of methanol by the THF pathway was most active during S-phase, causing rapid methanol depletion in culture. The font color of each VOC indicates its grouping into each cluster.
Isoprene, DMS, methanethiol, and their diel patterns
Isoprene (m/z 69.07; Cluster 1*); dimethyl sulfide (DMS, m/z 63.02; Cluster 3*), and methanethiol (m/z 49.01; Cluster 3*) were among the VOCs that were isostatically produced. These metabolites are commonly observed during cell stress, such as high light (Van Rijssel and Gieskes, 2002; Halsey et al., 2017; Steinke et al., 2018) and grazing (Ploug and Grossart, 2000; Fisher et al., 2020). Isoprene exhibits antioxidant properties (Pollastri et al., 2021; Zhao et al., 2024), but the balanced growth conditions in our experiments likely did not induce photic stress, which could explain low isoprene, DMS, and methanethiol concentrations. Daytime maxima in the concentrations of these compounds have been observed in other diatom cultures (Kameyama et al., 2011), coral reefs and the open ocean (Galí et al., 2013; Davie-Martin et al., 2020; Masdeu-Navarro et al., 2024). However, the isostatic production of isoprene, DMS, and methanethiol observed in this study suggests that other factors besides daily changes in light intensity, such as bacterial heterotrophic consumption (Boden et al., 2011; Moore et al., 2022), daytime grazing (Ng and Liu, 2016), photooxidation (Kieber et al., 1990; Masdeu-Navarro et al., 2024), and variable mixing depth (Halsey and Giovannoni, 2023) are likely to underlie daily fluctuations in VOC concentrations.
VOCs measured in this study exhibited diel variations associated with changes in the light and dark reactions of photosynthesis that shift over hourly timescales (Halsey et al., 2013). The culture experiment was limited to 24 h to ensure that cells remained in exponential growth. The observed VOC dynamics are expected to repeat if cell physiology is maintained via dilution with fresh media and light regime replicated. Subtle differences in the temporal dynamics of VOCs produced from metabolic pathways common across phytoplankton are hypothesized to be associated with species-specific regulation of their metabolic pathways. An important goal will be to understand taxonomic variation in VOC production patterns, which could help describe surface ocean VOC accumulation depending on community composition and may be predictable given sufficient understanding of the genetic components involved.
Summary
Relating in-water biological VOC production to atmospheric chemistry processes, such as ozone formation and depletion and condensation reactions that form SOA, is challenging because of the complexity of compounds involved, diversity of biological sources and sinks, and a lack of knowledge about the conditions that influence VOC accumulation in the surface ocean. Here, we showed that VOCs were produced throughout the diel cycle in a cultured diatom, with the highest rates of VOC production occurring during the daytime and involved in amino acid and fatty acid biosynthesis. Light intensities were saturating at midday, causing slight photoinhibition, but the patterns of VOC production were more strongly associated with diel shifts in cell physiology than stress responses. We documented strong VOC accumulation during the morning hours, consistent with the requirement for photosynthesis to produce sufficient organic carbon to support cell needs for growth and division. Nighttime metabolism includes cell division, protein synthesis, and pre-dawn preparation for photosynthesis. Thus, nighttime VOC production was caused by the ongoing flux of carbon through glycolysis, the TCA cycle, and fatty acid beta-oxidation. Positive air-sea emissions (into the atmosphere) of benzene, toluene, and xylenes (ethylbenzene + xylene) can be at least as large as DMS emissions, indicating these benzenoids should be recognized as major contributors to atmospheric chemistry and SOA formation. Our results show these benzenoids primarily originate from daytime photosynthetic processes in the model diatom P. tricornutum.
Diel VOC production dynamics give a fundamental understanding of the first steps in VOC accumulation in the surface ocean. The complex biological and abiotic processes that next act upon in-water VOCs will impact their concentrations in ways that are specific to each VOC. To understand the chemical composition of the surface ocean and overlying atmosphere will require an interdisciplinary approach to incorporate into next generation ocean-atmosphere models the dynamic roles of biological VOC metabolism, chemical transformations, and physical transport.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
VP: Writing – review & editing, Validation, Formal analysis, Software, Investigation, Data curation, Conceptualization, Methodology, Visualization, Writing – original draft. EP: Writing – review & editing, Data curation, Formal analysis, Writing – original draft, Methodology. YJ: Formal analysis, Writing – review & editing, Writing – original draft. HB: Writing – original draft, Writing – review & editing, Methodology. JK: Writing – review & editing, Data curation, Writing – original draft. KH: Validation, Writing – review & editing, Conceptualization, Investigation, Funding acquisition, Supervision, Methodology, Project administration, Resources, Writing – original draft.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by award ID OCE1948163 from the National Science Foundation, Biological Oceanography Program to KH and award ID 80NSSC24K0480 from the National Aeronautics and Space Administration to KH. Publication of this manuscript was supported by an award to VP from Applied and Environmental Microbiology Gordon Research Conference. This work was also supported by the US Department of Energy's (DOE) Genomic Science Program through the LLNL μBioSpheres Science Focus Area grant # SCW1039 and carried out at Lawrence Livermore National Laboratory (LLNL) under Contract DE-AC52-07NA27344.
Acknowledgments
We thank Kaylene M. Abraham and Kerri Jacobs for experimental assistance and Nicholas Baetge and James Fox for help with data analysis. We also thank Xavier Mayali and Peter Weber at Lawrence Livermore National Laboratory for the helpful discussions and interest in this VOC research that improved manuscript presentation.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1620542/full#supplementary-material
References
Achyuthan, K., Harper, J., Manginell, R., and Moorman, M. (2017). Volatile metabolites emission by in vivo microalgae-an overlooked opportunity? Metabolites 7:39. doi: 10.3390/metabo7030039
Beale, R., Dixon, J. L., Arnold, S. R., Liss, P. S., and Nightingale, P. D. (2013). Methanol, acetaldehyde, and acetone in the surface waters of the Atlantic ocean: OVOCs in the atlantic ocean. J. Geophys. Res. Oceans 118, 5412–5425. doi: 10.1002/jgrc.20322
Boden, R., Murrell, J. C., and Schäfer, H. (2011). Dimethylsulfide is an energy source for the heterotrophic marine bacterium Sagittula stellata. FEMS Microbiol. Lett. 322, 188–193. doi: 10.1111/j.1574-6968.2011.02349.x
Bonsang, B., Gros, V., Peeken, I., Yassaa, N., Bluhm, K., Zoellner, E., et al. (2010). Isoprene emission from phytoplankton monocultures: the relationship with chlorophyll-a, cell volume and carbon content. Environ. Chem. 7:554. doi: 10.1071/EN09156
Bou-Nader, C., Stull, F. W., Pecqueur, L., Simon, P., Guérineau, V., Royant, A., et al. (2021). An enzymatic activation of formaldehyde for nucleotide methylation. Nat. Commun. 12:4542. doi: 10.1038/s41467-021-24756-8
Broadgate, W. J., Liss, P. S., and Penkett, S. A. (1997). Seasonal emissions of isoprene and other reactive hydrocarbon gases from the ocean. Geophys. Res. Lett. 24, 2675–2678. doi: 10.1029/97GL02736
Broddrick, J. T., Du, N., Smith, S. R., Tsuji, Y., Jallet, D., Ware, M. A., et al. (2019). Cross-compartment metabolic coupling enables flexible photoprotective mechanisms in the diatom Phaeodactylum tricornutum. New Phytol. 222, 1364–1379. doi: 10.1111/nph.15685
Broz, A. K., Broeckling, C. D., De-la-Peña, C., Lewis, M. R., Greene, E., Callaway, R. M., et al. (2010). Plant neighbor identity influences plant biochemistry and physiology related to defense. BMC Plant Biol. 10:115. doi: 10.1186/1471-2229-10-115
Bruyant, F., Babin, M., Genty, B., Prasil, O., Behrenfeld, M. J., Claustre, H., et al. (2005). Diel variations in the photosynthetic parameters of Prochlorococcus strain PCC 9511: combined effects of light and cell cycle. Limnol. Oceanogr. 50, 850–863. doi: 10.4319/lo.2005.50.3.0850
Brzezinski, M. A. (1992). Cell-cycle effects on the kinetics of silicic acid uptake and resource competition among diatoms. J. Plankton Res. 14, 1511–1539. doi: 10.1093/plankt/14.11.1511
Catalanotti, C., Yang, W., Posewitz, M. C., and Grossman, A. R. (2013). Fermentation metabolism and its evolution in algae. Front. Plant Sci. 4:150. doi: 10.3389/fpls.2013.00150
Charlson, R. J., Lovelock, J. E., Andreae, M. O., and Warren, S. G. (1987). Oceanic phytoplankton, atmospheric sulphur, cloud albedo and climate. Nature 326, 655–661. doi: 10.1038/326655a0
Chauton, M. S., Winge, P., Brembu, T., Vadstein, O., and Bones, A. M. (2013). Gene regulation of carbon fixation, storage, and utilization in the diatom phaeodactylum tricornutum acclimated to light/dark cycles. Plant Physiol. 161, 1034–1048. doi: 10.1104/pp.112.206177
Collart, L., Jiang, D., and Halsey, K. H. (2023). The volatilome reveals microcystin concentration, microbial composition, and oxidative stress in a critical Oregon freshwater lake. mSystems 8, e00379–23. doi: 10.1128/msystems.00379-23
Colomb, A., Yassaa, N., Williams, J., Peeken, I., and Lochte, K. (2008). Screening volatile organic compounds (VOCs) emissions from five marine phytoplankton species by head space gas chromatography/mass spectrometry (HS-GC/MS). J. Environ. Monit. 10:325. doi: 10.1039/b715312k
Dani, K. G. S., Torzillo, G., Michelozzi, M., Baraldi, R., and Loreto, F. (2020). Isoprene emission in darkness by a facultative heterotrophic green alga. Front. Plant Sci. 11:598786. doi: 10.3389/fpls.2020.598786
Davie-Martin, C. L., Giovannoni, S. J., Behrenfeld, M. J., Penta, W. B., and Halsey, K. H. (2020). Seasonal and spatial variability in the biogenic production and consumption of volatile organic compounds (VOCs) by marine plankton in the north Atlantic ocean. Front. Mar. Sci. 7:611870. doi: 10.3389/fmars.2020.611870
Dawson, R., Crombie, A., Pichon, P., Steinke, M., McGenity, T., and Murrell, J. (2021). The microbiology of isoprene cycling in aquatic ecosystems. Aquat. Microb. Ecol. 87, 79–98. doi: 10.3354/ame01972
De Bruyn, W. J., Clark, C. D., and Senstad, M. (2020). Production of acetaldehyde from ethanol in coastal waters. Environ. Sci. Pollut. Res. 27, 12673–12682. doi: 10.1007/s11356-020-07880-8
Domingues, N., Matos, A. R., Marques Da Silva, J., and Cartaxana, P. (2012). Response of the diatom Phaeodactylum tricornutum to photooxidative stress resulting from high light exposure. PLoS ONE 7:e38162. doi: 10.1371/journal.pone.0038162
Dudareva, N., Klempien, A., Muhlemann, J. K., and Kaplan, I. (2013). Biosynthesis, function and metabolic engineering of plant volatile organic compounds. New Phytol. 198, 16–32. doi: 10.1111/nph.12145
Fall, R., and Benson, A. A. (1996). Leaf methanol—the simplest natural product from plants. Trends Plant Sci. 1, 296–301. doi: 10.1016/S1360-1385(96)88175-0
Field, C. B., Behrenfeld, M. J., Randerson, J. T., and Falkowski, P. (1998). Primary production of the biosphere: integrating terrestrial and oceanic components. Science 281, 237–240. doi: 10.1126/science.281.5374.237
Fink, P. (2007). Ecological functions of volatile organic compounds in aquatic systems. Mar. Freshw. Behav. Physiol. 40, 155–168. doi: 10.1080/10236240701602218
Fisher, C. L., Lane, P. D., Russell, M., Maddalena, R., and Lane, T. W. (2020). Low molecular weight volatile organic compounds indicate grazing by the marine rotifer Brachionus plicatilis on the microalgae microchloropsis salina. Metabolites 10:361. doi: 10.3390/metabo10090361
Forlani, G., Mantelli, M., and Nielsen, E. (1999). Biochemical evidence for multiple acetoin-forming enzymes in cultured plant cells - Part 1 in the series ‘'Acetoin Synthesis in Higher Plants‘'. 1. Phytochem. 50, 255–262. doi: 10.1016/S0031-9422(98)00550-0
Galí, M., Simó, R., Vila-Costa, M., Ruiz-González, C., Gasol, J. M., and Matrai, P. (2013). Diel patterns of oceanic dimethylsulfide (DMS) cycling: microbial and physical drivers. Glob. Biogeochem. Cycles 27, 620–636. doi: 10.1002/gbc.20047
Gu, Z., Eils, R., and Schlesner, M. (2016). Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32, 2847–2849. doi: 10.1093/bioinformatics/btw313
Guillard, R. R. L., and Ryther, J. H. (1962). Studies of marine planktonic diatoms: I. Cyclotella nana hustedt, and detonula confervacea (Cleve) gran. Can. J. Microbiol. 8, 229–239. doi: 10.1139/m62-029
Hackenberg, S. C., Andrews, S. J., Airs, R. L., Arnold, S. R., Bouman, H. A., Cummings, D., et al. (2017). Basin-scale observations of monoterpenes in the arctic and atlantic oceans. Environ. Sci. Technol. 51, 10449–10458. doi: 10.1021/acs.est.7b02240
Halsey, K. H., and Giovannoni, S. J. (2023). Biological controls on marine volatile organic compound emissions: a balancing act at the sea-air interface. Earth-Sci. Rev. 240:104360. doi: 10.1016/j.earscirev.2023.104360
Halsey, K. H., Giovannoni, S. J., Graus, M., Zhao, Y., Landry, Z., Thrash, J. C., et al. (2017). Biological cycling of volatile organic carbon by phytoplankton and bacterioplankton. Limnol. Oceanogr. 62, 2650–2661. doi: 10.1002/lno.10596
Halsey, K. H., O'Malley, R. T., Graff, J. R., Milligan, A. J., and Behrenfeld, M. J. (2013). A common partitioning strategy for photosynthetic products in evolutionarily distinct phytoplankton species. New Phytol. 198, 1030–1038. doi: 10.1111/nph.12209
Hansell, D. A. (2013). Recalcitrant dissolved organic carbon fractions. Ann. Rev. Mar. Sci. 5, 421–445. doi: 10.1146/annurev-marine-120710-100757
Hastings, J., De Matos, P., Dekker, A., Ennis, M., Harsha, B., Kale, N., et al. (2012). The ChEBI reference database and ontology for biologically relevant chemistry: enhancements for 2013. Nucleic Acids Res. 41, D456–D463. doi: 10.1093/nar/gks1146
He, H., Höper, R., Dodenhöft, M., Marlière, P., and Bar-Even, A. (2020). An optimized methanol assimilation pathway relying on promiscuous formaldehyde-condensing aldolases in E. coli. Metab. Eng. 60, 1–13. doi: 10.1016/j.ymben.2020.03.002
Holzinger, R., Acton, W. J. F., Bloss, W. J., Breitenlechner, M., Crilley, L. R., Dusanter, S., et al. (2019). Validity and limitations of simple reaction kinetics to calculate concentrations of organic compounds from ion counts in PTR-MS. Atmos. Meas. Tech. 12, 6193–6208. doi: 10.5194/amt-12-6193-2019
Huysman, M. J. J., Vyverman, W., and De Veylder, L. (2014). Molecular regulation of the diatom cell cycle. J. Exp. Botany 65, 2573–2584. doi: 10.1093/jxb/ert387
Jackson, R. L., Gabric, A. J., Cropp, R., and Woodhouse, M. T. (2020). Dimethylsulfide (DMS), marine biogenic aerosols and the ecophysiology of coral reefs. Biogeosciences 17, 2181–2204. doi: 10.5194/bg-17-2181-2020
Jallet, D., Caballero, M. A., Gallina, A. A., Youngblood, M., and Peers, G. (2016). Photosynthetic physiology and biomass partitioning in the model diatom Phaeodactylum tricornutum grown in a sinusoidal light regime. Algal Res. 18, 51–60. doi: 10.1016/j.algal.2016.05.014
Jardine, K. J., and McDowell, N. (2023). Fermentation-mediated growth, signaling, and defense in plants. New Phytol. 239, 839–851. doi: 10.1111/nph.19015
Kameyama, S., Tanimoto, H., Inomata, S., Suzuki, K., Komatsu, D. D., Hirota, A., Konno, U., and Tsunogai, U. (2011). Application of PTR-MS to an incubation experiment of the marine diatom Thalassiosira pseudonana. Geochem J. 45, 355–363. doi: 10.2343/geochemj.1.0127
Kameyama, S., Tanimoto, H., Inomata, S., Tsunogai, U., Ooki, A., Yokouchi, Y., et al. (2009). Equilibrator inlet-proton transfer reaction-mass spectrometry (EI-PTR-MS) for sensitive, high-resolution measurement of dimethyl sulfide dissolved in seawater. Anal. Chem. 81, 9021–9026. doi: 10.1021/ac901630h
Kanehisa, M. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457-D462. doi: 10.1093/nar/gkv1070
Kelsey, R. G., and Westlind, D. J. (2017). Physiological stress and ethanol accumulation in tree stems and woody tissues at sublethal temperatures from fire. BioScience 67, 443–451. doi: 10.1093/biosci/bix037
Kieber, R. J., Zhou, X., and Mopper, K. (1990). Formation of carbonyl compounds from UV-induced photodegradation of humic substances in natural waters: fate of riverine carbon in the sea. Limnol. Oceanogr. 35, 1503–1515. doi: 10.4319/lo.1990.35.7.1503
Kim, J., Brown, C. M., Kim, M. K., Burrows, E. H., Bach, S., Lun, D. S., et al. (2017). Effect of cell cycle arrest on intermediate metabolism in the marine diatom Phaeodactylum tricornutum. Proc. Natl. Acad. Sci. U.S.A. 114, E8007–E8016. doi: 10.1073/pnas.1711642114
Koeduka, T., Baiga, T. J., Noel, J. P., and Pichersky, E. (2009). Biosynthesis of t - anethole in anise: characterization of t - anol/isoeugenol synthase and an O - methyltransferase specific for a C7–C8 propenyl side chain. Plant Physiol. 149, 384–394. doi: 10.1104/pp.108.128066
Kolber, Z. S., Prášil, O., and Falkowski, P. G. (1998). Measurements of variable chlorophyll fluorescence using fast repetition rate techniques: defining methodology and experimental protocols. Biochim. Biophys. Acta - Bioenerg. 1367, 88–106. doi: 10.1016/S0005-2728(98)00135-2
Kroth, P. G., Chiovitti, A., Gruber, A., Martin-Jezequel, V., Mock, T., Parker, M. S., et al. (2008). A model for carbohydrate metabolism in the diatom Phaeodactylum tricornutum deduced from comparative whole genome analysis. PLoS ONE 3:e1426. doi: 10.1371/journal.pone.0001426
Ledford, S. M., and Meredith, L. K. (2024). Volatile organic compound metabolism on early earth. J. Mol. Evol. 92, 605–617. doi: 10.1007/s00239-024-10184-x
Lees, G. J., and Jago, G. R. (1978). Role of acetaldehyde in metabolism: a review 2. The metabolism of acetaldehyde in cultured dairy products. J. Dairy Sci. 61, 1216–1224. doi: 10.3168/jds.S0022-0302(78)83709-6
Lelieveld, J., Hoor, P., Jöckel, P., Pozzer, A., Hadjinicolaou, P., Cammas, J.-P., et al. (2009). Severe ozone air pollution in the Persian Gulf region. Atmos. Chem. Phys. 9, 1393–1406. doi: 10.5194/acp-9-1393-2009
Liebelt, D. J., Jordan, J. T., and Doherty, C. J. (2019). Only a matter of time: the impact of daily and seasonal rhythms on phytochemicals. Phytochem. Rev. 18, 1409–1433. doi: 10.1007/s11101-019-09617-z
Liu, M., and Matsui, H. (2022). Secondary organic aerosol formation regulates cloud condensation nuclei in the global remote troposphere. Geophys. Res. Lett. 49:e2022GL100543. doi: 10.1029/2022GL100543
Liu, Z., Wang, M., Wu, M., Li, X., Liu, H., Niu, N., et al. (2023). Volatile organic compounds (VOCs) from plants: from release to detection. Trends Analyt. Chem. 158:116872. doi: 10.1016/j.trac.2022.116872
Maeda, H., and Dudareva, N. (2012). The shikimate pathway and aromatic amino acid biosynthesis in plants. Annu. Rev. Plant Biol. 63, 73–105. doi: 10.1146/annurev-arplant-042811-105439
Masdeu-Navarro, M., Mangot, J.-F., Xue, L., Cabrera-Brufau, M., Kieber, D. J., Rodríguez-Ros, P., et al. (2024). Diel variation of seawater volatile organic compounds, DMSP-related compounds, and microbial plankton inside and outside a tropical coral reef ecosystem. Front. Mar. Sci. 11:1341619. doi: 10.3389/fmars.2024.1341619
Mayer, K. J., Wang, X., Santander, M. V., Mitts, B. A., Sauer, J. S., Sultana, C. M., et al. (2020). Secondary marine aerosol plays a dominant role over primary sea spray aerosol in cloud formation. ACS Cent. Sci. 6, 2259–2266. doi: 10.1021/acscentsci.0c00793
Mazeed, M., Singh, R., Kumar, P., Roy, A., Raman, B., Kruparani, S. P., et al. (2021). Recruitment of archaeal DTD is a key event toward the emergence of land plants. Sci. Adv. 7:eabe8890. doi: 10.1126/sciadv.abe8890
Mccormick, A. C., Gershenzon, J., and Unsicker, S. B. (2014). Little peaks with big effects: establishing the role of minor plant volatiles in plant-insect interactions. Plant Cell Environ. 37, 1836–1844. doi: 10.1111/pce.12357
McKay, W. A., Turner, M. F., Jones, B. M. R., and Halliwell, C. M. (1996). Emissions of hydrocarbons from marine phytoplankton-Some results from controlled laboratory experiments. Atmos. Environ. 30, 2583–2593. doi: 10.1016/1352-2310(95)00433-5
Mella-Flores, D., Six, C., Ratin, M., Partensky, F., Boutte, C., Le Corguillé, G., et al. (2012). Prochlorococcus and Synechococcus have evolved different adaptive mechanisms to cope with light and UV stress. Front. Microbiol. 3:285. doi: 10.3389/fmicb.2012.00285
Meredith, L. K., Ledford, S. M., Riemer, K., Geffre, P., Graves, K., Honeker, L. K., et al. (2023). Automating methods for estimating metabolite volatility. Front. Microbiol. 14:1267234. doi: 10.3389/fmicb.2023.1267234
Meskhidze, N., and Nenes, A. (2006). Phytoplankton and cloudiness in the Southern ocean. Science 314, 1419–1423. doi: 10.1126/science.1131779
Meskhidze, N., Sabolis, A., Reed, R., and Kamykowski, D. (2015). Quantifying environmental stress-induced emissions of algal isoprene and monoterpenes using laboratory measurements. Biogeosciences 12, 637–651. doi: 10.5194/bg-12-637-2015
Mincer, T. J., and Aicher, A. C. (2016). Methanol production by a broad phylogenetic array of marine phytoplankton. PLoS ONE 11:e0150820. doi: 10.1371/journal.pone.0150820
Misztal, P. K., Hewitt, C. N., Wildt, J., Blande, J. D., Eller, A. S. D., Fares, S., et al. (2015). Atmospheric benzenoid emissions from plants rival those from fossil fuels. Sci. Rep. 5:12064. doi: 10.1038/srep12064
Moore, E. R., Davie-Martin, C. L., Giovannoni, S. J., and Halsey, K. H. (2020). Pelagibacter metabolism of diatom-derived volatile organic compounds imposes an energetic tax on photosynthetic carbon fixation. Environ. Microbiol. 22, 1720–1733. doi: 10.1111/1462-2920.14861
Moore, E. R., Weaver, A. J., Davis, E. W., Giovannoni, S. J., and Halsey, K. H. (2022). Metabolism of key atmospheric volatile organic compounds by the marine heterotrophic bacterium Pelagibacter HTCC1062 (SAR11). Environ. Microbiol. 24, 212–222. doi: 10.1111/1462-2920.15837
Moran, M. A., Ferrer-González, F. X., Fu, H., Nowinski, B., Olofsson, M., Powers, M. A., et al. (2022). The ocean's labile DOC supply chain. Limnol. Oceanogr. 67, 1007–1021. doi: 10.1002/lno.12053
Ng, W. H. A., and Liu, H. (2016). Diel periodicity of grazing by heterotrophic nanoflagellates influenced by prey cell properties and intrinsic grazing rhythm. J. Plankton Res. 38, 636–651. doi: 10.1093/plankt/fbw014
Nuccio, J., Seaton, P. J., and Kieber, R. J. (1995). Biological production of form aldehyde in the marine environment. Limnol. Oceanogr. 40, 521–527. doi: 10.4319/lo.1995.40.3.0521
Orlova, I., Marshall-Colón, A., Schnepp, J., Wood, B., Varbanova, M., Fridman, E., et al. (2007). Reduction of benzenoid synthesis in petunia flowers reveals multiple pathways to benzoic acid and enhancement in auxin transport. Plant Cell 18, 3458–3475. doi: 10.1105/tpc.106.046227
Padaki, V. G., Mayali, X., Weber, P. K., Giovannoni, S. J., Abraham, K., Jacobs, K., et al. (2024). Bacterial volatile organic compound specialists in the phycosphere. BioRxiv. doi: 10.1101/2024.06.14.599129
Ploug, H., and Grossart, H. (2000). Bacterial growth and grazing on diatom aggregates: respiratory carbon turnover as a function of aggregate size and sinking velocity. Limnol. Oceanogr. 45, 1467–1475. doi: 10.4319/lo.2000.45.7.1467
Pollastri, S., Baccelli, I., and Loreto, F. (2021). Isoprene: an antioxidant itself or a molecule with multiple regulatory functions in plants? Antioxidants 10:684. doi: 10.3390/antiox10050684
Pozzer, A. C., Gómez, P. A., and Weiss, J. (2022). Volatile organic compounds in aquatic ecosystems - detection, origin, significance and applications. Sci. Total Environ. 838:156155. doi: 10.1016/j.scitotenv.2022.156155
Prézelin, B. B. (1992). “Diel periodicity in phytoplankton productivity,” in The Daily Growth Cycle of Phytoplankton, eds. T. Berman, H. J. Gons, and L. R. Mur (Dordrecht: Springer Netherlands), 1–35. doi: 10.1007/978-94-011-2805-6_1
Quinn, P. K., and Bates, T. S. (2011). The case against climate regulation via oceanic phytoplankton sulphur emissions. Nature 480, 51–56. doi: 10.1038/nature10580
Ratte, M., Plass-Dülmer, C., Koppmann, R., Rudolph, J., and Denga, J. (1993). Production mechanism of C2-C4 hydrocarbons in seawater: field measurements and experiments. Glob. Biogeochem. Cycles 7, 369–378. doi: 10.1029/93GB00054
Rico, M., Santiago-Díaz, P., Rivero, A., and Santana-Casiano, J. M. (2024). Characterization of polyphenols and carbohydrates exuded by Phaeodactylum tricornutum diatom grown under Cu stress. Sci. Rep. 14:9367. doi: 10.1038/s41598-024-60252-x
Ritchie, R. J. (2006). Consistent Sets of spectrophotometric chlorophyll equations for acetone, methanol and ethanol solvents. Photosynth. Res. 89, 27–41. doi: 10.1007/s11120-006-9065-9
Rocco, M., Dunne, E., Peltola, M., Barr, N., Williams, J., Colomb, A., et al. (2021). Oceanic phytoplankton are a potentially important source of benzenoids to the remote marine atmosphere. Commun. Earth Environ. 2:175. doi: 10.1038/s43247-021-00253-0
Ruiz-Halpern, S., Echeveste, P., AgustÃ, S., and Duarte, C. M. (2014). Size-dependence of volatile and semi-volatile organic carbon content in phytoplankton cells. Front. Mar. Sci. 1:25. doi: 10.3389/fmars.2014.00025
Sargeant, S. L., Murrell, J. C., Nightingale, P. D., and Dixon, J. L. (2018). Basin-scale variability of microbial methanol uptake in the Atlantic ocean. Biogeosciences 15, 5155–5167. doi: 10.5194/bg-15-5155-2018
Sartori, R. B., Siqueira, S. F., Maroneze, M. M., Fagundes, M. B., Wagner, R., Zepka, L. Q., et al. (2021). Microalgal secondary metabolites: effect of climatic variables, seasons, and photocycles on the biogeneration of volatile organic compounds (VOCs). J. Appl. Phycol. 33, 1457–1472. doi: 10.1007/s10811-021-02391-6
Shaw, S. L., Chisholm, S. W., and Prinn, R. G. (2003). Isoprene production by Prochlorococcus, a marine cyanobacterium, and other phytoplankton. Mar. Chem. 80, 227–245 doi: 10.1016/S0304-4203(02)00101-9
Shaw, S. L., Gantt, B., and Meskhidze, N. (2010). Production and emissions of marine isoprene and monoterpenes: a review. Adv. Meteorol. 2010:408696. doi: 10.1155/2010/408696
Shemi, A., Ben-Dor, S., Rotkopf, R., Dym, O., and Vardi, A. (2023). Phylogeny and biogeography of the algal DMS-releasing enzyme in the global ocean. ISME Commun. 3:72. doi: 10.1038/s43705-023-00280-2
Shinomura, Y., Iwata, T., and Suzuki, Y. (2005). Diel changes in dissolved organic carbon in the upper layer of Suruga Bay, Japan. Estuar. Coast. Shelf Sci. 62, 699–709. doi: 10.1016/j.ecss.2004.10.006
Simó, R., Cortés-Greus, P., Rodríguez-Ros, P., and Masdeu-Navarro, M. (2022). Substantial loss of isoprene in the surface ocean due to chemical and biological consumption. Commun. Earth Environ. 3:20. doi: 10.1038/s43247-022-00352-6
Smith, S. R., Gillard, J. T. F., Kustka, A. B., McCrow, J. P., Badger, J. H., Zheng, H., et al. (2016). Transcriptional orchestration of the global cellular response of a model pennate diatom to diel light cycling under iron limitation. PLoS Genet. 12:e1006490. doi: 10.1371/journal.pgen.1006490
Steinke, M., Hodapp, B., Subhan, R., Bell, T. G., and Martin-Creuzburg, D. (2018). Flux of the biogenic volatiles isoprene and dimethyl sulfide from an oligotrophic lake. Sci. Rep. 8:630. doi: 10.1038/s41598-017-18923-5
Su, Y., Hu, J., Xia, M., Chen, J., Meng, W., Qian, C., et al. (2024). An undiscovered circadian clock to regulate phytoplankton photosynthesis. PNAS 3:497. doi: 10.1093/pnasnexus/pgae497
Suganuma, N., Shimokawa, R., Katoh, T., and Nagai, T. (1993). Presence of acetone and acetoacetate decarboxylase in soybean root nodules. J. Soil Sci. Plant Nutr. 39, 653–660. doi: 10.1080/00380768.1993.10419182
Sun, J., Steindler, L., Thrash, J. C., Halsey, K. H., Smith, D. P., Carter, A. E., et al. (2011). One carbon metabolism in SAR11 pelagic marine bacteria. PLoS ONE 6:e23973. doi: 10.1371/journal.pone.0023973
Sunda, W., Kieber, D. J., Kiene, R. P., and Huntsman, S. (2002). An antioxidant function for DMSP and DMS in marine algae. Nature 418, 317–320. doi: 10.1038/nature00851
Sunda, W. G., Hardison, R., Kiene, R. P., Bucciarelli, E., and Harada, H. (2007). The effect of nitrogen limitation on cellular DMSP and DMS release in marine phytoplankton: climate feedback implications. Aquat. Sci. 69, 341–351. doi: 10.1007/s00027-007-0887-0
Uchimiya, M., Schroer, W., Olofsson, M., Edison, A. S., and Moran, M. A. (2022). Diel investments in metabolite production and consumption in a model microbial system. ISME J. 16, 1306–1317. doi: 10.1038/s41396-021-01172-w
Van Lis, R., Popek, M., Couté, Y., Kosta, A., Drapier, D., Nitschke, W., et al. (2017). Concerted up-regulation of aldehyde/alcohol dehydrogenase (ADHE) and starch in Chlamydomonas reinhardtii increases survival under dark anoxia. J. Biol. Chem. 292, 2395–2410. doi: 10.1074/jbc.M116.766048
Van Rijssel, M., and Gieskes, W. W. C. (2002). Temperature, light, and the dimethylsulfoniopropionate (DMSP) content of Emiliania huxleyi (Prymnesiophyceae). J. Sea Res. 48, 17–27. doi: 10.1016/S1385-1101(02)00134-X
Vivaldo, G., Masi, E., Taiti, C., Caldarelli, G., and Mancuso, S. (2017). The network of plants volatile organic compounds. Sci. Rep. 7:11050. doi: 10.1038/s41598-017-10975-x
Wang, Y., Zhao, Y., Bollas, A., Wang, Y., and Au, K. F. (2021). Nanopore sequencing technology, bioinformatics and applications. Nat. Biotechnol. 39, 1348–1365. doi: 10.1038/s41587-021-01108-x
Xiao, Z., and Lu, J. R. (2014). Generation of acetoin and its derivatives in foods. J. Agric. Food Chem. 62, 6487–6497. doi: 10.1021/jf5013902
Yáñez-Serrano, A. M., Filella, I., LLusià, J., Gargallo-Garriga, A., Granda, V., Bourtsoukidis, E., et al. (2021). GLOVOCS - master compound assignment guide for proton transfer reaction mass spectrometry users. Atmosph. Environ. 244:117929. doi: 10.1016/j.atmosenv.2020.117929
Yu, Z., and Li, Y. (2021). Marine volatile organic compounds and their impacts on marine aerosol-A review. Sci. Total Environ. 768:145054. doi: 10.1016/j.scitotenv.2021.145054
Zhao, D., Li, R., Fang, G., Tham, Y. J., Zou, S., Ling, Z., et al. (2024). Isoprene emission dynamics in chaetoceros curvisetus: insights from transcriptome analysis under light/dark cycles. Algal Res. 82:103641. doi: 10.1016/j.algal.2024.103641
Zuo, Z. (2019). Why algae release volatile organic compounds-the emission and roles. Front. Microbiol. 10:491. doi: 10.3389/fmicb.2019.00491
Keywords: volatile organic compounds, diel VOC metabolism, phytoplankton physiology, cell cycle, central carbon metabolism
Citation: Padaki VG, Palmer E, Jiang Y, Buchholz HH, Kimbrel JA and Halsey KH (2025) Diatom volatile organic compound production is driven by diel metabolism and the cell cycle. Front. Microbiol. 16:1620542. doi: 10.3389/fmicb.2025.1620542
Received: 29 April 2025; Accepted: 15 September 2025;
Published: 06 October 2025.
Edited by:
Taichi Yokokawa, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), JapanReviewed by:
Sohiko Kameyama, Hokkaido University, JapanSanjeeb Kumar Mandal, Chaitanya Bharathi Institute of Technology, India
Copyright © 2025 Padaki, Palmer, Jiang, Buchholz, Kimbrel and Halsey. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kimberly H. Halsey, a2ltYmVybHkuaGFsc2V5QG9yZWdvbnN0YXRlLmVkdQ==