Mechanism for microbial population collapse in a fluctuating resource environment

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing non-genetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupts acquired resistance in GBM.

Large-scale genetic characterization of the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough

Sulfate-reducing bacteria (SRB) are obligate anaerobes that can couple their growth to the reduction of sulfate. Despite the importance of SRB to global nutrient cycles and their damage to the petroleum industry, our molecular understanding of their physiology remains limited. To systematically provide new insights into SRB biology, we generated a randomly barcoded transposon mutant library in the model SRB Desulfovibrio vulgaris Hildenborough (DvH) and used this genome-wide resource to assay the importance of its genes under a range of metabolic and stress conditions. In addition to defining the essential gene set of DvH, we identified a conditional phenotype for 1,137 non-essential genes. Through examination of these conditional phenotypes, we were able to make a number of novel insights into our molecular understanding of DvH, including how this bacterium synthesizes vitamins. For example, we identified DVU0867 as an atypical L-aspartate decarboxylase required for the synthesis of pantothenic acid, provided the first experimental evidence that biotin synthesis in DvH occurs via a specialized acyl carrier protein and without methyl esters, and demonstrated that the uncharacterized dehydrogenase DVU0826:DVU0827 is necessary for the synthesis of pyridoxal phosphate. In addition, we used the mutant fitness data to identify genes involved in the assimilation of diverse nitrogen sources and gained insights into the mechanism of inhibition of chlorate and molybdate. Our large-scale fitness dataset and RB-TnSeq mutant library are community-wide resources that can be used to generate further testable hypotheses into the gene functions of this environmentally and industrially important group of bacteria.

Microscale dynamics of dark zone alterations in anthropized karstic cave shows abrupt microbial community switch

Strong anthropization of karstic caves may result in formation of various wall alterations including dark zones, whose microbial community differs from that of non-altered surfaces nearby. Dark zones grow quickly and without gradual visual changes, leading to the hypothesis of a simple process rather than complex microbial successions, but this is counter-intuitive as underground microbial changes are typically slow and dark zones are microbiologically very distinct from unmarked surfaces. We tested this hypothesis in Paleolithic Lascaux Cave, across two years of microscale sampling. Indeed, Illumina MiSeq metabarcoding evidenced only three community stages for bacteria, fungi and all microeukaryotes together (i.e. unmarked surfaces, newly-formed dark zones and intermediate/old dark zones) and just two stages for archaea (unmarked surfaces vs dark zones), indicating abrupt community changes. The onset of dark zone formation coincided with the development of Ochroconis fungi, Bacteroidota and the bacterial genera Labrys, Nonomuraea and Sphingomonas, in parallel to Pseudomonas counter-selection. Modeling of community assembly processes highlighted that the dynamics of rare taxa in unmarked surfaces adjacent to dark zones and in newly-formed dark zones were governed in part by deterministic processes. This suggests that cooperative relationships between these taxa might be important to promote dark zone formation. Taken together, these findings indicate an abrupt community switch as these new alterations form on Lascaux cave walls.

Microbial maintenance energy quantified and modeled with microcalorimetry

Refining the energetic costs of cellular maintenance is essential for predicting microbial growth and survival in the environment. Here, we evaluate a simple batch culture method to quantify energy partitioning between growth and maintenance using microcalorimetry and thermodynamic modeling. The constants derived from the batch culture system were comparable to those that have been reported from meta-analyses of data derived from chemostat studies. The model accurately predicted temperature-dependent biomass yield and the upper temperature limit of growth for Desulfovibrio alaskensis G20, suggesting the method may have broad application. An Arrhenius temperature dependence for the specific energy consumption rate, inferred from substrate consumption and heat evolution, was observed over the entire viable temperature range. By combining this relationship for specific energy consumption rates and observed specific growth rates, the model describes an increase in nongrowth associated maintenance at higher temperatures and the corresponding decrease in energy available for growth. This analytical and thermodynamic formulation suggests that simply monitoring heat evolution in batch culture could be a useful complement to the recognized limitations of estimating maintenance using extrapolation to zero growth in chemostats.

Adaptive laboratory evolution of microbial co-cultures for improved metabolite secretion

Adaptive laboratory evolution has proven highly effective for obtaining microorganisms with enhanced capabilities. Yet, this method is inherently restricted to the traits that are positively linked to cell fitness, such as nutrient utilization. Here, we introduce coevolution of obligatory mutualistic communities for improving secretion of fitness‐costly metabolites through natural selection. In this strategy, metabolic cross‐feeding connects secretion of the target metabolite, despite its cost to the secretor, to the survival and proliferation of the entire community. We thus co‐evolved wild‐type lactic acid bacteria and engineered auxotrophic Saccharomyces cerevisiae in a synthetic growth medium leading to bacterial isolates with enhanced secretion of two B‐group vitamins, viz., riboflavin and folate. The increased production was specific to the targeted vitamin, and evident also in milk, a more complex nutrient environment that naturally contains vitamins. Genomic, proteomic and metabolomic analyses of the evolved lactic acid bacteria, in combination with flux balance analysis, showed altered metabolic regulation towards increased supply of the vitamin precursors. Together, our findings demonstrate how microbial metabolism adapts to mutualistic lifestyle through enhanced metabolite exchange.

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

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.

Energetics of Acidianus ambivalens growth in response to oxygen availability

All life requires energy to drive metabolic reactions such as growth and cell maintenance; therefore, fluctuations in energy availability can alter microbial activity. There is a gap in our knowledge concerning how energy availability affects the growth of extreme chemolithoautotrophs. Toward this end, we investigated the growth of thermoacidophile Acidianus ambivalens during sulfur oxidation under aerobic to microaerophilic conditions. Calorimetry was used to measure enthalpy (ΔH_inc ) of microbial activity, and chemical changes in growth media were measured to calculate Gibbs energy change (ΔG_inc ) during incubation. In all experiments, Gibbs energy was primarily dissipated through the release of heat, which suggests enthalpy-driven growth. In microaerophilic conditions, growth was significantly more efficient in terms of biomass yield (defined as C-mol biomass per mole sulfur consumed) and resulted in lower ΔG_inc and ΔH_inc . ΔG_inc in oxygen-limited (OL) and oxygen- and CO₂ -limited (OCL) microaerophilic growth conditions resulted in averages of -1.44 × 10³  kJ/C-mol and -7.56 × 10²  kJ/C-mol, respectively, and average ΔH_inc values of -1.11 × 10⁵  kJ/C-mol and -4.43 × 10⁴  kJ/C-mol, respectively. High-oxygen experiments resulted in lower biomass yield values, an increase in ΔG_inc to -1.71 × 10⁴  kJ/C-mol, and more exothermic ΔH_inc values of -4.71 × 10⁵  kJ/C-mol. The observed inefficiency in high-oxygen conditions may suggest larger maintenance energy demands due to oxidative stresses and a preference for growth in microaerophilic environments.

Microbial communities display alternative stable states in a fluctuating environment

The effect of environmental fluctuations is a major question in ecology. While it is widely accepted that fluctuations and other types of disturbances can increase biodiversity, there are fewer examples of other types of outcomes in a fluctuating environment. Here we explore this question with laboratory microcosms, using cocultures of two bacterial species, P. putida and P. veronii. At low dilution rates we observe competitive exclusion of P. veronii, whereas at high dilution rates we observe competitive exclusion of P. putida. When the dilution rate alternates between high and low, we do not observe coexistence between the species, but rather alternative stable states, in which only one species survives and initial species’ fractions determine the identity of the surviving species. The Lotka-Volterra model with a fluctuating mortality rate predicts that this outcome is independent of the timing of the fluctuations, and that the time-averaged mortality would also lead to alternative stable states, a prediction that we confirm experimentally. Other pairs of species can coexist in a fluctuating environment, and again consistent with the model we observe coexistence in the time-averaged dilution rate. We find a similar time-averaging result holds in a three-species community, highlighting that simple linear models can in some cases provide powerful insight into how communities will respond to environmental fluctuations.

The effect of environmental fluctuations on community structure and function is a fundamental question in ecology. A significant body of work suggests that fluctuations increase diversity due to a variety of proposed mechanisms. In this study, we compare the effects of constant and fluctuating dilution regimes on simple microbial communities with two or three species. We find that in all cases, the outcome in a fluctuating environment is the same as that in a constant environment in which the fluctuations are time-averaged. This surprising result highlights that in some communities, ecological stable states may be predicted by averaging environmental parameters, rather than by the variation itself.

Microbial energy management-A product of three broad tradeoffs

Wherever thermodynamics allows, microbial life has evolved to transform and harness energy. Microbial life thus abounds in the most unexpected places, enabled by profound metabolic diversity. Within this diversity, energy is transformed primarily through variations on a few core mechanisms. Energy is further managed by the physiological processes of cell growth and maintenance that use energy. Some aspects of microbial physiology are streamlined for energetic efficiency while other aspects seem suboptimal or even wasteful. We propose that the energy that a microbe harnesses and devotes to growth and maintenance is a product of three broad tradeoffs: (i) economic, trading enzyme synthesis or operational cost for functional benefit, (ii) environmental, trading optimization for a single environment for adaptability to multiple environments, and (iii) thermodynamic, trading energetic yield for forward metabolic flux. Consideration of these tradeoffs allows one to reconcile features of microbial physiology that seem to opposingly promote either energetic efficiency or waste.

Legacy of prior host and soil selection on rhizobial fitness in planta

Measuring selection acting on microbial populations in natural or even seminatural environments is challenging because many microbial populations experience variable selection. The majority of rhizobial bacteria are found in the soil. However, they also live symbiotically inside nodules of legume hosts and each nodule can release thousands of daughter cells back into the soil. We tested how past selection (i.e., legacies) by two plant genotypes and by the soil alone affected selection and genetic diversity within a population of 101 strains of Ensifer meliloti. We also identified allelic variants most strongly associated with soil- and host-dependent fitness. In addition to imposing direct selection on rhizobia populations, soil and host environments had lasting effects across host generations. Host presence and genotype during the legacy period explained 22% and 12% of the variance in the strain composition of nodule communities in the second cohort, respectively. Although strains with high host fitness in the legacy cohort tended to be enriched in the second cohort, the diversity of the strain community was greater when the second cohort was preceded by host rather than soil legacies. Our results indicate the potential importance of soil selection driving the evolution of these plant-associated microbes.

The rate of environmental fluctuations shapes ecological dynamics in a two-species microbial system

Species interactions change when the external conditions change. How these changes affect microbial community properties is an open question. We address this question using a two-species consortium in which species interactions change from exploitation to competition depending on the carbon source provided. We built a mathematical model and calibrated it using single-species growth measurements. This model predicted that low frequencies of change between carbon sources lead to species loss, while intermediate and high frequencies of change maintained both species. We experimentally confirmed these predictions by growing co-cultures in fluctuating environments. These findings complement more established concepts of a diversity peak at intermediate disturbance frequencies. They also provide a mechanistic understanding for how the dynamics at the community level emerges from single-species behaviours and interspecific interactions. Our findings suggest that changes in species interactions can profoundly impact the ecological dynamics and properties of microbial systems.

An Escherichia coli Nitrogen Starvation Response Is Important for Mutualistic Coexistence with Rhodopseudomonas palustris

Microbial mutualistic cross-feeding interactions are ubiquitous and can drive important community functions. Engaging in cross-feeding undoubtedly affects the physiology and metabolism of individual species involved. However, the nature in which an individual species' physiology is influenced by cross-feeding and the importance of those physiological changes for the mutualism have received little attention. We previously developed a genetically tractable coculture to study bacterial mutualisms. The coculture consists of fermentative Escherichia coli and phototrophic Rhodopseudomonas palustris In this coculture, E. coli anaerobically ferments sugars into excreted organic acids as a carbon source for R. palustris In return, a genetically engineered R. palustris strain constitutively converts N₂ into NH₄⁺, providing E. coli with essential nitrogen. Using transcriptome sequencing (RNA-seq) and proteomics, we identified transcript and protein levels that differ in each partner when grown in coculture versus monoculture. When in coculture with R. palustris, E. coli gene expression changes resembled a nitrogen starvation response under the control of the transcriptional regulator NtrC. By genetically disrupting E. coli NtrC, we determined that a nitrogen starvation response is important for a stable coexistence, especially at low R. palustris NH₄⁺ excretion levels. Destabilization of the nitrogen starvation regulatory network resulted in variable growth trends and, in some cases, extinction. Our results highlight that alternative physiological states can be important for survival within cooperative cross-feeding relationships.IMPORTANCE Mutualistic cross-feeding between microbes within multispecies communities is widespread. Studying how mutualistic interactions influence the physiology of each species involved is important for understanding how mutualisms function and persist in both natural and applied settings. Using a bacterial mutualism consisting of Rhodopseudomonas palustris and Escherichia coli growing cooperatively through bidirectional nutrient exchange, we determined that an E. coli nitrogen starvation response is important for maintaining a stable coexistence. The lack of an E. coli nitrogen starvation response ultimately destabilized the mutualism and, in some cases, led to community collapse after serial transfers. Our findings thus inform on the potential necessity of an alternative physiological state for mutualistic coexistence with another species compared to the physiology of species grown in isolation.

Ocean acidification conditions increase resilience of marine diatoms

The fate of diatoms in future acidified oceans could have dramatic implications on marine ecosystems, because they account for ~40% of marine primary production. Here, we quantify resilience of Thalassiosira pseudonana in mid-20th century (300 ppm CO2) and future (1000 ppm CO2) conditions that cause ocean acidification, using a stress test that probes its ability to recover from incrementally higher amount of low-dose ultraviolet A (UVA) and B (UVB) radiation and re-initiate growth in day-night cycles, limited by nitrogen. While all cultures eventually collapse, those growing at 300 ppm CO2 succumb sooner. The underlying mechanism for collapse appears to be a system failure resulting from "loss of relational resilience," that is, inability to adopt physiological states matched to N-availability and phase of the diurnal cycle. Importantly, under elevated CO2 conditions diatoms sustain relational resilience over a longer timeframe, demonstrating increased resilience to future acidified ocean conditions. This stress test framework can be extended to evaluate and predict how various climate change associated stressors may impact microbial community resilience.

Evolution on the bright side of life: microorganisms and the evolution of mutualism

Mutualistic interactions, where two interacting species have a net beneficial effect on each other's fitness, play a crucial role in the survival and evolution of many species. Despite substantial empirical and theoretical work in past decades, the impact of these interactions on natural selection is not fully understood. In addition, mutualisms between microorganisms have been largely ignored, even though they are ecologically important and can be used as tools to bridge the gap between theory and empirical work. Here, I describe two problems with our current understanding of natural selection in mutualism and highlight the properties of microbial mutualisms that could help solve them. One problem is that bias and methodological problems have limited our understanding of the variety of mechanisms by which species may adapt to mutualism. Another problem is that it is rare for experiments testing coevolution in mutualism to address whether each species has adapted to evolutionary changes in its partner. These problems can be addressed with genome resequencing and time-shift experiments, techniques that are easier to perform in microorganisms. In addition, microbial mutualisms may inspire novel insights and hypotheses about natural selection in mutualism.

Tools for Genomic and Transcriptomic Analysis of Microbes at Single-Cell Level

Microbiologists traditionally study population rather than individual cells, as it is generally assumed that the status of individual cells will be similar to that observed in the population. However, the recent studies have shown that the individual behavior of each single cell could be quite different from that of the whole population, suggesting the importance of extending traditional microbiology studies to single-cell level. With recent technological advances, such as flow cytometry, next-generation sequencing (NGS), and microspectroscopy, single-cell microbiology has greatly enhanced the understanding of individuality and heterogeneity of microbes in many biological systems. Notably, the application of multiple ‘omics’ in single-cell analysis has shed light on how individual cells perceive, respond, and adapt to the environment, how heterogeneity arises under external stress and finally determines the fate of the whole population, and how microbes survive under natural conditions. As single-cell analysis involves no axenic cultivation of target microorganism, it has also been demonstrated as a valuable tool for dissecting the microbial ‘dark matter.’ In this review, current state-of-the-art tools and methods for genomic and transcriptomic analysis of microbes at single-cell level were critically summarized, including single-cell isolation methods and experimental strategies of single-cell analysis with NGS. In addition, perspectives on the future trends of technology development in the field of single-cell analysis was also presented.

Robustness of a model microbial community emerges from population structure among single cells of a clonal population

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.

Mechanism for microbial population collapse in a fluctuating resource environment

Managing trade-offs through gene regulation is believed to confer resilience to a microbial community in a fluctuating resource environment. To investigate this hypothesis, we imposed a fluctuating environment that required the sulfate-reducer Desulfovibrio vulgaris to undergo repeated ecologically relevant shifts between retaining metabolic independence (active capacity for sulfate respiration) and becoming metabolically specialized to a mutualistic association with the hydrogen-consuming Methanococcus maripaludis Strikingly, the microbial community became progressively less proficient at restoring the environmentally relevant physiological state after each perturbation and most cultures collapsed within 3-7 shifts. Counterintuitively, the collapse phenomenon was prevented by a single regulatory mutation. We have characterized the mechanism for collapse by conducting RNA-seq analysis, proteomics, microcalorimetry, and single-cell transcriptome analysis. We demonstrate that the collapse was caused by conditional gene regulation, which drove precipitous decline in intracellular abundance of essential transcripts and proteins, imposing greater energetic burden of regulation to restore function in a fluctuating environment.