Revisiting the notion of deleterious sweeps

Do chromosome rearrangements fix by genetic drift or natural selection? Insights from Brenthis butterflies

Large-scale chromosome rearrangements, such as fissions and fusions, are a common feature of eukaryote evolution. They can have considerable influence on the evolution of populations, yet it remains unclear exactly how rearrangements become established and eventually fix. Rearrangements could fix by genetic drift if they are weakly deleterious or neutral, or they may instead be favoured by positive natural selection. Here, we compare genome assemblies of three closely related Brenthis butterfly species and characterize a complex history of fission and fusion rearrangements. An inferred demographic history of these species suggests that rearrangements became fixed in populations with large long-term effective size (N e), consistent with rearrangements being selectively neutral or only very weakly underdominant. Using a recently developed analytic framework for characterizing hard selective sweeps, we find that chromosome fusions are not enriched for evidence of past sweeps compared to other regions of the genome. Nonetheless, we do infer a strong and recent selective sweep around one chromosome fusion in the B. daphne genome. Our results suggest that rearrangements in these species likely have weak absolute fitness effects and fix by genetic drift. However, one putative selective sweep raises the possibility that natural selection may sometimes play a role in the fixation of chromosome fusions.

Genetic diversity during selective sweeps in non-recombining populations

Selective sweeps, resulting from the spread of beneficial, neutral, or deleterious mutations through a population, shape patterns of genetic variation at linked neutral sites. While many theoretical, computational, and statistical advances have been made in understanding the genomic signatures of selective sweeps in recombining populations, substantially less is understood in populations with little/no recombination. We present a mathematical framework based on diffusion theory for obtaining the site frequency spectrum (SFS) at linked neutral sites immediately post and during the fixation of moderately or strongly beneficial mutations. We find that when a single hard sweep occurs, the SFS decays as 1/ x for low derived allele frequencies ( x ), similar to the neutral SFS at equilibrium, whereas at higher derived allele frequencies, it follows a 1/ x 2 power law. These power laws are universal in the sense that they are independent of the dominance and inbreeding coefficient, and also characterize the SFS during the sweep. Additionally, we find that the derived allele frequency where the SFS shifts from the 1/ x to 1/ x 2 law, is inversely proportional to the selection strength: thus under strong selection, the SFS follows the 1/ x 2 dependence for most allele frequencies, resembling a rapidly expanding neutral population. When clonal interference is pervasive, the SFS immediately post-fixation becomes U-shaped and is better explained by the equilibrium SFS of selected sites. Our results will be important in developing statistical methods to infer the timing and strength of recent selective sweeps in asexual populations, genomic regions that lack recombination, and clonally propagating tumor populations.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counterintuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here, we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify seven deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild nonhuman primates.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counter intuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify 7 deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild non-human primates.

Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modeled by a realistic mutation rate and as part of a realistic distribution of fitness effects, as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modeled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false-positive rates are in excess of true-positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

Genetic Drift Shapes the Evolution of a Highly Dynamic Metapopulation

The dynamics of extinction and (re)colonization in habitat patches are characterizing features of dynamic metapopulations, causing them to evolve differently than large, stable populations. The propagule model, which assumes genetic bottlenecks during colonization, posits that newly founded subpopulations have low genetic diversity and are genetically highly differentiated from each other. Immigration may then increase diversity and decrease differentiation between subpopulations. Thus, older and/or less isolated subpopulations are expected to have higher genetic diversity and less genetic differentiation. We tested this theory using whole-genome pool-sequencing to characterize nucleotide diversity and differentiation in 60 subpopulations of a natural metapopulation of the cyclical parthenogen Daphnia magna. For comparison, we characterized diversity in a single, large, and stable D. magna population. We found reduced (synonymous) genomic diversity, a proxy for effective population size, weak purifying selection, and low rates of adaptive evolution in the metapopulation compared with the large, stable population. These differences suggest that genetic bottlenecks during colonization reduce effective population sizes, which leads to strong genetic drift and reduced selection efficacy in the metapopulation. Consistent with the propagule model, we found lower diversity and increased differentiation in younger and also in more isolated subpopulations. Our study sheds light on the genomic consequences of extinction-(re)colonization dynamics to an unprecedented degree, giving strong support for the propagule model. We demonstrate that the metapopulation evolves differently from a large, stable population and that evolution is largely driven by genetic drift.

The effects of weak selection on neutral diversity at linked sites

The effects of selection on variability at linked sites have an important influence on levels and patterns of within-population variation across the genome. Most theoretical models of these effects have assumed that selection is sufficiently strong that allele frequency changes at the loci concerned are largely deterministic. These models have led to the conclusion that directional selection for selectively favorable mutations, or against recurrent deleterious mutations, reduces nucleotide site diversity at linked neutral sites. Recent work has shown, however, that fixations of weakly selected mutations, accompanied by significant stochastic changes in allele frequencies, can sometimes cause higher diversity at linked sites when compared with the effects of fixations of neutral mutations. This study extends this work by deriving approximate expressions for the mean conditional times to fixation and loss of mutations subject to selection, and analyzing the conditions under which selection increases rather than reduces these times. Simulations are used to examine the relations between diversity at a neutral site and the fixation and loss times of mutations at a linked site that is subject to selection. It is shown that the long-term level of neutral diversity can be increased over the purely neutral value by recurrent fixations and losses of linked, weakly selected dominant or partially dominant favorable mutations, or linked recessive or partially recessive deleterious mutations. The results are used to examine the conditions under which associative overdominance, as opposed to background selection, is likely to operate.

Strong neutral sweeps occurring during a population contraction

A strong reduction in diversity around a specific locus is often interpreted as a recent rapid fixation of a positively selected allele, a phenomenon called a selective sweep. Rapid fixation of neutral variants can however lead to a similar reduction in local diversity, especially when the population experiences changes in population size, e.g. bottlenecks or range expansions. The fact that demographic processes can lead to signals of nucleotide diversity very similar to signals of selective sweeps is at the core of an ongoing discussion about the roles of demography and natural selection in shaping patterns of neutral variation. Here, we quantitatively investigate the shape of such neutral valleys of diversity under a simple model of a single population size change, and we compare it to signals of a selective sweep. We analytically describe the expected shape of such "neutral sweeps" and show that selective sweep valleys of diversity are, for the same fixation time, wider than neutral valleys. On the other hand, it is always possible to parametrize our model to find a neutral valley that has the same width as a given selected valley. Our findings provide further insight into how simple demographic models can create valleys of genetic diversity similar to those attributed to positive selection.

Detecting Selection in Multiple Populations by Modeling Ancestral Admixture Components

One of the most powerful and commonly used approaches for detecting local adaptation in the genome is the identification of extreme allele frequency differences between populations. In this article, we present a new maximum likelihood method for finding regions under positive selection. It is based on a Gaussian approximation to allele frequency changes and it incorporates admixture between populations. The method can analyze multiple populations simultaneously and retains power to detect selection signatures specific to ancestry components that are not representative of any extant populations. Using simulated data, we compare our method to related approaches, and show that it is orders of magnitude faster than the state-of-the-art, while retaining similar or higher power for most simulation scenarios. We also apply it to human genomic data and identify loci with extreme genetic differentiation between major geographic groups. Many of the genes identified are previously known selected loci relating to hair pigmentation and morphology, skin, and eye pigmentation. We also identify new candidate regions, including various selected loci in the Native American component of admixed Mexican-Americans. These involve diverse biological functions, such as immunity, fat distribution, food intake, vision, and hair development.

Time to fixation in changing environments

Although many experimental and theoretical studies on natural selection have been carried out in a constant environment, as natural environments typically vary in time, it is important to ask if and how the results of these investigations are affected by a changing environment. Here, we study the properties of the conditional fixation time defined as the time to fixation of a new mutant that is destined to fix in a finite, randomly mating diploid population with intermediate dominance that is evolving in a periodically changing environment. It is known that in a static environment, the conditional mean fixation time of a co-dominant beneficial mutant is equal to that of a deleterious mutant with the same magnitude of selection coefficient. We find that this symmetry is not preserved, even when the environment is changing slowly. More generally, we find that the conditional mean fixation time of an initially beneficial mutant in a slowly changing environment depends weakly on the dominance coefficient and remains close to the corresponding result in the static environment. However, for an initially deleterious mutant under moderate and slowly varying selection, the fixation time differs substantially from that in a constant environment when the mutant is recessive. As fixation times are intimately related to the levels and patterns of genetic diversity, our results suggest that for beneficial sweeps, these quantities are only mildly affected by temporal variation in environment. In contrast, environmental change is likely to impact the patterns due to recessive deleterious sweeps strongly.

Effects of Selection at Linked Sites on Patterns of Genetic Variability

Patterns of variation and evolution at a given site in a genome can be strongly influenced by the effects of selection at genetically linked sites. In particular, the recombination rates of genomic regions correlate with their amount of within-population genetic variability, the degree to which the frequency distributions of DNA sequence variants differ from their neutral expectations, and the levels of adaptation of their functional components. We review the major population genetic processes that are thought to lead to these patterns, focusing on their effects on patterns of variability: selective sweeps, background selection, associative overdominance, and Hill–Robertson interference among deleterious mutations. We emphasize the difficulties in distinguishing among the footprints of these processes and disentangling them from the effects of purely demographic factors such as population size changes. We also discuss how interactions between selective and demographic processes can significantly affect patterns of variability within genomes.

Decoding 'Unnecessary Complexity': A Law of Complexity and a Concept of Hidden Variation Behind "Missing Heritability" in Precision Medicine

The high hopes for the Human Genome Project and personalized medicine were not met because the relationship between genotypes and phenotypes turned out to be more complex than expected. In a previous study we laid the foundation of a theory of complexity and showed that because of the blind nature of evolution, and molecular and historical contingency, cells have accumulated unnecessary complexity, complexity beyond what is necessary and sufficient to describe an organism. Here we provide empirical evidence and show that unnecessary complexity has become integrated into the genome in the form of redundancy and is relevant to molecular evolution of phenotypic complexity. Unnecessary complexity creates uncertainty between molecular and phenotypic complexity, such that phenotypic complexity (CP) is higher than molecular complexity (CM), which is higher than DNA complexity (CD). The qualitative inequality in complexity is based on the following hierarchy: CP > CM > CD. This law-like relationship holds true for all complex traits, including complex diseases. We present a hypothesis of two types of variation, namely open and closed (hidden) systems, show that hidden variation provides a hitherto undiscovered "third source" of phenotypic variation, beside genotype and environment, and argue that "missing heritability" for some complex diseases is likely to be a case of "diluted heritability". There is a need for radically new ways of thinking about the principles of genotype-phenotype relationship. Understanding how cells use hidden, pathway variation to respond to stress can shed light on why two individuals who share the same risk factors may not develop the same disease, or how cancer cells escape death.