The population genetics of beneficial mutations

Dynamic expedition of leading mutations in SARS-CoV-2 spike glycoproteins

The continuous evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the recent pandemic, has generated countless new variants with varying fitness. Mutations of the spike glycoprotein play a particularly vital role in shaping its evolutionary trajectory, as they have the capability to alter its infectivity and antigenicity. We present a time-resolved statistical method, Dynamic Expedition of Leading Mutations (deLemus), to analyze the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein. The proposed L -index of the deLemus method is effective in quantifying the mutation strength of each amino acid site and outlining evolutionarily significant sites, allowing the comprehensive characterization of the evolutionary mutation pattern of the spike glycoprotein.

Fixation times of de novo and standing beneficial variants in subdivided populations

The rate at which beneficial alleles fix in a population depends on the probability of and time to fixation of such alleles. Both of these quantities can be significantly impacted by population subdivision and limited gene flow. Here, we investigate how limited dispersal influences the rate of fixation of beneficial de novo mutations, as well as fixation time from standing genetic variation. We investigate this for a population structured according to the island model of dispersal allowing us to use the diffusion approximation, which we complement with simulations. We find that fixation may take on average fewer generations under limited dispersal than under panmixia when selection is moderate. This is especially the case if adaptation occurs from de novo recessive mutations, and dispersal is not too limited (such that approximately FST<0.2). The reason is that mildly limited dispersal leads to only a moderate increase in effective population size (which slows down fixation), but is sufficient to cause a relative excess of homozygosity due to inbreeding, thereby exposing rare recessive alleles to selection (which accelerates fixation). We also explore the effect of metapopulation dynamics through local extinction followed by recolonization, finding that such dynamics always accelerate fixation from standing genetic variation, while de novo mutations show faster fixation interspersed with longer waiting times. Finally, we discuss the implications of our results for the detection of sweeps, suggesting that limited dispersal mitigates the expected differences between the genetic signatures of sweeps involving recessive and dominant alleles.

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.

The distribution of fitness effects during adaptive walks using a simple genetic network

The tempo and mode of adaptation depends on the availability of beneficial alleles. Genetic interactions arising from gene networks can restrict this availability. However, the extent to which networks affect adaptation remains largely unknown. Current models of evolution consider additive genotype-phenotype relationships while often ignoring the contribution of gene interactions to phenotypic variance. In this study, we model a quantitative trait as the product of a simple gene regulatory network, the negative autoregulation motif. Using forward-time genetic simulations, we measure adaptive walks towards a phenotypic optimum in both additive and network models. A key expectation from adaptive walk theory is that the distribution of fitness effects of new beneficial mutations is exponential. We found that both models instead harbored distributions with fewer large-effect beneficial alleles than expected. The network model also had a complex and bimodal distribution of fitness effects among all mutations, with a considerable density at deleterious selection coefficients. This behavior is reminiscent of the cost of complexity, where correlations among traits constrain adaptation. Our results suggest that the interactions emerging from genetic networks can generate complex and multimodal distributions of fitness effects.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Origins and evolution of biological novelty

Understanding the origins and impacts of novel traits has been a perennial interest in many realms of ecology and evolutionary biology. Here, we build on previous evolutionary and philosophical treatments of this subject to encompass novelties across biological scales and eco-evolutionary perspectives. By defining novelties as new features at one biological scale that have emergent effects at other biological scales, we incorporate many forms of novelty that have previously been treated in isolation (such as novelty from genetic mutations, new developmental pathways, new morphological features, and new species). Our perspective is based on the fundamental idea that the emergence of a novelty, at any biological scale, depends on its environmental and genetic context. Through this lens, we outline a broad array of generative mechanisms underlying novelty and highlight how genomic tools are transforming our understanding of the origins of novelty. Lastly, we present several case studies to illustrate how novelties across biological scales and systems can be understood based on common mechanisms of change and their environmental and genetic contexts. Specifically, we highlight how gene duplication contributes to the evolution of new complex structures in visual systems; how genetic exchange in symbiosis alters functions of both host and symbiont, resulting in a novel organism; and how hybridisation between species can generate new species with new niches.

Eco-evolutionary extinction and recolonization dynamics reduce genetic load and increase time to extinction in highly inbred populations

Understanding how genetic and ecological effects can interact to shape genetic loads within and across local populations is key to understanding ongoing persistence of systems that should otherwise be susceptible to extinction through mutational meltdown. Classic theory predicts short persistence times for metapopulations comprising small local populations with low connectivity, due to accumulation of deleterious mutations. Yet, some such systems have persisted over evolutionary time, implying the existence of mechanisms that allow metapopulations to avoid mutational meltdown. We first hypothesize a mechanism by which the combination of stochasticity in the numbers and types of mutations arising locally (genetic stochasticity), resulting local extinction, and recolonization through evolving dispersal facilitates metapopulation persistence. We then test this mechanism using a spatially and genetically explicit individual-based model. We show that genetic stochasticity in highly structured metapopulations can result in local extinctions, which can favor increased dispersal, thus allowing recolonization of empty habitat patches. This causes fluctuations in metapopulation size and transient gene flow, which reduces genetic load and increases metapopulation persistence over evolutionary time. Our suggested mechanism and simulation results provide an explanation for the conundrum presented by the continued persistence of highly structured populations with inbreeding mating systems that occur in diverse taxa.

Massive-scale genomic analysis reveals SARS-CoV-2 mutation characteristics and evolutionary trends

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic resulted in significant societal costs. Hence, an in-depth understanding of SARS-CoV-2 virus mutation and its evolution will help determine the direction of the COVID-19 pandemic. In this study, we identified 296,728 de novo mutations in more than 2,800,000 high-quality SARS-CoV-2 genomes. All possible factors affecting the mutation frequency of SARS-CoV-2 in human hosts were analyzed, including zinc finger antiviral proteins, sequence context, amino acid change, and translation efficiency. As a result, we proposed that when adenine (A) and tyrosine (T) bases are in the context of AM (M stands for adenine or cytosine) or TA motif, A or T base has lower mutation frequency. Furthermore, we hypothesized that translation efficiency can affect the mutation frequency of the third position of the codon by the selection, which explains why SARS-CoV-2 prefers AT3 codons usage. In addition, we found a host-specific asymmetric dinucleotide mutation frequency in the SARS-CoV-2 genome, which provides a new basis for determining the origin of the SARS-CoV-2. Finally, we summarize all possible factors affecting mutation frequency and provide insights into the mutation characteristics and evolutionary trends of SARS-CoV-2.

The genetic basis of selfing rate evolution

Evolution of selfing is common in plant populations, but the genetic basis of selfing rate evolution remains unclear. Although the effects of genetic properties on fixation for mating-unrelated alleles have been investigated, loci that modify the selfing rate (selfing modifiers) differ from mating-unrelated loci in several aspects. Using population genetic models, I investigate the genetic basis of selfing rate evolution. For mating-unrelated alleles, selfing promotes fixation only for recessive mutations, but for selfing modifiers, because the selection coefficient depends on the background selfing rate, selfing can promote fixation even for dominant modifiers. For mating-unrelated alleles, the fixation probability from standing variation is independent of dominance and decreases with an increased background selfing rate. However, for selfing modifiers, the fixation probability peaks at an intermediate selfing rate and when alleles are recessive, because a change of its selection coefficient necessarily involves a change of the inbreeding coefficient, because both depend on the level of inbreeding depression. Furthermore, evolution of selfing involving multiple modifier loci is more likely when selfing is controlled by few large-effect rather than many slight-effect modifiers. I discuss how these characteristics of selfing modifiers have implications for the unidirectional transition from outcrossing to selfing and other empirical patterns.

Population size mediates the contribution of high-rate and large-benefit mutations to parallel evolution

Mutations with large fitness benefits and mutations occurring at high rates may both cause parallel evolution, but their contribution is predicted to depend on population size. Moreover, high-rate and large-benefit mutations may have different long-term adaptive consequences. We show that small and 100-fold larger bacterial populations evolve resistance to a β-lactam antibiotic by using similar numbers, but different types of mutations. Small populations frequently substitute similar high-rate structural variants and loss-of-function point mutations, including the deletion of a low-activity β-lactamase, and evolve modest resistance levels. Large populations more often use low-rate, large-benefit point mutations affecting the same targets, including mutations activating the β-lactamase and other gain-of-function mutations, leading to much higher resistance levels. Our results demonstrate the separation by clonal interference of mutation classes with divergent adaptive consequences, causing a shift from high-rate to large-benefit mutations with increases in population size.

An unbiased test reveals no enrichment of sexually antagonistic polymorphisms on the human X chromosome

Mutations with beneficial effects in one sex can have deleterious effects in the other. Such 'sexually antagonistic' (SA) variants contribute to variation in life-history traits and overall fitness, yet their genomic distribution is poorly resolved. Theory predicts that SA variants could be enriched on the X chromosome or autosomes, yet current empirical tests face two formidable challenges: (i) identifying SA selection in genomic data is difficult; and (ii) metrics of SA variation show persistent biases towards the X, even when SA variants are randomly distributed across the genome. Here, we present an unbiased test of the theory that SA variants are enriched on the X. We first develop models for reproductive F_ST-a metric for quantifying sex-differential (including SA) effects of genetic variants on lifetime reproductive success-that control for X-linked biases. Comparing data from approximately 250 000 UK Biobank individuals to our models, we find F_ST elevations consistent with both X-linked and autosomal SA polymorphisms affecting reproductive success in humans. However, the extent of F_ST elevations does not differ from a model in which SA polymorphisms are randomly distributed across the genome. We argue that the polygenic nature of SA variation, along with sex asymmetries in SA effects, might render X-linked enrichment of SA polymorphisms unlikely.

Weak coupling among barrier loci and waves of neutral and adaptive introgression across an expanding hybrid zone

Hybridization can serve as an evolutionary stimulus, but we have little understanding of introgression at early stages of hybrid zone formation. We analyze reproductive isolation and introgression between a range‐limited and a widespread species. Reproductive barriers are estimated based on differences in flowering time, ecogeographic distributions, and seed set from crosses. We find an asymmetrical mating barrier due to cytonuclear incompatibility that is consistent with observed clusters of coincident and concordant tension zone clines (barrier loci) for mtDNA haplotypes and nuclear SNPs. These groups of concordant clines are spread across the hybrid zone, resulting in weak coupling among barrier loci and extensive introgression. Neutral clines had nearly equal introgression into both species’ ranges, whereas putative cases of adaptive introgression had exceptionally wide clines with centers shifted toward one species. Analyses of cline shape indicate that secondary contact was initiated within the last 800 generations with the per‐generation dispersal between 200 and 400 m, and provide some of the first estimates of the strength of selection required to account for observed levels of adaptive introgression. The weak species boundary between these species appears to be in early stages of dissolution, and ultimately will precipitate genetic swamping of the range‐limited species.

Allen Orr and the genetics of adaptation

Over most of the 20th century, evolutionary biologists predominantly subscribed to a strong form of "micro-mutationism," in which adaptive phenotypic divergence arises from allele frequency changes at many loci, each with a small effect on the phenotype. To be sure, there were well-known examples of large-effect alleles contributing to adaptation, yet such cases were generally regarded as atypical and unrepresentative of evolutionary change in general. In 1998, Allen Orr published a landmark theoretical paper in Evolution, which showed that both small- and large-effect mutations are likely to contribute to "adaptive walks" of a population to an optimum. Coupled with a growing set of empirical examples of large-effect alleles contributing to divergence (e.g., from QTL studies), Orr's paper provided a mathematical formalism that converted many evolutionary biologists from micro-mutationism to a more pluralistic perspective on the genetic basis of evolutionary change. We revisit the theoretical insights emerging from Orr's paper within the historical context leading up to 1998, and track the influence of this paper on the field of evolutionary biology through an examination of its citations over the last two decades and an analysis of the extensive body of theoretical and empirical research that Orr's pioneering paper inspired.

Extreme value theory as a framework for understanding mutation frequency distribution in cancer genomes

Currently, the population dynamics of preclonal cancer cells before clonal expansion of tumors has not been sufficiently addressed thus far. By focusing on preclonal cancer cell population as a Darwinian evolutionary system, we formulated and analyzed the observed mutation frequency among tumors (MFaT) as a proxy for the hypothesized sequence read frequency and beneficial fitness effect of a cancer driver mutation. Analogous to intestinal crypts, we assumed that sample donor patients are separate culture tanks where proliferating cells follow certain population dynamics described by extreme value theory (EVT). To validate this, we analyzed three large-scale cancer genome datasets, each harboring > 10000 tumor samples and in total involving > 177898 observed mutation sites. We clarified the necessary premises for the application of EVT in the strong selection and weak mutation (SSWM) regime in relation to cancer genome sequences at scale. We also confirmed that the stochastic distribution of MFaT is likely of the Fréchet type, which challenges the well-known Gumbel hypothesis of beneficial fitness effects. Based on statistical data analysis, we demonstrated the potential of EVT as a population genetics framework to understand and explain the stochastic behavior of driver-mutation frequency in cancer genomes as well as its applicability in real cancer genome sequence data.

The maintenance of polygenic sex determination depends on the dominance of fitness effects which are predictive of the role of sexual antagonism

In species with polygenic sex determination (PSD), multiple male- and female-determining loci on different proto-sex chromosomes segregate as polymorphisms within populations. The extent to which these polymorphisms are at stable equilibria is not yet resolved. Previous work demonstrated that PSD is most likely to be maintained as a stable polymorphism when the proto-sex chromosomes have opposite (sexually antagonistic) fitness effects in males and females. However, these models usually consider PSD systems with only two proto-sex chromosomes, or they do not broadly consider the dominance of the alleles under selection. To address these shortcomings, I used forward population genetic simulations to identify selection pressures that can maintain PSD under different dominance scenarios in a system with more than two proto-sex chromosomes (modeled after the house fly). I found that overdominant fitness effects of male-determining proto-Y chromosomes are more likely to maintain PSD than dominant, recessive, or additive fitness effects. The overdominant fitness effects that maintain PSD tend to have proto-Y chromosomes with sexually antagonistic effects (male-beneficial and female-detrimental). In contrast, dominant fitness effects that maintain PSD tend to have sexually antagonistic multi-chromosomal genotypes, but the individual proto-sex chromosomes do not have sexually antagonistic effects. These results demonstrate that sexual antagonism can be an emergent property of the multi-chromosome genotype without individual sexually antagonistic chromosomes. My results further illustrate how the dominance of fitness effects has consequences for both the likelihood that PSD will be maintained as well as the role sexually antagonistic selection is expected to play in maintaining the polymorphism.

On the waiting time until coordinated mutations get fixed in regulatory sequences

In this paper we consider the time evolution of a population of size N with overlapping generations, in the vicinity of m genes. We assume that this population is subject to point mutations, genetic drift, and selection. More specifically, we analyze the statistical distribution of the waiting time T_m until the expression of these genes have changed for all individuals, when transcription factors recognize and attach to short DNA-sequences (binding sites) within regulatory sequences in the neighborhoods of the m genes. The evolutionary dynamics is described by a multitype Moran process, where each individual is assigned a m×L regulatory array that consists of regulatory sequences with L nucleotides for all m genes. We study how the waiting time distribution depends on the number of genes, the mutation rate, the length of the binding sites, the length of the regulatory sequences, and the way in which the targeted binding sites are coordinated for different genes in terms of selection coefficients. These selection coefficients depend on how many binding sites have appeared so far, and possibly on their order of appearance. We also allow for back mutations, whereby some acquired binding sites may be lost over time. It is further assumed that the mutation rate is small enough to warrant a fixed state population, so that all individuals have the same regulatory array, at any given time point, until the next successful mutation arrives in some individual and spreads to the rest of the population. We further incorporate stochastic tunneling, whereby successful mutations get mutated before their fixation. A crucial part of our approach is to divide the huge state space of regulatory arrays into a small number of components, assuming that the array component varies as a Markov process over time. This implies that T_m is the time until this Markov process hits an absorbing state, with a phase-type distribution. A number of interesting results can be derived from our general setup, for instance that the expected waiting time increases exponentially with m, for a selectively neutral model, when back-mutations are possible.

Genomic Background Governs Opposing Responses to Nalidixic Acid upon Megaplasmid Acquisition in Pseudomonas

Horizontal gene transfer is a significant driver of evolutionary dynamics across microbial populations. Although the benefits of the acquisition of new genetic material are often quite clear, experiments across systems have demonstrated that gene transfer events can cause significant phenotypic changes and entail fitness costs in a way that is dependent on the genomic and environmental context. Here, we test for the generality of one previously identified cost, sensitization of cells to the antibiotic nalidixic acid after acquisition of an ∼1-Mb megaplasmid, across Pseudomonas strains and species. Overall, we find that the presence of this megaplasmid sensitizes many different Pseudomonas strains to nalidixic acid but that this same horizontal gene transfer event increases resistance of Pseudomonas putida KT2440 to nalidixic acid across assays as well as to ciprofloxacin under competitive conditions. These phenotypic results are not easily explained away as secondary consequences of overall fitness effects and appear to occur independently of another cost associated with this megaplasmid, sensitization to higher temperatures. Lastly, we draw parallels between these reported results and the phenomenon of sign epistasis for de novo mutations and explore how context dependence of effects of plasmid acquisition could impact overall evolutionary dynamics and the evolution of antimicrobial resistance.IMPORTANCE Numerous studies have demonstrated that gene transfer events (e.g., plasmid acquisition) can entail a variety of costs that arise as by-products of the incorporation of foreign DNA into established physiological and genetic systems. These costs can be ameliorated through evolutionary time by the occurrence of compensatory mutations, which stabilize the presence of a horizontally transferred region within the genome but which also may skew future adaptive possibilities for these lineages. Here, we demonstrate another possible outcome, that phenotypic changes arising as a consequence of the same horizontal gene transfer (HGT) event are costly to some strains but may actually be beneficial in other genomic backgrounds under the right conditions. These results provide a new viewpoint for considering conditions that promote plasmid maintenance and highlight the influence of genomic and environmental contexts when considering amelioration of fitness costs after HGT events.

Is the X chromosome a hot spot for sexually antagonistic polymorphisms? Biases in current empirical tests of classical theory

Females and males carry nearly identical genomes, which can constrain the evolution of sexual dimorphism and generate conditions that are favourable for maintaining sexually antagonistic (SA) polymorphisms, in which alleles beneficial for one sex are deleterious for the other. An influential theoretical prediction, by Rice (Rice 1984 Evolution38, 735-742), is that the X chromosome should be a 'hot spot' (i.e. enriched) for SA polymorphisms. While important caveats to Rice's theoretical prediction have since been highlighted (e.g. by Fry (2010) Evolution64, 1510-1516), several empirical studies appear to support it. Here, we show that current tests of Rice's theory-most of which are based on quantitative genetic measures of fitness (co)variance-are frequently biased towards detecting X-linked effects. We show that X-linked genes tend to contribute disproportionately to quantitative genetic patterns of SA fitness variation whether or not the X is enriched for SA polymorphisms. Population genomic approaches for detecting SA loci, including genome-wide association study of fitness and analyses of intersexual FST, are similarly biased towards detecting X-linked effects. In the light of our models, we critically re-evaluate empirical evidence for Rice's theory and discuss prospects for empirically testing it.

A chromosome-scale draft sequence of the Canada fleabane genome

BACKGROUND

Due to the accessibility of underlying technologies the 'Omics', in particular genomics, are becoming commonplace in several fields of research, including the study of agricultural pests. The weed community is starting to embrace these approaches; genome sequences have been made available in the past years, with several other sequencing projects underway, as promoted by the International Weed Genome Consortium. Chromosome-scale sequences are essential to fully exploit the power of genetics and genomics.

RESULTS

We report such an assembly for Conyza canadensis, an important agricultural weed. Third-generation sequencing technology was used to create a genome assembly of 426 megabases, of which nine chromosome-scale scaffolds cover more than 98% of the entire assembled sequence. As this weed was the first to be identified with glyphosate resistance, and since we do not have a firm handle on the genetic mechanisms responsible for several herbicide resistances in the species, the genome sequence was annotated with genes known to be associated with herbicide resistance. A high number of ABC-type transporters, cytochrome P450 and glycosyltransferases (159, 352 and 181, respectively) were identified among the list of ab initio predicted genes.

CONCLUSION

As C. canadensis has a small genome that is syntenic with other Asteraceaes, has a short life cycle and is relatively easy to cross, it has the potential to become a model weed species and, with the chromosome-scale genome sequence, contribute to a paradigm shift in the way non-target site resistance is studied. © 2020 Her Majesty the Queen in Right of CanadaPest Management Science © 2020 Society of Chemical Industry.

Identifying and Classifying Shared Selective Sweeps from Multilocus Data

Positive selection causes beneficial alleles to rise to high frequency, resulting in a selective sweep of the diversity surrounding the selected sites. Accordingly, the signature of a selective sweep in an ancestral population may still remain in its descendants. Identifying signatures of selection in the ancestor that are shared among its descendants is important to contextualize the timing of a sweep, but few methods exist for this purpose. We introduce the statistic SS-H12, which can identify genomic regions under shared positive selection across populations and is based on the theory of the expected haplotype homozygosity statistic H12, which detects recent hard and soft sweeps from the presence of high-frequency haplotypes. SS-H12 is distinct from comparable statistics because it requires a minimum of only two populations, and properly identifies and differentiates between independent convergent sweeps and true ancestral sweeps, with high power and robustness to a variety of demographic models. Furthermore, we can apply SS-H12 in conjunction with the ratio of statistics we term [Formula: see text] and [Formula: see text] to further classify identified shared sweeps as hard or soft. Finally, we identified both previously reported and novel shared sweep candidates from human whole-genome sequences. Previously reported candidates include the well-characterized ancestral sweeps at LCT and SLC24A5 in Indo-Europeans, as well as GPHN worldwide. Novel candidates include an ancestral sweep at RGS18 in sub-Saharan Africans involved in regulating the platelet response and implicated in sudden cardiac death, and a convergent sweep at C2CD5 between European and East Asian populations that may explain their different insulin responses.

Impact of In Vivo Protein Folding Probability on Local Fitness Landscapes

It is incompletely understood how biophysical properties like protein stability impact molecular evolution and epistasis. Epistasis is defined as specific when a mutation exclusively influences the phenotypic effect of another mutation, often at physically interacting residues. In contrast, nonspecific epistasis results when a mutation is influenced by a large number of nonlocal mutations. As most mutations are pleiotropic, the in vivo folding probability-governed by basal protein stability-is thought to determine activity-enhancing mutational tolerance, implying that nonspecific epistasis is dominant. However, evidence exists for both specific and nonspecific epistasis as the prevalent factor, with limited comprehensive data sets to support either claim. Here, we use deep mutational scanning to probe how in vivo enzyme folding probability impacts local fitness landscapes. We computationally designed two different variants of the amidase AmiE with statistically indistinguishable catalytic efficiencies but lower probabilities of folding in vivo compared with wild-type. Local fitness landscapes show slight alterations among variants, with essentially the same global distribution of fitness effects. However, specific epistasis was predominant for the subset of mutations exhibiting positive sign epistasis. These mutations mapped to spatially distinct locations on AmiE near the initial mutation or proximal to the active site. Intriguingly, the majority of specific epistatic mutations were codon dependent, with different synonymous codons resulting in fitness sign reversals. Together, these results offer a nuanced view of how protein folding probability impacts local fitness landscapes and suggest that transcriptional-translational effects are as important as stability in determining evolutionary outcomes.

Challenges and Prospects for Testing the Mother's Curse Hypothesis

Maternal inheritance of mitochondrial DNA (mtDNA) renders selection blind to mutations whose effects are limited to males. Evolutionary theory predicts this will lead to the accumulation of a male-specific genetic load within the mitochondrial genomes of populations; that is, a pool of mutations that negatively affects male, but not female, fitness components. This principle has been termed the Mother's Curse hypothesis. While the hypothesis has received some empirical support, its relevance to natural populations of metazoans remains unclear, and these ambiguities are compounded by the lack of a clear predictive framework for studies attempting to test Mother's Curse. Here, we seek to redress this by outlining the core predictions of the hypothesis, as well as the key features of the experimental designs that are required to enable direct testing of the predictions. Our goal is to provide a roadmap for future research seeking to elucidate the evolutionary significance of the Mother's Curse hypothesis.

The Impact of Genomic and Traditional Selection on the Contribution of Mutational Variance to Long-Term Selection Response and Genetic Variance

De novo mutations (DNM) create new genetic variance and are an important driver for long-term selection response. We hypothesized that genomic selection exploits mutational variance less than traditional selection methods such as mass selection or selection on pedigree-based breeding values, because DNM in selection candidates are not captured when the selection candidates' own phenotype is not used in genomic selection, DNM are not on SNP chips and DNM are not in linkage disequilibrium with the SNP on the chip. We tested this hypothesis with Monte Carlo simulation. From whole-genome sequence data, a subset of ∼300,000 variants was used that served as putative markers, quantitative trait loci or DNM. We simulated 20 generations with truncation selection based on breeding values from genomic best linear unbiased prediction without (GBLUP_no_OP) or with own phenotype (GBLUP_OP), pedigree-based BLUP without (BLUP_no_OP) or with own phenotype (BLUP_OP), or directly on phenotype. GBLUP_OP was the best strategy in exploiting mutational variance, while GBLUP_no_OP and BLUP_no_OP were the worst in exploiting mutational variance. The crucial element is that GBLUP_no_OP and BLUP_no_OP puts no selection pressure on DNM in selection candidates. Genetic variance decreased faster with GBLUP_no_OP and GBLUP_OP than with BLUP_no_OP, BLUP_OP or mass selection. The distribution of mutational effects, mutational variance, number of DNM per individual and nonadditivity had a large impact on mutational selection response and mutational genetic variance, but not on ranking of selection strategies. We advocate that more sustainable genomic selection strategies are required to optimize long-term selection response and to maintain genetic diversity.

Adaptation by Loss of Heterozygosity in Saccharomyces cerevisiae Clones Under Divergent Selection

Loss of heterozygosity (LOH) is observed during vegetative growth and reproduction of diploid genotypes through mitotic crossovers, aneuploidy caused by nondisjunction, and gene conversion. We aimed to test the role that LOH plays during adaptation of two highly heterozygous Saccharomyces cerevisiae genotypes to multiple environments over a short time span in the laboratory. We hypothesized that adaptation would be observed through parallel LOH events across replicate populations. Using genome resequencing of 70 clones, we found that LOH was widespread with 5.2 LOH events per clone after ∼500 generations. The most common mode of LOH was gene conversion (51%) followed by crossing over consistent with either break-induced replication or double Holliday junction resolution. There was no evidence that LOH involved nondisjunction of whole chromosomes. We observed parallel LOH in both an environment-specific and environment-independent manner. LOH largely involved recombining existing variation between the parental genotypes, but also was observed after de novo, presumably beneficial, mutations occurred in the presence of canavanine, a toxic analog of arginine. One highly parallel LOH event involved the ENA salt efflux pump locus on chromosome IV, which showed repeated LOH to the allele from the European parent, an allele originally derived by introgression from S. paradoxus Using CRISPR-engineered LOH we showed that the fitness advantage provided by this single LOH event was 27%. Overall, we found extensive evidence that LOH could be adaptive and is likely to be a greater source of initial variation than de novo mutation for rapid evolution of diploid genotypes.

The population genetics of crypsis in vertebrates: recent insights from mice, hares, and lizards

By combining well-established population genetic theory with high-throughput sequencing data from natural populations, major strides have recently been made in understanding how, why, and when vertebrate populations evolve crypsis. Here, we focus on background matching, a particular facet of crypsis that involves the ability of an organism to conceal itself through matching its color to the surrounding environment. While interesting in and of itself, the study of this phenotype has also provided fruitful population genetic insights into the interplay of strong positive selection with other evolutionary processes. Specifically, and predicated upon the findings of previous candidate gene association studies, a primary focus of this recent literature involves the realization that the inference of selection from DNA sequence data first requires a robust model of population demography in order to identify genomic regions which do not conform to neutral expectations. Moreover, these demographic estimates provide crucial information about the origin and timing of the onset of selective pressures associated with, for example, the colonization of a novel environment. Furthermore, such inference has revealed crypsis to be a particularly useful phenotype for investigating the interplay of migration and selection—with examples of gene flow constraining rates of adaptation, or alternatively providing the genetic variants that may ultimately sweep through the population. Here, we evaluate the underlying evidence, review the strengths and weaknesses of the many population genetic methodologies used in these studies, and discuss how these insights have aided our general understanding of the evolutionary process.

Molecular Origins of Complex Heritability in Natural Genotype-to-Phenotype Relationships

The statistical complexity of heredity has long been evident, but its molecular origins remain elusive. To investigate, we charted 90 comprehensive genotype-to-phenotype maps in a large population of wild diploid yeast. In contrast to long-standing assumptions, all types of genetic variation contributed similarly to phenotype. Causal synonymous and regulatory variants exhibited distinct molecular signatures, as did nonlinearities in heterozygote fitness that likely contribute to hybrid vigor. Highly pleiotropic variants altered disordered sequences within signaling hubs, and their effects correlated across environments-even when antagonistic-suggesting that large fitness gains bring concomitant costs. Natural genetic networks defined by the causal loci differed from those determined by precise gene deletions or protein-protein interactions. Finally, we found that traits that would appear omnigenic in less powered studies do in fact have finite genetic determinants. Integrating these molecular principles will be crucial as genome reading and writing become routine in research, industry, and medicine.

Dominance reversals and the maintenance of genetic variation for fitness

Antagonistic selection between different fitness components (e.g., survival versus fertility) or different types of individuals in a population (e.g., females versus males) can potentially maintain genetic diversity and thereby account for the high levels of fitness variation observed in natural populations. However, the degree to which antagonistic selection can maintain genetic variation critically depends on the dominance relations between antagonistically selected alleles in diploid individuals. Conditions for stable polymorphism of antagonistically selected alleles are narrow, particularly when selection is weak, unless the alleles exhibit "dominance reversals"-in which each allele is partially or completely dominant in selective contexts in which it is favored and recessive in contexts in which it is harmful. Although theory predicts that dominance reversals should emerge under biologically plausible conditions, evidence for dominance reversals is sparse. In this primer, we review theoretical arguments and data supporting a role for dominance reversals in the maintenance of genetic variation. We then highlight an illuminating new study by Grieshop and Arnqvist, which reports a genome-wide signal of dominance reversals between male and female fitness in seed beetles.

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast-Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.

Quantifying the relative contributions of the X chromosome, autosomes, and mitochondrial genome to local adaptation

During local adaptation with gene flow, some regions of the genome are inherently more responsive to selection than others. Recent theory predicts that X-linked genes should disproportionately contribute to local adaptation relative to other genomic regions, yet this prediction remains to be tested. We carried out a multigeneration crossing scheme, using two cline-end populations of Drosophila melanogaster, to estimate the relative contributions of the X chromosome, autosomes, and mitochondrial genome to divergence in four traits involved in local adaptation (wing size, resistance to heat, desiccation, and starvation stresses). We found that the mitochondrial genome and autosomes contributed significantly to clinal divergence in three of the four traits. In contrast, the X made no significant contribution to divergence in these traits. Given the small size of the mitochondrial genome, our results indicate that it plays a surprisingly large role in clinal adaptation. In contrast, the X, which represents roughly 20% of the Drosophila genome, contributes negligibly-a pattern that conflicts with theoretical predictions. These patterns reinforce recent work implying a central role of mitochondria in climatic adaptation, and suggest that different genomic regions may play fundamentally different roles in processes of divergence with gene flow.

The fitness landscape of the codon space across environments

Fitness landscapes map the relationship between genotypes and fitness. However, most fitness landscape studies ignore the genetic architecture imposed by the codon table and thereby neglect the potential role of synonymous mutations. To quantify the fitness effects of synonymous mutations and their potential impact on adaptation on a fitness landscape, we use a new software based on Bayesian Monte Carlo Markov Chain methods and re-estimate selection coefficients of all possible codon mutations across 9 amino acid positions in Saccharomyces cerevisiae Hsp90 across 6 environments. We quantify the distribution of fitness effects of synonymous mutations and show that it is dominated by many mutations of small or no effect and few mutations of larger effect. We then compare the shape of the codon fitness landscape across amino acid positions and environments, and quantify how the consideration of synonymous fitness effects changes the evolutionary dynamics on these fitness landscapes. Together these results highlight a possible role of synonymous mutations in adaptation and indicate the potential mis-inference when they are neglected in fitness landscape studies.

Insights into the Structure of the Spruce Budworm (Choristoneura fumiferana) Genome, as Revealed by Molecular Cytogenetic Analyses and a High-Density Linkage Map

Genome structure characterization can contribute to a better understanding of processes such as adaptation, speciation, and karyotype evolution, and can provide useful information for refining genome assemblies. We studied the genome of an important North American boreal forest pest, the spruce budworm, Choristoneura fumiferana, through a combination of molecular cytogenetic analyses and construction of a high-density linkage map based on single nucleotide polymorphism (SNP) markers obtained through a genotyping-by-sequencing (GBS) approach. Cytogenetic analyses using fluorescence in situ hybridization methods confirmed the haploid chromosome number of n = 30 in both sexes of C. fumiferana and showed, for the first time, that this species has a WZ/ZZ sex chromosome system. Synteny analysis based on a comparison of the Bombyx mori genome and the C. fumiferana linkage map revealed the presence of a neo-Z chromosome in the latter species, as previously reported for other tortricid moths. In this neo-Z chromosome, we detected an ABC transporter C2 (ABCC2) gene that has been associated with insecticide resistance. Sex-linkage of the ABCC2 gene provides a genomic context favorable to selection and rapid spread of resistance against Bacillus thuringiensis serotype kurstaki (Btk), the main insecticide used in Canada to control spruce budworm populations. Ultimately, the linkage map we developed, which comprises 3586 SNP markers distributed over 30 linkage groups for a total length of 1720.41 cM, will be a valuable tool for refining our draft assembly of the spruce budworm genome.

Strong disorder leads to scale invariance in complex biological systems

Despite the innate complexity of the cell, emergent scale-invariant behavior is observed in many biological systems. We investigate one example of this phenomenon: the dynamics of large complexes in the bacterial cytoplasm. The observed dynamics of these complexes is scale invariant in three measures of dynamics: mean-squared displacement (MSD), velocity autocorrelation function, and the step-size distribution. To investigate the physical mechanism for this emergent scale invariance, we explore minimal models in which mobility is modeled as diffusion on a rough free-energy landscape in one dimension. We discover that all three scale-invariant characteristics emerge generically in the strong disorder limit. (Strong disorder is defined by the divergence of the ensemble-averaged hop time between lattice sites.) In particular, we demonstrate how the scale invariance of the relative step-size distribution can be understood from the perspective of extreme-value theory in statistics (EVT). We show that the Gumbel scale parameter is simply related to the MSD scaling parameter. The EVT mechanism of scale invariance is expected to be generic to strongly disordered systems and therefore a powerful tool for the analysis of other systems in biology and beyond.

Temperature effects on life-history trade-offs, germline maintenance and mutation rate under simulated climate warming

Mutation has a fundamental influence over evolutionary processes, but how evolutionary processes shape mutation rate remains less clear. In asexual unicellular organism, increased mutation rates have been observed in stressful environments and the reigning paradigm ascribes this increase to selection for evolvability. However, this explanation does not apply in sexually reproducing species, where little is known about how the environment affects mutation rate. Here we challenged experimental lines of seed beetle, evolved at ancestral temperature or under simulated climate warming, to repair induced mutations at ancestral and stressful temperature. Results show that temperature stress causes individuals to pass on a greater mutation load to their grand-offspring. This suggests that stress-induced mutation rates, in unicellular and multicellular organisms alike, can result from compromised germline DNA repair in low condition individuals. Moreover, lines adapted to simulated climate warming had evolved increased longevity at the cost of reproduction, and this allocation decision improved germline repair. These results suggest that mutation rates can be modulated by resource allocation trade-offs encompassing life-history traits and the germline and have important implications for rates of adaptation and extinction as well as our understanding of genetic diversity in multicellular organisms.

Increase in viability due to the accumulation of X chromosome mutations in Drosophila melanogaster males

To increase our understanding of the role of new X-chromosome mutations in adaptive evolution, single-X Drosophila melanogaster males were mated with attached-X chromosome females, allowing the male X chromosome to accumulate mutations over 28 generations. Contrary to our hypothesis that male viability would decrease over time, due to the accumulation and expression of X-linked recessive deleterious mutations in hemizygous males, viability significantly increased. This increase may be attributed to germinal selection and to new X-linked beneficial or compensatory mutations, possibly supporting the faster-X hypothesis.

Selfing Can Facilitate Transitions between Pollination Syndromes

Pollinator-mediated selection on plants can favor transitions to a new pollinator depending on the relative abundances and efficiencies of pollinators present in the community. A frequently observed example is the transition from bee pollination to hummingbird pollination. We present a population genetic model that examines whether the ability to inbreed can influence evolutionary change in traits that underlie pollinator attraction. We find that a transition to a more efficient but less abundant pollinator is favored under a broadened set of ecological conditions if plants are capable of delayed selfing rather than obligately outcrossing. Delayed selfing allows plants carrying an allele that attracts the novel pollinator to reproduce even when this pollinator is rare, providing reproductive assurance. In addition, delayed selfing weakens the effects of Haldane's sieve by increasing the fixation probability for recessive alleles that confer adaptation to the new pollinator. Our model provides novel insight into the paradoxical abundance of recessive mutations in adaptation to hummingbird attraction. It further predicts that transitions to efficient but less abundant pollinators (such as hummingbirds in certain communities) should disproportionately occur in self-compatible lineages. Currently available mating system data sets are consistent with this prediction, and we suggest future areas of research that will enable a rigorous test of this theory.

Genetic constraints on adaptation: a theoretical primer for the genomics era

Genetic constraints are features of inheritance systems that slow or prohibit adaptation. Several population genetic mechanisms of constraint have received sustained attention within the field since they were first articulated in the early 20th century. This attention is now reflected in a rich, and still growing, theoretical literature on the genetic limits to adaptive change. In turn, empirical research on constraints has seen a rapid expansion over the last two decades in response to changing interests of evolutionary biologists, along with new technologies, expanding data sets, and creative analytical approaches that blend mathematical modeling with genomics. Indeed, one of the most notable and exciting features of recent progress in genetic constraints is the close connection between theoretical and empirical research. In this review, we discuss five major population genetic contexts of genetic constraint: genetic dominance, pleiotropy, fitness trade-offs between types of individuals of a population, sign epistasis, and genetic linkage between loci. For each, we outline historical antecedents of the theory, specific contexts where constraints manifest, and their quantitative consequences for adaptation. From each of these theoretical foundations, we discuss recent empirical approaches for identifying and characterizing genetic constraints, each grounded and motivated by this theory, and outline promising areas for future work.

Critical Mutation Rate has an Exponential Dependence on Population Size for Eukaryotic-length Genomes with Crossover

The critical mutation rate (CMR) determines the shift between survival-of-the-fittest and survival of individuals with greater mutational robustness ("flattest"). We identify an inverse relationship between CMR and sequence length in an in silico system with a two-peak fitness landscape; CMR decreases to no more than five orders of magnitude above estimates of eukaryotic per base mutation rate. We confirm the CMR reduces exponentially at low population sizes, irrespective of peak radius and distance, and increases with the number of genetic crossovers. We also identify an inverse relationship between CMR and the number of genes, confirming that, for a similar number of genes to that for the plant Arabidopsis thaliana (25,000), the CMR is close to its known wild-type mutation rate; mutation rates for additional organisms were also found to be within one order of magnitude of the CMR. This is the first time such a simulation model has been assigned input and produced output within range for a given biological organism. The decrease in CMR with population size previously observed is maintained; there is potential for the model to influence understanding of populations undergoing bottleneck, stress, and conservation strategy for populations near extinction.

Single-mutation fitness landscapes for an enzyme on multiple substrates reveal specificity is globally encoded

Our lack of total understanding of the intricacies of how enzymes behave has constrained our ability to robustly engineer substrate specificity. Furthermore, the mechanisms of natural evolution leading to improved or novel substrate specificities are not wholly defined. Here we generate near-comprehensive single-mutation fitness landscapes comprising >96.3% of all possible single nonsynonymous mutations for hydrolysis activity of an amidase expressed in E. coli with three different substrates. For all three selections, we find that the distribution of beneficial mutations can be described as exponential, supporting a current hypothesis for adaptive molecular evolution. Beneficial mutations in one selection have essentially no correlation with fitness for other selections and are dispersed throughout the protein sequence and structure. Our results further demonstrate the dependence of local fitness landscapes on substrate identity and provide an example of globally distributed sequence-specificity determinants for an enzyme.

Systematically understanding the sequence determinants to substrate specificity for enzymes has implications in areas from evolutionary biology to biocatalysis. Here, Whitehead and colleagues generate and analyse near-comprehensive single-mutation fitness landscapes for an amidase with three different substrates.

The Relative Contributions of the X Chromosome and Autosomes to Local Adaptation

Models of sex chromosome and autosome evolution yield key predictions about the genomic basis of adaptive divergence, and such models have been important in guiding empirical research in comparative genomics and studies of speciation. In addition to the adaptive differentiation that occurs between species over time, selection also favors genetic divergence across geographic space, with subpopulations of single species evolving conspicuous differences in traits involved in adaptation to local environmental conditions. The potential contribution of sex chromosomes (the X or Z) to local adaptation remains unclear, as we currently lack theory that directly links spatial variation in selection to local adaptation of X-linked and autosomal genes. Here, we develop population genetic models that explicitly consider the effects of genetic dominance, effective population size, and sex-specific migration and selection on the relative contributions of X-linked and autosomal genes to local adaptation. We show that X-linked genes should nearly always disproportionately contribute to local adaptation in the presence of gene flow. We also show that considerations of dominance and effective population size-which play pivotal roles in the theory of faster-X adaptation between species-have surprisingly little influence on the relative contribution of the X chromosome to local adaptation. Instead, sex-biased migration is the primary mediator of the strength of spatial large-X effects. Our results yield novel predictions about the role of sex chromosomes in local adaptation. We outline empirical approaches in evolutionary quantitative genetics and genomics that could build upon this new theory.

Genomic Trajectories to Desiccation Resistance: Convergence and Divergence Among Replicate Selected Drosophila Lines

Adaptation to environmental stress is critical for long-term species persistence. With climate change and other anthropogenic stressors compounding natural selective pressures, understanding the nature of adaptation is as important as ever in evolutionary biology. In particular, the number of alternative molecular trajectories available for an organism to reach the same adaptive phenotype remains poorly understood. Here, we investigate this issue in a set of replicated Drosophila melanogaster lines selected for increased desiccation resistance-a classical physiological trait that has been closely linked to Drosophila species distributions. We used pooled whole-genome sequencing (Pool-Seq) to compare the genetic basis of their selection responses, using a matching set of replicated control lines for characterizing laboratory (lab-)adaptation, as well as the original base population. The ratio of effective population size to census size was high over the 21 generations of the experiment at 0.52-0.88 for all selected and control lines. While selected SNPs in replicates of the same treatment (desiccation-selection or lab-adaptation) tended to change frequency in the same direction, suggesting some commonality in the selection response, candidate SNP and gene lists often differed among replicates. Three of the five desiccation-selection replicates showed significant overlap at the gene and network level. All five replicates showed enrichment for ovary-expressed genes, suggesting maternal effects on the selected trait. Divergence between pairs of replicate lines for desiccation-candidate SNPs was greater than between pairs of control lines. This difference also far exceeded the divergence between pairs of replicate lines for neutral SNPs. Overall, while there was overlap in the direction of allele frequency changes and the network and functional categories affected by desiccation selection, replicates showed unique responses at all levels, likely reflecting hitchhiking effects, and highlighting the challenges in identifying candidate genes from these types of experiments when traits are likely to be polygenic.

Adaptation prevents the extinction of Chlamydomonas reinhardtii under toxic beryllium

The current biodiversity crisis represents a historic challenge for natural communities: the environmental rate of change exceeds the population’s adaptation capability. Integrating both ecological and evolutionary responses is necessary to make reliable predictions regarding the loss of biodiversity. The race against extinction from an eco-evolutionary perspective is gaining importance in ecological risk assessment. Here, we performed a classical study of population dynamics—a fluctuation analysis—and evaluated the results from an adaption perspective. Fluctuation analysis, widely used with microorganisms, is an effective empirical procedure to study adaptation under strong selective pressure because it incorporates the factors that influence demographic, genetic and environmental changes. The adaptation of phytoplankton to beryllium (Be) is of interest because human activities are increasing the concentration of Be in freshwater reserves; therefore, predicting the effects of human-induced pollutants is necessary for proper risk assessment. The fluctuation analysis was performed with phytoplankton, specifically, the freshwater microalgae Chlamydomonas reinhardtii, under acute Be exposure. High doses of Be led to massive microalgae death; however, by conducting a fluctuation analysis experiment, we found that C. reinhardtii was able to adapt to 33 mg/l of Be due to pre-existing genetic variability. The rescuing adapting genotype presented a mutation rate of 9.61 × 10⁻⁶ and a frequency of 10.42 resistant cells per million wild-type cells. The genetic adaptation pathway that was experimentally obtained agreed with the theoretical models of evolutionary rescue (ER). Furthermore, the rescuing genotype presented phenotypic and physiologic differences from the wild-type genotype, was 25% smaller than the Be-resistant genotype and presented a lower fitness and quantum yield performance. The abrupt distinctions between the wild-type and the Be-resistant genotype suggest a pleiotropic effect mediated by an advantageous mutation; however, no sequencing confirmation was performed.

The implications of small stem cell niche sizes and the distribution of fitness effects of new mutations in aging and tumorigenesis

Somatic tissue evolves over a vertebrate's lifetime due to the accumulation of mutations in stem cell populations. Mutations may alter cellular fitness and contribute to tumorigenesis or aging. The distribution of mutational effects within somatic cells is not known. Given the unique regulatory regime of somatic cell division, we hypothesize that mutational effects in somatic tissue fall into a different framework than whole organisms; one in which there are more mutations of large effect. Through simulation analysis, we investigate the fit of tumor incidence curves generated using exponential and power‐law distributions of fitness effects (DFE) to known tumorigenesis incidence. Modeling considerations include the architecture of stem cell populations, that is, a large number of very small populations, and mutations that do and do not fix neutrally in the stem cell niche. We find that the typically quantified DFE in whole organisms is sufficient to explain tumorigenesis incidence. Further, deleterious mutations are predicted to accumulate via genetic drift, resulting in reduced tissue maintenance. Thus, despite there being a large number of stem cells throughout the intestine, its compartmental architecture leads to the accumulation of deleterious mutations and significant aging, making the intestinal stem cell niche a prime example of Muller's Ratchet.

Greedy adaptive walks on a correlated fitness landscape

We study adaptation of a haploid asexual population on a fitness landscape defined over binary genotype sequences of length L. We consider greedy adaptive walks in which the population moves to the fittest among all single mutant neighbors of the current genotype until a local fitness maximum is reached. The landscape is of the rough mount Fuji type, which means that the fitness value assigned to a sequence is the sum of a random and a deterministic component. The random components are independent and identically distributed random variables, and the deterministic component varies linearly with the distance to a reference sequence. The deterministic fitness gradient c is a parameter that interpolates between the limits of an uncorrelated random landscape (c=0) and an effectively additive landscape (c→∞). When the random fitness component is chosen from the Gumbel distribution, explicit expressions for the distribution of the number of steps taken by the greedy walk are obtained, and it is shown that the walk length varies non-monotonically with the strength of the fitness gradient when the starting point is sufficiently close to the reference sequence. Asymptotic results for general distributions of the random fitness component are obtained using extreme value theory, and it is found that the walk length attains a non-trivial limit for L→∞, different from its values for c=0 and c=∞, if c is scaled with L in an appropriate combination.

Standing genetic variation as a major contributor to adaptation in the Virginia chicken lines selection experiment

Background

Artificial selection provides a powerful approach to study the genetics of adaptation. Using selective-sweep mapping, it is possible to identify genomic regions where allele-frequencies have diverged during selection. To avoid false positive signatures of selection, it is necessary to show that a sweep affects a selected trait before it can be considered adaptive. Here, we confirm candidate, genome-wide distributed selective sweeps originating from the standing genetic variation in a long-term selection experiment on high and low body weight of chickens.

Results

Using an intercross between the two divergent chicken lines, 16 adaptive selective sweeps were confirmed based on their association with the body weight at 56 days of age. Although individual additive effects were small, the fixation for alternative alleles across the loci contributed at least 40 % of the phenotypic difference for the selected trait between these lines. The sweeps contributed about half of the additive genetic variance present within and between the lines after 40 generations of selection, corresponding to a considerable portion of the additive genetic variance of the base population.

Conclusions

Long-term, single-trait, bi-directional selection in the Virginia chicken lines has resulted in a gradual response to selection for extreme phenotypes without a drastic reduction in the genetic variation. We find that fixation of several standing genetic variants across a highly polygenic genetic architecture made a considerable contribution to long-term selection response. This provides new fundamental insights into the dynamics of standing genetic variation during long-term selection and adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0785-z) contains supplementary material, which is available to authorized users.

The dominance effect of the adaptive transposable element insertion Bari-Jheh depends on the genetic background

Although adaptive mutations are often considered to be dominant, it has been recently shown that a substantial proportion of adaptive mutations should display heterozygote advantage. In this work, we take advantage of a recently characterized transposable element insertion mediating oxidative stress response in Drosophila melanogaster to test the dominance effect of an adaptive mutation. The comparison of the survival curves of heterozygous and the two corresponding homozygous flies indicated that the dominance effect of Bari-Jheh depends on the genetic background. Both in homozygous and in heterozygous flies, Bari-Jheh was associated with upregulation of Jheh1 (Juvenile Hormone Epoxyde Hydrolase 1) and/or Jheh2 genes. Our results add to the limited number of studies in which the dominance effect of adaptive mutations has been empirically estimated and highlights the complexity of their inheritance.

The evolutionarily stable distribution of fitness effects

The distribution of fitness effects (DFE) of new mutations is a key parameter in determining the course of evolution. This fact has motivated extensive efforts to measure the DFE or to predict it from first principles. However, just as the DFE determines the course of evolution, the evolutionary process itself constrains the DFE. Here, we analyze a simple model of genome evolution in a constant environment in which natural selection drives the population toward a dynamic steady state where beneficial and deleterious substitutions balance. The distribution of fitness effects at this steady state is stable under further evolution and provides a natural null expectation for the DFE in a population that has evolved in a constant environment for a long time. We calculate how the shape of the evolutionarily stable DFE depends on the underlying population genetic parameters. We show that, in the absence of epistasis, the ratio of beneficial to deleterious mutations of a given fitness effect obeys a simple relationship independent of population genetic details. Finally, we analyze how the stable DFE changes in the presence of a simple form of diminishing-returns epistasis.

Sex-specific trans-regulatory variation on the Drosophila melanogaster X chromosome

The X chromosome constitutes a unique genomic environment because it is present in one copy in males, but two copies in females. This simple fact has motivated several theoretical predictions with respect to how standing genetic variation on the X chromosome should differ from the autosomes. Unmasked expression of deleterious mutations in males and a lower census size are expected to reduce variation, while allelic variants with sexually antagonistic effects, and potentially those with a sex-specific effect, could accumulate on the X chromosome and contribute to increased genetic variation. In addition, incomplete dosage compensation of the X chromosome could potentially dampen the male-specific effects of random mutations, and promote the accumulation of X-linked alleles with sexually dimorphic phenotypic effects. Here we test both the amount and the type of genetic variation on the X chromosome within a population of Drosophila melanogaster, by comparing the proportion of X linked and autosomal trans-regulatory SNPs with a sexually concordant and discordant effect on gene expression. We find that the X chromosome is depleted for SNPs with a sexually concordant effect, but hosts comparatively more SNPs with a sexually discordant effect. Interestingly, the contrasting results for SNPs with sexually concordant and discordant effects are driven by SNPs with a larger influence on expression in females than expression in males. Furthermore, the distribution of these SNPs is shifted towards regions where dosage compensation is predicted to be less complete. These results suggest that intrinsic properties of dosage compensation influence either the accumulation of different types of trans-factors and/or their propensity to accumulate mutations. Our findings document a potential mechanistic basis for sex-specific genetic variation, and identify the X as a reservoir for sexually dimorphic phenotypic variation. These results have general implications for X chromosome evolution, as well as the genetic basis of sex-specific evolutionary change.