Inferring selection in instances of long‐range colonization: The Aleppo pine (Pinus halepensis) in the Mediterranean Basin

Teasing apart the effects of natural selection and demography on current allele frequencies is challenging, due to both processes leaving a similar molecular footprint. In particular, when attempting to identify selection in species that have undergone a recent range expansion, the increase in genetic drift at the edges of range expansions (“allele surfing”) can be a confounding factor. To address this potential issue, we first assess the long‐range colonization history of the Aleppo pine across the Mediterranean Basin, using molecular markers. We then look for single nucleotide polymorphisms (SNPs) involved in local adaptation using: (a) environmental correlation methods (bayenv2 ), focusing on bioclimatic variables important for the species’ adaptation (i.e., temperature, precipitation and water availability); and (b) F_ST‐related methods (pcadapt ). To assess the rate of false positives caused by the allele surfing effect, these results are compared with results from simulated SNP data that mimics the species’ past range expansions and the effect of genetic drift, but with no selection. We find that the Aleppo pine shows a previously unsuspected complex genetic structure across its range, as well as evidence of selection acting on SNPs involved with the response to bioclimatic variables such as drought. This study uses an original approach to disentangle the confounding effects of drift and selection in range margin populations. It also contributes to the increased evidence that plant populations are able to adapt to new environments despite the expected accumulation of deleterious mutations that takes place during long‐range colonizations.     

Evolutionarily stable strategies and short-term selection in Mendelian populations re-visited

This note concerns a one locus, two allele, random mating diploid population, subject to frequency-dependent viability selection. It is already known that in such a population, any evolutionarily stable strategies (ESS), if only accessible by the genotype-to-phenotype mapping, is the phenotypic image of a stable genetic equilibrium (Eshel, I. 1982. Evolutionarily stable strategies and viability selection in Mendelian populations. Theor. Popul. Biol. 22(2), 204-217; Cressman et al. 1996. Evolutionary stability in strategic models of single-locus frequency-dependent viability selection. J. Math. Biol. 34, 707-733). The opposite is not true. We find necessary and sufficient parametric conditions for global convergence to the ESS, but we also demonstrate conditions under which, although a unique, genetically accessible ESS exists, there is another, "non-phenotypic" genetically stable equilibrium.     

MIGRATION AND GENETIC DRIFT IN HUMAN POPULATIONS

In humans and many other species, mortality is concentrated early in the life cycle, and is low during the ages of dispersal and reproduction. Yet precisely the opposite is assumed by classical population-genetics models of migration and genetic drift. We introduce a model in which population regulation occurs before migration. In contrast to the conventional model, our model implies that geographic variation in the allele frequencies of newborns should exceed that of adults. Thus, it is important to distinguish genetic variation of adults from that of newborns in species with human-like life cycles. Classical models deal with the variance of group allele frequencies about the allele frequency of a hypothetical “continent” or “foundation stock.” Empirical studies, however, can only measure “reduced” variance, i.e., variance about the current population mean. Our model deals with reduced variance, and should therefore be more relevant to field studies. We show that reduced variance converges faster, which implies that populations are more likely to be at equilibrium with respect to reduced than unreduced variance. To summarize the effect of migration on genetic population structure, we introduce a new parameter, the effective migration rate. Unlike most population structure statistics, it does not confound the effects of mobility and population size, and it should therefore be useful for comparisons between populations. Finally, we show that the difference between geographic variation of newborn and adult allele frequencies contains information about both effective population size and effective migration rate.     

Population Analysis of Large Copy Number Variants and Hotspots of Human Genetic Disease

Copy number variants (CNVs) contribute to human genetic and phenotypic diversity. However, the distribution of larger CNVs in the general population remains largely unexplored. We identify large variants in ~2500 individuals by using Illumina SNP data, with an emphasis on “hotspots” prone to recurrent mutations. We find variants larger than 500 kb in 5%–10% of individuals and variants greater than 1 Mb in 1%–2%. In contrast to previous studies, we find limited evidence for stratification of CNVs in geographically distinct human populations. Importantly, our sample size permits a robust distinction between truly rare and polymorphic but low-frequency copy number variation. We find that a significant fraction of individual CNVs larger than 100 kb are rare and that both gene density and size are strongly anticorrelated with allele frequency. Thus, although large CNVs commonly exist in normal individuals, which suggests that size alone can not be used as a predictor of pathogenicity, such variation is generally deleterious. Considering these observations, we combine our data with published CNVs from more than 12,000 individuals contrasting control and neurological disease collections. This analysis identifies known disease loci and highlights additional CNVs (e.g., 3q29, 16p12, and 15q25.2) for further investigation. This study provides one of the first analyses of large, rare (0.1%–1%) CNVs in the general population, with insights relevant to future analyses of genetic disease.     

COMPETITION AND GENOTYPE-BY-ENVIRONMENT INTERACTION IN NATURAL BREEDING SUBSTRATES OF DROSOPHILA

Although empirical studies frequently suggest that genotype-by-environment (G X E) interaction can maintain genetic variation, very few data are available to test for the specific conditions necessary for the existence of a protected polymorphism (i.e., the property of persistence of an allele even when initially rare). Drosophila species live in patchy environments and their local population structure may be characterized to some extent by Levene's migration pattern, namely by a single pool of individuals that presumably mate at random and breed on discrete and ephemeral resources. We present here a field experiment that links Drosophila ecology and population genetics, which used the alcohol dehydrogenase (Adh) and α-glycerophosphate dehydrogenase (αGpdh) polymorphic loci in D. melanogaster flies raised from Opuntia ficus-indica fruits (prickly pears). The results show that there is density-dependent mortality in those fruits with a relatively high number of larvae (i.e., selection is “soft”) and suggest that there is differential viability for αGpdh genotypes. Additionally, a pattern of G X E interaction for fitness values, which is fully compatible with the theoretical conditions required for the existence of a protected polymorphism, was found after weighting the fitness estimates by the relative contribution that each fruit makes to the total adult population. The strong association between Adh^S and αGpdh^F alleles suggests that the occurrence of the common cosmopolitan inversion In(2L)t in the population might be responsible for the negative frequency-dependent selection predicted by Levene's model when genetic variation persists in heterogeneous environments.     

Natural selection affects multiple aspects of genetic variation at putatively neutral sites across the human genome

A major question in evolutionary biology is how natural selection has shaped patterns of genetic variation across the human genome. Previous work has documented a reduction in genetic diversity in regions of the genome with low recombination rates. However, it is unclear whether other summaries of genetic variation, like allele frequencies, are also correlated with recombination rate and whether these correlations can be explained solely by negative selection against deleterious mutations or whether positive selection acting on favorable alleles is also required. Here we attempt to address these questions by analyzing three different genome-wide resequencing datasets from European individuals. We document several significant correlations between different genomic features. In particular, we find that average minor allele frequency and diversity are reduced in regions of low recombination and that human diversity, human-chimp divergence, and average minor allele frequency are reduced near genes. Population genetic simulations show that either positive natural selection acting on favorable mutations or negative natural selection acting against deleterious mutations can explain these correlations. However, models with strong positive selection on nonsynonymous mutations and little negative selection predict a stronger negative correlation between neutral diversity and nonsynonymous divergence than observed in the actual data, supporting the importance of negative, rather than positive, selection throughout the genome. Further, we show that the widespread presence of weakly deleterious alleles, rather than a small number of strongly positively selected mutations, is responsible for the correlation between neutral genetic diversity and recombination rate. This work suggests that natural selection has affected multiple aspects of linked neutral variation throughout the human genome and that positive selection is not required to explain these observations.     

EVOLUTION OF PHENOTYPE–ENVIRONMENT ASSOCIATIONS BY GENETIC RESPONSES TO SELECTION AND PHENOTYPIC PLASTICITY IN A TEMPORALLY AUTOCORRELATED ENVIRONMENT

Covariation between population‐mean phenotypes and environmental variables, sometimes termed a “phenotype–environment association” (PEA), can result from phenotypic plasticity, genetic responses to natural selection, or both. PEAs can potentially provide information on the evolutionary dynamics of a particular set of populations, but this requires a full theoretical characterization of PEAs and their evolution. Here, we derive formulas for the expected PEA in a temporally fluctuating environment for a quantitative trait with a linear reaction norm. We compare several biologically relevant scenarios, including constant versus evolving plasticity, and the situation in which an environment affects both development and selection but at different time periods. We find that PEAs are determined not only by biological factors (e.g., magnitude of plasticity, genetic variation), but also environmental factors, such as the association between the environments of development and of selection, and in some cases the level of temporal autocorrelation. We also describe how a PEA can be used to estimate the relationship between an optimum phenotype and an environmental variable (i.e., the environmental sensitivity of selection), an important parameter for determining the extinction risk of populations experiencing environmental change. We illustrate this ability using published data on the predator‐induced morphological responses of tadpoles to predation risk.     

Cryptic genetic variation can make “irreducible complexity” a common mode of adaptation in sexual populations

The existence of complex (multiple‐step) genetic adaptations that are “irreducible” (i.e., all partial combinations are less fit than the original genotype) is one of the longest standing problems in evolutionary biology. In standard genetics parlance, these adaptations require the crossing of a wide adaptive valley of deleterious intermediate stages. Here, we demonstrate, using a simple model, that evolution can cross wide valleys to produce “irreducibly complex” adaptations by making use of previously cryptic mutations. When revealed by an evolutionary capacitor, previously cryptic mutants have higher initial frequencies than do new mutations, bringing them closer to a valley‐crossing saddle in allele frequency space. Moreover, simple combinatorics implies an enormous number of candidate combinations exist within available cryptic genetic variation. We model the dynamics of crossing of a wide adaptive valley after a capacitance event using both numerical simulations and analytical approximations. Although individual valley crossing events become less likely as valleys widen, by taking the combinatorics of genotype space into account, we see that revealing cryptic variation can cause the frequent evolution of complex adaptations.     

Random copying,frequency-dependent copying and culture change

Previous evolutionary analyses of human culture have found that a simple model of random copying, analogous to neutral genetic drift, can generate the distinct power-law frequency distribution of cultural traits that is typical of various real-world cultural datasets, such as first names, patent citations and prehistoric pottery types. Here, we use agent-based simulations to explore the effects of frequency-dependent copying (e.g., conformity and anti-conformity) on this power-law distribution. We find that when traits are actively selected on the basis of their frequency, then the power-law distribution is severely disrupted. Conformity generates a “winner-takes-all” distribution in which popular traits dominate, while anti-conformity generates a “humped” distribution in which traits of intermediate frequency are favoured. However, a more passive frequency-dependent “trimming”, in which traits are selectively ignored on the basis of their frequency, generates reasonable approximations to the power-law distribution. This frequency-dependent trimming may therefore be difficult to distinguish from genuine random copying using population-level data alone. Implications for the study of both human and nonhuman culture are discussed.     

Roles of maternal effects in maintaining genetic variation: Maternal storage effect

Understanding the maintenance of genetic variation remains a central challenge in evolutionary biology. Recent empirical studies suggest the importance of temporally varying selection, as allele frequencies have been found to fluctuate substantially in the wild. However, previous theory suggests that the conditions for the maintenance of genetic variation under temporally fluctuating selection are quite restrictive. Using mathematical models, we demonstrate that maternal genetic effects, whereby maternal genotypes affect offspring phenotypes, can facilitate the maintenance of polymorphism in temporally varying environments. Maternal effects result in mismatches between genotypes and phenotypes, thereby buffering the influence of selection on allele frequency. This decreases the magnitude of allele‐frequency fluctuations and creates conditions for the maintenance of variation when selection causes fluctuations. Therefore, maternal effects may result in a temporal storage effect (“maternal storage effect”). On the other hand, when selection does not cause fluctuations (e.g., linear negative frequency‐dependent selection), maternal genetic effects moderate the relative importance of selection compared to genetic drift and promote stochastic allele extinction in finite populations. Thus, maternal effects can play an important role in the maintenance of polymorphism, but the direction of the effect depends on the nature of selection.     

THE MAINTENANCE OF GENETIC VARIATION FOR OVIPOSITION RATE IN TWO‐SPOTTED SPIDER MITES: INFERENCES FROM ARTIFICIAL SELECTION

Despite the directional selection acting on life‐history traits, substantial amounts of standing variation for these traits have frequently been found. This variation may result from balancing selection (e.g., through genetic trade‐offs) or from mutation‐selection balance. These mechanisms affect allele frequencies in different ways: Under balancing selection alleles are maintained at intermediate frequencies, whereas under mutation‐selection balance variation is generated by deleterious mutations and removed by directional selection, which leads to asymmetry in the distribution of allele frequencies. To investigate the importance of these two mechanisms in maintaining heritable variation in oviposition rate of the two‐spotted spider mite, we analyzed the response to artificial selection. In three replicate experiments, we selected for higher and lower oviposition rate, compared to control lines. A response to selection only occurred in the downward direction. Selection for lower oviposition rate did not lead to an increase in any other component of fitness, but led to a decline in female juvenile survival. The results suggest standing variation for oviposition rate in this population consists largely of deleterious alleles, as in a mutation‐selection balance. Consequently, the standing variation for this trait does not appear to be indicative of its adaptive potential.     

Deleterious mutations are characterized by higher genomic heterozygosity than other genic variants in plant genomes

Deleterious mutations can reduce the fitness of crop varieties, which limits the plant breeding efficacy. While crop deleterious mutations have been extensively examined, most studies focused on one specific crop with different analyzing methods, which hinders unveiling shared genomic characteristics of deleterious mutations across diverse crops. Here we used standardized approaches to characterize the deleterious mutations in genomes of domesticated inbreeding (i.e., rice, soybean, and tomato) and clonally propagated crops (i.e., grape and pineapple). We found that deleterious mutations are commonly targeted by purifying selection, and are over-presented in a nearly fixed derived allele frequency in the course of plant domestication. Further, a generally negative correlation between genetic load and the artificial selection strength is observed. Importantly, we consistently uncovered the higher derived genomic heterozygosity for deleterious mutations compared to other genic variants. This study broadens our understanding of the evolution of deleterious mutations in plant genomes.     

Frequency and density-dependent selection on life-history strategies--a field experiment

Negative frequency-dependence, which favors rare genotypes, promotes the maintenance of genetic variability and is of interest as a potential explanation for genetic differentiation. Density-dependent selection may also promote cyclic changes in frequencies of genotypes. Here we show evidence for both density-dependent and negative frequency-dependent selection on opposite life-history tactics (low or high reproductive effort, RE) in the bank vole (Myodes glareolus). Density-dependent selection was evident among the females with low RE, which were especially favored in low densities. Instead, both negative frequency-dependent and density-dependent selection were shown in females with high RE, which were most successful when they were rare in high densities. Furthermore, selection at the individual level affected the frequencies of tactics at the population level, so that the frequency of the rare high RE tactic increased significantly at high densities. We hypothesize that these two selection mechanisms (density- and negative frequency-dependent selection) may promote genetic variability in cyclic mammal populations. Nevertheless, it remains to be determined whether the origin of genetic variance in life-history traits is causally related to density variation (e.g. population cycles).     

The evolution of genetic architecture under frequency-dependent disruptive selection

We propose a model to analyze a quantitative trait under frequency-dependent disruptive selection. Selection on the trait is a combination of stabilizing selection and intraspecific competition, where competition is maximal between individuals with equal phenotypes. In addition, there is a density-dependent component induced by population regulation. The trait is determined additively by a number of biallelic loci, which can have different effects on the trait value. In contrast to most previous models, we assume that the allelic effects at the loci can evolve due to epistatic interactions with the genetic background. Using a modifier approach, we derive analytical results under the assumption of weak selection and constant population size, and we investigate the full model by numerical simulations. We find that frequency-dependent disruptive selection favors the evolution of a highly asymmetric genetic architecture, where most of the genetic variation is concentrated on a small number of loci. We show that the evolution of genetic architecture can be understood in terms of the ecological niches created by competition. The phenotypic distribution of a population with an adapted genetic architecture closely matches this niche structure. Thus, evolution of the genetic architecture seems to be a plausible way for populations to adapt to regimes of frequency-dependent disruptive selection. As such, it should be seen as a potential evolutionary pathway to discrete polymorphisms and as a potential alternative to other evolutionary responses, such as the evolution of sexual dimorphism or assortative mating.     

Sperm head morphology is associated with sperm swimming speed: A comparative study of songbirds using electron microscopy

Sperm exhibit extraordinary levels of morphological diversification across the animal kingdom. In songbirds, sperm have a helically shaped head incorporating a distinct acrosomal membrane or “helical keel,” the form and extent of which varies across species. The functional significance of this helical shape, however, remains unknown. Using scanning electron microscopy, we quantified inter‐ and intraspecific variation in sperm head morphology across 36 songbird species (Passeriformes: Passerida). Using phylogenetic comparative methods, we investigated the relationship between sperm head morphology and both sperm swimming speed and the frequency of extra‐pair young (EPY). We found that species whose sperm had a relatively more pronounced helical form (i.e., long acrosome, short nucleus, wide helical membrane, and a more pronounced waveform along the sperm head “core”) had faster‐swimming sperm. We found no evidence of a relationship between interspecific variation in sperm head morphology and EPY, although we did find that among‐ and within‐male variation in sperm head traits were negatively correlated with EPY. Applying principles of fluid mechanics, we discuss how the helical form of the sperm head may influence swimming speed, and suggest that further studies considering aspects of sperm morphology beyond sperm length are needed to improve our understanding of sperm structure‐function relationships.     

A POPULATION GENETIC THEORY OF CANALIZATION

Canalization is the suppression of phenotypic variation. Depending on the causes of phenotypic variation, one speaks either of genetic or environmental canalization. Genetic canalization describes insensitivity of a character to mutations, and the insensitivity to environmental factors is called environmental canalization. Genetic canalization is of interest because it influences the availability of heritable phenotypic variation to natural selection, and is thus potentially important in determining the pattern of phenotypic evolution. In this paper a number of population genetic models are considered of a quantitative character under stabilizing selection. The main purpose of this study is to define the population genetic conditions and constraints for the evolution of canalization. Environmental canalization is modeled as genotype specific environmental variance. It is shown that stabilizing selection favors genes that decrease environmental variance of quantitative characters. However, the theoretical limit of zero environmental variance has never been observed. Of the many ways to explain this fact, two are addressed by our model. It is shown that a “canalization limit” is reached if canalizing effects of mutations are correlated with direct effects on the same character. This canalization limit is predicted to be independent of the strength of stabilizing selection, which is inconsistent with recent experimental data (Sterns et al. 1995). The second model assumes that the canalizing genes have deleterious pleiotropic effects. If these deleterious effects are of the same magnitude as all the other mutations affecting fitness very strong stabilizing selection is required to allow the evolution of environmental canalization. Genetic canalization is modeled as an influence on the average effect of mutations at a locus of other genes. It is found that the selection for genetic canalization critically depends on the amount of genetic variation present in the population. The more genetic variation, the stronger the selection for canalizing effects. All factors that increase genetic variation favor the evolution of genetic canalization (large population size, high mutation rate, large number of genes). If genetic variation is maintained by mutation-selection balance, strong stabilizing selection can inhibit the evolution of genetic canalization. Strong stabilizing selection eliminates genetic variation to a level where selection for canalization does not work anymore. It is predicted that the most important characters (in terms of fitness) are not necessarily the most canalized ones, if they are under very strong stabilizing selection (k > 0.2V_e). The rate of decrease of mutational variance V_m is found to be less than 10% of the initial V_m. From this result it is concluded that characters with typical mutational variances of about 10^–3 V_e are in a metastable state where further evolution of genetic canalization is too slow to be of importance at a microevolutionary time scale. The implications for the explanation of macroevolutionary patterns are discussed.     

Phenotypic evolution in stochastic environments: The contribution of frequency- and density-dependent selection

Understanding how environmental variation affects phenotypic evolution requires models based on ecologically realistic assumptions that include variation in population size and specific mechanisms by which environmental fluctuations affect selection. Here we generalize quantitative genetic theory for environmentally induced stochastic selection to include general forms of frequency- and density-dependent selection. We show how the relevant fitness measure under stochastic selection relates to Fisher's fundamental theorem of natural selection, and present a general class of models in which density regulation acts through total use of resources rather than just population size. In this model, there is a constant adaptive topography for expected evolution, and the function maximized in the long run is the expected factor restricting population growth. This allows us to generalize several previous results and to explain why apparently “-selected” species with slow life histories often have low carrying capacities. Our joint analysis of density- and frequency-dependent selection reveals more clearly the relationship between population dynamics and phenotypic evolution, enabling a broader range of eco-evolutionary analyses of some of the most interesting problems in evolution in the face of environmental variation.     

A genome-wide comparison of the functional properties of rare and common genetic variants in humans

One of the longest running debates in evolutionary biology concerns the kind of genetic variation that is primarily responsible for phenotypic variation in species. Here, we address this question for humans specifically from the perspective of population allele frequency of variants across the complete genome, including both coding and noncoding regions. We establish simple criteria to assess the likelihood that variants are functional based on their genomic locations and then use whole-genome sequence data from 29 subjects of European origin to assess the relationship between the functional properties of variants and their population allele frequencies. We find that for all criteria used to assess the likelihood that a variant is functional, the rarer variants are significantly more likely to be functional than the more common variants. Strikingly, these patterns disappear when we focus on only those variants in which the major alleles are derived. These analyses indicate that the majority of the genetic variation in terms of phenotypic consequence may result from a mutation-selection balance, as opposed to balancing selection, and have direct relevance to the study of human disease.     

The maintenance of phenotypic and genetic variation in threshold traits by frequency-dependent selection

Many traits are phenotypically dimorphic but determined by the action of many loci, the phenotype being a result of a threshold of sensitivity. Quantitative genetic analysis has shown that generally there is considerable additive genetic variation for the trait, the average heritability being 0.52. In numerous cases threshold traits have been shown, or are assumed, to be under frequency-dependent selection; examples include satellite-territorial behaviour, sex-determination, wing dimorphism and trophic dimorphism. In this paper I investigate the potential for frequency-dependent selection to maintain both phenotypic and additive genetic variation in threshold traits. The qualitative results are robust to the particular form of the frequency-dependent selection function. The equilibrium proportion is more or less independent of population size but the heritability increases with population size, typically approaching its maximal value at a population size of 5000, when the mutation rate is 10^?4. A tenfold decrease in the mutation rate requires an approximate doubling of the population size before an asymptotic value is approached. Thus frequency-dependent selection can account for both the existence of two morphs in a population and the observed levels of heritability. It is also shown, both via simulation and theory, that the quantitative genetic model and a simple phenotypic analysis predict the same equilibrium morph proportion.