A population genetic test of selection at the molecular level

Prevalence of CCR5delta32 in Northeastern Iran

Background

A 32-base pair deletion (∆32) in the open reading frame (ORF) of C-C motif chemokine receptor 5 (CCR5) seems to be a protective variant against immune system diseases, especially human immunodeficiency virus type 1 (HIV-1). We aimed to assess the frequency of CCR5∆32 in the healthy Iranian population.

Methods

In this study, 400 normal samples from Khorasan, northeastern Iran, were randomly selected. The frequency of CCR5∆32 carriers was investigated using PCR analysis. Allele prevalence and the fit to the Hardy-Weinberg equilibrium were analyzed.

Results

The prevalence of CCR5∆32 in the northeastern population of Iran was 0.016. Four hundred samples were studied, among which one with CCR5^∆32/∆32 and 11 with CCR5^Wild/∆32 genotype were detected.

Conclusion

This study was the first investigation for an assessment of the prevalence of CCR5∆32 in northeastern Iran. The low prevalence of CCR5∆32 allele in the Iranian population may result in the increased susceptibility to HIV-1. In addition, this prevalence is the same as that of reported in East Asia, while is lower than that in the Europeans.

Early selection of bZIP73 facilitated adaptation of japonica rice to cold climates

Cold stress is a major factor limiting production and geographic distribution of rice (Oryza sativa). Although the growth range of japonica subspecies has expanded northward compared to modern wild rice (O. rufipogon), the molecular basis of the adaptation remains unclear. Here we report bZIP73, a bZIP transcription factor-coding gene with only one functional polymorphism (+511 G>A) between the two subspecies japonica and indica, may have facilitated japonica adaptation to cold climates. We show the japonica version of bZIP73 (bZIP73^Jap) interacts with bZIP71 and modulates ABA levels and ROS homeostasis. Evolutionary and population genetic analyses suggest bZIP73 has undergone balancing selection; the bZIP73^Jap allele has firstly selected from standing variations in wild rice and likely facilitated cold climate adaptation during initial japonica domestication, while the indica allele bZIP73^Ind was subsequently selected for reasons that remain unclear. Our findings reveal early selection of bZIP73^Jap may have facilitated climate adaptation of primitive rice germplasms.

Japonica rice can grow further north than wild or indica rice and is more tolerant of cold climates. Here, the authors show that bZIP73 likely underwent selection in the early phase of rice domestication to facilitate cold tolerance in japonica by modulating ABA and ROS homeostasis.

Contemporary evolution of a Lepidopteran species, Heliothis virescens, in response to modern agricultural practices

Adaptation to human-induced environmental change has the potential to profoundly influence the genomic architecture of affected species. This is particularly true in agricultural ecosystems, where anthropogenic selection pressure is strong. Heliothis virescens primarily feeds on cotton in its larval stages, and US populations have been declining since the widespread planting of transgenic cotton, which endogenously expresses proteins derived from Bacillus thuringiensis (Bt). No physiological adaptation to Bt toxin has been found in the field, so adaptation in this altered environment could involve (i) shifts in host plant selection mechanisms to avoid cotton, (ii) changes in detoxification mechanisms required for cotton-feeding vs. feeding on other hosts or (iii) loss of resistance to previously used management practices including insecticides. Here, we begin to address whether such changes occurred in H. virescens populations between 1997 and 2012, as Bt-cotton cultivation spread through the agricultural landscape. For our study, we produced an H. virescens genome assembly and used this in concert with a ddRAD-seq-enabled genome scan to identify loci with significant allele frequency changes over the 15-year period. Genetic changes at a previously described H. virescens insecticide target of selection were detectable in our genome scan and increased our confidence in this methodology. Additional loci were also detected as being under selection, and we quantified the selection strength required to elicit observed allele frequency changes at each locus. Potential contributions of genes near loci under selection to adaptive phenotypes in the H. virescens cotton system are discussed.

Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel

Insecticide resistance is considered a classic model of microevolution, where a strong selective agent is applied to a large natural population, resulting in a change in frequency of alleles that confer resistance. While many insecticide resistance variants have been characterized at the gene level, they are typically single genes of large effect identified in highly resistant pest species. In contrast, multiple variants have been implicated in DDT resistance in Drosophila melanogaster; however, only the Cyp6g1 locus has previously been shown to be relevant to field populations. Here we use genome-wide association studies (GWAS) to identify DDT-associated polygenes and use selective sweep analyses to assess their adaptive significance. We identify and verify two candidate DDT resistance loci. A largely uncharacterized gene, CG10737, has a function in muscles that ameliorates the effects of DDT, while a putative detoxifying P450, Cyp6w1, shows compelling evidence of positive selection.

Ecological release and venom evolution of a predatory marine snail at Easter Island

BACKGROUND

Ecological release is coupled with adaptive radiation and ecological diversification yet little is known about the molecular basis of phenotypic changes associated with this phenomenon. The venomous, predatory marine gastropod Conus miliaris has undergone ecological release and exhibits increased dietary breadth at Easter Island.

METHODOLOGY/PRINCIPAL FINDINGS

We examined the extent of genetic differentiation of two genes expressed in the venom of C. miliaris among samples from Easter Island, American Samoa and Guam. The population from Easter Island exhibits unique frequencies of alleles that encode distinct peptides at both loci. Levels of divergence at these loci exceed observed levels of divergence observed at a mitochondrial gene region at Easter Island.

CONCLUSIONS/SIGNIFICANCE

Patterns of genetic variation at two genes expressed in the venom of this C. miliaris suggest that selection has operated at these genes and contributed to the divergence of venom composition at Easter Island. These results show that ecological release is associated with strong selection pressures that promote the evolution of new phenotypes.

Introgression from modern hybrid varieties into landrace populations of maize (Zea mays ssp. mays L.) in central Italy

Landraces are domesticated local plant varieties that did not experience a deliberate and intensive selection during a formal breeding programme. In Europe, maize landraces are still cultivated, particularly in marginal areas where traditional farming is often practiced. Here, we have studied the evolution of flint maize landraces from central Italy over 50 years of on-farm cultivation, when dent hybrid varieties were introduced and their use was widespread. We have compared an 'old' collection, obtained during the 1950s, before the introduction of hybrids, and a recent collection of maize landraces. For comparison, a sample of maize landraces from north Italy, and of improved germplasm, including hybrids and inbred lines were also used. A total of 296 genotypes were analysed using 21 microsatellites. Our results show that the maize landraces collected in the last 5-10 years have evolved directly from the flint landrace gene pool cultivated in central Italy before the introduction of modern hybrids. The population structure, diversity and linkage disequilibrium analyses indicate a significant amount of introgression from hybrid varieties into the recent landrace populations. No evidence of genetic erosion of the maize landraces was seen, suggesting that in situ conservation of landraces is an efficient strategy for preserving genetic diversity. Finally, the level of introgression detected was very variable among recent landraces, with most of them showing a low level of introgression; this suggests that coexistence between different types of agriculture is possible, with the adoption of correct practices that are aimed at avoiding introgression from undesired genetic sources.

Signatures of selection in natural populations adapted to chronic pollution

Background

Populations of the teleost fish Fundulus heteroclitus appear to flourish in heavily polluted and geographically separated Superfund sites. Populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) have independently evolved adaptive resistance to chemical pollutants. In these polluted populations, natural selection likely has altered allele frequencies of loci that affect fitness or that are linked to these loci. The aim of this study was to identify loci that exhibit non-neutral behavior in the F. heteroclitus genome in polluted populations versus clean reference populations.

Results

To detect signatures of natural selection and thus identify genetic bases for adaptation to anthropogenic stressors, we examined allele frequencies for many hundreds of amplified fragment length polymorphism markers among populations of F. heteroclitus. Specifically, we contrasted populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) to clean reference populations flanking the polluted sites. When empirical F_ST values were compared to a simulated distribution of F_ST values, 24 distinct outlier loci were identified among pairwise comparisons of pollutant impacted F. heteroclitus populations and both surrounding reference populations. Upon removal of all outlier loci, there was a strong correlation (R² = 0.79, p < 0.0001) between genetic and geographical distance. This apparently neutral evolutionary pattern was not evident when outlier loci were included (R² = 0.092, p = 0.0721). Two outlier loci were shared between New Bedford Harbor and Elizabeth River populations, and two different loci were shared between Newark Bay and Elizabeth River populations.

Conclusion

In total, 1% to 6% of loci are implicated as being under selection or linked to areas of the genome under selection in three F. heteroclitus populations that reside in polluted estuaries. Shared loci among polluted sites indicate that selection may be acting on multiple loci involved in adaptation, and loci shared between polluted sites potentially are involved in a generalized adaptive response.

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.

Identifying selected regions from heterozygosity and divergence using a light-coverage genomic dataset from two human populations

When a selective sweep occurs in the chromosomal region around a target gene in two populations that have recently separated, it produces three dramatic genomic consequences: 1) decreased multi-locus heterozygosity in the region; 2) elevated or diminished genetic divergence (F_ST) of multiple polymorphic variants adjacent to the selected locus between the divergent populations, due to the alternative fixation of alleles; and 3) a consequent regional increase in the variance of F_ST (S²F_ST) for the same clustered variants, due to the increased alternative fixation of alleles in the loci surrounding the selection target. In the first part of our study, to search for potential targets of directional selection, we developed and validated a resampling-based computational approach; we then scanned an array of 31 different-sized moving windows of SNP variants (5–65 SNPs) across the human genome in a set of European and African American population samples with 183,997 SNP loci after correcting for the recombination rate variation. The analysis revealed 180 regions of recent selection with very strong evidence in either population or both. In the second part of our study, we compared the newly discovered putative regions to those sites previously postulated in the literature, using methods based on inspecting patterns of linkage disequilibrium, population divergence and other methodologies. The newly found regions were cross-validated with those found in nine other studies that have searched for selection signals. Our study was replicated especially well in those regions confirmed by three or more studies. These validated regions were independently verified, using a combination of different methods and different databases in other studies, and should include fewer false positives. The main strength of our analysis method compared to others is that it does not require dense genotyping and therefore can be used with data from population-based genome SNP scans from smaller studies of humans or other species.

Spatial variation in colour morph, spotting and allozyme frequencies in the candy-stripe spider, Enoplognatha ovata (Theridiidae) on two Swedish archipelagos

The selective significance, if any, of many invertebrate visible polymorphisms is still not fully understood. Here we examine colour- and black spotting-morph frequencies in the spider Enoplognatha ovata in populations on two Swedish archipelagos with respect to different spatial scales and, in one archipelago, against the background of variation at four putative neutral allozyme marker loci. Every population studied was polymorphic for colour and 28 out of 30 contained all three colour morphs--lineata, redimita and ovata. We found no evidence for a breakdown in the traditional colour morph designation previously suggested for other northern European populations of this species. For colour there is no significant heterogeneity at spatial scales greater than between local populations within islands. Black spotting frequencies show a similar lack of pattern over larger spatial scales except that there are significant differences between the Stockholm and Göteborg archipelagos. Measures of population differentiation (theta) within the Stockholm islands for the two visible systems show them to be significantly more differentiated than the neutral markers, suggesting local selection acting on them in a population-specific manner. On the basis of previous observations and the distribution of spotting phenotypes on a European scale, it is argued that thermal selection might operate on black spotting during the juvenile stages favouring more spots in continental climates. It is not clear what selective forces act on colour.

Linkage disequilibrium as a tool for detecting signatures of natural selection

Natural selection has been theoretically and empirically proven to alter patterns of linkage disequilibrium (LD). Reciprocally, recombination, the driving force behind LD, modifies the signature of natural selection by reintroducing variation in a punctuate manner across the genome. To date, efforts to identify genes that have been subjected to historical selective pressure by examining polymorphic variation and allelic association have frequently fallen short of unambiguously distinguishing selection from other biological mechanisms. Contemporary genetic maps that describe LD in fine detail represent a much needed tool that can be exploited by researchers aiming to tease apart these opposing signals.

Which came first: insecticides or resistance?

Mutations that confer resistance to insecticides are well documented. However, so far, we have been unable to determine whether these mutations arose before or after the introduction of insecticides. Recently, a landmark study showed that resistance to Malathion can be detected in pinned specimens of Australian sheep blowflies that were collected before the introduction of the insecticide. This finding has numerous implications for our understanding of the prevalence of resistance to new compounds. It also indicates that pre-existing resistance alleles might not carry the fitness cost that is associated with new mutations.

Genetic susceptibility to infectious diseases: big is beautiful, but will bigger be even better?

Genetic epidemiology, including twin studies, provides robust evidence that genetic variation in human populations contributes to susceptibility to infectious disease. One of the major limitations of studies that attempt to identify the genes and mechanisms that underlie this susceptibility has been lack of power caused by small sample size. With the development of novel technologies, burgeoning information on the human genome, the HapMap project, and human genetic diversity, we are at the beginning of a new era in the study of the genetics of complex diseases. This review looks afresh at the epidemiological evidence that supports a role for genetics in susceptibility to infectious disease, examines the somewhat limited achievements to date, and discusses current advances in methodology and technology that will potentially lead to translational data in the future.

Detecting natural selection on cis-regulatory DNA

Changes in transcriptional regulation play an important role in the genetic basis for evolutionary change. Here I review a growing body of literature that seeks to determine the forces governing the non-coding regulatory sequences underlying these changes. I address the challenges present in studying natural selection without the familiar structure and regularity of protein-coding sequences, but show that most tests of neutrality that have been used for coding regions are applicable to non-coding regions, albeit with some caveats. While some experimental investment is necessary to identify heritable regulatory variation, the most basic inferences about selection require very little functional information. A growing body of research on cis-regulatory variation has uncovered all the forms of selection common to coding regions, in addition to novel forms of selection. An emerging pattern seems to be the ubiquity of local adaptation and balancing selection, possibly due to the greater freedom organisms have to fine-tune gene expression without changing protein function. It is clear from multiple single locus and whole genome studies of non-coding regulatory DNA that the effects of natural selection reach far beyond the start and stop codons.

An analysis of selection on a colour polymorphism in the northern leopard frog

In this study, we investigated the role of selection in the maintenance of a dorsal colour polymorphism in natural populations of the northern leopard frog, Rana pipiens. We determined genetic structure both spatially and temporally from a suite of putatively neutral molecular markers and tested whether or not the colour locus exhibited patterns of genetic variation that differed from those of the neutral loci. Spatial genetic structure at the colour locus was indistinguishable from structure at neutral loci [95% confidence intervals of F(ST) (neutral) = (0.07, 0.35), F(ST) (colour locus) = 0.114]. In the temporal analysis, we found that the variance among populations in the change in allele frequency at the colour locus (equal to 0.004) lies within the 95% confidence intervals for the variance among populations in changes in allele frequencies at neutral loci. In light of our inability to show evidence for the selective maintenance of the colour polymorphism, we used computer simulations to infer the power of our spatial and temporal techniques to detect selection. The computer simulations showed that although the strength of selection (s) would need to be relatively strong to have been detected by the temporal approach (s = 0.1-0.4, depending on the model tested), the spatial analysis would have detected all but weak selection (s = 0.01-0.04, depending on the model tested). This study illustrates the importance of using a locus comparison approach to detect evidence for selective maintenance before conducting studies to measure the selective mechanisms maintaining a polymorphism.

A nonparametric test reveals selection for rapid flowering in the Arabidopsis genome

The detection of footprints of natural selection in genetic polymorphism data is fundamental to understanding the genetic basis of adaptation, and has important implications for human health. The standard approach has been to reject neutrality in favor of selection if the pattern of variation at a candidate locus was significantly different from the predictions of the standard neutral model. The problem is that the standard neutral model assumes more than just neutrality, and it is almost always possible to explain the data using an alternative neutral model with more complex demography. Today's wealth of genomic polymorphism data, however, makes it possible to dispense with models altogether by simply comparing the pattern observed at a candidate locus to the genomic pattern, and rejecting neutrality if the pattern is extreme. Here, we utilize this approach on a truly genomic scale, comparing a candidate locus to thousands of alleles throughout the Arabidopsis thaliana genome. We demonstrate that selection has acted to increase the frequency of early-flowering alleles at the vernalization requirement locus FRIGIDA. Selection seems to have occurred during the last several thousand years, possibly in response to the spread of agriculture. We introduce a novel test statistic based on haplotype sharing that embraces the problem of population structure, and so should be widely applicable.

Selective sweep mapping of genes with large phenotypic effects

Many domestic dog breeds have originated through fixation of discrete mutations by intense artificial selection. As a result of this process, markers in the proximity of genes influencing breed-defining traits will have reduced variation (a selective sweep) and will show divergence in allele frequency. Consequently, low-resolution genomic scans can potentially be used to identify regions containing genes that have a major influence on breed-defining traits. We model the process of breed formation and show that the probability of two or three adjacent marker loci showing a spurious signal of selection within at least one breed (i.e., Type I error or false-positive rate) is low if highly variable and moderately spaced markers are utilized. We also use simulations with selection to demonstrate that even a moderately spaced set of highly polymorphic markers (e.g., one every 0.8 cM) has high power to detect regions targeted by strong artificial selection in dogs. Further, we show that a gene responsible for black coat color in the Large Munsterlander has a 40-Mb region surrounding the gene that is very low in heterozygosity for microsatellite markers. Similarly, we survey 302 microsatellite markers in the Dachshund and find three linked monomorphic microsatellite markers all within a 10-Mb region on chromosome 3. This region contains the FGFR3 gene, which is responsible for achondroplasia in humans, but not in dogs. Consequently, our results suggest that the causative mutation is a gene or regulatory region closely linked to FGFR3.

Searching for signals of evolutionary selection in 168 genes related to immune function

Pathogens have played a substantial role in human evolution, with past infections shaping genetic variation at loci influencing immune function. We selected 168 genes known to be involved in the immune response, genotyped common single nucleotide polymorphisms across each gene in three population samples (CEPH Europeans from Utah, Han Chinese from Guangxi, and Yoruba Nigerians from Southwest Nigeria) and searched for evidence of selection based on four tests for non-neutral evolution: minor allele frequency (MAF), derived allele frequency (DAF), Fst versus heterozygosity and extended haplotype homozygosity (EHH). Six of the 168 genes show some evidence for non-neutral evolution in this initial screen, with two showing similar signals in independent data from the International HapMap Project. These analyses identify two loci involved in immune function that are candidates for having been subject to evolutionary selection, and highlight a number of analytical challenges in searching for selection in genome-wide polymorphism data.

The case for selection at CCR5-Delta32

The C-C chemokine receptor 5, 32 base-pair deletion (CCR5-Δ32) allele confers strong resistance to infection by the AIDS virus HIV. Previous studies have suggested that CCR5-Δ32 arose within the past 1,000 y and rose to its present high frequency (5%–14%) in Europe as a result of strong positive selection, perhaps by such selective agents as the bubonic plague or smallpox during the Middle Ages. This hypothesis was based on several lines of evidence, including the absence of the allele outside of Europe and long-range linkage disequilibrium at the locus. We reevaluated this evidence with the benefit of much denser genetic maps and extensive control data. We find that the pattern of genetic variation at CCR5-Δ32 does not stand out as exceptional relative to other loci across the genome. Moreover using newer genetic maps, we estimated that the CCR5-Δ32 allele is likely to have arisen more than 5,000 y ago. While such results can not rule out the possibility that some selection may have occurred at C-C chemokine receptor 5 (CCR5), they imply that the pattern of genetic variation seen atCCR5-Δ32 is consistent with neutral evolution. More broadly, the results have general implications for the design of future studies to detect the signs of positive selection in the human genome.

Sabeti and colleagues use dense genetic maps to show that the HIV-resistance CCR5-Δ32 allele is more than 5,000 years old and is likely to have been under mainly neutral selection.

Geographical distribution of selected and putatively neutral SNPs in Southeast Asian malaria parasites

Loci targeted by directional selection are expected to show elevated geographical population structure relative to neutral loci, and a flurry of recent papers have used this rationale to search for genome regions involved in adaptation. Studies of functional mutations that are known to be under selection are particularly useful for assessing the utility of this approach. Antimalarial drug treatment regimes vary considerably between countries in Southeast Asia selecting for local adaptation at parasite loci underlying resistance. We compared the population structure revealed by 10 nonsynonymous mutations (nonsynonymous single-nucleotide polymorphisms [nsSNPs]) in four loci that are known to be involved in antimalarial drug resistance, with patterns revealed by 10 synonymous mutations (synonymous single-nucleotide polymorphisms [sSNPs]) in housekeeping genes or genes of unknown function in 755 Plasmodium falciparum infections collected from 13 populations in six Southeast Asian countries. Allele frequencies at known nsSNPs underlying resistance varied markedly between locations (F(ST) = 0.18-0.66), with the highest frequencies on the Thailand-Burma border and the lowest frequencies in neighboring Lao PDR. In contrast, we found weak but significant geographic structure (F(ST) = 0-0.14) for 8 of 10 sSNPs. Importantly, all 10 nsSNPs showed significantly higher F(ST) (P < 8 x 10(-5)) than simulated neutral expectations based on observed F(ST) values in the putatively neutral sSNPs. This result was unaffected by the methods used to estimate allele frequencies or the number of populations used in the simulations. Given that dense single-nucleotide polymorphism (SNP) maps and rapid SNP assay methods are now available for P. falciparum, comparing genetic differentiation across the genome may provide a valuable aid to identifying parasite loci underlying local adaptation to drug treatment regimes or other selective forces. However, the high proportion of polymorphic sites that appear to be under balancing selection (or linked to selected sites) in the P. falciparum genome violates the central assumption that selected sites are rare, which complicates identification of outlier loci, and suggests that caution is needed when using this approach.

Population genetic and phylogenetic evidence for positive selection on regulatory mutations at the factor VII locus in humans

The abundance of cis-regulatory polymorphisms in humans suggests that many may have been important in human evolution, but evidence for their role is relatively rare. Four common polymorphisms in the 5' promoter region of factor VII (F7), a coagulation factor, have been shown to affect its transcription and protein abundance both in vitro and in vivo. Three of these polymorphisms have low-frequency alleles that decrease expression of F7 and may provide protection against myocardial infarction (heart attacks). The fourth polymorphism has a minor allele that increases the level of transcription. To look for evidence of natural selection on the cis-regulatory variants flanking F7, we genotyped three of the polymorphisms in six Old World populations for which we also have data from a group of putatively neutral SNPs. Our population genetic analysis shows evidence for selection within humans; surprisingly, the strongest evidence is due to a large increase in frequency of the high-expression variant in Singaporean Chinese. Further characterization of a Japanese population shows that at least part of the increase in frequency of the high-expression allele is found in other East Asian populations. In addition, to examine interspecific patterns of selection we sequenced the homologous 5' noncoding region in chimpanzees, bonobos, a gorilla, an orangutan, and a baboon. Analysis of these data reveals an excess of fixed differences within transcription factor binding sites along the human lineage. Our results thus further support the hypothesis that regulatory mutations have been important in human evolution.

Using genome scans of DNA polymorphism to infer adaptive population divergence

Elucidating the genetic basis of adaptive population divergence is a goal of central importance in evolutionary biology. In principle, it should be possible to identify chromosomal regions involved in adaptive divergence by screening genome-wide patterns of DNA polymorphism to detect the locus-specific signature of positive directional selection. In the case of spatially separated populations that inhabit different environments or sympatric populations that exploit different ecological niches, it is possible to identify loci that underlie divergently selected traits by comparing relative levels of differentiation among large numbers of unlinked markers. In this review I first address the question of whether diversifying selection on polygenic traits can be expected to produce predictable patterns of allelic variation at the underlying quantitative trait loci (QTL), and whether the locus-specific effects of selection can be reliably detected against the genome-wide backdrop of stochastic variability. I then review different approaches that have been developed to identify loci involved in adaptive population divergence and I discuss the relative merits of model-based approaches that rely on assumptions about population structure vs. model-free approaches that are based on empirical distributions of summary statistics. Finally, I consider the evolutionary and functional insights that might be gained by conducting genome scans for loci involved in adaptive population divergence.

Interaction among different heterocromatic variants in the grasshopper Dichroplus elongatus

Simultaneous chromosome polymorphisms for supernumerary elements allow us to analyse the relationships among different forms of heterochromatic variation in nature. We report simultaneous variation patterns for supernumerary segments in chromosomes S10 (SS10), S9 (SS9) and S6 (SS6) and B chromosomes in nine populations of the grasshopper Dichroplus elongatus from two biogeographic provinces from east Argentina. Our results show spatial chromosome differentiation for three out of four supernumerary heterochromatic variants (B chromosomes, SS6 and SS10). The incidence of B chromosomes was negatively correlated with the SS10 frequency. The distribution pattern analysis shows different degree of differentiation among populations for each supernumerary heterochromatic variant suggesting that the detected chromosome variation cannot be explained by interaction between migration and genetic drift. Moreover, the observed population chromosome differentiation was not in agreement with the hierarchical analysis of molecular of heterogeneity at mitochondrial DNA level (mtDNA). The present results point out the importance of the interaction among heterochromatic variants in the chromosome intraspecific variation in east Argentina natural populations of the grasshopper D. elongatus.

World-wide survey of an Accord insertion and its association with DDT resistance in Drosophila melanogaster

Previous work showed that insecticide resistance in Drosophila melanogaster is correlated with the insertion of an Accord-like element into the 5' region of the cytochrome P450 gene, Cyp6g1. Here, we study the distribution of the Accord-like element in 673 recently collected D. melanogaster lines from 34 world-wide populations. We also examine the extent of microsatellite variability along a 180-kilobase (kb) genomic region of chromosome II encompassing the resistance gene. We confirm a 100% correlation of the Accord insertion with insecticide resistance and a significant reduction in variability extending at least 20 kb downstream of the Cyp6g1 gene. The frequency of the Accord insertion differs significantly between East African (32-55%) and nonAfrican (85-100%) populations. This pattern is consistent with a selective sweep driving the Accord insertion close to fixation in nonAfrican populations as a result of the insecticide resistance phenotype it confers. This study confirms that hitchhiking mapping can be used to identify beneficial mutations in natural populations.

Assessing the effect of natural selection in malaria parasites

There are few concepts that have been used across disciplines; one of them is natural selection. The impact that this process has on parasite genetic diversity is reviewed here by discussing examples on drug resistance and vaccine antigens. Emphasis is made on how mechanisms need to be addressed rather than associations, and how such investigations were out of reach of biomedical researchers only a decade ago.

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.

Geographic and haplotype structure of candidate type 2 diabetes susceptibility variants at the calpain-10 locus

Recently, a positional cloning study proposed that haplotypes at the calpain-10 locus (CAPN10) are associated with increased risk of type 2 diabetes, or non-insulin-dependent diabetes mellitus, in Mexican Americans, Finns, and Germans. To inform the interpretation of the original mapping results and to look for evidence for the action of natural selection on CAPN10, we undertook a population-based genotyping survey of the candidate susceptibility variants. First, we genotyped sites 43, 19, and 63 (the haplotype-defining variants previously proposed) and four closely linked SNPs, in 561 individuals from 11 populations from five continents, and we examined the linkage disequilibrium among them. We then examined the ancestral state of these sites by sequencing orthologous portions of CAPN10 in chimpanzee and orangutan (the identity of sites 43 and 19 was further investigated in a limited sample of other great apes and Old World and New World monkeys). Our survey suggests larger-than-expected differences in the distribution of CAPN10 susceptibility variants between African and non-African populations, with common, derived haplotypes in European and Asian samples (including one of two proposed risk haplotypes) being rare or absent in African samples. These results suggest a history of positive natural selection at the locus, resulting in significant geographic differences in polymorphism frequencies. The relationship of these differences to disease risk is discussed.

Complex signatures of natural selection at the Duffy blood group locus

The Duffy blood group locus (FY) has long been considered a likely target of natural selection, because of the extreme pattern of geographic differentiation of its three major alleles (FY*B, FY*A, and FY*O). In the present study, we resequenced the FY region in samples of Hausa from Cameroon (fixed for FY*O), Han Chinese (fixed for FY*A), Italians, and Pakistanis. Our goals were to characterize the signature of directional selection on FY*O in sub-Saharan Africa and to understand the extent to which natural selection has also played a role in the extreme geographic differentiation of the other derived allele at this locus, FY*A. The data from the FY region are compared with the patterns of variation observed at 10 unlinked, putatively neutral loci from the same populations, as well as with theoretical expectations from the neutral-equilibrium model. The FY region in the Hausa shows evidence of directional selection in two independent properties of the data (i.e., level of sequence variation and frequency spectrum), observations that are consistent with the FY*O mutation being the target. The Italian and Chinese FY data show patterns of variation that are very unusual, particularly with regard to frequency spectrum and linkage disequilibrium, but do not fit the predictions of any simple model of selection. These patterns may represent a more complex and previously unrecognized signature of positive selection.

A single gene (yes) controls pigmentation of eyes and scales in Heliothis virescens

A yellow-eyed mutant was discovered in a strain of Heliothis virescens, the tobacco budworm, that already exhibited a mutation for yellow scale, y. We investigated the inheritance of these visible mutations as candidate markers for transgenesis. Yellow eye was controlled by a single, recessive, autosomal factor, the same type of inheritance previously known for y. Presence of the recombinant mutants with yellow scales and wild type eyes in test crosses indicated independent segregation of genes for these traits. The recombinant class with wild type scales and yellow eyes was completely absent and there was a corresponding increase of the double mutant parental class having yellow scales and yellow eyes. These results indicated that a single factor for yellow eye also controlled yellow scales independently of y. This gene was named yes, for yellow eye and scale. We hypothesize that yes controls both eye and scale color through a deficiency in transport of pigment precursors in both the ommochrome and melanin pathways. The unlinked gene y likely controls an enzyme affecting the melanin pathway only. Both y and yes segregated independently of AceIn, acetylcholinesterase insensitivity, and sodium channel hscp, which are genes related to insecticide resistance.

Population genomics: genome-wide sampling of insect populations

Modern population genetics underwent a major paradigm shift during the last decade of the 20th century with the discovery that thousands of genes of known function and position in a genome can be analyzed simultaneously in a single individual. The impact of this technology on insect population genetics is potentially profound. Sampling distributions of genetic statistics can now be derived from many individual loci or among many segregating sites within a gene. Inferences regarding random mating, gene flow, effective population sizes, disequilibrium, and relatedness among populations can now be based on patterns of variation at many loci. More importantly, genome-wide sampling enables population geneticists to distinguish effects that act on the whole genome from those that act on individual loci or nucleotides. We introduce the term "population genomics" to describe the process of simultaneous sampling of numerous variable loci within a genome and the inference of locus-specific effects from the sample distributions. The four critical assumptions implicit in the population genomics approach are explained in detail. Studies adopting this paradigm are reviewed, and the steps necessary to complete a population genomics study are outlined.

Genetic microstructure in two Spanish cat populations. I: genic diversity, gene flow and selection

The genetic microstructure of two Spanish cat populations (in Barcelona and Alicante) was studied. These populations were in Hardy-Weinberg equilibrium at the locus O . There was significant genetic heterogeneity for most of the loci studied at the colony level, especially, and the subpopulation level in Barcelona, and at the subpopulation level in Alicante, although the amount of heterogeneity was relatively small compared with that found in other mammal species. Therefore, the major part of the gene diversity found was at the level of the small population structure (colonies). This agrees quite well with high theoretical gene flow estimates. The Lewontin-Krakauer test showed, in some cases, significant F tests, which appears to indicate the presence of several selective events on some of the loci studied (diversifying and/or unifying), although some other causes, such as differential gene flow, some different historical and demographic parameters at the time of the introduction of the different alleles within the cities studied, as well as different mutation rates of the loci analyzed, could produce significant Lewontin-Krakauer test values, as well.

Natural selection mapping of the warfarin-resistance gene

In theory, genes under natural selection can be revealed by unique patterns of linkage disequilibrium (LD) and polymorphism at physically linked loci. However, given the effects of recombination and mutation, the physical extent and persistence of LD patterns in natural populations is uncertain. To assess the LD signature of selection, we survey variation in 26 microsatellite loci spanning an approximately 32-cM region that includes the warfarin-resistance gene (Rw) in five wild rat populations having resistance levels between 0 and 95%. We find a high frequency of heterozygote deficiency at microsatellite loci in resistant populations, and a negative association between gene diversity (H) and resistance. Contrary to previous studies, these data suggest that directional rather than overdominant selection may predominate during periods of intense anticoagulant treatment. In highly resistant populations, extensive LD was observed over a chromosome segment spanning approximately 14% of rat chromosome 1. In contrast, LD in a moderately resistant population was more localized and, in conjunction with likelihood ratios, allowed assignment of Rw to a 2. 2-cM interval. Within this genomic window, a diagnostic marker, D1Rat219, assigned 91% of rats to the correct resistance category. These results further demonstrate that "natural selection mapping" in field populations can detect and map major fitness-related genes, and question overdominance as the predominant mode of selection in anticoagulant-resistant rat populations.

A principal target of human immunity to malaria identified by molecular population genetic and immunological analyses

New strategies are required to identify the most important targets of protective immunity in complex eukaryotic pathogens. Natural selection maintains allelic variation in some antigens of the malaria parasite Plasmodium falciparum. Analysis of allele frequency distributions could identify the loci under most intense selection. The merozoite surface protein 1 (Msp1) is the most-abundant surface component on the erythrocyte-invading stage of P. falciparum. Immunization with whole Msp1 has protected monkeys completely against homologous and partially against non-homologous parasite strains. The single-copy msp1 gene, of about 5 kilobases, has highly divergent alleles with stable frequencies in endemic populations. To identify the region of msp1 under strongest selection to maintain alleles within populations, we studied multiple intragenic sequence loci in populations in different regions of Africa and Southeast Asia. On both continents, the locus with the lowest inter-population variance in allele frequencies was block 2, indicating selection in this part of the gene. To test the hypothesis of immune selection, we undertook a large prospective longitudinal cohort study. This demonstrated that serum IgG antibodies against each of the two most frequent allelic types of block 2 of the protein were strongly associated with protection from P. falciparum malaria.

Do discrepancies between microsatellite and allozyme variation reveal differential selection between sea and lagoon in the sea bass (Dicentrarchus labrax)?

In the present study the genetic structure of Dicentrarchus labrax (14 samples from the Mediterranean) was analysed at six microsatellite loci, in order to test the hypothesis that some enzymatic loci undergo selection between marine and lagoon habitat. Eight of the 14 samples were analysed at both microsatellite and allozyme markers. The analysis of the genetic variation among the Mediterranean samples showed that (i) &Fcirc;ST values obtained with the six microsatellite loci were much smaller than those obtained with the 28 allozymes and (ii) microsatellite loci seemed to reflect more the geographical proximity than an ecological one. Thirteen enzymatic loci exhibited moderate to high values compared with microsatellites. This was interpreted as evidence that these allozymes are non-neutral. However, only six loci seemed to be implicated in differentiation between marine and lagoon samples, the causes of selection being unknown for the others. A possible scenario of population dynamics of the sea bass between marine and lagoon habitat is suggested.

Different trajectories of parallel evolution during viral adaptation

The molecular basis of adaptation is a major focus of evolutionary biology, yet the dynamic process of adaptation has been explored only piecemeal. Experimental evolution of two bacteriophage lines under strong selection led to over a dozen nucleotide changes genomewide in each replicate. At least 96 percent of the amino acid substitutions appeared to be adaptive, and half the changes in one line also occurred in the other. However, the order of these changes differed between replicates, and parallel substitutions did not reflect the changes with the largest beneficial effects or indicate a common trajectory of adaptation.

From genes to individuals: developmental genes and the generation of the phenotype

The success of the genetic approach to developmental biology has provided us with a suite of genes that are involved in the regulation of ontogenetic pathways. It is therefore time to ask whether and how such genes might be involved in the generation of adaptive phenotypes. Unfortunately, the current results do not provide a clear answer. Most of the genes that have been studied by developmental biologists affect early embryonic traits with significant effects on the whole organism. These genes are often highly conserved which allows us to do comparative studies even across phyla. However, whether the same genes are also involved in short-term ecological adaptations remains unclear. The suggestion that early acting ontogenetic genes may also affect late phenotypes comes from the genetic analysis of quantitative traits like bristle numbers in Drosophila. A rough mapping of the major loci affecting these traits shows that these loci might correspond to well known early acting genes. On the other hand, there are also many minor effect loci that are as yet uncharacterized. We suggest that these minor loci might correspond to a different class of genes. In comparative studies of randomly drawn cDNAs from Drosophila we find that there is a large group of genes that evolve fast and that are significantly under-represented in normal genetic screens. We speculate that these genes might provide a large, as yet poorly understood, reservoir of genes that might be involved in the evolution of quantitative traits and short-term adaptations.

Variation at the b Mating Type Locus of Ustilago maydis

ABSTRACT Population level diversity at the Ustilago maydis b mating type locus was determined in samples from four Minnesota locations using a combination of plate mating techniques and a polymerase chain reaction (PCR)-based assay. The PCR method allows rapid identification of b types from samples of natural populations and utilizes the hypervariable regions of the b locus that determine mating type specificity. Results demonstrated high levels of b diversity within populations, with one population yielding 17 of the total 18 b types found in the study. Pairwise G(ST) values were in the range of 0.02 to 0.05, and common b mating types were found across broad geographic distances. These data demonstrated that very low levels of differentiation among U. maydis populations occur with respect to b locus variation. Consistent with frequency-dependent selection models, b types were represented at approximately equal frequencies within the entire Minnesota population. However, neutral evolutionary models for patterns of geographic distribution and variation at b cannot be entirely excluded. The importance to agricultural practices of understanding population genetic processes is discussed.

A valine421 to methionine mutation in IS6 of the hscp voltage-gated sodium channel associated with pyrethroid resistance in Heliothis virescens F

Multiple mutations in a locus encoding a voltage-gated sodium channel have been predicted for pyrethroid resistance in insects. Previously we reported a mutation associated with pyrethroid resistance, Leu1029 to His, in domain II transmembrane segment S6 (IIS6) of the Heliothis virescens F. sodium channel (para homologue) hscp locus. Sequence analysis of additional resistance haplotypes 5' to this mutation in the hscp locus has uncovered a G to A transition leading to a Val to Met mutation at amino acid position 421 in IS6 (V421M, numbering from Drosophila para). The V421M mutation is found only in a unique resistant haplotype, but not in two susceptible and a distinct resistant haplotype carrying the L1029H mutation. Implications of this finding in the evolution and mechanisms of pyrethroid resistance are discussed.

A novel mutation L1029H in sodium channel gene hscp associated with pyrethroid resistance for Heliothis virescens (Lepidoptera:Noctuidae)

Sequencing of the sodium channel gene hscp (para-homologous) revealed that a T to A mutation encoding a Leu to His change (L1029H) is associated with pyrethroid resistance in populations of Heliothis virescens F. This mutation occurs at a position homologous to that of the mutation Leu to Phe (L1014F) reportedly associated with knock-down resistance (kdr) to pyrethroids in a house fly strain. Another mutation reportedly associated with super-kdr in house fly was not found in H. virescens samples. In previous work, we found that one haplotype Hpy3 of a marker locus Hpy that lies approximately 5 kb away from the L1029H polymorphism in the gene, increased in frequency among survivors of selection with pyrethroids, presumably through "hitch-hiking" with a resistance-conferring mutation. Not all resistant individuals examined in this study carried the L1029H mutation, but only those that also carried Hpy3 haplotypes. Resistance has been genetically linked to hscp markers for the strain PEG87, but this strain was found not to carry the L1029H mutation. Assuming that the L1029H mutation does in fact confer resistance, this suggests that more than one sodium channel mutation may be contributing to pyrethroid resistance in field populations of H. virescens.