Generalized population models and the nature of genetic drift

Determinants of rapid adaptation in species with large variance in offspring production

The speed of population adaptation to changing biotic and abiotic environments is determined by the interaction between genetic drift, positive selection and linkage effects. Many marine species (fish, crustaceans), invertebrates and pathogens of humans and crops, exhibit sweepstakes reproduction characterized by the production of a very large amount of offspring (fecundity phase) from which only a small fraction may survive to the next generation (viability phase). Using stochastic simulations, we investigate whether the occurrence of sweepstakes reproduction affects the efficiency of a positively selected unlinked locus, and thus, the speed of adaptation since fecundity and/or viability have distinguishable consequences on mutation rate, probability and fixation time of advantageous alleles. We observe that the mean number of mutations at the next generation is always the function of the population size, but the variance increases with stronger sweepstakes reproduction when mutations occur in the parents. On the one hand, stronger sweepstakes reproduction magnifies the effect of genetic drift thus increasing the probability of fixation of neutral allele and decreasing that of selected alleles. On the other hand, the time to fixation of advantageous (as well as neutral) alleles is shortened by stronger sweepstakes reproduction. Importantly, fecundity and viability selection exhibit different probabilities and times to fixation of advantageous alleles under intermediate and weak sweepstakes reproduction. Finally, alleles under both strong fecundity and viability selection display a synergistic efficiency of selection. We conclude that measuring and modelling accurately fecundity and/or viability selection are crucial to predict the adaptive potential of species with sweepstakes reproduction.

Sweepstakes reproduction facilitates rapid adaptation in highly fecund populations

Adaptation enables natural populations to survive in a changing environment. Understanding the mechanics of adaptation is therefore crucial for learning about the evolution and ecology of natural populations. We focus on the impact of random sweepstakes on selection in highly fecund haploid and diploid populations partitioned into two genetic types, with one type conferring selective advantage. For the diploid populations, we incorporate various dominance mechanisms. We assume that the populations may experience recurrent bottlenecks. In random sweepstakes, the distribution of individual recruitment success is highly skewed, resulting in a huge variance in the number of offspring contributed by the individuals present in any given generation. Using computer simulations, we investigate the joint effects of random sweepstakes, recurrent bottlenecks and dominance mechanisms on selection. In our framework, bottlenecks allow random sweepstakes to have an effect on the time to fixation, and in diploid populations, the effect of random sweepstakes depends on the dominance mechanism. We describe selective sweepstakes that are approximated by recurrent sweeps of strongly beneficial allelic types arising by mutation. We demonstrate that both types of sweepstakes reproduction may facilitate rapid adaptation (as defined based on the average time to fixation of a type conferring selective advantage conditioned on fixation of the type). However, whether random sweepstakes cause rapid adaptation depends also on their interactions with bottlenecks and dominance mechanisms. Finally, we review a case study in which a model of recurrent sweeps is shown to essentially explain population genomic data from Atlantic cod.

Extinction scenarios in evolutionary processes: a multinomial Wright-Fisher approach

We study a discrete-time multi-type Wright-Fisher population process. The mean-field dynamics of the stochastic process is induced by a general replicator difference equation. We prove several results regarding the asymptotic behavior of the model, focusing on the impact of the mean-field dynamics on it. One of the results is a limit theorem that describes sufficient conditions for an almost certain path to extinction, first eliminating the type which is the least fit at the mean-field equilibrium. The effect is explained by the metastability of the stochastic system, which under the conditions of the theorem spends almost all time before the extinction event in a neighborhood of the equilibrium. In addition to the limit theorems, we propose a maximization principle for a general deterministic replicator dynamics and study its implications for the stochastic model.

Modelling genetic stability in engineered cell populations

Predicting the evolution of engineered cell populations is a highly sought-after goal in biotechnology. While models of evolutionary dynamics are far from new, their application to synthetic systems is scarce where the vast combination of genetic parts and regulatory elements creates a unique challenge. To address this gap, we here-in present a framework that allows one to connect the DNA design of varied genetic devices with mutation spread in a growing cell population. Users can specify the functional parts of their system and the degree of mutation heterogeneity to explore, after which our model generates host-aware transition dynamics between different mutation phenotypes over time. We show how our framework can be used to generate insightful hypotheses across broad applications, from how a device's components can be tweaked to optimise long-term protein yield and genetic shelf life, to generating new design paradigms for gene regulatory networks that improve their functionality.

Allele frequency divergence reveals ubiquitous influence of positive selection in Drosophila

Resolving the role of natural selection is a basic objective of evolutionary biology. It is generally difficult to detect the influence of selection because ubiquitous non-selective stochastic change in allele frequencies (genetic drift) degrades evidence of selection. As a result, selection scans typically only identify genomic regions that have undergone episodes of intense selection. Yet it seems likely such episodes are the exception; the norm is more likely to involve subtle, concurrent selective changes at a large number of loci. We develop a new theoretical approach that uncovers a previously undocumented genome-wide signature of selection in the collective divergence of allele frequencies over time. Applying our approach to temporally resolved allele frequency measurements from laboratory and wild Drosophila populations, we quantify the selective contribution to allele frequency divergence and find that selection has substantial effects on much of the genome. We further quantify the magnitude of the total selection coefficient (a measure of the combined effects of direct and linked selection) at a typical polymorphic locus, and find this to be large (of order 1%) even though most mutations are not directly under selection. We find that selective allele frequency divergence is substantially elevated at intermediate allele frequencies, which we argue is most parsimoniously explained by positive-not negative-selection. Thus, in these populations most mutations are far from evolving neutrally in the short term (tens of generations), including mutations with neutral fitness effects, and the result cannot be explained simply as an ongoing purging of deleterious mutations.

On inferring evolutionary stability in finite populations using infinite population models

Models of evolution by natural selection often make the simplifying assumption that populations are infinitely large. In this infinite population limit, rare mutations that are selected against always go extinct, whereas in finite populations they can persist and even reach fixation. Nevertheless, for mutations of arbitrarily small phenotypic effect, it is widely believed that in sufficiently large populations, if selection opposes the invasion of rare mutants, then it also opposes their fixation. Here, we identify circumstances under which infinite-population models do or do not accurately predict evolutionary outcomes in large, finite populations. We show that there is no population size above which considering only invasion generally suffices: for any finite population size, there are situations in which selection opposes the invasion of mutations of arbitrarily small effect, but favours their fixation. This is not an unlikely limiting case; it can occur when fitness is a smooth function of the evolving trait, and when the selection process is biologically sensible. Nevertheless, there are circumstances under which opposition of invasion does imply opposition of fixation: in fact, for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document}n-player snowdrift game (a common model of cooperation) we identify sufficient conditions under which selection against rare mutants of small effect precludes their fixation—in sufficiently large populations—for any selection process. We also find conditions under which—no matter how large the population—the trait that fixes depends on the selection process, which is important because any particular selection process is only an approximation of reality.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

Fixation probabilities in evolutionary dynamics under weak selection

In evolutionary dynamics, a key measure of a mutant trait's success is the probability that it takes over the population given some initial mutant-appearance distribution. This "fixation probability" is difficult to compute in general, as it depends on the mutation's effect on the organism as well as the population's spatial structure, mating patterns, and other factors. In this study, we consider weak selection, which means that the mutation's effect on the organism is small. We obtain a weak-selection perturbation expansion of a mutant's fixation probability, from an arbitrary initial configuration of mutant and resident types. Our results apply to a broad class of stochastic evolutionary models, in which the size and spatial structure are arbitrary (but fixed). The problem of whether selection favors a given trait is thereby reduced from exponential to polynomial complexity in the population size, when selection is weak. We conclude by applying these methods to obtain new results for evolutionary dynamics on graphs.

Evaluating the contributions of purifying selection and progeny-skew in dictating within-host Mycobacterium tuberculosis evolution

The within-host evolutionary dynamics of tuberculosis (TB) remain unclear, and underlying biological characteristics render standard population genetic approaches based upon the Wright-Fisher model largely inappropriate. In addition, the compact genome combined with an absence of recombination is expected to result in strong purifying selection effects. Thus, it is imperative to establish a biologically relevant evolutionary framework incorporating these factors in order to enable an accurate study of this important human pathogen. Further, such a model is critical for inferring fundamental evolutionary parameters related to patient treatment, including mutation rates and the severity of infection bottlenecks. We here implement such a model and infer the underlying evolutionary parameters governing within-patient evolutionary dynamics. Results demonstrate that the progeny skew associated with the clonal nature of TB severely reduces genetic diversity and that the neglect of this parameter in previous studies has led to significant mis-inference of mutation rates. As such, our results suggest an underlying de novo mutation rate that is considerably faster than previously inferred, and a progeny distribution differing significantly from Wright-Fisher assumptions. This inference represents a more appropriate evolutionary null model, against which the periodic effects of positive selection, associated with drug-resistance for example, may be better assessed.

Drivers of diversity in individual life courses: Sensitivity of the population entropy of a Markov chain

Individuals differ in their life courses, but how this diversity is generated, how it has evolved and how it is maintained is less understood. However, this understanding is crucial to comprehend evolutionary and ecological population dynamics. In structured populations, individual life courses represent sequences of stages that end in death. These life course trajectories or sequences can be described by a Markov chain and individuals diversify over the course of their lives by transitioning through diverse discrete stages. The rate at which stage sequences diversify with age can be quantified by the population entropy of a Markov chain. Here, we derive sensitivities of the population entropy of a Markov chain to identify which stage transitions generate - or contribute - most to diversification in stage sequences, i.e. life courses. We then use these sensitivities to reveal potential selective forces on the dynamics of life courses. To do so we correlated the sensitivity of each matrix element (stage transition) with respect to the population entropy, to its sensitivity with respect to fitness λ, the population growth rate. Positive correlation between the two sensitivities would suggest that the stage transitions that selection has acted most strongly on (high sensitivities with respect to λ) are also those that contributed most to the diversification of life courses. Using an illustrative example on a seabird population, the Thick-billed Murres on Coats Island, that is structured by reproductive stages, we show that the most influential stage transitions for diversification of life courses are not correlated with the most influential transitions for population growth. Our finding suggests that observed diversification in life courses is neutral rather than adaptive, note this does not imply that the life histories themselves are not adaptive. We are at an early stage of understanding how individual level dynamics shape ecological and evolutionary dynamics, and many discoveries await.

Ecological and Evolutionary Processes Shaping Viral Genetic Diversity

The contemporary genomic diversity of viruses is a result of the continuous and dynamic interaction of past ecological and evolutionary processes. Thus, genome sequences of viruses can be a valuable source of information about these processes. In this review, we first describe the relevant processes shaping viral genomic variation, with a focus on the role of host–virus coevolution and its potential to give rise to eco-evolutionary feedback loops. We further give a brief overview of available methodology designed to extract information about these processes from genomic data. Short generation times and small genomes make viruses ideal model systems to study the joint effect of complex coevolutionary and eco-evolutionary interactions on genetic evolution. This complexity, together with the diverse array of lifetime and reproductive strategies in viruses ask for extensions of existing inference methods, for example by integrating multiple information sources. Such integration can broaden the applicability of genetic inference methods and thus further improve our understanding of the role viruses play in biological communities.

Inferring Demography and Selection in Organisms Characterized by Skewed Offspring Distributions

The recent increase in time-series population genomic data from experimental, natural, and ancient populations has been accompanied by a promising growth in methodologies for inferring demographic and selective parameters from such data. However, these methods have largely presumed that the populations of interest are well-described by the Kingman coalescent. In reality, many groups of organisms, including viruses, marine organisms, and some plants, protists, and fungi, typified by high variance in progeny number, may be best characterized by multiple-merger coalescent models. Estimation of population genetic parameters under Wright-Fisher assumptions for these organisms may thus be prone to serious mis-inference. We propose a novel method for the joint inference of demography and selection under the Ψ-coalescent model, termed Multiple-Merger Coalescent Approximate Bayesian Computation, or MMC-ABC. We first demonstrate mis-inference under the Kingman, and then exhibit the superior performance of MMC-ABC under conditions of skewed offspring distributions. In order to highlight the utility of this approach, we reanalyzed previously published drug-selection lines of influenza A virus. We jointly inferred the extent of progeny-skew inherent to viral replication and identified putative drug-resistance mutations.

Stationary frequencies and mixing times for neutral drift processes with spatial structure

We study a general setting of neutral evolution in which the population is of finite, constant size and can have spatial structure. Mutation leads to different genetic types (traits), which can be discrete or continuous. Under minimal assumptions, we show that the marginal trait distributions of the evolutionary process, which specify the probability that any given individual has a certain trait, all converge to the stationary distribution of the mutation process. In particular, the stationary frequencies of traits in the population are independent of its size, spatial structure and evolutionary update rule, and these frequencies can be calculated by evaluating a simple stochastic process describing a population of size one (i.e. the mutation process itself). We conclude by analysing mixing times, which characterize rates of convergence of the mutation process along the lineages, in terms of demographic variables of the evolutionary process.

Effective models and the search for quantitative principles in microbial evolution

Microbes evolve rapidly. Yet they do so in idiosyncratic ways, which depend on the specific mutations that are beneficial or deleterious in a given situation. At the same time, some population-level patterns of adaptation are strikingly similar across different microbial systems, suggesting that there may also be simple, quantitative principles that unite these diverse scenarios. We review the search for simple principles in microbial evolution, ranging from the biophysical level to emergent evolutionary dynamics. A key theme has been the use of effective models, which coarse-grain over molecular and cellular details to obtain a simpler description in terms of a few effective parameters. Collectively, these theoretical approaches provide a set of quantitative principles that facilitate understanding, prediction, and potentially control of evolutionary phenomena, though formidable challenges remain due to the ecological complexity of natural populations.

Evaluating genetic drift in time-series evolutionary analysis

The Wright-Fisher model is the most popular population model for describing the behaviour of evolutionary systems with a finite population size. Approximations have commonly been used but the model itself has rarely been tested against time-resolved genomic data. Here, we evaluate the extent to which it can be inferred as the correct model under a likelihood framework. Given genome-wide data from an evolutionary experiment, we validate the Wright-Fisher drift model as the better option for describing evolutionary trajectories in a finite population. This was found by evaluating its performance against a Gaussian model of allele frequency propagation. However, we note a range of circumstances under which standard Wright-Fisher drift cannot be correctly identified.

Coalescent inferences in conservation genetics: should the exception become the rule?

Genetic estimates of effective population size (Ne) are an established means to develop informed conservation policies. Another key goal to pursue the conservation of endangered species is keeping the connectivity across fragmented environments, to which genetic inferences of gene flow and dispersal greatly contribute. Most current statistical tools for estimating such population demographic parameters are based on Kingman's coalescent (KC). However, KC is inappropriate for taxa displaying skewed reproductive variance, a property widely observed in natural species. Coalescent models that consider skewed reproductive success-called multiple merger coalescents (MMCs)-have been shown to substantially improve estimates of Ne when the distribution of offspring per capita is highly skewed. MMCs predictions of standard population genetic parameters, including the rate of loss of genetic variation and the fixation probability of strongly selected alleles, substantially depart from KC predictions. These extended models also allow studying gene genealogies in a spatial continuum, providing a novel theoretical framework to investigate spatial connectivity. Therefore, development of statistical tools based on MMCs should substantially improve estimates of population demographic parameters with major conservation implications.

Evolution of highly fecund haploid populations

We consider a model of viability selection in a highly fecund haploid population with sweepstakes reproduction. We use simulations to estimate the time until the allelic type with highest fitness has reached high frequency in a finite population. We compare the time between two reproduction modes of high and low fecundity. We also consider the probability that the allelic type with highest fitness is lost from the population before reaching high frequency. Our simulation results indicate that highly fecund populations can evolve faster (in some cases much faster) than populations of low fecundity. However, high fecundity and sweepstakes reproduction also confer much higher risk of losing the allelic type with highest fitness from the population by chance. The impact of selection on driving alleles to high frequency varies depending on the trait value conferring highest fitness; in some cases the effect of selection can hardly be detected.

Variability in fitness effects can preclude selection of the fittest

Evolutionary biologists often predict the outcome of natural selection on an allele by measuring its effects on lifetime survival and reproduction of individual carriers. However, alleles affecting traits like sex, evolvability, and cooperation can cause fitness effects that depend heavily on differences in the environmental, social, and genetic context of individuals carrying the allele. This variability makes it difficult to summarize the evolutionary fate of an allele based solely on its effects on any one individual. Attempts to average over this variability can sometimes salvage the concept of fitness. In other cases evolutionary outcomes can only be predicted by considering the entire genealogy of an allele, thus limiting the utility of individual fitness altogether. We describe a number of intriguing new evolutionary phenomena that have emerged in studies that explicitly model long-term lineage dynamics and discuss implications for the evolution of infectious diseases.

Genotype by random environmental interactions gives an advantage to non-favored minor alleles

Fixation probability, the probability that the frequency of a newly arising mutation in a population will eventually reach unity, is a fundamental quantity in evolutionary genetics. Here we use a number of models (several versions of the Moran model and the haploid Wright-Fisher model) to examine fixation probabilities for a constant size population where the fitness is a random function of both allelic state and spatial position, despite neither allele being favored on average. The concept of fitness varying with respect to both genotype and environment is important in models of cancer initiation and progression, bacterial dynamics, and drug resistance. Under our model spatial heterogeneity redefines the notion of neutrality for a newly arising mutation, as such mutations fix at a higher rate than that predicted under neutrality. The increased fixation probability appears to be due to rare alleles having an advantage. The magnitude of this effect can be large, and is an increasing function of the spatial variance and skew in fitness. The effect is largest when the fitness values of the mutants and wild types are anti-correlated across environments. We discuss results for both a spatial ring geometry of cells (such as that of a colonic crypt), a 2D lattice and a mass-action (complete graph) arrangement.

On selection in finite populations

Two major forces shaping evolution are drift and selection. The standard models of neutral drift-the Wright-Fisher (WF) and Moran processes-can be extended to include selection. However, these standard models are not always applicable in practice, and-even without selection-many other drift models make very different predictions. For example, "generalised Wright-Fisher" models (so-called because their first two conditional moments agree with those of the WF process) can yield wildly different absorption times from WF. Additionally, evolutionary stability in finite populations depends only on fixation probabilities, which can be evaluated under less restrictive assumptions than those required to estimate fixation times or more complex population-genetic quantities. We therefore distill the notion of a selection process into a broad class of finite-population, mutationless models of drift and selection (including the WF and Moran processes). We characterize when selection favours fixation of one strategy over another, for any selection process, which allows us to derive finite-population conditions for evolutionary stability independent of the selection process. In applications, the precise details of the selection process are seldom known, yet by exploiting these new theoretical results it is now possible to make rigorously justifiable inferences about fixation of traits.

On the importance of skewed offspring distributions and background selection in virus population genetics

Many features of virus populations make them excellent candidates for population genetic study, including a very high rate of mutation, high levels of nucleotide diversity, exceptionally large census population sizes, and frequent positive selection. However, these attributes also mean that special care must be taken in population genetic inference. For example, highly skewed offspring distributions, frequent and severe population bottleneck events associated with infection and compartmentalization, and strong purifying selection all affect the distribution of genetic variation but are often not taken into account. Here, we draw particular attention to multiple-merger coalescent events and background selection, discuss potential misinference associated with these processes, and highlight potential avenues for better incorporating them into future population genetic analyses.

The Λ-Fleming-Viot Process and a Connection with Wright-Fisher Diffusion

The d-dimensional Λ-Fleming-Viot generator acting on functions g(x), with x being a vector of d allele frequencies, can be written as a Wright-Fisher generator acting on functions g with a modified random linear argument of x induced by partitioning occurring in the Λ-Fleming-Viot process. The eigenvalues and right polynomial eigenvectors are easy to see from this representation. The two-dimensional process, which has a one-dimensional generator, is considered in detail. A nonlinear equation is found for the Green's function. In a model with genic selection a proof is given that there is a critical selection value such that if the selection coefficient is greater than or equal to the critical value then fixation, when the boundary 1 is hit, has probability 1 beginning from any nonzero frequency. This is an analytic proof different from the proofs of Der, Epstein and Plotkin (2011) and Foucart (2013). An application in the infinitely-many-alleles Λ-Fleming-Viot process is finding an interesting identity for the frequency spectrum of alleles that is based on size biasing. The moment dual process in the Fleming-Viot process is the usual Λ-coalescent tree back in time. The Wright-Fisher representation using a different set of polynomials gn(x) as test functions produces a dual death process which has a similarity to the Kingman coalescent and decreases by units of one. The eigenvalues of the process are analogous to the Jacobi polynomials when expressed in terms of gn(x), playing the role of xn. Under the stationary distribution when there is mutation, is analogous to the nth moment in a beta distribution. There is a d-dimensional version gn(X), and even an intriguing Ewens' sampling formula analogy when d → ∞.

The Λ-Fleming-Viot Process and a Connection with Wright-Fisher Diffusion

Thed-dimensional Λ-Fleming-Viot generator acting on functionsg(x), withxbeing a vector ofdallele frequencies, can be written as a Wright-Fisher generator acting on functionsgwith a modified random linear argument ofxinduced by partitioning occurring in the Λ-Fleming-Viot process. The eigenvalues and right polynomial eigenvectors are easy to see from this representation. The two-dimensional process, which has a one-dimensional generator, is considered in detail. A nonlinear equation is found for the Green's function. In a model with genic selection a proof is given that there is a critical selection value such that if the selection coefficient is greater than or equal to the critical value then fixation, when the boundary 1 is hit, has probability 1 beginning from any nonzero frequency. This is an analytic proof different from the proofs of Der, Epstein and Plotkin (2011) and Foucart (2013). An application in the infinitely-many-alleles Λ-Fleming-Viot process is finding an interesting identity for the frequency spectrum of alleles that is based on size biasing. The moment dual process in the Fleming-Viot process is the usual Λ-coalescent tree back in time. The Wright-Fisher representation using a different set of polynomialsgn(x) as test functions produces a dual death process which has a similarity to the Kingman coalescent and decreases by units of one. The eigenvalues of the process are analogous to the Jacobi polynomials when expressed in terms ofgn(x), playing the role ofxn. Under the stationary distribution when there is mutation,is analogous to thenth moment in a beta distribution. There is ad-dimensional versiongn(X), and even an intriguing Ewens' sampling formula analogy whend→ ∞.

Population genomic analyses based on 1 million SNPs in commercial egg layers

Identifying signatures of selection can provide valuable insight about the genes or genomic regions that are or have been under selective pressure, which can lead to a better understanding of genotype-phenotype relationships. A common strategy for selection signature detection is to compare samples from several populations and search for genomic regions with outstanding genetic differentiation. Wright's fixation index, F_ST, is a useful index for evaluation of genetic differentiation between populations. The aim of this study was to detect selective signatures between different chicken groups based on SNP-wise F_ST calculation. A total of 96 individuals of three commercial layer breeds and 14 non-commercial fancy breeds were genotyped with three different 600K SNP-chips. After filtering a total of 1 million SNPs were available for F_ST calculation. Averages of F_ST values were calculated for overlapping windows. Comparisons of these were then conducted between commercial egg layers and non-commercial fancy breeds, as well as between white egg layers and brown egg layers. Comparing non-commercial and commercial breeds resulted in the detection of 630 selective signatures, while 656 selective signatures were detected in the comparison between the commercial egg-layer breeds. Annotation of selection signature regions revealed various genes corresponding to productions traits, for which layer breeds were selected. Among them were NCOA1, SREBF2 and RALGAPA1 associated with reproductive traits, broodiness and egg production. Furthermore, several of the detected genes were associated with growth and carcass traits, including POMC, PRKAB2, SPP1, IGF2, CAPN1, TGFb2 and IGFBP2. Our approach demonstrates that including different populations with a specific breeding history can provide a unique opportunity for a better understanding of farm animal selection.

MODELING THE ROLE OF MUTATIONS AND DENSITY INDEPENDENT DISPERSAL IN EVOLUTIONARY RESCUE

Faced with a strong and sudden deterioration of environment, a population encounters two possible options — adapt or perish. In general, it is not known which of those outcomes the environmental changes will lead to. Building on experimental research, we introduce a discrete-space, discrete-time model for environmental rescue based on the influence of population dispersal, as well as, potentially beneficial mutations. Numerical results obtained by the model are shown to correspond well to experimentally obtained data.

The equilibrium allele frequency distribution for a population with reproductive skew

We study the population genetics of two neutral alleles under reversible mutation in a model that features a skewed offspring distribution, called the Λ-Fleming-Viot process. We describe the shape of the equilibrium allele frequency distribution as a function of the model parameters. We show that the mutation rates can be uniquely identified from this equilibrium distribution, but the form of the offspring distribution cannot itself always be so identified. We introduce an estimator for the mutation rate that is consistent, independent of the form of reproductive skew. We also introduce a two-allele infinite-sites version of the Λ-Fleming-Viot process, and we use it to study how reproductive skew influences standing genetic diversity in a population. We derive asymptotic formulas for the expected number of segregating sites as a function of sample size and offspring distribution. We find that the Wright-Fisher model minimizes the equilibrium genetic diversity, for a given mutation rate and variance effective population size, compared to all other Λ-processes.

Linkage Disequilibrium Estimation of Effective Population Size with Immigrants from Divergent Populations: A Case Study on Spanish Mackerel (Scomberomorus commerson)

Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. However, pitfalls are predicted when the effective size is large because estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes of approximately 5000. This number was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow-barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest SD of Ne estimates occurred when low-frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or nonbreeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were greater than 24,000.

Dynamics of neutral and selected alleles when the offspring distribution is skewed

We analyze the dynamics of two alternative alleles in a simple model of a population that allows for large family sizes in the distribution of offspring number. This population model was first introduced by Eldon and Wakeley, who described the backward-time genealogical relationships among sampled individuals, assuming neutrality. We study the corresponding forward-time dynamics of allele frequencies, with or without selection. We derive a continuum approximation, analogous to Kimura's diffusion approximation, and we describe three distinct regimes of behavior that correspond to distinct regimes in the coalescent processes of Eldon and Wakeley. We demonstrate that the effect of selection is strongly amplified in the Eldon-Wakeley model, compared to the Wright-Fisher model with the same variance effective population size. Remarkably, an advantageous allele can even be guaranteed to fix in the Eldon-Wakeley model, despite the presence of genetic drift. We compute the selection coefficient required for such behavior in populations of Pacific oysters, based on estimates of their family sizes. Our analysis underscores that populations with the same effective population size may nevertheless experience radically different forms of genetic drift, depending on the reproductive mechanism, with significant consequences for the resulting allele dynamics.