Neutral Models of Microbiome Evolution

Comparative gut microbiome research through the lens of ecology: theoretical considerations and best practices

Comparative approaches in animal gut microbiome research have revealed patterns of phylosymbiosis, dietary and physiological convergences, and environment-host interactions. However, most large-scale comparative studies, especially those that are highly cited, have focused on mammals, and efforts to integrate comparative approaches with existing ecological frameworks are lacking. While mammals serve as useful model organisms, developing generalised principles of how animal gut microbiomes are shaped and how these microbiomes interact bidirectionally with host ecology and evolution requires a more complete sampling of the animal kingdom. Here, we provide an overview of what past comparative studies have taught us about the gut microbiome, and how community ecology theory may help resolve certain contradictions in comparative gut microbiome research. We explore whether certain hypotheses are supported across clades, and how the disproportionate focus on mammals has introduced potential bias into gut microbiome theory. We then introduce a methodological solution by which public gut microbiome data of understudied hosts can be compiled and analysed in a comparative context. Our aggregation and analysis of 179 studies shows that generating data sets with rich host diversity is possible with public data and that key gut microbes associated with mammals are widespread across the animal kingdom. We also show the effects that sample size and taxonomic rank have on comparative gut microbiome studies and that results of multivariate analyses can vary significantly with these two parameters. While challenges remain in developing a universal model of the animal gut microbiome, we show that existing ecological frameworks can help bring us one step closer to integrating the gut microbiome into animal ecology and evolution.

Maternal transmission as a microbial symbiont sieve, and the absence of lactation in male mammals

Gut microbiomes of mammals carry a complex symbiotic assemblage of microorganisms. Feeding newborn infants milk from the mammary gland allows vertical transmission of the parental milk microbiome to the offspring's gut microbiome. This has benefits, but also has hazards for the host population. Using mathematical models, we demonstrate that biparental vertical transmission enables deleterious microbial elements to invade host populations. In contrast, uniparental vertical transmission acts as a sieve, preventing these invasions. Moreover, we show that deleterious symbionts generate selection on host modifier genes that keep uniparental transmission in place. Since microbial transmission occurs during birth in placental mammals, subsequent transmission of the milk microbiome needs to be maternal to avoid the spread of deleterious elements. This paper therefore argues that viviparity and the hazards from biparental transmission of the milk microbiome, together generate selection against male lactation in placental mammals.

A computational framework for resolving the microbiome diversity conundrum

Recent empirical studies offer conflicting findings regarding the relation between host fitness and the composition of its microbiome, a conflict which we term 'the microbial β- diversity conundrum'. The microbiome is crucial for host wellbeing and survival. Surprisingly, different healthy individuals' microbiome compositions, even in the same population, often differ dramatically, contrary to the notion that a vital trait should be highly conserved. Moreover, gnotobiotic individuals exhibit highly deleterious phenotypes, supporting the view that the microbiome is paramount to host fitness. However, the introduction of almost arbitrarily selected microbiota into the system often achieves a significant rescue effect of the deleterious phenotypes. This is true even for microbiota from soil or phylogenetically distant host species, highlighting an apparent paradox. We suggest several solutions to the paradox using a computational framework, simulating the population dynamics of hosts and their microbiomes over multiple generations. The answers invoke factors such as host population size, the specific mode of microbial contribution to host fitness, and typical microbiome richness, offering solutions to the conundrum by highlighting scenarios where even when a host's fitness is determined in full by its microbiome composition, this composition has little effect on the natural selection dynamics of the population.

Gut microbes predominantly act as living beneficial partners rather than raw nutrients

Animals and their gut microbes mutually benefit their health. Nutrition plays a central role in this, directly influencing both host and microbial fitness and the nature of their interactions. This makes nutritional symbioses a complex and dynamic tri-system of diet-microbiota-host. Despite recent discoveries on this field, full control over the interplay among these partners is challenging and hinders the resolution of fundamental questions, such as how to parse the gut microbes' effect as raw nutrition or as symbiotic partners? To tackle this, we made use of the well-characterized Drosophila melanogaster/Lactiplantibacillus plantarum experimental model of nutritional symbiosis to generate a quantitative framework of gut microbes' effect on the host. By coupling experimental assays and Random Forest analysis, we show that the beneficial effect of L. plantarum strains primarily results from the active relationship as symbionts rather than raw nutrients, regardless of the bacterial strain. Metabolomic analysis of both active and inactive bacterial cells further demonstrated the crucial role of the production of beneficial bacterial metabolites, such as N-acetylated-amino-acids, as result of active bacterial growth and function. Altogether, our results provide a ranking and quantification of the main bacterial features contributing to sustain animal growth. We demonstrate that bacterial activity is the predominant and necessary variable involved in bacteria-mediated benefit, followed by strain-specific properties and the nutritional potential of the bacterial cells. This contributes to elucidate the role of beneficial bacteria and probiotics, creating a broad quantitative framework for host-gut microbiome that can be expanded to other model systems.

On the effect of inheritance of microbes in commensal microbiomes

Background

Our current view of nature depicts a world where macroorganisms dwell in a landscape full of microbes. Some of these microbes not only transit but establish themselves in or on hosts. Although hosts might be occupied by microbes for most of their lives, a microbe-free stage during their prenatal development seems to be the rule for many hosts. The questions of who the first colonizers of a newborn host are and to what extent these are obtained from the parents follow naturally.

Results

We have developed a mathematical model to study the effect of the transfer of microbes from parents to offspring. Even without selection, we observe that microbial inheritance is particularly effective in modifying the microbiome of hosts with a short lifespan or limited colonization from the environment, for example by favouring the acquisition of rare microbes.

Conclusion

By modelling the inheritance of commensal microbes to newborns, our results suggest that, in an eco-evolutionary context, the impact of microbial inheritance is of particular importance for some specific life histories.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02029-2.

Gut microbiota plasticity in insular lizards under reversed island syndrome

Animals living on small islands are more drastically exposed to environmental changes, such as food or water starvation, and rapid temperature shifts. Facing such conditions, and probably thank to adaptive plasticity mechanisms, some animals display a Reversed Island Syndrome (RIS), a suite of traits, including skin pigmentation, voracity, sexual dimorphism, showed differently from mainland relatives. Here, we analyse a so far poorly explored aspect of RIS: the effect of this on the microbiota composition of host Italian wall lizard (Podarcis siculus), strongly influenced by the animal’s lifestyle, and conditioning the same. We compare mainland and island populations, assessing the difference between their microbial communities and their response under unexpected food, experimentally provided. Our observations showed a significant difference in microbiota communities between island and mainland groups, depended mainly from changes in relative abundance of the shared genera (difference due to decrease/increase). Exposure to experimental diet regimes resulted into significative reshaping of bacterial composition of microbiota and a greater variation in body mass only in the island population. Our results could be an evidence that gut microbial community contributes to adaptive plasticity mechanisms of island lizards under RIS to efficiently respond to unexpected changes.

Host genomic influence on bacterial composition in the switchgrass rhizosphere

Host genetic variation can shape the diversity and composition of associated microbiomes, which may reciprocally influence host traits and performance. While the genetic basis of phenotypic diversity of plant populations in nature has been studied, comparatively little research has investigated the genetics of host effects on their associated microbiomes. Switchgrass (Panicum virgatum) is a highly outcrossing, perennial, grass species with substantial locally adaptive diversity across its native North American range. Here, we compared 383 switchgrass accessions in a common garden to determine the host genotypic influence on rhizosphere bacterial composition. We hypothesized that the composition and diversity of rhizosphere bacterial assemblages would differentiate due to genotypic differences between hosts (potentially due to root phenotypes and associated life history variation). We observed higher alpha diversity of bacteria associated with upland ecotypes and tetraploids, compared to lowland ecotypes and octoploids, respectively. Alpha diversity correlated negatively with flowering time and plant height, indicating that bacterial composition varies along switchgrass life history axes. Narrow-sense heritability (h² ) of the relative abundance of twenty-one core bacterial families was observed. Overall compositional differences among tetraploids, due to genetic variation, supports wide-spread genotypic influence on the rhizosphere microbiome. Tetraploids were only considered due to complexities associated with the octoploid genomes. Lastly, a genome-wide association study identified 1,861 single-nucleotide polymorphisms associated with 110 families and genes containing them related to potential regulatory functions. Our findings suggest that switchgrass genomic and life-history variation influences bacterial composition in the rhizosphere, potentially due to host adaptation to local environments.

Species Identity Dominates over Environment in Driving Bacterial Community Assembly in Wild Invasive Leaf Miners

The microbiota of invasive animal species may be pivotal to their adaptation and spread, yet the processes driving the assembly and potential sources of host-microbiota remain poorly understood. Here, we characterized microbiota of four Liriomyza leaf miner fly species totaling 310 individuals across 43 geographical populations in China and assessed whether the microbiota of the wild leaf miner was acquired from the soil microbiota or the host plant microbiota, using high-throughput 16S rRNA sequencing. Bacterial communities differed significantly among four leaf miner species but did not mirror host phylogeny. Microbiota diversity in the native L. chinensis was significantly higher than in three invasive leaf miners (i.e., L. trifolii, L. huidobrensis, and L. sativae), yet the microbial community of the invasive species exhibited a more connected and complex network structure. Structural equation models revealed that host species identity was more important than environmental factors (e.g., geography, climate, or plants) in shaping microbiota composition. Using neutral and null model analyses, we found that deterministic processes like variable selection played a primary role in driving microbial community assembly, with some influence by stochastic processes like drift. The relative degree of these processes governing microbiota was likely correlated with host species but independent of either geographical or climatic factors. Finally, source tracking analysis showed that leaf miners might acquire microbes from their host plant rather than the soil. Our results provide a robust assessment of the ecological processes governing bacterial community assembly and potential sources of microbes in invasive leaf miners. IMPORTANCE The invasion of foreign species, including leaf miners, is a major threat to world biota. Host-associated microbiota may facilitate host adaption and expansion in a variety of ways. Thus, understanding the processes that drive leaf miner microbiota assembly is imperative for better management of invasive species. However, how microbial communities assemble during the leaf miner invasions and how predictable the processes remain unexplored. This work quantitatively deciphers the relative importance of deterministic process and stochastic process in governing the assembly of four leaf miner microbiotas and identifies potential sources of leaf miner-colonizing microbes from the soil-plant-leaf miner continuum. Our study provides new insights into the mechanisms underlying the drive of leaf miner microbiota assembly.

Host, microbiome, and complex space: applying population and landscape genetic approaches to gut microbiome research in wild populations

In recent years, emerging sequencing technologies and computational tools have driven a tidal wave of research on host-associated microbiomes, particularly the gut microbiome. These studies demonstrate numerous connections between the gut microbiome and vital host functions, primarily in humans, model organisms, and domestic animals. As the adaptive importance of the gut microbiome becomes clearer, interest in studying the gut microbiomes of wild populations has increased, in part due to the potential for discovering conservation applications. The study of wildlife gut microbiomes holds many new challenges and opportunities due to the complex genetic, spatial, and environmental structure of wild host populations, and the potential for these factors to interact with the microbiome. The emerging picture of adaptive coevolution in host-microbiome relationships highlights the importance of understanding microbiome variation in the context of host population genetics and landscape heterogeneity across a wide range of host populations. We propose a conceptual framework for understanding wildlife gut microbiomes in relation to landscape variables and host population genetics, including the potential of approaches derived from landscape genetics. We use this framework to review current research, synthesize important trends, highlight implications for conservation, and recommend future directions for research. Specifically, we focus on how spatial structure and environmental variation interact with host population genetics and microbiome variation in natural populations, and what we can learn from how these patterns of covariation differ depending on host ecological and evolutionary traits.

Emergent spatial patterns of coexistence in species-rich plant communities

Statistical physics has proved essential to analyze multiagent environments. Motivated by the empirical observation of various nonequilibrium features in Barro Colorado and other ecological systems, we analyze a plant-species abundance model of neutral competition, presenting analytical evidence of scale-invariant plant clusters and nontrivial emergent modular correlations. Such first theoretical confirmation of a scale-invariant region, based on percolation processes, reproduces the key features in natural rainforest ecosystems and can confer the most stable equilibrium for ecosystems with vast biodiversity.

Cross-Scale Analyses of Animal and Human Gut Microbiome Assemblies from Metacommunity to Global Landscape

Animal (human) gut microbiomes have been coevolving with their hosts for many millions of years. Understanding how the coevolution shapes the processes of microbiome assembly and diversity maintenance is important but rather challenging. An effort may start with the understanding of how and why animals and humans may differ in their microbiome neutrality (stochasticity) levels. Here, we attempted to perform layered comparative stochasticity analyses across animal species (including humans), class, and kingdom scales, corresponding to microbial metacommunity, landscape, and global-landscape scales. By analyzing 4,903 microbiome samples from 274 animal species covering 4 major invertebrate classes and all 6 vertebrate classes and including 1,787 human gut microbiome samples, we discovered the following: (i) at the microbial metacommunity (animal species) scale, although the general trend of stochasticity (measured in the relationships between fundamental biodiversity/dispersal numbers of Hubbell’s neutral theory and host species phylogenetic timeline) seems continuous, there seems to be a turning point from animals to humans in the passing rate of neutrality tests (12% to 45% versus 100%). We postulate that it should be the human experiences from agricultural/industrial activities (e.g., diet effects) and frequent social/familial contacts that are responsible for the dramatically rising stochastic neutrality in human gut microbiomes. (ii) At the microbial landscape (animal class) and global landscape (animal kingdom) scales, neutrality is not detectable, suggesting that the landscape is niche differentiated—animal species may possess “home niches” for their coadapted microbiomes. We further analyze the reliabilities of our findings by using variable P value thresholds (type I error) and performing power analysis (type II error) of neutrality tests.

IMPORTANCE Understanding how the coevolution (evolutionary time scale) and/or the interactions (ecological time scale) between animal (human) gut microbiomes and their hosts shape the processes of the microbiome assembly and diversity maintenance is important but rather challenging. An effort may start with the understanding of how and why animals and humans may differ in their microbiome neutrality (stochasticity) levels. Here, we attempted to perform layered comparative stochasticity analyses across animal species (including humans), class, and kingdom scales, corresponding to microbial metacommunity, landscape, and global-landscape scales by analyzing 4,903 microbiome samples from 274 animal species covering 4 major invertebrate classes and all 6 vertebrate classes, and including 1,787 human gut microbiome samples. The analyses were implemented by fitting the multisite neutral model and further augmented by checking false-positive and false-negative errors, respectively. It appears that there is a turning (tipping) point in the neutrality level from animal to human microbiomes.

Towards an ecosystem model of infectious disease

Increasingly intimate associations between human society and the natural environment are driving the emergence of novel pathogens, with devastating consequences for humans and animals alike. Prior to emergence, these pathogens exist within complex ecological systems that are characterized by trophic interactions between parasites, their hosts and the environment. Predicting how disturbance to these ecological systems places people and animals at risk from emerging pathogens-and the best ways to manage this-remains a significant challenge. Predictive systems ecology models are powerful tools for the reconstruction of ecosystem function but have yet to be considered for modelling infectious disease. Part of this stems from a mistaken tendency to forget about the role that pathogens play in structuring the abundance and interactions of the free-living species favoured by systems ecologists. Here, we explore how developing and applying these more complete systems ecology models at a landscape scale would greatly enhance our understanding of the reciprocal interactions between parasites, pathogens and the environment, placing zoonoses in an ecological context, while identifying key variables and simplifying assumptions that underly pathogen host switching and animal-to-human spillover risk. As well as transforming our understanding of disease ecology, this would also allow us to better direct resources in preparation for future pandemics.

Global climate change, diet, and the complex relationship between human host and microbiome: Towards an integrated picture

Dietary changes can alter the human microbiome with potential detrimental consequences for health. Given that environment, health, and evolution are interconnected, we ask: Could diet-driven microbiome perturbations have consequences that extend beyond their immediate impact on human health? We address this question in the context of the urgent health challenges posed by global climate change. Drawing on recent studies, we propose that not only can diet-driven microbiome changes lead to dysbiosis, they can also shape life-history traits and fuel human evolution. We posit that dietary shifts prompt mismatched microbiome-host genetics configurations that modulate human longevity and reproductive success. These mismatches can also induce a heritable intra-holobiont stress response, which encourages the holobiont to re-establish equilibrium within the changed nutritional environment. Thus, while mismatches between climate change-related genetic and epigenetic configurations within the holobiont increase the risk and severity of diseases, they may also affect life-history traits and facilitate adaptive responses. These propositions form a framework that can help systematize and address climate-related dietary challenges for policy and health interventions.

A theoretic approach to the mode of gut microbiome translocation in SIV-infected Asian macaques

Human gut microbiome could translocate to other tissues, and the relocation triggered by HIV/SIV infection has received increasing attention. However, the underlying mode of this translocation, whether it is deterministic or random (passive) process, is not clear, not to mention quantitative estimation of the relocation probability and rates. Using multi-tissue microbiome datasets collected from SIV-infected macaques, originally reported by Klase et al. (2015), we apply Hubbell's unified neutral theory of biodiversity (UNTB) implemented by Harris et al. (2017) in the form of multi-site neutral (MSN) model to explore the translocation mode and rates of the gut microbiome. We found that (i) The translocation from gastrointestinal tract to tissues was driven by stochastic (neutral) forces as revealed by 100% neutrality-passing rates with MSN testing; (ii) The translocation probability from gastrointestinal tract to tissues is significantly larger than the baseline dispersal rates occurring within gastrointestinal tract (0.234 vs. 0.006 at the phylum level, P< 0.001). (iii) Approximately, 23% of phyla and 55% of genera were migrated from gastrointestinal tract to the tissues (liver and mesenteric lymph nodes). Our findings offer the first interpretation of the microbial translocation mode from gastrointestinal tract to tissues, and the first estimates of the translocation probability and level.

Agent Based Models of Polymicrobial Biofilms and the Microbiome-A Review

The human microbiome has been a focus of intense study in recent years. Most of the living organisms comprising the microbiome exist in the form of biofilms on mucosal surfaces lining our digestive, respiratory, and genito-urinary tracts. While health-associated microbiota contribute to digestion, provide essential nutrients, and protect us from pathogens, disturbances due to illness or medical interventions contribute to infections, some that can be fatal. Myriad biological processes influence the make-up of the microbiota, for example: growth, division, death, and production of extracellular polymers (EPS), and metabolites. Inter-species interactions include competition, inhibition, and symbiosis. Computational models are becoming widely used to better understand these interactions. Agent-based modeling is a particularly useful computational approach to implement the various complex interactions in microbial communities when appropriately combined with an experimental approach. In these models, each cell is represented as an autonomous agent with its own set of rules, with different rules for each species. In this review, we will discuss innovations in agent-based modeling of biofilms and the microbiota in the past five years from the biological and mathematical perspectives and discuss how agent-based models can be further utilized to enhance our comprehension of the complex world of polymicrobial biofilms and the microbiome.

The yellow perch (Perca flavescens) microbiome revealed resistance to colonisation mostly associated with neutralism driven by rare taxa under cadmium disturbance

Background

Disentangling the dynamics of microbial interactions within communities improves our comprehension of metacommunity assembly of microbiota during host development and under perturbations. To assess the impact of stochastic variation of neutral processes on microbiota structure and composition under disturbance, two types of microbial habitats, free-living (water), and host-associated (skin and gut) were experimentally exposed to either a constant or gradual selection regime exerted by two sublethal cadmium chloride dosages (CdCl₂). Yellow Perch (Perca flavescens) was used as a piscivorous ecotoxicological model. Using 16S rDNA gene based metataxonomics, quantitative diversity metrics of water, skin and gut microbial communities were characterized along with development and across experimental conditions.

Results

After 30 days, constant and gradual selection regimes drove a significant alpha diversity increase for both skin and gut microbiota. In the skin, pervasive negative correlations between taxa in both selection regimes in addition to the taxonomic convergence with the environmental bacterial community, suggest a loss of colonisation resistance resulting in the dysbiosis of yellow perch microbiota. Furthermore, the network connectivity in gut microbiome was exclusively maintained by rare (low abundance) OTUs, while most abundant OTUs were mainly composed of opportunistic invaders such as Mycoplasma and other genera related to fish pathogens such as Flavobacterium. Finally, the mathematical modelling of community assembly using both non-linear least squares models (NLS) based estimates of migration rates and normalized stochasticity ratios (NST) based beta-diversity distances suggested neutral processes drove by taxonomic drift in host and water communities for almost all treatments. The NLS models predicted higher demographic stochasticity in the cadmium-free host and water microbiomes, however, NST models suggested higher ecological stochasticity under perturbations.

Conclusions

Neutral models agree that water and host-microbiota assembly promoted by rare taxa have evolved predominantly under neutral processes with potential involvement of deterministic forces sourced from host filtering and cadmium selection. The early signals of perturbations in the skin microbiome revealed antagonistic interactions by a preponderance of negative correlations in the co-abundance networks. Our findings enhance our understanding of community assembly host-associated and free-living under anthropogenic selective pressure.

Stochastic colonization of hosts with a finite lifespan can drive individual host microbes out of equilibrium

Macroorganisms are inhabited by microbial communities that often change through the lifespan of an individual. One of the factors contributing to this change is colonization from the environment. The colonization of initially microbe-free hosts is particularly interesting, as their microbiome depends entirely on microbes of external origin. We present a mathematical model of this process with a particular emphasis on the effect of ecological drift and a finite host lifespan. Our results indicate the host lifespan becomes especially relevant for short-living organisms (e.g. Caenorhabditis elegans, Drosophila melanogaster, and Danio rerio). In this case, alternative microbiome states (often called enterotypes), the coexistence of microbe-free and colonized hosts, and a reduced probability of colonization can be observed in our model. These results unify multiple reported observations around colonization and suggest that no selective or deterministic drivers are necessary to explain them.

Microbial communities are prevalent not only in the environment but also in hosts. Although the drivers of environmental microbiomes have been studied extensively, less is known about the drivers distinguishing a host environment. Recent experimental observations have highlighted the influence of ecological drift in hosts with short lifespan, including model organisms like C. elegans, D. melanogaster and D. rerio. We have developed a theoretical model to study the effect of a finite host lifespan on relevant observables of the microbiome, including the microbial load, probability of colonization of a microbial taxon, and distribution of microbiome composition in a host population. Although we focus on a case free of any selection, our results indicate the possible coexistence of hosts with alternative microbiome composition, and to a larger extent the coexistence of colonized and microbe-free hosts. A quantitative description is provided.

Species Sorting and Neutral Theory Analyses Reveal Archaeal and Bacterial Communities Are Assembled Differently in Hot Springs

Although the recognition of archaea as one of the three kingdoms in the tree of life has been nearly a half-century long, the comparative investigations on their ecological adaptations with bacteria have been limited. The mechanisms of their community assembly and diversity maintenance in hot springs have not been addressed. The mechanistic study is critical not only for understanding the hot-spring microbiome structure and dynamics, but also for shedding light on their evolutionary adaptations. We applied the neutral theory model and species sorting paradigm of metacommunity theory to investigate how hot-spring microbial communities were assembled, how their diversities were maintained, and how the temperature and pH influence these mechanisms. Through rigorous statistical tests based on the neutral theory and species sorting paradigm, we found (i) According to the neutral theory, archaeal and bacterial communities are assembled differently, with stochastic neutral force playing a more significant role in archaeal communities than in bacterial communities (neutrality-rate = 52.9 vs. 15.8%, p-value < 0.05). (ii) Temperature and pH account for rather limited (<10%) variations in hot-spring microbiomes based on the species sorting paradigm. The pH has more significant influences than temperature on archaeal communities, and both pH and temperature have similarly low influences on bacterial community structure. (iii) We postulate that the differences between archaea and bacteria are likely due to the longer evolutionary history and better adaptation of archaea to host spring environments.

Emerging Priorities for Microbiome Research

Microbiome research has increased dramatically in recent years, driven by advances in technology and significant reductions in the cost of analysis. Such research has unlocked a wealth of data, which has yielded tremendous insight into the nature of the microbial communities, including their interactions and effects, both within a host and in an external environment as part of an ecological community. Understanding the role of microbiota, including their dynamic interactions with their hosts and other microbes, can enable the engineering of new diagnostic techniques and interventional strategies that can be used in a diverse spectrum of fields, spanning from ecology and agriculture to medicine and from forensics to exobiology. From June 19-23 in 2017, the NIH and NSF jointly held an Innovation Lab on Quantitative Approaches to Biomedical Data Science Challenges in our Understanding of the Microbiome. This review is inspired by some of the topics that arose as priority areas from this unique, interactive workshop. The goal of this review is to summarize the Innovation Lab's findings by introducing the reader to emerging challenges, exciting potential, and current directions in microbiome research. The review is broken into five key topic areas: (1) interactions between microbes and the human body, (2) evolution and ecology of microbes, including the role played by the environment and microbe-microbe interactions, (3) analytical and mathematical methods currently used in microbiome research, (4) leveraging knowledge of microbial composition and interactions to develop engineering solutions, and (5) interventional approaches and engineered microbiota that may be enabled by selectively altering microbial composition. As such, this review seeks to arm the reader with a broad understanding of the priorities and challenges in microbiome research today and provide inspiration for future investigation and multi-disciplinary collaboration.

Influences of a Prolific Gut Fungus (Zancudomyces culisetae) on Larval and Adult Mosquito (Aedes aegypti)-Associated Microbiota

Adult mosquitoes inherit a bacterial community from larvae via transstadial transmission, an understudied process that may influence host-microbe interactions. Microbes contribute to important host life history traits, and analyzing transmitted microbial communities, the interrelationship between larval and adult-associated microbiota, and factors influencing host-microbe relationships provides targets for research. During its larval stage, the yellow fever mosquito (Aedes aegypti) hosts the trichomycete gut fungus Zancudomyces culisetae, and fungal colonization coincides with environmental perturbations in the digestive tract microecosystem. Natural populations are differentially exposed to fungi, thereby potentially harboring distinct microbiota and experiencing disparate host-microbe interactions. This study's objectives were to characterize larval and initial adult microbiomes, investigate variation in diversity and distribution of microbial communities across individuals, and assess whether larval fungal colonization impacted microbiomes at these developmental stages. Laboratory-based fungal infestation assays, sequencing of 16S rRNA gene amplicons, and bacterial load quantification protocols revealed that initial adult microbiomes varied in diversity and distribution. Larval fungal colonization had downstream effects on initial adult microbiomes, significantly reducing microbial community variation, shifting relative abundances of certain bacterial families, and influencing transstadial transmission outcomes of particular genera. Further, abundances of several families consistently decreased in adults relative to levels in larvae, possibly reflecting impacts of host development on specific bacterial taxa. These findings demonstrated that a prolific gut fungus impacted mosquito-associated microbiota at two developmental stages in an insect connected with global human health.IMPORTANCE Mosquitoes are widespread vectors of numerous human pathogens and harbor microbiota known to affect host phenotypic traits. However, little research has directly investigated how bacterial communities associated with larvae and adults are connected. We characterized whole-body bacterial communities in mosquito larvae preceding pupation and in newly emerged adults, and investigated whether a significant biotic factor, fungal colonization of the larval hindgut, impacted these microbiomes. Results showed that fungal colonization reduced microbial community variation across individuals and differentially impacted the outcomes of transstadial transmission for certain bacterial genera, revealing downstream effects of the fungus on initial adult microbiomes. The importance of our research is in providing a thorough comparative analysis of whole-body microbiota harbored in larvae and adults of the yellow fever mosquito (Aedes aegypti) and in demonstrating the important role a widespread gut fungus played in a host-associated microbiome.

The microbiome of the Melitaea cinxia butterfly shows marked variation but is only little explained by the traits of the butterfly or its host plant

Understanding of the ecological factors that shape intraspecific variation of insect microbiota in natural populations is relatively poor. In Lepidopteran caterpillars, microbiota is assumed to be mainly composed of transient bacterial symbionts acquired from the host plant. We sampled Glanville fritillary (Melitaea cinxia) caterpillars from natural populations to describe their gut microbiome and to identify potential ecological factors that determine its structure. Our results demonstrate high variability of microbiota composition even among caterpillars that shared the same host plant individual and most likely the same genetic background. We observed that the caterpillars harboured microbial classes that varied among individuals and alternated between two distinct communities (one composed of mainly Enterobacteriaceae and another with more variable microbiota community). Even though the general structure of the microbiota was not attributed to the measured ecological factors, we found that phylogenetically similar microbiota showed corresponding responses to the sex and the parasitoid infection of the caterpillar and to those of the host plant's microbial and chemical composition. Our results indicate high among-individual variability in the microbiota of the M. cinxia caterpillar and contradict previous findings that the host plant is the major driver of the microbiota communities of insect herbivores.

Motion, fixation probability and the choice of an evolutionary process

Seemingly minor details of mathematical and computational models of evolution are known to change the effect of population structure on the outcome of evolutionary processes. For example, birth-death dynamics often result in amplification of selection, while death-birth processes have been associated with suppression. In many biological populations the interaction structure is not static. Instead, members of the population are in motion and can interact with different individuals at different times. In this work we study populations embedded in a flowing medium; the interaction network is then time dependent. We use computer simulations to investigate how this dynamic structure affects the success of invading mutants, and compare these effects for different coupled birth and death processes. Specifically, we show how the speed of the motion impacts the fixation probability of an invading mutant. Flows of different speeds interpolate between evolutionary dynamics on fixed heterogeneous graphs and well-stirred populations; this allows us to systematically compare against known results for static structured populations. We find that motion has an active role in amplifying or suppressing selection by fragmenting and reconnecting the interaction graph. While increasing flow speeds suppress selection for most evolutionary models, we identify characteristic responses to flow for the different update rules we test. In particular we find that selection can be maximally enhanced or suppressed at intermediate flow speeds.

Whether a mutation spreads in a population or not is one of the most important questions in biology. The evolution of cancer and antibiotic resistance, for example, are mediated by invading mutants. Recent work has shown that population structure can have important consequences for the outcome of evolution. For instance, a mutant can have a higher or a lower chance of invasion than in unstructured populations. These effects can depend on seemingly minor details of the evolutionary model, such as the order of birth and death events. Many biological populations are in motion, for example due to external stirring. Experimentally this is known to be important; the performance of mutants in E. coli populations, for example, depends on the rate of mixing. Here, we focus on simulations of populations in a flowing medium, and compare the success of a mutant for different flow speeds. We contrast different evolutionary models, and identify what features of the evolutionary model affect mutant success for different speeds of the flow. We find that the chance of mutant invasion can be at its highest (or lowest) at intermediate flow speeds, depending on the order in which birth and death events occur in the evolutionary process.

Processes shaping gut microbiota diversity in allopatric populations of the endemic lizard Podarcis lilfordi from Menorcan islets (Balearic Islands)

Compositional variation of the gut microbiota across host allopatric populations can reflect both adaptation and stochasticity since the time of separation. Major factors shaping this variation include the host phylogeographic and demographic history, the microbiota inheritance, environmental inputs and dispersal of bacteria. Here we explored the impact of these factors in driving gut community diversity in seven allopatric populations of the omnivorous lizard Podarcis lilfordi from the Menorcan coastal islets, all descending from an ancestral mainland population. Using 16S rRNA Illumina sequencing, we showed that 'islet' and 'age' (time since islet separation from mainland) were the only significant variables in microbial community clustering, suggesting a partial islet-restricted diversification following these lizards phylogeography. Despite a significant variation, islets/populations were characterized by a remarkably low bacterial uniqueness (2.4% of total OTUs) and a minor differential enrichment of taxa, indicating a negligible impact of local inputs and important host common constraints. Overall, the extant pattern of similarity/dissimilarity among islets is compatible with partial retention of the ancestral mainland microbial pool, with differences among islets potentially explained by a differential loss of bacteria following population fragmentation and bottlenecks (i.e. ecological drift). While more quantitative data are needed to validate this hypothesis, this study unveils the importance of considering both neutral and niche-driven processes in driving contemporary patterns of gut metacommunity diversity.

Synthetic Mutualism and the Intervention Dilemma

Ecosystems are complex networks of interacting individuals co-evolving with their environment. As such, changes to an interaction can influence the whole ecosystem. However, to predict the outcome of these changes, considerable understanding of processes driving the system is required. Synthetic biology provides powerful tools to aid this understanding, but these developments also allow us to change specific interactions. Of particular interest is the ecological importance of mutualism, a subset of cooperative interactions. Mutualism occurs when individuals of different species provide a reciprocal fitness benefit. We review available experimental techniques of synthetic biology focused on engineered synthetic mutualistic systems. Components of these systems have defined interactions that can be altered to model naturally occurring relationships. Integrations between experimental systems and theoretical models, each informing the use or development of the other, allow predictions to be made about the nature of complex relationships. The predictions range from stability of microbial communities in extreme environments to the collapse of ecosystems due to dangerous levels of human intervention. With such caveats, we evaluate the promise of synthetic biology from the perspective of ethics and laws regarding biological alterations, whether on Earth or beyond. Just because we are able to change something, should we?

Competitive lottery-based assembly of selected clades in the human gut microbiome

BACKGROUND

While the composition of the gut microbiome has now been well described by several large-scale studies, models that can account for the range of microbiome compositions that have been observed are still lacking. One model that has been well studied in macro communities and that could be useful for understanding microbiome assembly is the competitive lottery model. This model posits that groups of organisms from a regional pool of species are able to colonize the same niche and that the first species to arrive will take over the entire niche, excluding other group members.

RESULTS

Here, we examined whether this model also plays a role in the assembly of the human gut microbiome, defining measures to identify groups of organisms whose distribution across samples conforms to the competitive lottery schema. Applying this model to multiple datasets with thousands of human gut microbiome samples, we identified several taxonomic groups that exhibit a lottery-like distribution, including the Akkermansia, Dialister, and Phascolarctobacterium genera. We validated that these groups exhibit lottery-like assembly in multiple independent microbiome datasets confirming that this assembly schema is universal and not cohort specific. Examining the distribution of species from these groups in the gut microbiome of developing infants, we found that the initial lottery winner can be replaced by a different member of the group. We further found that species from lottery-like groups tend to have fewer genes in their genomes, suggesting more specialized species that are less able to engage in niche differentiation.

CONCLUSIONS

Combined, our findings highlight the complex and dynamic process through which microbial communities assemble and suggest that different phylogenetic groups may follow different models during this process.

Is Host Filtering the Main Driver of Phylosymbiosis across the Tree of Life?

Host-associated microbiota composition can be conserved over evolutionary time scales. Indeed, closely related species often host similar microbiota; i.e., the composition of their microbiota harbors a phylogenetic signal, a pattern sometimes referred to as "phylosymbiosis." Elucidating the origins of this pattern is important to better understand microbiota ecology and evolution. However, this is hampered by our lack of theoretical expectations and a comprehensive overview of phylosymbiosis prevalence in nature. Here, we use simulations to provide a simple expectation for when we should expect this pattern to occur and then review the literature to document the prevalence and strength of phylosymbiosis across the host tree of life. We demonstrate that phylosymbiosis can readily emerge from a simple ecological filtering process, whereby a given host trait (e.g., gut pH) that varies with host phylogeny (i.e., harbors a phylogenetic signal) filters preadapted microbes. We found marked differences between methods used to detect phylosymbiosis, so we proposed a series of practical recommendations based on using multiple best-performing approaches. Importantly, we found that, while the prevalence of phylosymbiosis is mixed in nature, it appears to be stronger for microbiotas living in internal host compartments (e.g., the gut) than those living in external compartments (e.g., the rhizosphere). We show that phylosymbiosis can theoretically emerge without any intimate, long-term coevolutionary mechanisms and that most phylosymbiosis patterns observed in nature are compatible with a simple ecological process. Deviations from baseline ecological expectations might be used to further explore more complex hypotheses, such as codiversification. IMPORTANCE Phylosymbiosis is a pattern defined as the tendency of closely related species to host microbiota whose compositions resemble each other more than host species drawn at random from the same tree. Understanding the mechanisms behind phylosymbiosis is important because it can shed light on rules governing the assembly of host-associated microbiotas and, potentially, their coevolutionary dynamics with hosts. For example, is phylosymbiosis a result of coevolution, or can it be generated by simple ecological filtering processes? Beyond qualitative theoretical models, quantitative theoretical expectations can provide new insights. For example, deviations from a simple baseline of ecological filtering may be used to test more-complex hypotheses (e.g., coevolution). Here, we use simulations to provide evidence that simple host-related ecological filtering can readily generate phylosymbiosis, and we contrast these predictions with real-world data. We find that while phylosymbiosis is widespread in nature, phylosymbiosis patterns are compatible with a simple ecological model in the majority of taxa. Internal compartments of hosts, such as the animal gut, often display stronger phylosymbiosis than expected from a purely ecological filtering process, suggesting that other mechanisms are also involved.

Testing the neutral theory of biodiversity with the microbiome dataset from cystic fibrosis patients

Cystic fibrosis (CF) is a hereditary disease that is characterized by defective mucociliary clearance, airway obstruction, chronic infection, and persistent inflammation. Cystic fibrosis pulmonary exacerbation (CFPE) majorly causes the morbidity of CF patients. Although CF has been demonstrated to change the composition of lung microbial community, previous studies have not made efforts to study the differences in the mechanism of assembly and diversity maintenance of lung microbial community in CF patients. In this study, we applied the neutral theory of biodiversity to comparatively investigate the assembly and diversity maintenance of the lung microbial community before and after the antibiotic treatment by reanalyzing the dataset from Fodor et al's study. We found that no one sample in the lung microbial communities of the sputum samples of Exacerbation group, nor those of End-of-treatment group satisfied the predictions of neutral model, suggesting that the neutral-process does not dominate in CF patients before and after antibiotic treatments. By comparing the biodiversity parameter between Exacerbation and End-of-treatment group, we found that the former had the significantly higher biodiversity, but the change in diversity parameter is slight and the P value is close to.05 (P value = .41). Therefore, our second finding is that although CFPE may increase the biodiversity of lung microbial community, the change is not essential.

Neutral models of short-term microbiome dynamics with host subpopulation structure and migration limitation

BACKGROUND

Most empirical studies tend to focus on microbiome dynamics within hosts or microbiome compositional differences between hosts over short periods. However, there is still a dearth of formal models that allow us to investigate the observed short-term dynamics of microbiomes under a unified ecological and evolutionary framework. In our previous study, we developed a computational agent-based neutral framework that simulates microbiome dynamics spanning many host generations with the added dimension of a genealogy of hosts. Although this long-term framework revealed interesting microbial diversity patterns under a simple but plausible evolutionary process and provided a platform for future elaboration of more complex systems, it does not allow us to explore microbiome dynamics within a single host generation.

METHODS

In this paper, we developed a computational, agent-based, forward-time framework of microbiome dynamics within a single host generation. As we have done under our neutral long-term models, we incorporate neutral processes of environmental microbiome assembly and microbe acquisition from parents and environment. We also incorporate a Moran genealogical model of hosts, so that the dynamics of microbiome evolution can be studied within a single host generation. Furthermore, we allow host subpopulation structure and host migration to affect microbiome recruitment.

RESULTS

We show that microbiome diversity within hosts increases monotonically with increases in environmental contribution, while microbiome diversity between hosts increases with increasing parental inheritance. Host population division and dispersal limitation under high host contribution further shaped the patterns by elevating microbiome differences between hosts and depressing microbial diversity within hosts. Microbiome diversity within the whole population showed strong temporal stability regardless of the modes of microbiome acquisition and subpopulation structures.

CONCLUSIONS

We present a computational framework that integrates various processes including host genealogy, microbe recruitment, and host dispersal limitation acting on the short-term dynamics of microbiomes. Our framework demonstrates that the neutral dynamics of microbiomes within a population of hosts is strongly influenced by transmission mode and shared environment.

Allometry and Ecology of the Bilaterian Gut Microbiome

Classical ecology provides principles for construction and function of biological communities, but to what extent these apply to the animal-associated microbiota is just beginning to be assessed. Here, we investigated the influence of several well-known ecological principles on animal-associated microbiota by characterizing gut microbial specimens from bilaterally symmetrical animals (Bilateria) ranging from flies to whales. A rigorously vetted sample set containing 265 specimens from 64 species was assembled. Bacterial lineages were characterized by 16S rRNA gene sequencing. Previously published samples were also compared, allowing analysis of over 1,098 samples in total. A restricted number of bacterial phyla was found to account for the great majority of gut colonists. Gut microbial composition was associated with host phylogeny and diet. We identified numerous gut bacterial 16S rRNA gene sequences that diverged deeply from previously studied taxa, identifying opportunities to discover new bacterial types. The number of bacterial lineages per gut sample was positively associated with animal mass, paralleling known species-area relationships from island biogeography and implicating body size as a determinant of community stability and niche complexity. Samples from larger animals harbored greater numbers of anaerobic communities, specifying a mechanism for generating more-complex microbial environments. Predictions for species/abundance relationships from models of neutral colonization did not match the data set, pointing to alternative mechanisms such as selection of specific colonists by environmental niche. Taken together, the data suggest that niche complexity increases with gut size and that niche selection forces dominate gut community construction.IMPORTANCE The intestinal microbiome of animals is essential for health, contributing to digestion of foods, proper immune development, inhibition of pathogen colonization, and catabolism of xenobiotic compounds. How these communities assemble and persist is just beginning to be investigated. Here we interrogated a set of gut samples from a wide range of animals to investigate the roles of selection and random processes in microbial community construction. We show that the numbers of bacterial species increased with the weight of host organisms, paralleling findings from studies of island biogeography. Communities in larger organisms tended to be more anaerobic, suggesting one mechanism for niche diversification. Nonselective processes enable specific predictions for community structure, but our samples did not match the predictions of the neutral model. Thus, these findings highlight the importance of niche selection in community construction and suggest mechanisms of niche diversification.

The human microbiome in evolution

The trillions of microbes living in the gut-the gut microbiota-play an important role in human biology and disease. While much has been done to explore its diversity, a full understanding of our microbiomes demands an evolutionary perspective. In this review, we compare microbiomes from human populations, placing them in the context of microbes from humanity's near and distant animal relatives. We discuss potential mechanisms to generate host-specific microbiome configurations and the consequences of disrupting those configurations. Finally, we propose that this broader phylogenetic perspective is useful for understanding the mechanisms underlying human-microbiome interactions.

Models of microbiome evolution incorporating host and microbial selection

BACKGROUND

Numerous empirical studies suggest that hosts and microbes exert reciprocal selective effects on their ecological partners. Nonetheless, we still lack an explicit framework to model the dynamics of both hosts and microbes under selection. In a previous study, we developed an agent-based forward-time computational framework to simulate the neutral evolution of host-associated microbial communities in a constant-sized, unstructured population of hosts. These neutral models allowed offspring to sample microbes randomly from parents and/or from the environment. Additionally, the environmental pool of available microbes was constituted by fixed and persistent microbial OTUs and by contributions from host individuals in the preceding generation.

METHODS

In this paper, we extend our neutral models to allow selection to operate on both hosts and microbes. We do this by constructing a phenome for each microbial OTU consisting of a sample of traits that influence host and microbial fitnesses independently. Microbial traits can influence the fitness of hosts ("host selection") and the fitness of microbes ("trait-mediated microbial selection"). Additionally, the fitness effects of traits on microbes can be modified by their hosts ("host-mediated microbial selection"). We simulate the effects of these three types of selection, individually or in combination, on microbiome diversities and the fitnesses of hosts and microbes over several thousand generations of hosts.

RESULTS

We show that microbiome diversity is strongly influenced by selection acting on microbes. Selection acting on hosts only influences microbiome diversity when there is near-complete direct or indirect parental contribution to the microbiomes of offspring. Unsurprisingly, microbial fitness increases under microbial selection. Interestingly, when host selection operates, host fitness only increases under two conditions: (1) when there is a strong parental contribution to microbial communities or (2) in the absence of a strong parental contribution, when host-mediated selection acts on microbes concomitantly.

CONCLUSIONS

We present a computational framework that integrates different selective processes acting on the evolution of microbiomes. Our framework demonstrates that selection acting on microbes can have a strong effect on microbial diversities and fitnesses, whereas selection on hosts can have weaker outcomes.

Dynamic microbiome evolution in social bees

The highly social (eusocial) corbiculate bees, comprising the honey bees, bumble bees, and stingless bees, are ubiquitous insect pollinators that fulfill critical roles in ecosystem services and human agriculture. Here, we conduct wide sampling across the phylogeny of these corbiculate bees and reveal a dynamic evolutionary history behind their microbiota, marked by multiple gains and losses of gut associates, the presence of generalist as well as host-specific strains, and patterns of diversification driven, in part, by host ecology (for example, colony size). Across four continents, we found that different host species have distinct gut communities, largely independent of geography or sympatry. Nonetheless, their microbiota has a shared heritage: The emergence of the eusocial corbiculate bees from solitary ancestors appears to coincide with the acquisition of five core gut bacterial lineages, supporting the hypothesis that host sociality facilitates the development and maintenance of specialized microbiomes.

Making a microbiome: the many determinants of host-associated microbial community composition

The composition of many host-associated microbial communities is characterized by seemingly contradictory features: strong selection for specific taxa by the host, but substantial variability among hosts and over time within one host. Recent advances have revealed that both deterministic and stochastic processes operating across multiple spatial scales shape the composition of host-associated microbial communities. Although most research has focused on deterministic processes within individual hosts, the microbiota within each host is increasingly recognized to contribute to a wider metacommunity maintained by transmission between individual hosts and dispersal between host-associated and free-living microbial communities. By applying a community ecology perspective encompassing the microbial metacommunity, the many determinants of host-associated microbial community composition can be identified, guiding the directions of future research.

Testing the Neutral Theory of Biodiversity with Human Microbiome Datasets

The human microbiome project (HMP) has made it possible to test important ecological theories for arguably the most important ecosystem to human health-the human microbiome. Existing limited number of studies have reported conflicting evidence in the case of the neutral theory; the present study aims to comprehensively test the neutral theory with extensive HMP datasets covering all five major body sites inhabited by the human microbiome. Utilizing 7437 datasets of bacterial community samples, we discovered that only 49 communities (less than 1%) satisfied the neutral theory, and concluded that human microbial communities are not neutral in general. The 49 positive cases, although only a tiny minority, do demonstrate the existence of neutral processes. We realize that the traditional doctrine of microbial biogeography "Everything is everywhere, but the environment selects" first proposed by Baas-Becking resolves the apparent contradiction. The first part of Baas-Becking doctrine states that microbes are not dispersal-limited and therefore are neutral prone, and the second part reiterates that the freely dispersed microbes must endure selection by the environment. Therefore, in most cases, it is the host environment that ultimately shapes the community assembly and tip the human microbiome to niche regime.

Effects of dispersal limitation in the face of intense selection via dietary intervention on the faecal microbiota of rats

Four ecological processes mediate microbial diversity: selection, whereby host factors favour the persistence of particular microbial taxa; historical contingency, where differences in the timing and order of microbial acquisition results in different communities; stochastic factors that impact community assemblages; and dispersal limitation, where the presence of particular taxa is restricted by host population structure and local environmental factors. However, few studies have explored the impact of selection through diet modification on the differences in microbial community composition arising from dispersal limitation. At weaning, the faecal microbiota of 45 rats originating from six litters, as assessed by 16S rRNA gene sequencing, was strongly correlated with an animal's litter membership. Following 14 weeks on one of three the experimental diets, the faecal microbiota of the rats were again characterized. Clear effects of diet on microbial community composition were observed. However, after 14 weeks of dietary intervention, the effect of litter membership on microbial community composition was still significant. Our results demonstrate that intense selection (diet manipulation) cannot completely eliminate differences in microbial communities that occurred as a consequence of dispersal limitation (litter membership) prior to selection. These results have clear implications for efforts that attempt to achieve positive health outcomes through diet manipulation.

Glycomics: revealing the dynamic ecology and evolution of sugar molecules

Sugars are the most functionally and structurally diverse molecules in the biological world. Glycan structures range from tiny single monosaccharide units to giant chains thousands of units long. Some glycans are branched, their monosaccharides linked together in many different combinations and orientations. Some exist as solitary molecules; others are conjugated to proteins and lipids and alter their collective functional properties. In addition to structural and storage roles, glycan molecules participate in and actively regulate physiological and developmental processes. Glycans also mediate cellular interactions within and between individuals. Their roles in ecology and evolution are pivotal, but not well studied because glycan biochemistry requires different methods than standard molecular biology practice. The properties of glycans are in some ways convenient, and in others challenging. Glycans vary on organismal timescales, and in direct response to physiological and ecological conditions. Their mature structures are physical records of both genetic and environmental influences during maturation. We describe the scope of natural glycan variation and discuss how studying glycans will allow researchers to further integrate the fields of ecology and evolution.