Population Dynamics of Metastable Growth-Rate Phenotypes

Neo-Darwinian evolution has presented a paradigm for population dynamics built on random mutations and selection with a clear separation of time-scales between single-cell mutation rates and the rate of reproduction. Laboratory experiments on evolving populations until now have concentrated on the fixation of beneficial mutations. Following the Darwinian paradigm, these experiments probed populations at low temporal resolution dictated by the rate of rare mutations, ignoring the intermediate evolving phenotypes. Selection however, works on phenotypes rather than genotypes. Research in recent years has uncovered the complexity of genotype-to-phenotype transformation and a wealth of intracellular processes including epigenetic inheritance, which operate on a wide range of time-scales. Here, by studying the adaptation dynamics of genetically rewired yeast cells, we show a novel type of population dynamics in which the intracellular processes intervene in shaping the population structure. Under constant environmental conditions, we measure a wide distribution of growth rates that coexist in the population for very long durations (>100 generations). Remarkably, the fastest growing cells do not take over the population on the time-scale dictated by the width of the growth-rate distributions and simple selection. Additionally, we measure significant fluctuations in the population distribution of various phenotypes: the fraction of exponentially-growing cells, the distributions of single-cell growth-rates and protein content. The observed fluctuations relax on time-scales of many generations and thus do not reflect noisy processes. Rather, our data show that the phenotypic state of the cells, including the growth-rate, for large populations in a constant environment is metastable and varies on time-scales that reflect the importance of long-term intracellular processes in shaping the population structure. This lack of time-scale separation between the intracellular and population processes calls for a new framework for population dynamics which is likely to be significant in a wide range of biological contexts, from evolution to cancer.     

Ecological opportunity and the origin of adaptive radiations

Ecological opportunity – through entry into a new environment, the origin of a key innovation or extinction of antagonists – is widely thought to link ecological population dynamics to evolutionary diversification. The population‐level processes arising from ecological opportunity are well documented under the concept of ecological release. However, there is little consensus as to how these processes promote phenotypic diversification, rapid speciation and adaptive radiation. We propose that ecological opportunity could promote adaptive radiation by generating specific changes to the selective regimes acting on natural populations, both by relaxing effective stabilizing selection and by creating conditions that ultimately generate diversifying selection. We assess theoretical and empirical evidence for these effects of ecological opportunity and review emerging phylogenetic approaches that attempt to detect the signature of ecological opportunity across geological time. Finally, we evaluate the evidence for the evolutionary effects of ecological opportunity in the diversification of Caribbean Anolis lizards. Some of the processes that could link ecological opportunity to adaptive radiation are well documented, but others remain unsupported. We suggest that more study is required to characterize the form of natural selection acting on natural populations and to better describe the relationship between ecological opportunity and speciation rates.     

Experimental evolution with a multicellular host causes diversification within and between microbial parasite populations—Differences in emerging phenotypes of two different parasite strains

Host‐parasite coevolution is predicted to have complex evolutionary consequences, potentially leading to the emergence of genetic and phenotypic diversity for both antagonists. However, little is known about variation in phenotypic responses to coevolution between different parasite strains exposed to the same experimental conditions. We infected Caenorhabditis elegans with one of two strains of Bacillus thuringiensis and either allowed the host and the parasite to experimentally coevolve (coevolution treatment) or allowed only the parasite to adapt to the host (one‐sided parasite adaptation). By isolating single parasite clones from evolved populations, we found phenotypic diversification of the ancestral strain into distinct clones, which varied in virulence toward ancestral hosts and competitive ability against other parasite genotypes. Parasite phenotypes differed remarkably not only between the two strains, but also between and within different replicate populations, indicating diversification of the clonal population caused by selection. This study highlights that the evolutionary selection pressure mediated by a multicellular host causes phenotypic diversification, but not necessarily with the same phenotypic outcome for different parasite strains.     

Linking divergent selection on vegetative traits to environmental variation and phenotypic diversification in the Iberian columbines (Aquilegia)

Divergent selection is a key in the ecological theory of adaptive radiation. Most evidence on its causes and consequences relies on studies of pairs of populations or closely related taxa. However, adaptive radiation involves multiple taxa adapted to different environmental factors. We propose an operational definition of divergent selection to explore the continuum between divergent and convergent selection in multiple populations and taxa, and its links with environmental variation and phenotypic and taxonomic differentiation. We apply this approach to explore phenotypic differentiation of vegetative traits between 15 populations of four taxa of Iberian columbines (Gen. Aquilegia). Differences in soil rockiness impose divergent selection on inflorescence height and the number of flowers per inflorescence, likely affecting the processes of phenotypic and, in the case of inflorescence height, taxonomic diversification between taxa. Elevational variation imposes divergent selection on the number of leaves; however, the current pattern of divergent selection on this trait seems related to ecotypic differentiation within taxa but not to their taxonomic diversification.     

DETERMINISM IN THE DIVERSIFICATION OF HISPANIOLAN TRUNK‐GROUND ANOLES (ANOLIS CYBOTES SPECIES COMPLEX)

The evolutionary processes that produce adaptive radiations are enigmatic. They can only be studied after the fact, once a radiation has occurred and been recognized, rather than while the processes are ongoing. One way to connect pattern to process is to study the processes driving divergence today among populations of species that belong to an adaptive radiation, and compare the results to patterns observed at a deeper, macroevolutionary level. We tested whether evolution is a deterministic process with similar outcomes during different stages of the adaptive radiation of Anolis lizards. Using a clade of terrestrial–scansorial lizards in the genus Anolis, we inferred the adaptive basis of spatial variation among contemporary populations and tested whether axes of phenotypic differentiation among them mirror known axes of diversification at deeper levels of the anole radiation. Nonparallel change associated with genetic divergence explains the vast majority of geographic variation. However, we found phenotypic variation to be adaptive as confirmed by convergence in populations occurring in similar habitats in different mountain ranges. Morphological diversification among populations recurs deterministically along two axes of diversification previously identified in the anole radiation, but the characters involved differ from those involved in adaptation at higher levels of anole phylogeny.     

Genetic and phenotypic variation in a colourful treefrog across five geographic barriers

Aim Understanding the patterns and processes underlying phenotype in a polytypic species provides key insights into microevolutionary mechanisms of diversification. The red‐eyed treefrog, Agalychnis callidryas, exhibits strong regional differentiation in colour pattern, corresponding to five admixed mitochondrial DNA clades. We evaluated spatial diversity patterns across multiple, putative barriers to examine the fine‐scale processes that mediate phenotypic divergence between some regions while maintaining homogeneity between others. Location We examined patterns of phenotypic diversification among 17 sites that span five putative biogeographic barriers in lower Central America (Costa Rica and Panama). Methods We tested the extent to which genetic distance (F_ST) derived from six multilocus nuclear genotypes covaried with measures of phenotypic distance (leg coloration) within and between biogeographic regions. We used linear regression analyses to determine the role of geographic and genetic factors in structuring spatial patterns of phenotypic diversity. Results The factors that best explained patterns of phenotypic diversity varied among biogeographic regions. We identified one geographic barrier that impeded gene exchange and resulted in concordant phenotypic divergence across the Continental Divide, isolating Caribbean and Pacific populations. Across Caribbean Costa Rican populations, one barrier structured phenotypic but not genetic diversity patterns, indicating a role for selection. In other regions, the putative barriers had no determining effect on either genetic or leg colour structure. Main conclusions The processes mediating the distribution and diversification of colour pattern in this polytypic, wide‐ranging treefrog varied among biogeographic regions. Spatially varying selection combined with the isolating effects of geographic factors probably resulted in the patchy distribution of colour diversity across Costa Rican and Panamanian populations.     

Continuous cultivation of fission yeast: Analysis of single-cell protein synthesis kinetics

A fundamental problem in microbial reactor analysis is identification of the relationship between environment and individual cell metabolic activity. Population balance equations provide a link between experimental measurements of composition frequency functions in microbial populations on the one hand and macromolecular synthesis kinetics and cell division control parameters for single cells on the other. Flow microfluorometry measurements of frequency functions for single-cell protein content in Schizosaccharomyces pombe in balanced exponential growth have been analyzed by two different methods. One approach utilizes the integrated form of the population balance equation known as the Collins-Richmond equation, and the other method involves optimization of parameters in assumed kinetic and cell division functional forms in order to best fit measured frequency functions with corresponding model solutions. Both data interpretation techniques indicate that rates of protein synthesis increase most in small protein content cells as the population specific growth rate increases, leading to parabolic single-cell protein synthesis kinetics at large specific growth rates. Utilization of frequency function data for an asynchronous population is shown in this case to be a far more sensitive method for determination of single-cell kinetics than is monitoring the metabolic dynamics of a single cell or, equivalently, synchronous culture analyses.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

Applications of single-cell technology on bacterial analysis

Background: Traditionally, scientists studied microbiology through the manner of batch cultures, to conclude the dynamics or outputs by averaging all individuals. However, as the researches go further, the heterogeneities among the individuals have been proven to be crucial for the population dynamics and fates. Results: Due to the limit of technology, single-cell analysis methods were not widely used to decipher the inherent connections between individual cells and populations. Since the early decades of this century, the rapid development of microfluidics, fluorescent labelling, next-generation sequencing, and high-resolution microscopy have speeded up the development of single-cell technologies and further facilitated the applications of these technologies on bacterial analysis. Conclusions: In this review, we summarized the recent processes of single-cell technologies applied in bacterial analysis in terms of intracellular characteristics, cell physiology dynamics, and group behaviors, and discussed how single-cell technologies could be more applicable for future bacterial researches.     

Discordant patterns of evolutionary differentiation in two Neotropical treefrogs

Comparative studies of codistributed taxa test the degree to which historical processes have shaped contemporary population structure. Discordant patterns of lineage divergence among taxa indicate that species differ in their response to common historical processes. The complex geologic landscape of the Isthmus of Central America provides an ideal setting to test the effects of vicariance and other biogeographic factors on population history. We compared divergence patterns between two codistributed Neotropical frogs ( Dendropsophus ebraccatus and Agalychnis callidryas ) that exhibit colour pattern polymorphisms among populations, and found significant differences between them in phenotypic and genetic divergence among populations. Colour pattern in D. ebraccatus did not vary with genetic or geographic distance, while colour pattern co-varied with patterns of gene flow in A. callidryas . In addition, we detected significant species differences in the phylogenetic history of populations, gene flow among them, and the extent to which historical diversification and recent gene flow have been restricted by five biogeographic barriers in Costa Rica and Panama. We inferred that alternate microevolutionary processes explain the unique patterns of diversification in each taxon. Our study underscores how differences in selective regimes and species-typical ecological and life-history traits maintain spatial patterns of diversification.     

Untangling the evolutionary history of a highly polymorphic species: introgressive hybridization and high genetic structure in the desert cichlid fish Herichtys minckleyi

Understanding the origin of biodiversity requires knowledge on the evolutionary processes that drive divergence and speciation, as well as on the processes constraining it. Intraspecific polymorphisms can provide insight into the mechanisms that generate and maintain phenotypic, behavioural and life history diversification, and can help us understand not only the processes that lead to speciation but also the processes that prevent local fixation of morphs. The ‘desert cichlid’ Herichtys minckleyi is a highly polymorphic species endemic to a biodiversity hotspot in northern Mexico, the Cuatro Ciénegas valley. This species is polymorphic in body shape and trophic apparatus, and eco‐morphotypes coexist in small spring‐fed lagoons across the valley. We investigated the genetic structure of these polymorphisms and their phylogeographic history by analysing the entire control region of the mitochondrial DNA and 10 nuclear microsatellite markers in several populations from different sites and morphs. We found two very divergent mitochondrial lineages that most likely predate the closing of the valley and are not associated with morphotypes or sites. One of these lineages is also found in the sister species Herichthys cyanoguttatus. Data from neutral microsatellite markers suggest that most lagoons or drainages constitute their own genetic cluster with sympatric eco‐morphotypes forming panmictic populations. Alternative mechanisms such as phenotypic plasticity and a few loci controlled traits provide possible explanations for the sympatric coexistence of discrete nonoverlapping eco‐morphotypes with apparent lack of barriers to gene flow within multiple lagoons and drainages.     

Phenotypic variation in invasive and biocontrol populations of the harlequin ladybird, Harmonia axyridis

Despite numerous releases for biological control purposes during more than 20 years in Europe, Harmonia axyridis failed to become established until the beginning of the 21st century. Its status as invasive alien species is now widely recognised. Theory suggests that invasive populations should evolve toward greater phenotypic plasticity because they encounter differing environments during the invasion process. On the contrary, populations used for biological control have been maintained under artificial rearing conditions for many generations; they are hence expected to become specialised on a narrow range of environments and show lower phenotypic plasticity. Here we compared phenotypic traits and the extent of adaptive phenotypic plasticity in two invasive populations and two populations commercialized for biological control by (i) measuring six phenotypic traits related to fitness (eggs hatching rate, larval survival rate, development time, sex ratio, fecundity over 6 weeks and survival time of starving adults) at three temperatures (18, 24 and 30°C), (ii) recording the survival rate and quiescence aggregation behaviour when exposed to low temperatures (5, 10 and 15°C), and (iii) studying the cannibalistic behaviour of populations in the absence of food. Invasive and biocontrol populations displayed significantly different responses to temperature variation for a composite fitness index computed from the traits measured at 18, 24 and 30°C, but not for any of those traits considered independently. The plasticity measured on the same fitness index was higher in the two invasive populations, but this difference was not statistically significant. On the other hand, invasive populations displayed significantly higher survival and higher phenotypic plasticity when entering into quiescence at low temperatures. In addition, one invasive population displayed a singular cannibalistic behaviour. Our results hence only partly support the expectation of increased adaptive phenotypic plasticity of European invasive populations of H. axyridis, and stress the importance of the choice of the environmental parameters to be manipulated for assessing phenotypic plasticity variation among populations.     

Single-cell behavior and population heterogeneity: solving an inverse problem to compute the intrinsic physiological state functions

The dynamics of isogenic cell populations can be described by cell population balance models that account for phenotypic heterogeneity. To utilize the predictive power of these models, however, we must know the rates of single-cell reaction and division and the bivariate partition probability density function. These three intrinsic physiological state (IPS) functions can be obtained by solving an inverse problem that requires knowledge of the phenotypic distributions for the overall cell population, the dividing cell subpopulation and the newborn cell subpopulation. We present here a robust computational procedure that can accurately estimate the IPS functions for heterogeneous cell populations. A detailed parametric analysis shows how the accuracy of the inverse solution is affected by discretization parameters, the type of non-parametric estimators used, the qualitative characteristics of phenotypic distributions and the unknown partitioning probability density function. The effect of finite sampling and measurement errors on the accuracy of the recovered IPS functions is also assessed. Finally, we apply the procedure to estimate the IPS functions of an E. coli population carrying an IPTG-inducible genetic toggle network. This study completes the development of an integrated experimental and computational framework that can become a powerful tool for quantifying single-cell behavior using measurements from heterogeneous cell populations.     

Stochastic state transitions give rise to phenotypic equilibrium in populations of cancer cells

Cancer cells within individual tumors often exist in distinct phenotypic states that differ in functional attributes. While cancer cell populations typically display distinctive equilibria in the proportion of cells in various states, the mechanisms by which this occurs are poorly understood. Here, we study the dynamics of phenotypic proportions in human breast cancer cell lines. We show that subpopulations of cells purified for a given phenotypic state return towards equilibrium proportions over time. These observations can be explained by a Markov model in which cells transition stochastically between states. A prediction of this model is that, given certain conditions, any subpopulation of cells will return to equilibrium phenotypic proportions over time. A second prediction is that breast cancer stem-like cells arise de novo from non-stem-like cells. These findings contribute to our understanding of cancer heterogeneity and reveal how stochasticity in single-cell behaviors promotes phenotypic equilibrium in populations of cancer cells.     

Direct Correlation between Motile Behavior and Protein Abundance in Single Cells

Understanding how stochastic molecular fluctuations affect cell behavior requires the quantification of both behavior and protein numbers in the same cells. Here, we combine automated microscopy with in situ hydrogel polymerization to measure single-cell protein expression after tracking swimming behavior. We characterized the distribution of non-genetic phenotypic diversity in Escherichia coli motility, which affects single-cell exploration. By expressing fluorescently tagged chemotaxis proteins (CheR and CheB) at different levels, we quantitatively mapped motile phenotype (tumble bias) to protein numbers using thousands of single-cell measurements. Our results disagreed with established models until we incorporated the role of CheB in receptor deamidation and the slow fluctuations in receptor methylation. Beyond refining models, our central finding is that changes in numbers of CheR and CheB affect the population mean tumble bias and its variance independently. Therefore, it is possible to adjust the degree of phenotypic diversity of a population by adjusting the global level of expression of CheR and CheB while keeping their ratio constant, which, as shown in previous studies, confers functional robustness to the system. Since genetic control of protein expression is heritable, our results suggest that non-genetic diversity in motile behavior is selectable, supporting earlier hypotheses that such diversity confers a selective advantage.     

Morphological and genetic differentiation in populations of the dispersal‐limited coco de mer (Lodoicea maldivica): implications for management and conservation

Aims Developing plant conservation strategies requires knowledge of ecological and genetic processes underlying population dynamics. We aimed to quantify morphological and genetic differentiation among remnant populations of the iconic coco‐de‐mer palm Lodoicea maldivica. We hypothesized that limited gene flow among widely spaced populations would result in high genetic variation and large phenotypic differences among populations. Location Islands of Praslin and Curieuse (CU), Seychelles, Indian Ocean. Methods We conducted an extensive population survey and recorded morphological parameters for 447 Lodoicea in the main populations at Vallée de Mai (VM) and Fond Ferdinand (FF) on Praslin, and on CU. We collected leaf material from 180 trees in these populations for DNA genotyping using amplified fragment length polymorphisms. Results A total of 16,766 Lodoicea trees were recorded in the three populations (72.6% of Lodoicea on both islands). Lodoicea trees at VM and FF showed similar morphology, but differed in most parameters from those at CU, which were shorter, grew more slowly and produced fewer seeds. Mean overall genetic diversity was 0.337, and percentage of polymorphic loci was 91.1. Genetic diversity of the CU population was lower than that at VM and FF. There was weak genetic differentiation between CU and Praslin populations, but 99% of all genetic diversity was within populations. Main conclusions Trees on CU differed in growth and morphology from those of the two Praslin populations. These phenotypic differences, however, were not mirrored in the genetic structure of the populations. All populations were relatively genetically diverse with remarkably little differentiation among populations. This suggests that the capacity of Lodoicea to dominate across a range of habitats may be because of high phenotypic plasticity. High genetic connectivity may be maintained through long‐distance wind pollination. Given the uncertainty about the extent of underlying adaptive variation, we recommend that restoration projects avoid transferring seeds between island populations.     

Automated Characterization and Parameter-Free Classification of Cell Tracks Based on Local Migration Behavior

Cell migration is the driving force behind the dynamics of many diverse biological processes. Even though microscopy experiments are routinely performed today by which populations of cells are visualized in space and time, valuable information contained in image data is often disregarded because statistical analyses are performed at the level of cell populations rather than at the single-cell level. Image-based systems biology is a modern approach that aims at quantitatively analyzing and modeling biological processes by developing novel strategies and tools for the interpretation of image data. In this study, we take first steps towards a fully automated characterization and parameter-free classification of cell track data that can be generally applied to tracked objects as obtained from image data. The requirements to achieve this aim include: (i) combination of different measures for single cell tracks, such as the confinement ratio and the asphericity of the track volume, and (ii) computation of these measures in a staggered fashion to retrieve local information from all possible combinations of track segments. We demonstrate for a population of synthetic cell tracks as well as for in vitro neutrophil tracks obtained from microscopy experiment that the information contained in the track data is fully exploited in this way and does not require any prior knowledge, which keeps the analysis unbiased and general. The identification of cells that show the same type of migration behavior within the population of all cells is achieved via agglomerative hierarchical clustering of cell tracks in the parameter space of the staggered measures. The recognition of characteristic patterns is highly desired to advance our knowledge about the dynamics of biological processes.