DOM was collected in the Altamaha-Doboy-Sapelo Estuary off the southeastern U.S. coast in November 2012 and May 2013. Samples were collected in triplicate at three locations: the Altamaha River (salinity S = 0 to 11; strong terrigenous inputs), the mouth of Doboy Sound (S = 22.9 to 34.6; large marine influence) and near the head of Sapelo Sound (S = 9.1 to 31.5; mixture of terrigenous, marine and salt marsh inputs; Medeiros et al., 2015a) (Figure S1). All samples were collected during high tide. Immediately after collection, triplicate water samples from each site were filtered sequentially through Whatman GF/D filters (pre-combusted at 450°C for 5 h; nominal 2.7 μm pore size) and 0.2 μm Pall Supor membrane filters into pre-combusted amber glass bottles. Fifty milliliters aliquots of the filtrate from each sampling site were set aside in separated beakers for the preparation of microbial inocula. For that, the 0.2 μm filters from all locations were aseptically cut into pieces and equal areas of the filters were pooled into 50 mL aliquots of 0.2 μm filtrate from each of the three sites and stirred for 30 min. The resulting filtrates, now containing microbes from all sites, were added back to the appropriate set of bottles. This ensured that functional capabilities of the microbes were similar during all incubations and changes in DOM composition could therefore be attributed predominantly to the different initial composition of the DOM pools. We note that fine particles that passed through the GF/D filter but were retained in the 0.2 μm membrane were likely also introduced to the bottles. After the addition of inorganic nutrients (20 μM Na₂PO₄ and 50 μM NH₄Cl), a triplicate set from each sampling site was immediately filtered (0.2 μm) in order to characterize the initial condition for each experiment. Aliquots were stored frozen for DOC analysis, and the remaining filtrates were acidified to pH 2 (using HCl) and DOM was extracted using solid phase extraction (SPE) cartridges (Agilent Bond Elut PPL) as in Dittmar et al. (2008). We refer to those samples as T0.

The remaining triplicate sets were incubated in the dark at the temperature measured at the time of collection for up to 70 days. At days 35 and 70, samples were filtered, acidified, and DOM was extracted using PPL cartridges as described previously. These samples are referred to as T35 and T70. Shorter-term incubations lasting 24 h were also conducted with water collected at the mouth of Doboy Sound in February and in December of 2011 (referred to as T1). Aliquots for DOC analysis were collected for all incubations as described above (i.e., prior to acidification, stored frozen).

In all cases, DOC concentration in water samples and extracts (i.e., dried and redissolved in ultrapure water) was measured with a Shimadzu TOC-V_CPH analyzer. Analytical accuracy and precision were tested against the Consensus Reference Material (Hansell, 2005) and were better than 5%. SPE extraction efficiency across all samples was 75 ± 5% of the DOC. Bulk δ¹³C of extracted DOC (SPE-DOC) was measured with a Finnigan MAT 251 isotope ratio mass spectrometer after complete drying. Procedural blanks did not yield detectable amounts of carbon isotope contamination. The molecular composition of DOM extracts (15 mg C L⁻¹ in 1:1 methanol/water) was analyzed using a 15T Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS; Bruker Daltonics) with electrospray ionization (negative mode) as in Seidel et al. (2014). Spectra were internally calibrated with >100 known C_xH_yO_z molecular formulae over the mass range in the samples. With this calibration procedure, a mass error of <0.1 ppm was achieved. Before each sample set, blank checks with methanol and ultrapure water were measured. Molecular formulae were calculated in the mass range between 150 and 850 Da by applying the following restrictions: ¹²C₁₋₁₃₀¹H₁₋₂₀₀O₁₋₁₅₀¹⁴N₀₋₄S₀₋₂P₀₋₂. Assignment of molecular formulae was done considering a maximum mass error of 0.5 ppm and using the criteria described by Rossel et al. (2013). Only compounds with a signal-to-noise ratio of 4 or higher were used in the analysis to eliminate inter-sample variability based on peaks that were close to the limit of detection. Additional details of molecular formulae assignment are given in Seidel et al. (2014). A sample from the North Equatorial Pacific Intermediate Water collected at station NELHA off Big Island in Hawaii at 674 m water depth (Osterholz et al., 2015; Riedel et al., 2016) was also analyzed repeatedly (41 times) in order to characterize differences in DOM composition as identified by FT-ICR MS that arise solely due to instrument variability.

Patterns of variability in FT-ICR MS-derived DOM composition were distinguished by principal component (PC) analysis using all replicates simultaneously as in Medeiros et al. (2015a,b,c). All PCs shown here are significantly different (95% confidence level) from results of PC analysis of spatially and temporally uncorrelated random processes (Overland and Preisendorfer, 1982). The Wilcoxon rank-sum test was used for comparisons between groups of compounds, as in Osterholz et al. (2016).

A detailed characterization of spatial and temporal differences in DOM composition in this system has been presented in Medeiros et al. (2015a). We briefly summarize results here for comparison with microbially-mediated transformations described in the next section. DOC concentrations were substantially lower near the ocean at Doboy Sound than at the other sites (Table S1) and increased at all locations in May 2013 when river discharge was high. The pre-incubation samples (T0; solid symbols) were spread along the first PC (PC1) axis (Figure 1A), which separated samples with high salinity (e.g., Doboy in November 2012) from those with low salinity (e.g., Altamaha in May 2013) (Table S1) (see caption of Figure 1 for example of interpretation of the analysis). The correlation coefficient between PC1 and salinity for the pre-incubation samples was −0.87 (p < 0.01). The loading of PC1 (Figure 1B) shows a pattern typical of river-to-ocean transects (Medeiros et al., 2015c), with high salinity samples enriched with compounds with high H/C ratios compared to low-salinity samples, which is consistent with the higher abundance of aliphatic structures in marine DOM (Sleighter and Hatcher, 2008; Osterholz et al., 2016). This is also consistent with the strong gradient in DOM sources in the system associated with terrigenous vs. oceanic inputs (Medeiros et al., 2015a).

For any given location, samples from May 2013 had higher DOC concentrations and DOM with a more terrigenous character than samples from November 2012 (e.g., samples from Altamaha River in May 2013 had higher PC1 scores than samples from Altamaha River in November 2012; the same is true for samples from Sapelo to Doboy Sounds; Figure 1A), a pattern consistent with their lighter δ¹³C values (Table S1). This difference probably resulted from a severe drought in fall 2012, when the Altamaha River discharge was anomalously low (Medeiros et al., 2015a). During that time, salinity at all sampling sites was higher than in May 2013 (Table S1), which is consistent with the more marine character of the DOM observed in November 2012.

In addition to capturing the gradient in DOM composition in the estuary (section Spatial and Temporal Variations in DOM Composition), the PC analysis shown in Figure 1 also captured microbially-mediated transformations of the DOM. The largest transformations in DOM as captured by the first 2 principal components [i.e., the largest distances in PC space between T0 (solid symbols) and T70 (open symbols) for a given incubation in Figure 1A] occurred in samples with a more marine character (i.e., Doboy Sound in November 2012, black circles), decreasing progressively in samples with a more terrigenous signature. Incubations of Altamaha River water in May 2013 (green squares), when discharge was near the seasonal maximum (Medeiros et al., 2015a), yielded no significant changes in DOM composition as detected by the first 2 PCs (a different pattern of DOM transformation was observed in that case; see Figure S2).

It is interesting to note that microbially-mediated transformations resulted in changes on both the first and second PCs. The PC1 scores at the end of the incubations (T70) were generally higher than those at the start of the incubations (T0). Since large PC1 scores are associated with terrigenous DOM (see section Spatial and Temporal Variations in DOM Composition), this indicates that microbial degradation resulted in DOM with a more terrigenous character. This presumably resulted from preferential degradation of marine DOM, which is indeed thought to be more labile than freshwater DOM (D'Andrilli et al., 2015). As marine DOM (i.e., molecular formulae shown in blue in Figure 1B) was degraded during the experiment, terrigenous DOM (i.e., molecular formulae shown in red in Figure 1B) became relatively more important compared to the initial state. The shift toward a terrigenous signature (i.e., a shift toward high values along the PC1 axis) was stronger in the samples that had a stronger marine signature initially (Figure 1A). For the incubations of Doboy Sound water in November 2012, in particular, PC1 changed substantially during the incubation, representing 53% of the distance along the PC1 axis between marine (solid black circle, PC1 = −2.45 ± 0.03) and terrigenous (solid green square, PC1 = +0.86 ± 0.05) DOM. In other words, PC1 (Figures 1A,B) captured the fraction of the DOM that is statistically different between the sites and is related to the gradient in DOM composition due to terrigenous vs. marine inputs (as shown in section Spatial and Temporal Variations in DOM Composition). Biodegradation during the incubation of Doboy Sound water from November 2012 resulted in large enough changes in that fraction captured by PC1 that the DOM composition as identified by FT-ICR MS at the end of the experiment (open black circles in Figure 1A) was nearly as similar to riverine DOM (solid green squares) as it was to the original pre-incubated marine DOM (solid black circles). We note that photo-degradation generally produces the opposite effect, with riverine DOM resembling marine DOM in its broad molecular composition after extensive photo-degradation mainly due to removal of aromatics (Riedel et al., 2016).

Pre- and post-incubation samples were also characterized by different PC2 scores (Figure 1A). Microbially-mediated transformations therefore also resulted in a shift in DOM composition that did not fall along the terrigenous-marine gradient captured by PC1. The loading of PC2 reveals that molecules characterized by low O/C and high H/C ratios (molecules in red in Figure 1C) were relatively enriched at the end of the incubations. A detailed investigation of those changes is presented in the next section.

These analyses show that microbial degradation resulted in shifts in DOM composition so that the pool as a whole lost some of its marine characteristics and became more similar to terrigenous DOM. However, the analyses included samples from different locations and times. To examine the specific chemical formula changes occurring during biodegradation, the PC analysis was repeated considering only incubations of DOM from Doboy Sound in November 2012. This removed effects of spatial and temporal variability in DOM composition from the analysis, so that transformations could be primarily attributed to biodegradation.

DOC concentrations decreased progressively during the course of the incubation (Table S1), indicating remineralization of a fraction of the DOM pool. The score of the first principal component (referred to as PC1_DOB−long) increased with time (Figure 2A), and compounds associated with molecular formulae shown in blue in Figure 2B (negative loading) decreased in relative abundance during microbial processing while those associated with formulae shown in red (positive loading) increased with processing. As such, compounds that decreased or increased in relative abundance had distinct chemical characteristics. There was a relative enrichment of molecules that occupy the upper left corner of the van Krevelen diagram after the incubation (Figure 2B), i.e., a shift toward low O/C and high H/C ratios, which is consistent with molecular transformations observed in dark incubations of DOM from the Amazon River mouth (Medeiros et al., 2015c; Seidel et al., 2015). Additionally, while formulae whose relative abundances strongly decreased after the incubation (i.e., more negative loadings of PC1_DOB−long in Figure 2B) were characterized by high masses (553.5 ± 78.7 Da), formulae whose relative abundances were enriched after the incubation (i.e., more positive PC1_DOB−long loading in Figure 2B) were characterized by significantly lower masses (286.0 ± 93.9 Da) (Figure 3A and Table 1; Wilcoxon rank-sum test, p < 0.01). Molecular formulae associated with compounds that were preferentially depleted during the incubation were also characterized by a significantly higher numbers of carbon, hydrogen and oxygen atoms compared to compounds whose relative abundance enriched due to microbially-mediated transformations in DOM composition (Figures 3B–D; Wilcoxon rank-sum test, p < 0.01). They were also characterized by lower aromaticity indices (AI^*; a measure of the aromaticity of the molecules; Koch and Dittmar, 2006, 2016; Figure 3H), although the differences are small (but still significant; Wilcoxon rank-sum test, p < 0.01). Lastly, the transformations in DOM composition resulting from biodegradation (Figure 2B) were compared with the DOM composition variability associated with differences in terrigenous and marine inputs (Figure 1B) on a molecular basis. Molecular formulae associated with compounds whose relative abundances decreased after the microbial transformation (blue dots in Figure 2B) were generally also associated with compounds enriched in marine DOM (blue dots in Figure 1B), although there is substantial scatter around the bin averages (Figure 3E). Molecular formulae associated with compounds whose relative abundances increased after the incubation, on the other hand, were associated with compounds enriched in terrigenous DOM, with much smaller scatter around the bin averages. Thus, the marine DOM appeared to be preferentially transformed during the incubations, leaving behind relatively more terrigenous DOM, which is consistent with results shown in Figure 1. The large scatter around the formulae associated with the compounds depleted during the incubation (Figure 3E) indicated that, although there is a preference for biotransformation of marine DOM, a fraction of the compounds enriched in terrigenous DOM were also transformed during the incubation.

Repeating the PC analysis using only samples from days 0 to 35 (T0 and T35), T0 and T70, or T35 and T70 led to results similar to those obtained using all time points together. This indicates that the dominant pattern of DOM compositional change associated with biodegradation remained similar during the course of the incubation, at least at the monthly time scale captured by the sampling. Moreover, when the T0, T35 and T70 samples from Doboy Sound were analyzed together, approximately 60% of the total change in PC1_DOB−long occurred in the first 35 days (i.e., from −1.29 ± 0.13 to 0.15 ± 0.27 along PC1_DOB−long axis; Figure 2A), with the remaining 40% of the change (i.e., from 0.15 ± 0.27 to 1.10 ± 0.09; Figure 2A) occurring in the second half of the incubation. This indicates that active DOM transformations occurred throughout the experiment.

We have also investigated changes in DOM composition during the 24-h long incubations (T0-T1; Table 1). At those short time scales, the overall change in DOC concentrations were quite small (Table S1). Molecular formulae which were enriched after the 1-day experiments had chemical characteristics that were similar to those whose relative abundances decreased, occupying similar regions in van Krevelen space (Figure 4). For example, molecular formulae that were preferentially depleted after the 1-day incubations had masses of 485–500 Da, while formulae that were preferentially enriched had masses of 460–475 Da (Table 1; those differences were statistically significant; Wilcoxon rank-sum test, p < 0.01). Thus, the difference in the average mass of molecular formulae whose relative abundance decreased vs. increased during the short incubations (~Δ mass = 25–30 Da) was substantially smaller than the difference observed in the incubation lasting 70 days (~Δ mass = 270 Da; Table 1). Likewise, the differences in the number of carbon and oxygen atoms between compounds whose relative abundances decreased vs. increased in the 1-day incubations were also small (23 vs. 20–22 carbon atoms, respectively; unaltered for oxygen atoms; Table 1).

We also pursued a separate PC analysis for the other sets of replicates from each long incubation (e.g., from Sapelo Sound in May 2013, from Altamaha River in November 2012), similarly to the analysis presented in Figures 2, 3. The analysis revealed a progressive shift in the attributes of transformed DOM (Figure 5). In the 70-day incubation of DOM from the Altamaha River in May 2013 (i.e., water rich in terrigenous material; Figure 1), the compounds whose relative abundance enriched and decreased during the incubation had similar molecular masses (Figure 5). In that regard, results were qualitatively similar to the 1-day incubation (T0-T1) of the more marine-influenced Doboy Sound DOM (Table 1). As the marine character of the sample increased (as identified in Figure 1A), the mass difference between decreasing vs. increasing compounds also increased (Figure 5). For the incubation of Doboy Sound water from November 2012, increasing the duration of the experiment from 35 to 70 days led to a slight (but statistically significant) increase in the mass difference between compounds depleted vs. enriched.

Since the average mass difference between compounds whose relative abundance significantly increased vs. decreased in the short incubation of Doboy Sound water is small (average of 20–30 Da; Table 1) and they often occupy similar regions in van Krevelen space (Figure 4), it is possible that those differences are reflecting transformations of functional groups. Side chain oxidation (e.g., Filley et al., 2002), demethylation, or hydroxylation (e.g., Crawford and Crawford, 1980) of compounds present in the initial DOM pool (i.e., at T0) would lead to a reduction in the abundance of those compounds at the end of the incubation and to an increase in the abundance of compounds with slightly lower molecular masses. For example, vanillic acid (~168 Da) can be converted via demethylation to protocatechuic acid (~154 Da) (Merkens et al., 2005). In this case, the difference in the mass of the compounds that were consumed and produced is small at 14.0157 Da (-CH₂). If multiple progressive transformations occur, however, then the difference in the mass of the compounds consumed and produced will be larger (Bugg et al., 2011). As an example, the lignin model biphenyl compound 5,5′-dehydrodivanillate (~332 Da) can be transformed by bacterial degradation via an O-demethylation followed by a subsequent ring cleavage (Masai et al., 2007). If the resulting product is further hydrolyzed, decarboxylated, and then demethylated, protocatechuic acid is generated. In this case, the change in the mass between compounds consumed and produced progressively increases with each transformation (i.e., O-demethylation, ring cleavage, etc), reaching ~178 Da for the end result (i.e., comparing 5,5′-dehydrodivanillate and protocatechuic acid). Although, those examples come from investigations of microbial degradation of terrigenous material in soils, they provide illustrations of how biodegradation of DOM may result in small shifts in masses initially, and to large shifts in masses later on associated with additional, progressive transformations of the intermediate products. This scenario is consistent with the progression observed for microbially-mediated transformations in DOM composition in the incubations shown here. For Doboy Sound water, for example, the short term incubations yielded small differences in the average mass of the compounds whose relative abundances significantly decreased vs. increased (Table 1), which could be associated with the removal of a functional group (e.g., demethylation). If those resulting products are further transformed progressively via modification of additional functional groups (e.g., demethylation followed by hydroxylation and decarboxylation), then a large difference in the average mass of the compounds whose relative abundances significantly decreased vs. increased would be expected, which is consistent with results from the long-term incubations (Table 1). Alternatively, the large difference in masses between consumed and produced compounds in the long-term incubation could also be a result of different reactions. It will be interesting to see if emerging tools in analytical chemistry, microbiology and informatics (Moran et al., 2016) will allow for future studies to identify the specific pathways that occur in the transformation of DOM.

In the previous section, we hypothesized that progressive modifications of functional groups could explain why the differences in the masses of compounds whose relative abundance decreased vs. increased was small in short-term incubations, and comparatively larger in long-term experiments. Here, we compare the elemental composition of the compounds that changed significantly during the incubations to identify if those changes are consistent with that hypothesis.

The molecular formulae associated with the compounds whose relative abundance changed during the 1-day incubations have masses that are similar to each other. This does not necessarily mean, however, that they can be related to each other via transformations of side functional groups. Compounds with similar masses can have very different elemental compositions, and FT-ICR MS analysis does not give direct information about structure (e.g., Longnecker et al., 2015). We therefore investigated the possibility of transformations of functional groups using a strictly probabilistic approach. The analysis began by first selecting all molecular formulae associated with compounds whose relative abundance significantly increased after a given incubation. For the 1-day incubation (triplicate) from February 2011, for example, there were 496 such molecular formulae. For each of these, we computed the difference between its elemental composition and the elemental composition of all molecular formulae associated with compounds whose relative abundance decreased during the incubation (there were 408 such formulae in the February 2011 incubation). Therefore, the February 2011 calculation yielded a total of 202,368 (496 × 408) possible composition differences of potential substrates and products. Lastly, we sorted those composition differences to identify how often a potential change in elemental composition appeared (see Figure S3 for an example of the calculation). The rationale behind this calculation is that, if the microbial community is transforming the DOM via modification of a specific functional group, then it is more likely that a difference in elemental composition between compounds whose relative abundance decreased vs. increased equivalent to that modification will appear more often in the difference matrix. For example, if the microbial community is actively demethylating (i.e., targeting –CH₃), then it is more likely that a difference in elemental composition of a CH₂ (e.g., replacement of a CH₃ group by an H) between compounds whose relative abundance changed will be common. If, on the other hand, the differences in molecular formulae between relatively depleted and enriched compounds appear random, then it seems less likely that the microbial community is biasing their activity toward specific molecular modifications. If the most frequent composition differences in the 1-day incubations involve additions and removals of large numbers of carbon and heteroatoms, then the similarity in masses (Table 1) is likely simply coincidental and unlikely to have resulted from direct transformation of side functional groups. For example, although C₁₆H₆O₉ and C₁₃H₁₃O₇NS have a small mass difference of 14.9599 Da, the former cannot be converted into the latter via simple transformations of side functional groups given the large difference in their elemental composition.

The calculations revealed that, for the 1-day incubations, the most common differences in elemental composition were those involving a small number of C, H, and/or O atoms (e.g., removal or addition of a CH₂; removal of an O) (Table S2) and a small mass (Figure S4). This indicates that compounds whose relative abundances were significantly decreased vs. increased in the short incubations not only had similar masses (Table 1) and occupied similar regions in van Krevelen space (Figure 4), but many of them also had similar elemental compositions (i.e., small differences in the number of atoms in Table S2). Of the 200 most frequently observed composition differences shown in Table S2, 82 appeared in both the February and the December 2011 short incubations. If the analysis is expanded to include the 800 most frequently observed composition differences, then 268 appear in both incubations. Some of the most common composition differences that appeared in both the February and the December 2011 short incubations are listed in Table 2. For the 70-day incubation, on the other hand, the most frequent composition differences involve a net removal of a much larger number of C, H, and O (e.g., net removal of C₁₆H₈O₁₅, Table S2) and a larger mass (Figure S4). These could be related to changes in functional groups occurring sequentially (e.g., demethylation followed by hydroxylation and decarboxylation, whose cumulative result would be a large change in elemental composition) and/or other reactions. There is no overlap between the composition differences appearing in the short and long incubations. It is likely that similar changes do occur in the long incubations, but because of their cumulative effect (i.e., changes occurring on intermediate products) the net change in elemental composition is different than in the short incubations. These results are consistent with our previous hypothesis that changes in DOM composition in the 1-day incubations may be related to transformations in side functional groups (with a small change in elemental composition and mass), while changes in the long incubations may be related to progressive transformations (and thus with a larger change in elemental composition and mass). We note that previous incubations of coastal water lasting for up to 9 days have not revealed a significant impact of incubation time on DOM composition (Kujawinski et al., 2016). Our analysis of short- (1 day) and long-term (35 or 70 days) incubations indicate that significant differences can occur over these longer time scales.

The patterns just described for short- and long-term incubations are not due to instrument variability or noise. First, the principal component analysis separated the pre-incubation samples from those collected at the end of the experiments in both cases (Figures 2A, 4A,C). If the transformation in DOM composition shown in the van Krevelen diagrams (Figures 2B, 4B,D) were simply due to instrument noise, there would be no reason to expect separation between pre- and post-incubation samples in the principal component analysis. Second, we repeated the analysis using all 41 injections of the sample from the North Equatorial Pacific Intermediate Water collected at station NELHA off Hawaii. The differences in elemental composition between those samples (which in this case were entirely due to instrument variability) were mostly random, with small and large changes in molecular mass occurring frequently.

It is important to note some of the limitations of the analysis. The calculation provides a check for consistency only, since it is possible that even when the difference in elemental composition is small (small number of C, H, and/or O atoms), the structure of the compounds that underlies the predicted formula may not allow for that particular transformation. In addition, the analytical methods used here are focused on the investigation of compounds with molecular mass larger than about 150–160 Da. As reviewed by Kujawinski (2011), previous studies have shown that small molecules that are ubiquitous in cellular biomass can be easily assimilated into cells through dedicated transport systems. Since those compounds fall outside our analytical window, however, the transformations in DOM composition reported here are only representative of the fraction in the 150–850 Da range that can be detected by FT-ICR MS analysis. The analyses are also only representative of the fraction of the DOM that is retained in the solid phase extraction resin. Saccharides, for example, are known to be consumed by microorganisms (e.g., Rich et al., 1996), but are generally not well captured by the resin used here. As such, some of the DOC decrease observed in the incubations is likely related to a fraction of the DOM that is not captured by the FT-ICR MS analysis. The statistically significant transformations in composition reported in the short-term incubations reveal, however, that at least a portion of the labile DOM was captured by the solid phase extraction. Lastly, we note that we have focused here on transformations of estuarine DOM due to microbial degradation. In situ transformations are influenced by additional processes, including photochemistry and flocculation.

In summary, our analyses of microbial-mediated changes in DOM composition in an estuarine system characterized by strong terrigenous and marine inputs identified significant transformations. Analyses at the molecular level revealed a general preference for degradation of compounds that are enriched in marine DOM, which resulted in the remaining DOM becoming relatively more terrigenous. Transformations observed during short-term incubations resulted in a small difference in the average mass of compounds whose relative abundance decreased vs. increased, and a novel probabilistic approach revealed that changes in elemental composition were consistent with modifications of functional groups. Over longer time scales, the differences in mass were significantly larger, possibly a result of additional changes in function groups of intermediate products. Despite limitations, and given the challenges of actually identifying transformations in the complex DOM mixture, the analysis presented here serves as a step toward shedding light into these elusive processes. Future studies combining this probabilistic approach, NMR analysis and other techniques that provide direct information about structure, and biological gene expression at a highly-resolved scale will be particularly useful to tackle the biological complexity of carbon turnover.

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.