AbstractThe functioning of microbial ecosystems has important consequences from global climate to human health, but quantitative mechanistic understanding remains elusive. The components of microbial ecosystems can now be observed at high resolution, but interactions still have to be inferred e.g., a time-series may show a bloom of bacteria X followed by virus Y suggesting they interact. Existing inference approaches are mostly empirical, like correlation networks, which are not mechanistically constrained and do not provide quantitative mass fluxes, and thus have limited utility. We developed an inference method, where a mechanistic model with hundreds of species and thousands of parameters is calibrated to time series data. The large scale, nonlinearity and feedbacks pose a challenging optimization problem, which is overcome using a novel procedure that mimics natural speciation or diversification e.g., stepwise increase of bacteria species. The method allows for curation using species-level information from e.g., physiological experiments or genome sequences. The product is a mass-balancing, mechanistically-constrained, quantitative representation of the ecosystem. We apply the method to characterize phytoplankton—heterotrophic bacteria interactions via dissolved organic matter in a marine system. The resulting model predicts quantitative fluxes for each interaction and time point (e.g., 0.16 µmolC/L/d of chrysolaminarin to Polaribacter on April 16, 2009). At the system level, the flux network shows a strong correlation between the abundance of bacteria species and their carbon flux during blooms, with copiotrophs being relatively more important than oligotrophs. However, oligotrophs, like SAR11, are unexpectedly high carbon processors for weeks into blooms, due to their higher biomass. The fraction of exudates (vs. grazing/death products) in the DOM pool decreases during blooms, and they are preferentially consumed by oligotrophs. In addition, functional similarity of phytoplankton i.e., what they produce, decouples their association with heterotrophs. The methodology is applicable to other microbial ecosystems, like human microbiome or wastewater treatment plants.
Insights from Time Series of Atmospheric Carbon Dioxide and Related Tracers
