Powerful detection of polygenic selection and evidence of environmental adaptation in US beef cattle

Selection on complex traits can rapidly drive evolution, especially in stressful environments. This polygenic selection does not leave intense sweep signatures on the genome, rather many loci experience small allele frequency shifts, resulting in large cumulative phenotypic changes. Directional selection and local adaptation are changing populations; but, identifying loci underlying polygenic or environmental selection has been difficult. We use genomic data on tens of thousands of cattle from three populations, distributed over time and landscapes, in linear mixed models with novel dependent variables to map signatures of selection on complex traits and local adaptation. We identify 207 genomic loci associated with an animal’s birth date, representing ongoing selection for monogenic and polygenic traits. Additionally, hundreds of additional loci are associated with continuous and discrete environments, providing evidence for historical local adaptation. These candidate loci highlight the nervous system’s possible role in local adaptation. While advanced technologies have increased the rate of directional selection in cattle, it has likely been at the expense of historically generated local adaptation, which is especially problematic in changing climates. When applied to large, diverse cattle datasets, these selection mapping methods provide an insight into how selection on complex traits continually shapes the genome. Further, understanding the genomic loci implicated in adaptation may help us breed more adapted and efficient cattle, and begin to understand the basis for mammalian adaptation, especially in changing climates. These selection mapping approaches help clarify selective forces and loci in evolutionary, model, and agricultural contexts.     

Molecular markers to study genetic drift and selection in wheat populations

Studying the heterogeneity in variation of gene frequency among populations or between generations may be a possible way to detect genomic regions experiencing selection. In order to evaluate this approach, RFLP markers were used to compare the allelic frequencies in wheat populations that had been submitted to natural selection. In 1984, samples of two composite cross populations were distributed in the French network for dynamic management of genetic resources. Since then, all the sub-populations have been cultivated in the same sites with no human selection. The strong differentiation between populations found for agro-morphological traits (earliness, resistance to pathogens, ...) provided evidence of their adaptation to local conditions. The two initial populations and six derived sub-populations cultivated for 10 years in four contrasted sites were studied with RFLP markers. Differentiation between sub-populations based on RFLP diversity was highly significant. Variations on allelic frequencies of the 30 loci scored were found to be much greater than expected under genetic drift only. This led us to conclude that selection greatly influenced the evolution of the populations. Some of the loci clearly presented a higher differentiation than the others. This might indicate that they were genetically linked to other loci polymorphic in the populations and involved in adaptation. However, the effect of one selected gene on a marker, even located very close to the gene, could not be predicted with certainty. Hence, though the populations were predominantly selfing, it seems that initial linkage disequilibriums between markers and selected genes were not strong enough to control closely the evolution of allelic frequencies at the markers.     

Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird

Investigating the extent (or the existence) of local adaptation is crucial to understanding how populations adapt. When experiments or fitness measurements are difficult or impossible to perform in natural populations, genomic techniques allow us to investigate local adaptation through the comparison of allele frequencies and outlier loci along environmental clines. The thick‐billed murre (Uria lomvia) is a highly philopatric colonial arctic seabird that occupies a significant environmental gradient, shows marked phenotypic differences among colonies, and has large effective population sizes. To test whether thick‐billed murres from five colonies along the eastern Canadian Arctic coast show genomic signatures of local adaptation to their breeding grounds, we analyzed geographic variation in genome‐wide markers mapped to a newly assembled thick‐billed murre reference genome. We used outlier analyses to detect loci putatively under selection, and clustering analyses to investigate patterns of differentiation based on 2220 genomewide single nucleotide polymorphisms (SNPs) and 137 outlier SNPs. We found no evidence of population structure among colonies using all loci but found population structure based on outliers only, where birds from the two northernmost colonies (Minarets and Prince Leopold) grouped with birds from the southernmost colony (Gannet), and birds from Coats and Akpatok were distinct from all other colonies. Although results from our analyses did not support local adaptation along the latitudinal cline of breeding colonies, outlier loci grouped birds from different colonies according to their non‐breeding distributions, suggesting that outliers may be informative about adaptation and/or demographic connectivity associated with their migration patterns or nonbreeding grounds.     

MHC-mediated local adaptation in reciprocally translocated Chinook salmon

Most Pacific salmonid populations have faced significant population declines over the past 30 years. In order to effectively conserve and manage these populations, knowledge of the evolutionary adaptive state of individuals and the scale of adaptation across populations is needed. The vertebrate major histocompatibility complex (MHC) represents an important adaptation to parasites, and genes encoding for the MHC are widely held to be undergoing balancing selection. However, the generality of balancing selection across populations at MHC loci is not well documented. Using Chinook salmon (Oncorhynchus tshawytscha) from two populations, we follow the survival of full-sib family replicates reared in their natal river and reciprocally transplanted to a foreign river to examine selection and local adaptation at the MHC class I and II loci. In both populations, we found evidence of a survivorship advantage associated with nucleotide diversity at the MHC class I locus. In contrast, we found evidence that MHC class II diversity was disadvantageous in one population. There was no evidence that these effects occurred in translocated families, suggesting some degree of local adaptation at the MHC loci. Thus, our results implicate balancing selection at the MHC class I but potentially differing selection across populations at the class II locus.     

Comparing methods for detecting multilocus adaptation with multivariate genotype–environment associations

Identifying adaptive loci can provide insight into the mechanisms underlying local adaptation. Genotype–environment association (GEA) methods, which identify these loci based on correlations between genetic and environmental data, are particularly promising. Univariate methods have dominated GEA, despite the high dimensional nature of genotype and environment. Multivariate methods, which analyse many loci simultaneously, may be better suited to these data as they consider how sets of markers covary in response to environment. These methods may also be more effective at detecting adaptive processes that result in weak, multilocus signatures. Here, we evaluate four multivariate methods and five univariate and differentiation‐based approaches, using published simulations of multilocus selection. We found that Random Forest performed poorly for GEA. Univariate GEAs performed better, but had low detection rates for loci under weak selection. Constrained ordinations, particularly redundancy analysis (RDA), showed a superior combination of low false‐positive and high true‐positive rates across all levels of selection. These results were robust across the demographic histories, sampling designs, sample sizes and weak population structure tested here. The value of combining detections from different methods was variable and depended on the study goals and knowledge of the drivers of selection. Re‐analysis of genomic data from grey wolves highlighted the unique, covarying sets of adaptive loci that could be identified using RDA. Although additional testing is needed, this study indicates that RDA is an effective means of detecting adaptation, including signatures of weak, multilocus selection, providing a powerful tool for investigating the genetic basis of local adaptation.     

Locus-specific genetic differentiation at Rw among warfarin-resistant rat (Rattus norvegicus) populations

Populations may diverge at fitness-related genes as a result of adaptation to local conditions. The ability to detect this divergence by marker-based genomic scans depends on the relative magnitudes of selection, recombination, and migration. We survey rat (Rattus norvegicus) populations to assess the effect that local selection with anticoagulant rodenticides has had on microsatellite marker variation and differentiation at the warfarin resistance gene (Rw) relative to the effect on the genomic background. Initially, using a small sample of 16 rats, we demonstrate tight linkage of microsatellite D1Rat219 to Rw by association mapping of genotypes expressing an anticoagulant-rodenticide-insensitive vitamin K 2,3-epoxide reductase (VKOR). Then, using allele frequencies at D1Rat219, we show that predicted and observed resistance levels in 27 populations correspond, suggesting intense and recent selection for resistance. A contrast of F(ST) values between D1Rat219 and the genomic background revealed that rodenticide selection has overwhelmed drift-mediated population structure only at Rw. A case-controlled design distinguished these locus-specific effects of selection at Rw from background levels of differentiation more effectively than a population-controlled approach. Our results support the notion that an analysis of locus-specific population genetic structure may assist the discovery and mapping of novel candidate loci that are the object of selection or may provide supporting evidence for previously identified loci.     

A tight balance between natural selection and gene flow in a southern African arid-zone endemic bird

Gene flow is traditionally thought to be antagonistic to population differentiation and local adaptation. However, recent studies have demonstrated that local adaptation can proceed provided that selection is greater than the homogenizing effects of gene flow. We extend these initial studies by combining ecology (climate), phenotype (body size), physiological genetics (oxidative phosphorylation genes), and neutral loci (nuclear microsatellites and introns) to test whether selection can counter-balance gene flow and hence promote local adaptation in a bird whose distribution spans an aridity gradient. Our results show that the Karoo scrub-robin's climatic niche is spatially structured, providing the potential for local adaptation to develop. We found remarkably discordant patterns of divergence among mtDNA, morphology, and neutral loci. For the mitochondrial genes, two amino acid replacements, strong population structure and reduced gene flow were associated with the environmental gradient separating western coastal sites from the interior of southern Africa. In contrast, morphology and the neutral loci exhibited variation independent of environmental variables, and revealed extensive levels of gene flow across the aridity gradient, 50 times larger than the estimates for mitochondrial genes. Together, our results suggest that selective pressures on physiology, mediated by the mitochondrial genome, may well be a common mechanism for facilitating local adaptation to new climatic conditions.     

Population Genetic Studies Revealed Local Adaptation in a High Gene-Flow Marine Fish,the Small Yellow Croaker (Larimichthys polyactis)

The genetic differentiation of many marine fish species is low. Yet local adaptation may be common in marine fish species as the vast and changing marine environment provides more chances for natural selection. Here, we used anonymous as well as known protein gene linked microsatellites and mitochondrial DNA to detect the population structure of the small yellow croaker (Larimichthys polyactis) in the Northwest Pacific marginal seas. Among these loci, we detected at least two microsatellites, anonymous H16 and HSP27 to be clearly under diversifying selection in outlier tests. Sequence cloning and analysis revealed that H16 was located in the intron of BAHCC1 gene. Landscape genetic analysis showed that H16 mutations were significantly associated with temperature, which further supported the diversifying selection at this locus. These marker types presented different patterns of population structure: (i) mitochondrial DNA phylogeny showed no evidence of genetic divergence and demonstrated only one glacial linage; (ii) population differentiation using putatively neutral microsatellites presented a pattern of high gene flow in the L. polyactis. In addition, several genetic barriers were identified; (iii) the population differentiation pattern revealed by loci under diversifying selection was rather different from that revealed by putatively neutral loci. The results above suggest local adaptation in the small yellow croaker. In summary, population genetic studies based on different marker types disentangle the effects of demographic history, migration, genetic drift and local adaptation on population structure and also provide valuable new insights for the design of management strategies in L. polyactis.     

Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria)

Today, with the rapid development of population genomics, the genetic basis of adaptation can be unraveled directly at the genome level, without any prerequisites about the selectively advantageous genes or traits. For nonmodel species, it is now possible to screen many markers randomly scattered across the genome and to distinguish between the neutral genetic background and outlier loci displaying an atypical behavior (e.g., a higher differentiation between populations). This study investigated the genetic frame of adaptation to a gradient of altitude in the common frog (Rana temporaria) by means of a genome scan based on 392 amplified fragment length polymorphism markers. Using two outlier detection methods never applied to dominant data so far, we sought for loci with a genetic differentiation diverging from neutral expectations when comparing populations from different altitudes. All the detected loci were sorted out according to their most probable cause for outlier behavior and classified as false positives, outliers due to local effects, or outliers associated with altitude. Altogether, eight good candidate loci were identified as potentially involved in adaptation to altitude because they were picked out in several independent interaltitude comparisons. This result illustrated the potential of genome-wide surveys to reveal selection signatures along selection gradients, where the association between environmental variables and fitness-related traits may be complex and/or cryptic. In this article, we also underlined the need for confirmation of the selection footprints for the outlier loci. Finally, we provided some preliminary insights into the genetic basis of adaptation along an altitudinal cline in the common frog.     

Genomic and chemical evidence for local adaptation in resistance to different herbivores in Datura stramonium

Because most species are collections of genetically variable populations distributed to habitats differing in their abiotic/biotic environmental factors and community composition, the pattern and strength of natural selection imposed by species on each other's traits are also expected to be highly spatially variable. Here, we used genomic and quantitative genetic approaches to understand how spatially variable selection operates on the genetic basis of plant defenses to herbivores. To this end, an F₂ progeny was generated by crossing Datura stramonium (Solanaceae) parents from two populations differing in their level of chemical defense. This F₂ progeny was reciprocally transplanted into the parental plants’ habitats and by measuring the identity by descent (IBD) relationship of each F₂ plant to each parent, we were able to elucidate how spatially variable selection imposed by herbivores operated on the genetic background (IBD) of resistance to herbivory, promoting local adaptation. The results highlight that plants possessing the highest total alkaloid concentrations (sum of all alkaloid classes) were not the most well-defended or fit. Instead, specific alkaloids and their linked loci/alleles were favored by selection imposed by different herbivores. This has led to population differentiation in plant defenses and thus, to local adaptation driven by plant-herbivore interactions.     

Evolution of migration under kin selection and local adaptation

We present here a stochastic two-locus, two-habitat model for the evolution of migration with local adaptation and kin selection. One locus determines the migration rate while the other causes local adaptation. We show that the opposing forces of kin competition and local adaptation can lead to the existence of one or two convergence stable migration rates, notably depending on the recombination rate between the two loci. We show that linkage between migration and local adaptation loci has two antagonist effects: when linkage is tight, cost of local adaptation increases, leading to smaller equilibrium migration rates. However, when linkage is tighter, the population structure at the migration locus tends to be very high because of the indirect selection, and thus equilibrium migration rates increases. This result, qualitatively different from results obtained with other models of migration evolution, indicates that ignoring drift or the detail of the genetic architecture may lead to incorrect conclusions.     

Structured Coalescent Processes on Different Time Scales

It is demonstrated that the structured coalescent model can readily be extended to include phenomena such as partial selfing and background selection through the use of an approximation based on separation of time scales. A model that includes these phenomena, as well as geographic subdivision and linkage to a polymorphism maintained either by local adaptation or by balancing selection, is derived, and the expected coalescence time for a pair of genes is calculated. It is found that background selection reduces coalescence times within subpopulations and allelic classes, leading to a high degree of apparent differentiation. Extremely high levels of subpopulation differentiation are also expected for regions of the genome surrounding loci important in local adaptation. These regions will be wider the stronger the local selection, and the higher the selfing rate.     

A genome scan to detect candidate regions influenced by local natural selection in human populations

As human populations dispersed throughout the world, they were subjected to new selective forces, which must have led to local adaptation via natural selection and hence altered patterns of genetic variation. Yet, there are very few examples known in which such local selection has clearly influenced human genetic variation. A potential approach for detecting local selection is to screen random loci across the genome; those loci that exhibit unusually large genetic distances between human populations are then potential markers of genomic regions under local selection. We investigated this approach by genotyping 332 short tandem repeat (STR) loci in Africans and Europeans and calculating the genetic differentiation for each locus. Patterns of genetic diversity at these loci were consistent with greater variation in Africa and with local selection operating on populations as they moved out of Africa. For 11 loci exhibiting the largest genetic differences, we genotyped an additional STR locus located nearby; the genetic distances for these nearby loci were significantly larger than average. These genomic regions therefore reproducibly exhibit larger genetic distances between populations than the "average" genomic region, consistent with local selection. Our results demonstrate that genome scans are a promising means of identifying candidate regions that have been subjected to local selection.     

The coupling hypothesis: why genome scans may fail to map local adaptation genes

Genomic scans of multiple populations often reveal marker loci with greatly increased differentiation between populations. Often this differentiation coincides in space with contrasts in ecological factors, forming a genetic-environment association (GEA). GEAs imply a role for local adaptation, and so it is tempting to conclude that the strongly differentiated markers are themselves under ecologically based divergent selection, or are closely linked to loci under such selection. Here, we highlight an alternative and neglected explanation: intrinsic (i.e. environment-independent) pre- or post-zygotic genetic incompatibilities rather than local adaptation can be responsible for increased differentiation. Intrinsic genetic incompatibilities create endogenous barriers to gene flow, also known as tension zones, whose location can shift over time. However, tension zones have a tendency to become trapped by, and therefore to coincide with, exogenous barriers due to ecological selection. This coupling of endogenous and exogenous barriers can occur easily in spatially subdivided populations, even if the loci involved are unlinked. The result is that local adaptation explains where genetic breaks are positioned, but not necessarily their existence, which can be best explained by endogenous incompatibilities. More precisely, we show that (i) the coupling of endogenous and exogenous barriers can easily occur even when ecological selection is weak; (ii) when environmental heterogeneity is fine-grained, GEAs can emerge at incompatibility loci, but only locally, in places where habitats and gene pools are sufficiently intermingled to maintain linkage disequilibria between genetic incompatibilities, local-adaptation genes and neutral loci. Furthermore, the association between the locally adapted and intrinsically incompatible alleles (i.e. the sign of linkage disequilibrium between endogenous and exogenous loci) is arbitrary and can form in either direction. Reviewing results from the literature, we find that many predictions of our model are supported, including endogenous genetic barriers that coincide with environmental boundaries, local GEA in mosaic hybrid zones, and inverted or modified GEAs at distant locations. We argue that endogenous genetic barriers are often more likely than local adaptation to explain the majority of Fst-outlying loci observed in genome scan approaches - even when these are correlated to environmental variables.     

A Population Genetic Signal of Polygenic Adaptation

Adaptation in response to selection on polygenic phenotypes may occur via subtle allele frequencies shifts at many loci. Current population genomic techniques are not well posed to identify such signals. In the past decade, detailed knowledge about the specific loci underlying polygenic traits has begun to emerge from genome-wide association studies (GWAS). Here we combine this knowledge from GWAS with robust population genetic modeling to identify traits that may have been influenced by local adaptation. We exploit the fact that GWAS provide an estimate of the additive effect size of many loci to estimate the mean additive genetic value for a given phenotype across many populations as simple weighted sums of allele frequencies. We use a general model of neutral genetic value drift for an arbitrary number of populations with an arbitrary relatedness structure. Based on this model, we develop methods for detecting unusually strong correlations between genetic values and specific environmental variables, as well as a generalization of comparisons to test for over-dispersion of genetic values among populations. Finally we lay out a framework to identify the individual populations or groups of populations that contribute to the signal of overdispersion. These tests have considerably greater power than their single locus equivalents due to the fact that they look for positive covariance between like effect alleles, and also significantly outperform methods that do not account for population structure. We apply our tests to the Human Genome Diversity Panel (HGDP) dataset using GWAS data for height, skin pigmentation, type 2 diabetes, body mass index, and two inflammatory bowel disease datasets. This analysis uncovers a number of putative signals of local adaptation, and we discuss the biological interpretation and caveats of these results.     

Fitness costs of insecticide resistance in natural breeding sites of the mosquito Culex pipiens

Genetic changes conferring adaptation to a new environment may induce a fitness cost in the previous environment. Although this prediction has been verified in laboratory conditions, few studies have tried to document this cost directly in natural populations. Here, we evaluated the pleiotropic effects of insecticide resistance on putative fitness components of the mosquito Culex pipiens. Experiments using different larval densities were performed during the summer in two natural breeding sites. Two loci that possess alleles conferring organophosphate (OP) resistance were considered: ace-1 coding for an acetylcholinesterase (AChE1, the OP target) and Ester, a 'super locus" including two closely linked loci coding for esterases A and B. Resistance ace-1 alleles coding for a modified AChE1 were associated with a longer development time and shorter wing length. The pleiotropic effects of two resistance alleles Ester1 and Ester4 coding for the overproduced esterases A1 and A4-B4, respectively, were more variable. Both A1 and A4-B4 reduced wing length, although only A1 was associated with a longer preimaginal stage. The fluctuating asymmetry (FA) of the wing did not respond to the presence or to the interaction of resistance alleles at the two loci at any of the density levels tested. Conversely, the FA of one wing section decreased when larval density increased. This may be the consequence of selection against less developmentally stable individuals. The results are discussed in relation to the local evolution of insecticide resistance genes.     

Signatures of selection and sex-specific expression variation of a novel duplicate during the evolution of the Drosophila desaturase gene family

The tempo and mode of evolution of loci with a large effect on adaptation and reproductive isolation will influence the rate of evolutionary divergence and speciation. Desaturase loci are involved in key biochemical changes in long-chain fatty acids. In insects, these have been shown to influence adaptation to starvation or desiccation resistance and in some cases act as important pheromones. The desaturase gene family of Drosophila is known to have evolved by gene duplication and diversification, and at least one locus shows rapid evolution of sex-specific expression variation. Here, we examine the evolution of the gene family in species representing the Drosophila phylogeny. We find that the family includes more loci than have been previously described. Most are represented as single-copy loci, but we also find additional examples of duplications in loci which influence pheromone blends. Most loci show patterns of variation associated with purifying selection, but there are strong signatures of diversifying selection in new duplicates. In the case of a new duplicate of desat1 in the obscura group species, we show that strong selection on the coding sequence is associated with the evolution of sex-specific expression variation. It seems likely that both sexual selection and ecological adaptation have influenced the evolution of this gene family in Drosophila.     

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of Q(ST) and F(ST) at two spatial scales

Genetic differentiation in 20 hierarchically sampled populations of wild barley was analyzed with quantitative traits, allozymes and Random Amplified Polymorphic DNAs (RAPDs), and compared for three marker types at two hierarchical levels. Regional subdivision for both molecular markers was much lower than for quantitative traits. For both allozymes and RAPDs, most loci exhibited minor or no regional differentiation, and the relatively high overall estimates of the latter were due to several loci with exceptionally high regional differentiation. The allozyme- and RAPD-specific patterns of differentiation were concordant in general with one another, but not with quantitative trait differentiation. Divergent selection on quantitative traits inferred from very high regional Q(ST) was in full agreement with our previous results obtained from a test of local adaptation and multilevel selection analysis. In contrast, most variation in allozyme and RAPD variation was neutral, although several allozyme loci and RAPD markers were exceptional in their levels of regional differentiation. However, it is not possible to answer the question whether these exceptional loci are directly involved in the response to selection pressure or merely linked to the selected loci. The fact that Q(ST) and F(ST) did not differ at the population scale, that is, within regions, but differed at the regional scale, for which local adaptation has been previously shown, implies that comparison of the level of subdivision in quantitative traits, as compared with molecular markers, is indicative of adaptive population differentiation only when sampling is carried out at the appropriate scale.     

Habitat-specific natural selection at a flowering-time QTL is a main driver of local adaptation in two wild barley populations

Understanding the genetic basis of local adaptation requires insight in the fitness effects of individual loci under natural field conditions. While rapid progress is made in the search for genes that control differences between plant populations, it is typically unknown whether the genes under study are in fact key targets of habitat-specific natural selection. Using a quantitative trait loci (QTL) approach, we show that a QTL associated with flowering-time variation between two locally adapted wild barley populations is an important determinant of fitness in one, but not in the other population's native habitat. The QTL mapped to the same position as a habitat-specific QTL for field fitness that affected plant reproductive output in only one of the parental habitats, indicating that the genomic region is under differential selection between the native habitats. Consistent with the QTL results, phenotypic selection of flowering time differed between the two environments, whereas other traits (growth rate and seed weight) were under selection but experienced no habitat-specific differential selection. This implies the flowering-time QTL as a driver of adaptive population divergence. Our results from phenotypic selection and QTL analysis are consistent with local adaptation without genetic trade-offs in performance across environments, i.e. without alleles or traits having opposing fitness effects in contrasting environments.     

Genetic explorations of recent human metabolic adaptations: hypotheses and evidence

Since humans and chimpanzees split from a common ancestor over 6 million years ago, human metabolism has changed dramatically. This change includes adaptations to a high-quality diet, the evolution of an energetically expensive brain, dramatic increases in endurance abilities, and capacity for energy storage in white adipose tissue. Human metabolism continues to evolve in modern human populations in response to local environmental and cultural selective forces. Understanding the nature of these selective forces and the physiological responses during human evolution is a compelling challenge for evolutionary biologists. The complex genetic architecture surrounding metabolic phenotypes indicates that selection probably altered allelic frequencies across many loci in populations experiencing adaptive metabolic change to fit their environment. A recent analysis supports this hypothesis, finding that classic selective sweeps at single loci were rare during the past 250 000 years of human evolution. Detection of selective signatures at multiple loci, as well as exploration of physiological adaptation to environment in humans, will require cross-disciplinary collaboration, including the incorporation of biological pathway analysis. This review explores the Thrifty Genotype Hypothesis, high-altitude adaptation, cold-resistance adaptation, and genetic evidence surrounding these proposed metabolic adaptations in an attempt to clarify current challenges and avenues for future progress.