New sector morphologies emerge from anisotropic colony growth.

Competition during range expansions is of great interest from both practical and theoretical view points. Experimentally, range expansions are often studied in homogeneous Petri dishes, which lack spatial anisotropy that might be present in realistic populations. Here, we analyze a model of anisotropic growth, based on coupled Kardar-Parisi-Zhang and Fisher-Kolmogorov-Petrovsky-Piskunov equations that describe surface growth and lateral competition. Compared to a previous study of isotropic growth, anisotropy relaxes a constraint between parameters of the model. We completely characterize spatial patterns and invasion velocities in this generalized model. In particular, we find that strong anisotropy results in a distinct morphology of spatial invasion with a kink in the displaced strain ahead of the boundary between the strains. This morphology of the out-competed strain is similar to a shock wave and serves as a signature of anisotropic growth.

Mechanically-driven growth and competition in a Voronoi model of tissues.

The mechanisms leading cells to acquire a fitness advantage and establish themselves in a population are paramount to understanding the development and growth of cancer. Although there are many works that study separately either the evolutionary dynamics or the mechanics of cancer, little has been done to couple evolutionary dynamics to mechanics. To address this question, we study a confluent model of tissue using a Self-Propelled Voronoi (SPV) model with stochastic growth rates that depend on the mechanical variables of the system. The SPV model is an out-of-equilibrium model of tissue derived from an energy functional that has a jamming/unjamming transition between solid-like and liquid-like states. By considering several scenarios of mutants invading a resident population in both phases, we determine the range of parameters that confer a fitness advantage and show that the preferred area and perimeter are the most relevant ones. We find that the liquid-like state is more resistant to invasion and show that the outcome of the competition can be determined from the simulation of a non-growing mixture. Moreover, a mean-field approximation can accurately predict the fate of a mutation affecting mechanical properties of a cell. Our results can be used to infer evolutionary dynamics from tissue images, understand cancer-suppressing effects of tissue mechanics, and even search for mechanics-based therapies.

A mechanistic statistical approach to infer invasion characteristics of human-dispersed species with complex life cycle.

The rising introduction of invasive species through trade networks threatens biodiversity and ecosystem services. Yet, we have a limited understanding of how transportation networks determine patterns of range expansion. This is partly because current analytical models fail to integrate the invader's life-history dynamics with heterogeneity in human-mediated dispersal patterns. And partly because classical statistical methods often fail to provide reliable estimates of model parameters due to spatial biases in the presence-only records and lack of informative demographic data. To address these gaps, we first formulate an age-structured metapopulation model that uses a probability matrix to emulate human-mediated dispersal patterns. The model reveals that an invader spreads along the shortest network path, such that the inter-patch network distances decrease with increasing traffic volume and reproductive value of hitchhikers. Next, we propose a Bayesian statistical method to estimate model parameters using presence-only data and prior demographic knowledge. To show the utility of the statistical approach, we analyze zebra mussel (Dreissena polymorpha) expansion in North America through the commercial shipping network. Our analysis underscores the importance of correcting spatial biases and leveraging priors to answer questions, such as where and when the zebra mussels were introduced and what life-history characteristics make these mollusks successful invaders.

Interplay between morphology and competition in two-dimensional colony expansion.

In growing populations, the fate of mutations depends on their competitive ability against the ancestor and their ability to colonize new territory. Here we present a theory that integrates both aspects of mutant fitness by coupling the classic description of one-dimensional competition (Fisher equation) to the minimal model of front shape (Kardar-Parisi-Zhang equation). We solve these equations and find three regimes, which are controlled solely by the expansion rates, solely by the competitive abilities, or by both. Collectively, our results provide a simple framework to study spatial competition.

Inferring microbial co-occurrence networks from amplicon data: a systematic evaluation.

Microbes commonly organize into communities consisting of hundreds of species involved in complex interactions with each other. 16S ribosomal RNA (16S rRNA) amplicon profiling provides snapshots that reveal the phylogenies and abundance profiles of these microbial communities. These snapshots, when collected from multiple samples, can reveal the co-occurrence of microbes, providing a glimpse into the network of associations in these communities. However, the inference of networks from 16S data involves numerous steps, each requiring specific tools and parameter choices. Moreover, the extent to which these steps affect the final network is still unclear. In this study, we perform a meticulous analysis of each step of a pipeline that can convert 16S sequencing data into a network of microbial associations. Through this process, we map how different choices of algorithms and parameters affect the co-occurrence network and identify the steps that contribute substantially to the variance. We further determine the tools and parameters that generate robust co-occurrence networks and develop consensus network algorithms based on benchmarks with mock and synthetic data sets. The Microbial Co-occurrence Network Explorer, or MiCoNE (available at https://github.com/segrelab/MiCoNE) follows these default tools and parameters and can help explore the outcome of these combinations of choices on the inferred networks. We envisage that this pipeline could be used for integrating multiple data sets and generating comparative analyses and consensus networks that can guide our understanding of microbial community assembly in different biomes. IMPORTANCE Mapping the interrelationships between different species in a microbial community is important for understanding and controlling their structure and function. The surge in the high-throughput sequencing of microbial communities has led to the creation of thousands of data sets containing information about microbial abundances. These abundances can be transformed into co-occurrence networks, providing a glimpse into the associations within microbiomes. However, processing these data sets to obtain co-occurrence information relies on several complex steps, each of which involves numerous choices of tools and corresponding parameters. These multiple options pose questions about the robustness and uniqueness of the inferred networks. In this study, we address this workflow and provide a systematic analysis of how these choices of tools affect the final network and guidelines on appropriate tool selection for a particular data set. We also develop a consensus network algorithm that helps generate more robust co-occurrence networks based on benchmark synthetic data sets.

Competition on the edge of an expanding population.

In growing populations, the fate of mutations depends on their competitive ability against the ancestor and their ability to colonize new territory. Here we present a theory that integrates both aspects of mutant fitness by coupling the classic description of one-dimensional competition (Fisher equation) to the minimal model of front shape (KPZ equation). We solved these equations and found three regimes, which are controlled solely by the expansion rates, solely by the competitive abilities, or by both. Collectively, our results provide a simple framework to study spatial competition.

Ecological landscapes guide the assembly of optimal microbial communities.

Assembling optimal microbial communities is key for various applications in biofuel production, agriculture, and human health. Finding the optimal community is challenging because the number of possible communities grows exponentially with the number of species, and so an exhaustive search cannot be performed even for a dozen species. A heuristic search that improves community function by adding or removing one species at a time is more practical, but it is unknown whether this strategy can discover an optimal or nearly optimal community. Using consumer-resource models with and without cross-feeding, we investigate how the efficacy of search depends on the distribution of resources, niche overlap, cross-feeding, and other aspects of community ecology. We show that search efficacy is determined by the ruggedness of the appropriately-defined ecological landscape. We identify specific ruggedness measures that are both predictive of search performance and robust to noise and low sampling density. The feasibility of our approach is demonstrated using experimental data from a soil microbial community. Overall, our results establish the conditions necessary for the success of the heuristic search and provide concrete design principles for building high-performing microbial consortia.

Growth instabilities shape morphology and genetic diversity of microbial colonies.

Cellular populations assume an incredible variety of shapes ranging from circular molds to irregular tumors. While we understand many of the mechanisms responsible for these spatial patterns, little is known about how the shape of a population influences its ecology and evolution. Here, we investigate this relationship in the context of microbial colonies grown on hard agar plates. This a well-studied system that exhibits a transition from smooth circular disks to more irregular and rugged shapes as either the nutrient concentration or cellular motility is decreased. Starting from a mechanistic model of colony growth, we identify two dimensionless quantities that determine how morphology and genetic diversity of the population depend on the model parameters. Our simulations further reveal that population dynamics cannot be accurately described by the commonly-used surface growth models. Instead, one has to explicitly account for the emergent growth instabilities and demographic fluctuations. Overall, our work links together environmental conditions, colony morphology, and evolution. This link is essential for a rational design of concrete, biophysical perturbations to steer evolution in the desired direction.

Slow expanders invade by forming dented fronts in microbial colonies.

Most organisms grow in space, whether they are viruses spreading within a host tissue or invasive species colonizing a new continent. Evolution typically selects for higher expansion rates during spatial growth, but it has been suggested that slower expanders can take over under certain conditions. Here, we report an experimental observation of such population dynamics. We demonstrate that mutants that grow slower in isolation nevertheless win in competition, not only when the two types are intermixed, but also when they are spatially segregated into sectors. The latter was thought to be impossible because previous studies focused exclusively on the global competitions mediated by expansion velocities, but overlooked the local competitions at sector boundaries. Local competition, however, can enhance the velocity of either type at the sector boundary and thus alter expansion dynamics. We developed a theory that accounts for both local and global competitions and describes all possible sector shapes. In particular, the theory predicted that a slower on its own, but more competitive, mutant forms a dented V-shaped sector as it takes over the expansion front. Such sectors were indeed observed experimentally, and their shapes matched quantitatively with the theory. In simulations, we further explored several mechanisms that could provide slow expanders with a local competitive advantage and showed that they are all well-described by our theory. Taken together, our results shed light on previously unexplored outcomes of spatial competition and establish a universal framework to understand evolutionary and ecological dynamics in expanding populations.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS).

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

Genealogical structure changes as range expansions transition from pushed to pulled.

Range expansions accelerate evolution through multiple mechanisms, including gene surfing and genetic drift. The inference and control of these evolutionary processes ultimately rely on the information contained in genealogical trees. Currently, there are two opposing views on how range expansions shape genealogies. In invasion biology, expansions are typically approximated by a series of population bottlenecks producing genealogies with only pairwise mergers between lineages-a process known as the Kingman coalescent. Conversely, traveling wave models predict a coalescent with multiple mergers, known as the Bolthausen-Sznitman coalescent. Here, we unify these two approaches and show that expansions can generate an entire spectrum of coalescent topologies. Specifically, we show that tree topology is controlled by growth dynamics at the front and exhibits large differences between pulled and pushed expansions. These differences are explained by the fluctuations in the total number of descendants left by the early founders. High growth cooperativity leads to a narrow distribution of reproductive values and the Kingman coalescent. Conversely, low growth cooperativity results in a broad distribution, whose exponent controls the merger sizes in the genealogies. These broad distribution and non-Kingman tree topologies emerge due to the fluctuations in the front shape and position and do not occur in quasi-deterministic simulations. Overall, our results show that range expansions provide a robust mechanism for generating different types of multiple mergers, which could be similar to those observed in populations with strong selection or high fecundity. Thus, caution should be exercised in making inferences about the origin of non-Kingman genealogies.

Cooperation mitigates diversity loss in a spatially expanding microbial population.

The evolution and potentially even the survival of a spatially expanding population depends on its genetic diversity, which can decrease rapidly due to a serial founder effect. The strength of the founder effect is predicted to depend strongly on the details of the growth dynamics. Here, we probe this dependence experimentally using a single microbial species, Saccharomyces cerevisiae, expanding in multiple environments that induce varying levels of cooperativity during growth. We observe a drastic reduction in diversity during expansions when yeast grows noncooperatively on simple sugars, but almost no loss of diversity when cooperation is required to digest complex metabolites. These results are consistent with theoretical expectations: When cells grow independently from each other, the expansion proceeds as a pulled wave driven by growth at the low-density tip of the expansion front. Such populations lose diversity rapidly because of the strong genetic drift at the expansion edge. In contrast, diversity loss is substantially reduced in pushed waves that arise due to cooperative growth. In such expansions, the low-density tip of the front grows much more slowly and is often reseeded from the genetically diverse population core. Additionally, in both pulled and pushed expansions, we observe a few instances of abrupt changes in allele fractions due to rare fluctuations of the expansion front and show how to distinguish such rapid genetic drift from selective sweeps.

Genetic drift in range expansions is very sensitive to density dependence in dispersal and growth.

Theory predicts rapid genetic drift during invasions, yet many expanding populations maintain high genetic diversity. We find that genetic drift is dramatically suppressed when dispersal rates increase with the population density because many more migrants from the diverse, high-density regions arrive at the expansion edge. When density dependence is weak or negative, the effective population size of the front scales only logarithmically with the carrying capacity. The dependence, however, switches to a sublinear power law and then to a linear increase as the density dependence becomes strongly positive. We develop a unified framework revealing that the transitions between different regimes of diversity loss are controlled by a single, universal quantity: the ratio of the expansion velocity to the geometric mean of dispersal and growth rates at expansion edge. Our results suggest that positive density dependence could dramatically alter evolution in expanding populations even when its contribution to the expansion velocity is small.

Evolution at the Edge of Expanding Populations.

Predicting the evolution of expanding populations is critical to controlling biological threats such as invasive species and cancer metastasis. Expansion is primarily driven by reproduction and dispersal, but nature abounds with examples of evolution where organisms pay a reproductive cost to disperse faster. When does selection favor this "survival of the fastest"? We searched for a simple rule, motivated by evolution experiments where swarming bacteria evolved into a hyperswarmer mutant that disperses ∼100% faster but pays a growth cost of ∼10% to make many copies of its flagellum. We analyzed a two-species model based on the Fisher equation to explain this observation: the population expansion rate (v) results from an interplay of growth (r) and dispersal (D) and is independent of the carrying capacity: v=2(rD)1/2 . A mutant can take over the edge only if its expansion rate (v2) exceeds the expansion rate of the established species (v1); this simple condition ( v2>v1 ) determines the maximum cost in slower growth that a faster mutant can pay and still be able to take over. Numerical simulations and time-course experiments where we tracked evolution by imaging bacteria suggest that our findings are general: less favorable conditions delay but do not entirely prevent the success of the fastest. Thus, the expansion rate defines a traveling wave fitness, which could be combined with trade-offs to predict evolution of expanding populations.

Pinned, locked, pushed, and pulled traveling waves in structured environments.

Traveling fronts describe the transition between two alternative states in a great number of physical and biological systems. Examples include the spread of beneficial mutations, chemical reactions, and the invasions by foreign species. In homogeneous environments, the alternative states are separated by a smooth front moving at a constant velocity. This simple picture can break down in structured environments such as tissues, patchy landscapes, and microfluidic devices. Habitat fragmentation can pin the front at a particular location or lock invasion velocities into specific values. Locked velocities are not sensitive to moderate changes in dispersal or growth and are determined by the spatial and temporal periodicity of the environment. The synchronization with the environment results in discontinuous fronts that propagate as periodic pulses. We characterize the transition from continuous to locked invasions and show that it is controlled by positive density-dependence in dispersal or growth. We also demonstrate that velocity locking is robust to demographic and environmental fluctuations and examine stochastic dynamics and evolution in locked invasions.

Microbiota profile in new-onset pediatric Crohn's disease: data from a non-Western population.

BACKGROUND: The role of microbiota in Crohn's disease (CD) is increasingly recognized. However, most of the reports are from Western populations. Considering the possible variation from other populations, the aim of this study was to describe the microbiota profile in children with CD in Saudi Arabia, a non-Western developing country population. RESULTS: Significantly more abundant genera in children with CD included Fusobacterium, Peptostreptococcus, Psychrobacter, and Acinetobacter; whereas the most significantly-depleted genera included Roseburia, Clostridium, Ruminococcus, Ruminoclostridium, Intestinibacter, Mitsuokella, Megasphaera, Streptococcus, Lactobacillus, Turicibacter, and Paludibacter. Alpha diversity was significantly reduced in stool (p = 0.03) but not in mucosa (p = 0.31). Beta diversity showed significant difference in community composition between control and CD samples (p = 0.03). CONCLUSION: In this developing country, we found a pattern of microbiota in children with CD similar to Western literature, suggesting a role of recent dietary lifestyle changes in this population on microbiota structure.

Chirality provides a direct fitness advantage and facilitates intermixing in cellular aggregates.

Chirality in shape and motility can evolve rapidly in microbes and cancer cells. To determine how chirality affects cell fitness, we developed a model of chiral growth in compact aggregates such as microbial colonies and solid tumors. Our model recapitulates previous experimental findings and shows that mutant cells can invade by increasing their chirality or switching their handedness. The invasion results either in a takeover or stable coexistence between the mutant and the ancestor depending on their relative chirality. For large chiralities, the coexistence is accompanied by strong intermixing between the cells, while spatial segregation occurs otherwise. We show that the competition within the aggregate is mediated by bulges in regions where the cells with different chiralities meet. The two-way coupling between aggregate shape and natural selection is described by the chiral Kardar-Parisi-Zhang equation coupled to the Burgers' equation with multiplicative noise. We solve for the key features of this theory to explain the origin of selection on chirality. Overall, our work suggests that chirality could be an important ecological trait that mediates competition, invasion, and spatial structure in cellular populations.

Fungal dysbiosis predicts the diagnosis of pediatric Crohn's disease.

AIM: To investigate the accuracy of fungal dysbiosis in mucosa and stool for predicting the diagnosis of Crohn's disease (CD). METHODS: Children were prospectively enrolled in two medical centers: one university hospital and one private gastroenterology clinic in the city of Riyadh, Kingdom of Saudi Arabia. The children with confirmed diagnosis of CD by standard guidelines were considered cases, and the others were considered non-inflammatory bowel disease controls. Mucosal and stool samples were sequenced utilizing Illumina MiSeq chemistry following the manufacturer's protocols, and abundance and diversity of fungal taxa in mucosa and stool were analyzed. Sparse logistic regression was used to predict the diagnosis of CD. The accuracy of the classifier was tested by computing the receiver operating characteristic curves with 5-fold stratified cross-validation under 100 permutations of the training data partition and the mean area under the curve (AUC) was calculated. RESULTS: All the children were Saudi nationals. There were 15 children with CD and 20 controls. The mean age was 13.9 (range: 6.7-17.8) years for CD children and 13.9 (3.25-18.6) years for controls, and 10/15 (67%) of the CD and 13/20 (65%) of the control subjects were boys. CD locations at diagnosis were ileal (L1) in 4 and colonic (L3) in 11 children, while CD behavior was non-stricturing and non-penetrating (B1) in 12 and stricturing (B2) in 3 children. The mean AUC for the fungal dysbiosis classifier was significantly higher in stools (AUC = 0.85 ± 0.057) than in mucosa (AUC = 0.71 ± 0.067) (P < 0.001). Most fungal species were significantly more depleted in stools than mucosal samples, except for Saccharomyces cerevisiae and S. bayanus, which were significantly more abundant. Diversity was significantly more reduced in stools than in mucosa. CONCLUSION: We found high AUC of fungal dysbiosis in fecal samples of children with CD, suggesting high accuracy in predicting diagnosis of CD.

Fluctuations uncover a distinct class of traveling waves.

Epidemics, flame propagation, and cardiac rhythms are classic examples of reaction-diffusion waves that describe a switch from one alternative state to another. Only two types of waves are known: pulled, driven by the leading edge, and pushed, driven by the bulk of the wave. Here, we report a distinct class of semipushed waves for which both the bulk and the leading edge contribute to the dynamics. These hybrid waves have the kinetics of pushed waves, but exhibit giant fluctuations similar to pulled waves. The transitions between pulled, semipushed, and fully pushed waves occur at universal ratios of the wave velocity to the Fisher velocity. We derive these results in the context of a species invading a new habitat by examining front diffusion, rate of diversity loss, and fluctuation-induced corrections to the expansion velocity. All three quantities decrease as a power law of the population density with the same exponent. We analytically calculate this exponent, taking into account the fluctuations in the shape of the wave front. For fully pushed waves, the exponent is -1, consistent with the central limit theorem. In semipushed waves, however, the fluctuations average out much more slowly, and the exponent approaches 0 toward the transition to pulled waves. As a result, a rapid loss of genetic diversity and large fluctuations in the position of the front occur, even for populations with cooperative growth and other forms of an Allee effect. The evolutionary outcome of spatial spreading in such populations could therefore be less predictable than previously thought.

Interactions between species introduce spurious associations in microbiome studies.

Microbiota contribute to many dimensions of host phenotype, including disease. To link specific microbes to specific phenotypes, microbiome-wide association studies compare microbial abundances between two groups of samples. Abundance differences, however, reflect not only direct associations with the phenotype, but also indirect effects due to microbial interactions. We found that microbial interactions could easily generate a large number of spurious associations that provide no mechanistic insight. Using techniques from statistical physics, we developed a method to remove indirect associations and applied it to the largest dataset on pediatric inflammatory bowel disease. Our method corrected the inflation of p-values in standard association tests and showed that only a small subset of associations is directly linked to the disease. Direct associations had a much higher accuracy in separating cases from controls and pointed to immunomodulation, butyrate production, and the brain-gut axis as important factors in the inflammatory bowel disease.