Eco-evolutionary dynamics of complex social strategies in microbial communities

Microbial communities thrive in virtually every habitat on Earth and are essential to the function of diverse ecosystems. Most microbial habitats are not spatially continuous and well-mixed, but rather composed, at the microscale, of many isolated or semi-isolated local patches of different sizes, resulting in partitioning of microbial populations into discrete local populations. The impact of this spatial fragmentation on population dynamics is not well-understood. Here, we study how such variably sized micro-habitat patches affect the growth dynamics of clonal microbial populations and how dynamics in individual patches dictate those of the metapopulation. To investigate this, we developed the μ-SPLASH, an ecology-on-a-chip platform, enabling the culture of microbes in microscopic landscapes comprised of thousands of microdroplets, with a wide range of sizes. Using the μ-SPLASH, we cultured the model bacteria Escherichia coli and based on time-lapse microscopy, analyzed the population dynamics within thousands of individual droplets. Our results reveal that growth curves substantially vary with droplet size. Although growth rates generally increase with drop size, reproductive success and the time to approach carrying capacity, display non-monotonic patterns. Combining μ-SPLASH experiments with computational modeling, we show that these patterns result from both stochastic and deterministic processes, and demonstrate the roles of initial population density, patchiness, and patch size distribution in dictating the local and metapopulation dynamics. This study reveals basic principles that elucidate the effects of habitat fragmentation and population partitioning on microbial population dynamics. These insights deepen our understanding of natural microbial communities and have significant implications for microbiome engineering.

Rapid environmental change is forcing populations into environments where plasticity will no longer maintain fitness. When populations are exposed to novel environments, evolutionary theory predicts that genetic variation in fitness will increase and should be associated with genetic differences in plasticity. If true, then genetic variation in plasticity can increase adaptive potential in novel environments, and population persistence via evolutionary rescue is more likely. To test whether genetic variation in fitness increases in novel environments and is associated with plasticity, we transplanted 8149 clones of 314 genotypes of a Sicilian daisy (Senecio chrysanthemifolius) within and outside its native range, and quantified genetic variation in fitness, and plasticity in leaf traits and gene expression. Although mean fitness declined by 87% in the novel environment, genetic variance in fitness increased threefold and was correlated with plasticity in leaf traits. High fitness genotypes showed greater plasticity in gene expression, but lower plasticity in most leaf traits. Interestingly, genotypes with the highestfitness in the novel environment had the lowest fitness at the native site. These results suggest that standing genetic variation in plasticity could help populations to persist and adapt to novel environments, despite remaining hidden in native environments.

Microbial mat communities consist of dense populations of microorganisms embedded in exopolymers and/or biomineralized solid phases, and are often found in mm-cm thick assemblages, which can be stratified due to environmental gradients such as light, oxygen or sulfide. Microbial mat communities are commonly observed under extreme environmental conditions, deriving energy primarily from light and/or reduced chemicals to drive autotrophic fixation of carbon dioxide. Microbial mat ecosystems are regarded as living analogues of primordial systems on Earth, and they often form perennial structures with conspicuous stratifications of microbial populations that can be studied in situ under stable conditions for many years. Consequently, microbial mat communities are ideal natural laboratories and represent excellent model systems for studying microbial community structure and function, microbial dynamics and interactions, and discovery of new microorganisms with novel metabolic pathways potentially useful in future industrial and/or medical applications. Due to their relative simplicity and organization, microbial mat communities are often excellent testing grounds for new technologies in microbiology including micro-sensor analysis, stable isotope methodology and modern genomics. Integrative studies of microbial mat communities that combine modern biogeochemical and molecular biological methods with traditional microbiology, macro-ecological approaches, and community network modeling will provide new and detailed insights regarding the systems biology of microbial mats and the complex interplay among individual populations and their physicochemical environment. These processes ultimately control the biogeochemical cycling of energy and/or nutrients in microbial systems. Similarities in microbial community function across different types of communities from highly disparate environments may provide a deeper basis for understanding microbial community dynamics and the ecological role of specific microbial populations. Approaches and concepts developed in highly-constrained, relatively stable natural communities may also provide insights useful for studying and understanding more complex microbial communities.

Gene expression responds to the environment and can also evolve rapidly in response to altered selection regimes. Little is known, however, about the extent to which evolutionary adaptation to a particular type of stress involves changes in the within-generation ('plastic') responses of gene expression to the stress. We used microarrays to quantify gene expression plasticity in response to ethanol in laboratory populations of Drosophila melanogaster differing in their history of ethanol exposure. Two populations ('R' populations) were maintained on regular medium, two ('E') were maintained on medium supplemented with ethanol, and two ('M') were maintained in a mixed regime in which half of the population was reared on one medium type, and half on the other, each generation. After more than 300 generations, embryos from each population were collected and exposed to either ethanol or water as a control, and RNA was extracted from the larvae shortly after hatching. Nearly 2000 transcripts showed significant within-generation responses to ethanol exposure. Evolutionary history also affected gene expression: the E and M populations were largely indistinguishable in expression, but differed significantly in expression from the R populations for over 100 transcripts, the majority of which did not show plastic responses. Notably, in no case was the interaction between selection regime and ethanol exposure significant after controlling for multiple comparisons, indicating that adaptation to ethanol in the E and M populations did not involve substantial changes in gene expression plasticity. The results give evidence that expression plasticity evolves considerably more slowly than mean expression.

The dynamics of individual microbial populations and their gene functions in agricultural soils, especially after major activities such as nitrogen (N) fertilization, remain elusive but are important for a better understanding of nutrient cycling. Here, we analyzed 20 short-read metagenomes collected at four time points during 1 year from two depths (0 to 5 and 20 to 30 cm) in two Midwestern agricultural sites representing contrasting soil textures (sandy versus silty loam) with similar cropping histories. Although the microbial community taxonomic and functional compositions differed between the two locations and depths, they were more stable within a depth/site throughout the year than communities in natural aquatic ecosystems. For example, among the 69 population genomes assembled from the metagenomes, 75% showed a less than 2-fold change in abundance between any two sampling points. Interestingly, six deep-branching Thaumarchaeota and three complete ammonia oxidizer (comammox) Nitrospira populations increased up to 5-fold in abundance upon the addition of N fertilizer. These results indicated that indigenous archaeal ammonia oxidizers may respond faster (are more copiotrophic) to N fertilization than previously thought. None of 29 recovered putative denitrifier genomes encoded the complete denitrification pathway, suggesting that denitrification is carried out by a collection of different populations. Altogether, our study identified novel microbial populations and genes responding to seasonal and human-induced perturbations in agricultural soils that should facilitate future monitoring efforts and N-related studies.IMPORTANCE Even though the impact of agricultural management on the microbial community structure has been recognized, an understanding of the dynamics of individual microbial populations and what functions each population carries are limited. Yet, this information is important for a better understanding of nutrient cycling, with potentially important implications for preserving nitrogen in soils and sustainability. Here, we show that reconstructed metagenome-assembled genomes (MAGs) are relatively stable in their abundance and functional gene content year round, and seasonal nitrogen fertilization has selected for novel Thaumarchaeota and comammox Nitrospira nitrifiers that are potentially less oligotrophic than their marine counterparts previously studied.

Theories to explain the diversity of species have required that individual species occupy unique niches and/or vary in their response to environmental factors. Positive interactions within a species, although common in communities, have not been thought to maintain species diversity because in non-spatial models the more abundant species always outcompetes the rarer species. Here, we show, using a stochastic spatial model, that positive intraspecific interactions such as those caused by positive frequency dependence and/or priority effects, can maintain species diversity if interactions between individuals are primarily local and the habitat contains areas that cannot be colonized by any species, such as boulders or other physical obstructions. When intraspecific interactions are primarily neutral, species diversity will eventually erode to a single species. When the landscape is homogeneous (i.e. does not contain areas that cannot be colonized by any species), the presence of strong intraspecific interactions will not maintain diversity.

Microbial population dynamics on leaves in time and space are a function of immigration, emigration, growth, and death. Insight into the relative significance of each population process to the generation of specific dynamics for individual microorganisms is necessary to understanding the ecology and life history strategy of the microorganism and to developing effective control strategies. Additionally, information on the significance of within-leaf versus extra-leaf processes to the generation of phyllosphere dynamics is important to determining the range of spatial scales over which a population should be studied. Unfortunately, such information is difficult to obtain due to the lack of effective methodologies for distinguishing these processes within phyllosphere populations. Future research efforts should focus on the quantification of immigration, emigration, growth, and death relative to the population dynamics of phyllosphere microorganisms.

EDITORIAL article Front. Microbiol., 26 September 2014Sec. Terrestrial Microbiology Volume 5 - 2014 | https://doi.org/10.3389/fmicb.2014.00497

The difficulty involved in quantifying biogeochemically significant microbes in marine sediments limits our ability to assess interspecific interactions, population turnover times, and niches of uncultured taxa. We incubated surface sediments from Cape Lookout Bight, North Carolina, USA, anoxically at 21°C for 122 days. Sulfate decreased until day 68, after which methane increased, with hydrogen concentrations consistent with the predicted values of an electron donor exerting thermodynamic control. We measured turnover times using two relative quantification methods, quantitative PCR (qPCR) and the product of 16S gene read abundance and total cell abundance (FRAxC, which stands for "fraction of read abundance times cells"), to estimate the population turnover rates of uncultured clades. Most 16S rRNA reads were from deeply branching uncultured groups, and ∼98% of 16S rRNA genes did not abruptly shift in relative abundance when sulfate reduction gave way to methanogenesis. Uncultured Methanomicrobiales and Methanosarcinales increased at the onset of methanogenesis with population turnover times estimated from qPCR at 9.7 ± 3.9 and 12.6 ± 4.1 days, respectively. These were consistent with FRAxC turnover times of 9.4 ± 5.8 and 9.2 ± 3.5 days, respectively. Uncultured Syntrophaceae, which are possibly fermentative syntrophs of methanogens, and uncultured Kazan-3A-21 archaea also increased at the onset of methanogenesis, with FRAxC turnover times of 14.7 ± 6.9 and 10.6 ± 3.6 days. Kazan-3A-21 may therefore either perform methanogenesis or form a fermentative syntrophy with methanogens. Three genera of sulfate-reducing bacteria, Desulfovibrio, Desulfobacter, and Desulfobacterium, increased in the first 19 days before declining rapidly during sulfate reduction. We conclude that population turnover times on the order of days can be measured robustly in organic-rich marine sediment, and the transition from sulfate-reducing to methanogenic conditions stimulates growth only in a few clades directly involved in methanogenesis, rather than in the whole microbial community.IMPORTANCE Many microbes cannot be isolated in pure culture to determine their preferential growth conditions and predict their response to changing environmental conditions. We created a microcosm of marine sediments that allowed us to simulate a diagenetic profile using a temporal analog for depth. This allowed for the observation of the microbial community population dynamics caused by the natural shift from sulfate reduction to methanogenesis. Our research provides evidence for the population dynamics of uncultured microbes as well as the application of a novel method of turnover rate analysis for individual taxa within a mixed incubation, FRAxC, which stands for "fraction of read abundance times cells," which was verified by quantitative PCR. This allows for the calculation of population turnover times for microbes in a natural setting and the identification of uncultured clades involved in geochemical processes.

Understanding how populations adapt to changing environmental conditions is a long-standing theme in evolutionary biology. Gene expression changes have been recognized as an important driver of local adaptation, but relatively little is known regarding the direction of change and in particular, about the interplay between plastic and evolutionary gene expression. We have previously shown that the gene expression profiles of European grayling (Thymallus thymallus) populations inhabiting different thermal environments include both plastic and evolutionary components. However, whether the plastic and evolutionary responses were in the same direction was not investigated in detail, nor was the identity of the specific genes involved. In this study, we show that the plastic changes in protein expression in response to different temperatures are highly correlated with the evolutionary response in grayling subpopulations adapted to different thermal environments. This finding provides preliminary evidence that the plastic response most likely facilitates adaptation during the early phases of colonization of thermal environments. The proteins that showed significant changes in expression level between warm and cold temperature treatments were mostly related to muscle development, which is consistent with earlier findings demonstrating muscle mass differentiation between cold and warm grayling populations.

Although there is a long history of conjecture regarding the role and significance of bacterial chemotaxis in microbial ecology, only recently has a significant body of work appeared attempting to address this issue. The purpose of this paper is to provide a concise overview of this work, which combined mathematical modeling of bacterial population migration and experimental measurement of the model parameters with modeling of competitive microbial population dynamics in a nonmixed environment. Predictions from the population dynamics models, based on experimental estimates of the various motility and growth parameter values, are related to the small number of experimental observations available to date dealing with the effects of bacterial motility on competition in a nonmixed environment. Current results indicate that cell motility and chemotaxis properties can be as important to population dynamics as cell growth kinetic properties, so that greater attention to this aspect of microbial behavior is warranted in future studies of microbial ecology.

Absolute quantification and genome-centric analyses elucidate the dynamics of microbial populations in anaerobic digesters

Eco-evolutionary feedback loops, where rapid evolution influences the ecology of an organism and subsequently alters the evolutionary trajectory of the population, are intriguing possibilities, but evidence for or against them is scarce. Timema cristinae stick insects express variation within and among populations in the expression of death-feigning behaviour, but the causes of variation in this species is not known. Here, I test the hypothesis that variation in death feigning across populations stems from an eco-evolutionary feedback loop, whereby gene flow causes poor camouflage, which increases the strength of bird predation, and ultimately drives the evolution of increased death-feigning. By conducting behavioural trials on eight T. cristinae populations that differ in the degree of maladaptive gene flow incurred, I falsify the eco-evolutionary feedback hypothesis for the evolution of death-feigning. Instead, I show that smaller individuals are more likely to feign death than larger individuals. By rearing individuals in the lab, I offer suggestive evidence that the body size effect is explained by the age of wild-caught individuals: younger individuals feign death more than older individuals. These findings add an example to the literature where no eco-evolutionary feedback exists in a system for which other similar feedbacks have been found, and provide evidence that death-feigning behaviour depends on body size.

Expression analysis of synaptic plasticity genes in curcumin-treated amnesic mice model

Transgenic or genetically modified plants possess novel genes that impart beneficial characteristics such as herbicide resistance. One of the least understood areas in the environmental risk assessment of genetically modified crops is their impact on soil- and plant-associated microbial communities. The potential for interaction between transgenic plants and plant residues and the soil microbial community is not well understood. The recognition that these interactions could change microbial biodiversity and affect ecosystem functioning has initiated a limited number of studies in the area. At this time, studies have shown the possibility that transgenes can be transferred to native soil microorganisms through horizontal gene transfer, although there is not evidence of this occurring in the soil. Furthermore, novel proteins have been shown to be released from transgenic plants into the soil ecosystem, and their presence can influence the biodiversity of the microbial community by selectively stimulating the growth of organisms that can use them. Microbial diversity can be altered when associated with transgenic plants; however, these effects are both variable and transient. Soil- and plant-associated microbial communities are influenced not only by plant species and transgene insertion but also by environmental factors such as field site and sampling date. Minor alterations in the diversity of the microbial community could affect soil health and ecosystem functioning, and therefore, the impact that plant variety may have on the dynamics of the rhizosphere microbial populations and in turn plant growth and health and ecosystem sustainability, requires further study.