Microbial populations can withstand, overcome and persist in the face of environmental fluctuation. Previously, we demonstrated how conditional gene regulation in a fluctuating environment drives dilution of condition-specific transcripts, causing a population of Desulfovibrio vulgaris Hildenborough (DvH) to collapse after repeatedly transitioning from sulfate respiration to syntrophic conditions with the methanogen Methanococcus maripaludis. Failure of the DvH to successfully transition contributed to the collapse of this model community. We investigated the mechanistic basis for loss of robustness by examining whether conditional gene regulation altered heterogeneity in gene expression across individual DvH cells. We discovered that robustness of a microbial population across environmental transitions was attributable to the retention of cells in two states that exhibited different condition-specific gene expression patterns. In our experiments, a population with disrupted conditional regulation successfully alternated between cell states. Meanwhile, a population with intact conditional regulation successfully switched between cell states initially, but collapsed after repeated transitions, possibly due to the high energy requirements of regulation. These results demonstrate that the survival of this entire model microbial community is dependent on the regulatory system's influence on the distribution of distinct cell states among individual cells within a clonal population.

Background

Algae accumulate lipids to endure different kinds of environmental stresses including macronutrient starvation. Although this response has been extensively studied, an in depth understanding of the transcriptional regulatory network (TRN) that controls the transition into lipid accumulation remains elusive. In this study, we used a systems biology approach to elucidate the transcriptional program that coordinates the nitrogen starvation-induced metabolic readjustments that drive lipid accumulation in Chlamydomonas reinhardtii.

Results

We demonstrate that nitrogen starvation triggered differential regulation of 2147 transcripts, which were co-regulated in 215 distinct modules and temporally ordered as 31 transcriptional waves. An early-stage response was triggered within 12 min that initiated growth arrest through activation of key signaling pathways, while simultaneously preparing the intracellular environment for later stages by modulating transport processes and ubiquitin-mediated protein degradation. Subsequently, central metabolism and carbon fixation were remodeled to trigger the accumulation of triacylglycerols. Further analysis revealed that these waves of genome-wide transcriptional events were coordinated by a regulatory program orchestrated by at least 17 transcriptional regulators, many of which had not been previously implicated in this process. We demonstrate that the TRN coordinates transcriptional downregulation of 57 metabolic enzymes across a period of nearly 4 h to drive an increase in lipid content per unit biomass. Notably, this TRN appears to also drive lipid accumulation during sulfur starvation, while phosphorus starvation induces a different regulatory program. The TRN model described here is available as a community-wide web-resource at http://networks.systemsbiology.net/chlamy-portal.

Conclusions

In this work, we have uncovered a comprehensive mechanistic model of the TRN controlling the transition from N starvation to lipid accumulation. The program coordinates sequentially ordered transcriptional waves that simultaneously arrest growth and lead to lipid accumulation. This study has generated predictive tools that will aid in devising strategies for the rational manipulation of regulatory and metabolic networks for better biofuel and biomass production.

Resource partitioning is central to the incredible productivity of microbial communities, including gigatons in annual methane emissions through syntrophic interactions. Previous work revealed how a sulfate reducer (Desulfovibrio vulgaris, Dv) and a methanogen (Methanococcus maripaludis, Mm) underwent evolutionary diversification in a planktonic context, improving stability, cooperativity, and productivity within 300-1000 generations. Here, we show that mutations in just 15 Dv and 7 Mm genes within a minimal assemblage of this evolved community gave rise to co-existing ecotypes that were spatially enriched within a few days of culturing in a fluidized bed reactor. The spatially segregated communities partitioned resources in the simulated subsurface environment, with greater lactate utilization by attached Dv but partial utilization of resulting H₂ by low affinity hydrogenases of Mm in the same phase. The unutilized H₂ was scavenged by high affinity hydrogenases of planktonic Mm, producing copious amounts of methane. Our findings show how a few mutations can drive resource partitioning amongst niche-differentiated ecotypes, whose interplay synergistically improves productivity of the entire mutualistic community.

Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR- and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.