Population admixture,biological invasions and the balance between local adaptation and inbreeding depression

When previously isolated populations meet and mix, the resulting admixed population can benefit from several genetic advantages, including increased genetic variation, the creation of novel genotypes and the masking of deleterious mutations. These admixture benefits are thought to play an important role in biological invasions. In contrast, populations in their native range often remain differentiated and frequently suffer from inbreeding depression owing to isolation. While the advantages of admixture are evident for introduced populations that experienced recent bottlenecks or that face novel selection pressures, it is less obvious why native range populations do not similarly benefit from admixture. Here we argue that a temporary loss of local adaptation in recent invaders fundamentally alters the fitness consequences of admixture. In native populations, selection against dilution of the locally adapted gene pool inhibits unconstrained admixture and reinforces population isolation, with some level of inbreeding depression as an expected consequence. We show that admixture is selected against despite significant inbreeding depression because the benefits of local adaptation are greater than the cost of inbreeding. In contrast, introduced populations that have not yet established a pattern of local adaptation can freely reap the benefits of admixture. There can be strong selection for admixture because it instantly lifts the inbreeding depression that had built up in isolated parental populations. Recent work in Silene suggests that reduced inbreeding depression associated with post-introduction admixture may contribute to enhanced fitness of invasive populations. We hypothesize that in locally adapted populations, the benefits of local adaptation are balanced against an inbreeding cost that could develop in part owing to the isolating effect of local adaptation itself. The inbreeding cost can be revealed in admixing populations during recent invasions.     

Adaptive Role of Inversion Polymorphism of Drosophila subobscura in Lead Stressed Environment

Local adaptation to environmental stress at different levels of genetic polymorphism in various plants and animals has been documented through evolution of heavy metal tolerance. We used samples of Drosophila subobscura populations from two differently polluted environments to analyze the change of chromosomal inversion polymorphism as genetic marker during laboratory exposure to lead. Exposure to environmental contamination can affect the genetic content within a particular inversion and produce targets for selection in populations from different environments. The aims were to discover whether the inversion polymorphism is shaped by the local natural environments, and if lead as a selection pressure would cause adaptive divergence of two populations during the multigenerational laboratory experiment. The results showed that populations retain signatures from past contamination events, and that heavy metal pollution can cause adaptive changes in population. Differences in inversion polymorphism between the two populations increased over generations under lead contamination in the laboratory. The inversion polymorphism of population originating from the more polluted natural environment was more stable during the experiment, both under conditions with and without lead. Therefore, results showed that inversion polymorphism as a genetic marker reflects a strong signature of adaptation to the local environment, and that historical demographic events and selection are important for both prediction of evolutionary potential and long-term viability of natural populations.     

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.     

Selection in heterogeneous environments maintains the gene arrangement polymorphism of Drosophila pseudoobscura

Chromosomal rearrangements may play an important role in how populations adapt to a local environment. The gene arrangement polymorphism on the third chromosome of Drosophila pseudoobscura is a model system to help determine the role that inversions play in the evolution of this species. The gene arrangements are the likely target of strong selection because they form classical clines across diverse geographic habitats, they cycle in frequency over seasons, and they form stable equilibria in population cages. A numerical approach was developed to estimate the fitness sets for 15 gene arrangement karyotypes in six niches based on a model of selection-migration balance. Gene arrangement frequencies in the six different niches were able to reach a stable meta-population equilibrium that matched the observed gene arrangement frequencies when recursions used the estimated fitnesses with a variety of initial inversion frequencies. These analyses show that a complex pattern of selection is operating in the six niches to maintain the D. pseudoobscura gene arrangement polymorphism. Models of local adaptation predict that the new inversion mutations were able to invade populations because they held combinations of two to 13 local adaptation loci together.     

Comparison of genetic models for analysing the effects of a PvuII polymorphism in the oestrogen receptor 1 (ESR1) gene on prolificacy in an Iberian x Meishan pig population

The effect of a previously reported PvuII polymorphism in oestrogen receptor 1 (ESR1) was analysed in an F(2) population of Iberian x Meishan pigs. We tested three hypotheses: (1) that a causal mutation was fixed in the parental populations, (2) that a causal mutation existed that was in complete linkage disequilibrium with the alleles of the PvuII polymorphism and (3) that a causal mutation existed in linkage disequilibrium within each parental population. The third model was the most plausible based on the available data. ESR1 alleles displayed different patterns of linkage disequilibrium with the causal mutation in each of the parental populations and the PvuII polymorphism was clearly not the causal mutation. As a consequence, the use of the ESR1 mutation for selection must be evaluated for a particular pig population before it is applied.     

Balancing selection is the main force shaping the evolution of innate immunity genes

The evolutionarily recent geographic expansion of humans, and the even more recent development of large, relatively dense human settlements, has exposed our species to new pathogenic environments. Potentially lethal pathogens are likely to have exerted important selective pressures on our genome, so immunity genes can be expected to show molecular signatures of the adaptation of human populations to these recent conditions. While genes related to the acquired immunity system have indeed been reported to show traces of local adaptation, little is known about the response of the innate immunity system. In this study, we analyze the variability patterns in different human populations of fifteen genes related to innate immunity. We have used both single nucleotide polymorphism and sequence data, and through the analysis of interpopulation differentiation, the linkage disequilibrium pattern, and intrapopulation diversity, we have discovered some signatures of positive and especially balancing selection in these genes, thus confirming the importance of the immune system genetic plasticity in the evolutionary adaptive process. Interestingly, the strongest evidence is found in three TLR genes and CD14. These innate immunity genes play a pivotal role, being involved in the primary recognition of pathogens. In general, more evidences of selection appear in the European populations, in some case possibly related to severe population specific pressures. However, we also describe evidence from African populations, which may reflect parallel or long-term selective forces acting in different geographic areas.     

Inferences of Demography and Selection in an African Population of Drosophila melanogaster

It remains a central problem in population genetics to infer the past action of natural selection, and these inferences pose a challenge because demographic events will also substantially affect patterns of polymorphism and divergence. Thus it is imperative to explicitly model the underlying demographic history of the population whenever making inferences about natural selection. In light of the considerable interest in adaptation in African populations of Drosophila melanogaster, which are considered ancestral to the species, we generated a large polymorphism data set representing 2.1 Mb from each of 20 individuals from a Ugandan population of D. melanogaster. In contrast to previous inferences of a simple population expansion in eastern Africa, our demographic modeling of this ancestral population reveals a strong signature of a population bottleneck followed by population expansion, which has significant implications for future demographic modeling of derived populations of this species. Taking this more complex underlying demographic history into account, we also estimate a mean X-linked region-wide rate of adaptation of 6 × 10⁻¹¹/site/generation and a mean selection coefficient of beneficial mutations of 0.0009. These inferences regarding the rate and strength of selection are largely consistent with most other estimates from D. melanogaster and indicate a relatively high rate of adaptation driven by weakly beneficial mutations.     

Conditional neutrality at two adjacent NBS-LRR disease resistance loci in natural populations of Arabidopsis lyrata

We examined patterns of nucleotide diversity at a genomic region containing two linked candidate disease resistance (NBS-LRR) genes in seven populations of the outcrossing plant Arabidopsis lyrata. In comparison with two adjacent control genes and neutral reference genes across the genome, the NBS-LRR genes exhibited elevated nonsynonymous variation and a large number of major-effect polymorphisms causing early stop codons and/or frameshift mutations. In contrast, analysis of synonymous diversity provided no evidence that the region was subject to long-term balancing selection or recent selective sweeps in any of the seven populations surveyed. Also in contrast with earlier surveys of one of these R genes, there was no evidence that the resistance genes or the major-effect mutations were subject to elevated differentiation between populations. We suggest that conditional neutrality in the absence of the corresponding pathogen, rather than long-term balancing selection or local adaptation, may in some circumstances be a significant cause of elevated functional polymorphism at R genes. In contrast with the R genes, analysis of diversity and differentiation at the flanking FERONIA locus showed high population divergence, suggesting local adaptation on this locus controlling male-female signalling during fertilization.     

Influence of genetic architecture on contemporary local evolution in the soapberry bug,Jadera haematoloma: artificial selection on beak length

Little is known about the influence of genetic architecture on local adaptation. We investigated the genetic architecture of the rapid contemporary evolution of mouthparts, the flight polymorphism and life history traits in the soapberry bug Jadera haematoloma (Hemiptera) using laboratory selection. The mouthparts of these seed‐feeding bugs have adapted in 40–50 years by decreasing in length following novel natural selection induced by a host switch to the seeds of an introduced tree with smaller fruits than those of the native host vine. Laboratory selection on beak length in both an ancestral population feeding on the native host and a derived population feeding on the introduced host reveals genetic variance allowing a rapid response (heritabilities of 0.51–0.87) to selection for either longer or shorter beaks. This selection resulted in reverse evolution by restoring long beaks in the derived population and forward evolution by re‐creating short beaks in the ancestral bugs. There were strong genetic correlations (0.68–0.84) in both populations between beak lengths and the frequency of flight morphs, with short beaks associated with short wings. The results reveal a genetically interrelated set of adaptive multivariate traits including both beak length and flight morph. This suite of traits reflects host plant patchiness and seeding phenology. Weaker evidence suggests that egg mass and early egg production may be elements of the same suite. Reversible or forward evolution thus may occur in a broad set of genetically correlated multivariate traits undergoing rapid contemporary adaptation to altered local environments.     

Genome scans of DNA variability in humans reveal evidence for selective sweeps outside of Africa

The last 50,000-150,000 years of human history have been characterized by rapid demographic expansions and the colonization of novel environments outside of sub-Saharan Africa. Mass migrations outside the ancestral species range likely entailed many new selection pressures, suggesting that genetic adaptation to local environmental conditions may have been more prevalent in colonizing populations outside of sub-Saharan Africa. Here we report a test of this hypothesis using genome-wide patterns of DNA polymorphism. We conducted a multilocus scan of microsatellite variability to identify regions of the human genome that may have been subject to continent-specific hitchhiking events. Using published polymorphism data for a total of 624 autosomal loci in multiple populations of humans, we used coalescent simulations to identify candidate loci for geographically restricted selective sweeps. We identified a total of 13 loci that appeared as outliers in replicated population comparisons involving different reference samples for Africa. A disproportionate number of these loci exhibited reduced levels of relative variability in non-African populations alone, suggesting that recent episodes of positive selection have been more prevalent outside of sub-Saharan Africa.     

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.     

Theoretical Study of near Neutrality. I. Heterozygosity and Rate of Mutant Substitution

In order to clarify the nature of "near neutrality" in molecular evolution and polymorphism, extensive simulation studies were performed. Selection coefficients of new mutations are assumed to be small so that both random genetic drift and selection contribute to determining the behavior of mutants. The model also incorporates normally distributed spatial fluctuation of selection coefficients. If the system starts from "average neutrality," it will move to a better adapted state, and most new mutations will become "slightly deleterious." Monte Carlo simulations have indicated that such adaptation is attained, but that the rate of such "progress" is very low for weak selection. In general, the larger the population size, the more effective the selection becomes. Also, as selection becomes weaker, the behavior of the mutants approaches that of completely neutral genes. Thus, the weaker the selection, the smaller is the effect of population size on mutant dynamics. Increase of heterozygosity with population size is very pronounced for subdivided populations. The significance of these results is discussed in relation to various observed facts on molecular evolution and polymorphism, such as generation-time dependency and overdispersion of the molecular clock, or contrasting patterns of DNA and protein polymorphism among some closely related species.     

Limits to local adaptation in six populations of the annual plant Diodia teres

* Local adaptation is common, but tests for adaptive differentiation frequently compare populations from strongly divergent environments, making it unlikely that any influence of stochastic processes such as drift or mutation on local adaptation will be detected. Here, the hypothesis that local adaptation is more likely to develop when the native environments of populations are more distinct than when they are similar was tested. * A reciprocal transplant experiment including two populations from each of three habitats was conducted to determine the pattern of local adaptation. In addition to testing for local adaptation at the population level, the hypothesis was tested that local adaptation is more common between populations from different habitats than between populations from the same habitat. * Local adaptation was not common, but more evidence was found of local adaptation between populations from different habitats than between populations from the same habitat. Two instances of foreign genotype fitness advantage confirm that stochastic processes such as drift can limit local adaptation. * These results are consistent with the hypothesis that stochastic processes can inhibit local adaptation but are more likely to be overwhelmed by natural selection when populations occur in divergent environments.     

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.     

Evolutionary rescue and the limits of adaptation

Populations subject to severe stress may be rescued by natural selection, but its operation is restricted by ecological and genetic constraints. The cost of natural selection expresses the limited capacity of a population to sustain the load of mortality or sterility required for effective selection. Genostasis expresses the lack of variation that prevents many populations from adapting to stress. While the role of relative fitness in adaptation is well understood, evolutionary rescue emphasizes the need to recognize explicitly the importance of absolute fitness. Permanent adaptation requires a range of genetic variation in absolute fitness that is broad enough to provide a few extreme types capable of sustained growth under a stress that would cause extinction if they were not present. This principle implies that population size is an important determinant of rescue. The overall number of individuals exposed to selection will be greater when the population declines gradually under a constant stress, or is progressively challenged by gradually increasing stress. In gradually deteriorating environments, survival at lethal stress may be procured by prior adaptation to sublethal stress through genetic correlation. Neither the standing genetic variation of small populations nor the mutation supply of large populations, however, may be sufficient to provide evolutionary rescue for most populations.     

A new approach for using genome scans to detect recent positive selection in the human genome

Genome-wide scanning for signals of recent positive selection is essential for a comprehensive and systematic understanding of human adaptation. Here, we present a genomic survey of recent local selective sweeps, especially aimed at those nearly or recently completed. A novel approach was developed for such signals, based on contrasting the extended haplotype homozygosity (EHH) profiles between populations. We applied this method to the genome single nucleotide polymorphism (SNP) data of both the International HapMap Project and Perlegen Sciences, and detected widespread signals of recent local selection across the genome, consisting of both complete and partial sweeps. A challenging problem of genomic scans of recent positive selection is to clearly distinguish selection from neutral effects, given the high sensitivity of the test statistics to departures from neutral demographic assumptions and the lack of a single, accurate neutral model of human history. We therefore developed a new procedure that is robust across a wide range of demographic and ascertainment models, one that indicates that certain portions of the genome clearly depart from neutrality. Simulations of positive selection showed that our tests have high power towards strong selection sweeps that have undergone fixation. Gene ontology analysis of the candidate regions revealed several new functional groups that might help explain some important interpopulation differences in phenotypic traits.     

Local Adaptation in the Flowering-Time Gene Network of Balsam Poplar, Populus balsamifera L

Identifying the signature and targets of local adaptation is an increasingly important goal in empirical population genetics. Using data from 443 balsam poplar Populus balsamifera trees sampled from 31 populations, we tested for evidence of geographically variable selection shaping diversity at 27 homologues of the Arabidopsis flowering-time network. These genes are implicated in the control of seasonal phenology, an important determinant of fitness. Using 335 candidate and 412 reference single nucleotide polymorphisms (SNPs), we tested for evidence of local adaptation by searching for elevated population differentiation using F(ST)-based outlier analyses implemented in BayeScan or a Hierarchical Model in Arelquin and by testing for significant associations between allele frequency and environmental variables using BAYENV. A total of 46 SNPs from 14 candidate genes had signatures of local adaptation-either significantly greater population differentiation or significant covariance with one or more environmental variable relative to reference SNP distributions. Only 11 SNPs from two genes exhibited both elevated population differentiation and covariance with one or more environmental variables. Several genes including the abscisic acid gene ABI1B and the circadian clock genes ELF3 and GI5 harbored a large number of SNPs with signatures of local adaptation-with SNPs in GI5 strongly covarying with both latitude and precipitation and SNPs in ABI1B strongly covarying with temperature. In contrast to several other systems, we find little evidence that photoreceptors, including phytochromes, play an important role in local adaptation. Our results additionally show that detecting local adaptation is sensitive to the analytical approaches used and that model-based significance thresholds should be viewed with caution.     

Spatiotemporal variation in local adaptation of a specialist insect herbivore to its long‐lived host plant

Local adaptation of interacting species to one another indicates geographically variable reciprocal selection. This process of adaptation is central in the organization and maintenance of genetic variation across populations. Given that the strength of selection and responses to it often vary in time and space, the strength of local adaptation should in theory vary between generations and among populations. However, such spatiotemporal variation has rarely been explicitly demonstrated in nature and local adaptation is commonly considered to be relatively static. We report persistent local adaptation of the short‐lived herbivore Abrostola asclepiadis to its long‐lived host plant Vincetoxicum hirundinaria over three successive generations in two studied populations and considerable temporal variation in local adaptation in six populations supporting the geographic mosaic theory. The observed variation in local adaptation among populations was best explained by geographic distance and population isolation, suggesting that gene flow reduces local adaptation. Changes in herbivore population size did not conclusively explain temporal variation in local adaptation. Our results also imply that short‐term studies are likely to capture only a part of the existing variation in local adaptation.     

Multilocus coalescent analysis of haemoglobin differentiation between low- and high-altitude populations of crested ducks (Lophonetta specularioides)

Hypoxia is a key factor determining survival, and haemoglobins are targets of selection in species native to high-altitude regions. We studied population genetic structure and evaluated evidence for local adaptation in the crested duck (Lophonetta specularioides). Differentiation, gene flow and time since divergence between highland and lowland populations were assessed for three haemoglobin genes (α(A) , α(D) , β(A) ) and compared to seven reference loci (six autosomal introns and mtDNA). Four derived amino acid replacements were found in the globin genes that had elevated Φ(ST) values between the Andean highlands and Patagonian lowlands. A single β(A) -globin polymorphism at a site known to influence O(2) affinity was fixed for different alleles in the two populations, whereas three α(A) - and α(D) -globin polymorphisms exhibited high heterozygosity in the highlands but not in the lowlands. Coalescent analyses supported restricted gene flow for haemoglobin alleles and mitochondrial DNA but nonzero gene flow for the introns. Simulating genetic data under a drift-migration model of selective neutrality, the β(A) -globin fell outside the 95% confidence limit of simulated data, suggesting that directional selection is maintaining different variants in the contrasting elevational environments, thereby restricting migration of β(A) -globin alleles. The α(A) - and α(D) -globins, by contrast, did not differ from the simulated values, suggesting that variants in these genes are either selectively neutral, or that the effects of selection could not be differentiated from background levels of population structure and linkage disequilibrium. This study illustrates the combined effects of selection and population history on inferring levels of population divergence for a species distributed across an altitudinal gradient in which selection for hypoxia resistance has likely played an important role.