Detection of the signature of natural selection in humans: evidence from the Duffy blood group locus

The Duffy blood group locus, which encodes a chemokine receptor, is characterized by three alleles-FY*A, FY*B, and FY*O. The frequency of the FY*O allele, which corresponds to the absence of Fy antigen on red blood cells, is at or near fixation in most sub-Saharan African populations but is very rare outside Africa. The FST value for the FY*O allele is the highest observed for any allele in humans, providing strong evidence for the action of natural selection at this locus. Homozygosity for the FY*O allele confers complete resistance to vivax malaria, suggesting that this allele has been the target of selection by Plasmodium vivax or some other infectious agent. To characterize the signature of directional selection at this locus, we surveyed DNA sequence variation, both in a 1.9-kb region centered on the FY*O mutation site and in a 1-kb region 5-6 kb away from it, in 17 Italians and in a total of 24 individuals from five sub-Saharan African populations. The level of variation across both regions is two- to threefold lower in the Africans than in the Italians. As a result, the pooled African sample shows a significant departure from the neutral expectation for the number of segregating sites, whereas the Italian sample does not. The FY*O allele occurs on two major haplotypes in three of the five African populations. This finding could be due to recombination, recurrent mutation, population structure, and/or mutation accumulation and drift. Although we are unable to distinguish among these alternative hypotheses, it is likely that the two major haplotypes originated prior to selection on the FY*O mutation.     

The effect of natural selection on estimates of genetic divergence among populations of the Atlantic salmon

The effect of natural selection on the _mMEP-2 * locus on measures of genetic divergence among Atlantic salmon populations was investigated by examining the pattern of change in the level of genetic differentiation (F_ST) averaged over loci when data on the _mMEP-2 * locus were either included or excluded. The level of F_ST among populations at various geographic scales was estimated from allele frequencies at up to four loci (_s AAT-4 *, IDDH-1 *, IDHP-3 *, and _mMEP-2 *). At smaller geographic scales (within river systems or limited geographic regions) levels of variance in _mMEP-2 * allele frequencies were reduced relative to mean levels. At larger geographic scales (across continents or the species range) variation in _mMEP-2 * allele frequencies was greater than mean levels. These results suggest an a priori hypothesis for the effect of selection on the _mMEP-2 * locus which may be applied in future studies on variation in protein coding or other (e.g. mini- and microsatellite) loci in the Atlantic salmon. It is recommended that estimates of gene flow among populations of the Atlantic salmon based on mean F _ST estimates which include data on the _mMEP-2 * locus should be viewed with caution.     

Simulated climate change provokes rapid genetic change in the Mediterranean shrub Fumana thymifolia

Rapid climate change will impose strong directional selection pressures on natural plant populations. Climate-linked genetic variation in natural populations indicates that an evolutionary response is possible. We investigated such a response by comparing individuals subjected to elevated drought and warming treatments with individuals establishing in an unmanipulated climate within the same population. We report that reduction in seedling establishment in response to climate manipulations is nonrandom and results from the selection pressure imposed by artificially warmed and droughted conditions. When compared against control samples, high single-locus genetic divergence occurred in drought and warming treatment samples, with genetic differentiation up to 37 times higher than background (mean neutral locus) genetic differentiation. These loci violate assumptions of selective neutrality, indicating the signature of natural selection by drought. Our results demonstrate that rapid evolution in response to climate change may be widespread in natural populations, based on genetic variation already present within the population.     

Duffy blood group genotypes among African-Brazilian communities of the Amazon region

Duffy blood group genotype was studied in 95 unrelated subjects from four African-Brazilian communities of the Amazon region: Trombetas, Pitimandeua, Curiaú, and Mazag?o Velho. Genotyping was performed using an allele-specific primer polymerase chain reaction technique for determining the three major alleles at FY blood group, and as expected, FY*O allele was the most common one, with frequencies ranging from 56.4% in Mazag?o Velho to 72.2% in Pitimandeua, whereas the FY*O/FY*O genotype was found with frequencies between 32.3% in Mazag?o Velho and 58.8% in Curiaú. Genotype and allele distributions in the four Amazonian communities are consistent with a predominantly African origin with some degree of local differentiation and admixture with people of Caucasian ancestry and/or Amerindians. These results reveal that the impact of the FY*O/FY*O genotype on the transmission and endemicity of the vivax malaria deserves to be investigated in full detail in an attempt to identify the contribution of host biological factors and explain the non-homogeneous prevalence of malaria in the region expressed by its different levels of exposure.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

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.     

Common garden and natural selection experiments support ecotypic differentiation in the Dominican anole (Anolis oculatus)

The theory behind ecotypic differentiation and ecological speciation assumes a predominant role for natural selection working on characteristics with genetic variance, but experimental support for these assumptions is limited. Lesser Antillean anoles show marked ecotypic variation within islands and the potential for ecological speciation. Common garden rearing experiments on the Dominican anole (Anolis oculatus) suggest that the characters showing geographic variation have genetic variance and are not primarily determined by environmental plasticity. Replicated natural selection experiments using large-scale enclosures show that translocated montane samples experience significant (multivariate) directional selection in both wet and dry seasons in both males and females. The targets of selection appear to be spread among the various character systems. An experiment on 12 geographically segregated populations along a coastal xeric-montane rainforest gradient (four replicate enclosures) clearly showed that the magnitude of the directional selection intensity is positively related to the position along this gradient. The results of the common garden and natural selection experiments support the interpretation that the geographic differentiation is primarily driven by natural selection and are compatible with the potential for ecological speciation in this system.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Hitchhiking mapping reveals a candidate genomic region for natural selection in three-spined stickleback chromosome VIII

Identification of genes and genomic regions under directional natural selection has become one of the major goals in evolutionary genetics, but relatively little work to this end has been done by applying hitchhiking mapping to wild populations. Hitchhiking mapping starts from a genome scan using a randomly spaced set of molecular markers followed by a fine-scale analysis in the flanking regions of the candidate regions under selection. We used the hitchhiking mapping approach to narrow down a selective sweep in the genomic region flanking a candidate locus (Stn90) in chromosome VIII in the three-spined stickleback (Gasterosteus aculeatus). Twenty-four microsatellite markers were screened in an approximately 800-kb region around the candidate locus in three marine and four freshwater populations. The patterns of genetic diversity and differentiation in the candidate region were compared to those of a putatively neutral set of markers. The Bayesian FST-test indicated an elevated genetic differentiation, deviating significantly from neutral expectations, at a continuous region of approximately 20 kb upstream from the candidate locus. Furthermore, a method developed for an array of microsatellite markers rejected neutrality in a region of approximately 90 kb flanking the candidate locus supporting the selective sweep hypothesis. Likewise, the genomewide pattern of genetic diversity differed from the candidate region in a bottleneck analysis suggesting that selection, rather than demography, explains the reduced genetic diversity at the candidate interval. The neutrality tests suggest that the selective sweep had occurred mainly in the Lake Pulmanki population, but the results from bottleneck analyses indicate that selection might have operated in other populations as well. These results suggest that the narrow interval around locus Stn90 has likely been under directional selection, but the region contains several predicted genes, each of which can be the actual targets of selection. Understanding of the functional significance of this genomic region in an ecological context will require a more detailed sequence analysis.     

Haplotype structure and evidence for positive selection at the human IL13 locus

Interleukin-13 (IL13) is believed to play an important role in the pathogenesis of atopy and allergic asthma. To better understand genetic variation at the IL13 locus, we resequenced a 5.1-kb genomic region spanning the entire locus and identified 26 single-nucleotide polymorphisms (SNPs) in 74 individuals from three major populations-Chinese, Caucasian, and African. Our survey suggests exceptionally high and significant geographic structure at the IL13 locus between African and outside Africa populations. This unusual pattern suggests that positive selection that acts in some local populations may have played a role on the IL13 locus. In support of this suggestion, we found a significant excess of high frequency-derived SNPs in the Chinese population and Caucasian population, respectively, as expected after a recent episode of positive selection. Further, the unusual haplotype structure indicates that different scenarios of the action of positive selection on the IL13 locus in different populations may exist. In the Caucasian population, the skewed haplotype distribution dominated by one common haplotype supports the hypothesis of simple directional selection. Whereas, in the Chinese population, the two-round hitchhiking hypothesis may explain the skewed haplotype structure with three dominant ones. These findings may provide insight into the likely relative roles of selection and population history in establishing present-day variation at the IL13 locus, and, motivate further studies of this locus as an important candidate in common diseases association studies.     

Strong positive selection and habitat-specific amino acid substitution patterns in MHC from an estuarine fish under intense pollution stress

Population-level studies using the major histocompatibility complex (Mhc) have linked specific alleles with specific diseases, but data requirements are high and the power to detect disease association is low. A novel use of Mhc population surveys involves mapping allelic substitutions onto the inferred structural molecular model to show functional differentiation related to local selective pressures. In the estuarine fish Fundulus heteroclitus, populations experiencing strong differences in antigenic challenges show significant differences in amino acid substitution patterns that are reflected as variation in the structural location of changes between populations. Fish from a population genetically adapted to severe chemical pollution also show novel patterns of DNA substitution at a highly variable Mhc class II B locus including strong signals of positive selection at inferred antigen-binding sites and population-specific signatures of amino acid substitution. Heavily parasitized fish from an extreme PCB-contaminated (U.S. Environmental Protection Agency Superfund) site show enhanced population-specific substitutions in the a-helix portion of the inferred antigen-binding region. In contrast, fish from an unpolluted site show a significantly different pattern focused on the first strand of the B-pleated sheet. Whether Mhc population profile differences represent the direct effects of chemical toxicants or indirect parasite-mediated selection, the result is a composite habitat-specific signature of strong selection and evolution affecting the genetic repertoire of the major histocompatibility complex.     

Genetic evidence for depth- and spatially separated stocks of the deep-water spikey oreo in Australasian waters

Samples of Neocyttus rhomboidalis from eight areas off southern Australia and eastern New Zealand were examined to assess the stock structure of the species. No spatial heterogeneity was detected for mitochondrial DNA variation after restriction analysis by two enzymes detecting polymorphic cut sites. Twenty mtDNA haplotypes were recorded, with an overall haplotype diversity of 0·683. Twenty-eight allozyme loci were examined; the mean heterozygosity per locus was 13·0%. Significant sample heterogeneity was detected at three of twelve polymorphic loci: MPI*, PGM-1 * and, most strikingly (G_ST=0·43), sSOD *. Most of the sSOD * differentiation was associated with depth rather than geographical separation; sSOD * allele frequencies of shallow and deep samples being very different. It is suggested that the correlation between the sSOD * allele frequencies and depth is more likely to reflect natural selection than reproductive isolation. The spatial patterning of the MPI * and PGM-1 * data suggest there are at least three geographic stocks of N. rhomboidalis in Australasian waters: (1) New Zealand, (2) Western Australia and South Australia, (3) west Tasmania to New South Wales and Lord Howe Rise.     

Genetic analysis of population differentiation and adaptation in Leuciscus waleckii

Demographic events and natural selection both influence animal phenotypic and genetic variation; exploring the effects of demography and selection on population divergence is of great significance in evolutionary biology. To uncover the causes behind the patterns of genetic differentiation and adaptation among six populations of Leuciscus waleckii from Dali Basin (two populations, alkaline vs. freshwater) and Amur Basin (four populations, freshwater rivers vs. alkaline lake), a set of 21 unlinked polymorphic microsatellite markers and two mitochondrial DNA sequences (Cytb and D-loop) were applied to examine whether populations from different environments or habitats have distinct genetic differentiation and whether alkalinity is the major factor that caused population divergence. Bayesian analysis and principal component analysis as well as haplotype network analysis showed that these populations are primarily divided into two groups, which are congruent with geographic separation but not inconsistent with the habitat environment (alkalinity). Using three different approaches, outlier detection indicated that one locus, HLJYL017, may be under directional selection and involved in local adaptation processes. Overall, this study suggested that demographic events and selection of local environmental conditions including of alkalinity are jointly responsible for population divergence. These findings constitute an important step towards the understanding of the genetic basis of differentiation and adaptation, as well as towards the conservation of L. waleckii.     

Geographical delimitation of a partial selective sweep in African Drosophila melanogaster

Positive selection leaves characteristic footprints on DNA variation but detecting such patterns is challenging as the age, the intensity and the mode of selection as well as demography and evolutionary parameters (mutation and recombination rates) all play roles and these are difficult to disentangle. We recorded nucleotide variation in a sample of isogenic chromosomes from a western African population of Drosophila melanogaster at a locus (Fbp2) for which a partial selective sweep had previously been reported. We compared this locus to four other genes from the same chromosomes and from a European and an East African population. Then, we assessed Fbp2 variation in a sample of 370 chromosomes covering a comprehensive geographic sampling of 16 African localities. The signature of selection was tested while accounting for the demographic history of the populations. We found a significant signal of selection in two West African localities including Ivory Coast. Variation at Fpb2 would thus represent a case of an ongoing selective sweep in the range of this species. A weaker, nonsignificant, signal of selection was, however, apparent in some other populations, thus leaving open several possibilities: (i) the selective sweep originated in Ivory Coast and has spread to the rest of the continent; (ii) several African populations report the signature of a selective event having occurred in an ancestral population; (iii) this genome region is subject to independent selective events in African populations; and (iv) A neutral scenario with population subdivision and local bottleneck cannot be fully excluded to explain the molecular patterns observed in some populations.     

Genetic variation in the acorn barnacle from allozymes to population genomics

Understanding the patterns of genetic variation within and among populations is a central problem in population and evolutionary genetics. We examine this question in the acorn barnacle, Semibalanus balanoides, in which the allozyme loci Mpi and Gpi have been implicated in balancing selection due to varying selective pressures at different spatial scales. We review the patterns of genetic variation at the Mpi locus, compare this to levels of population differentiation at mtDNA and microsatellites, and place these data in the context of genome-wide variation from high-throughput sequencing of population samples spanning the North Atlantic. Despite considerable geographic variation in the patterns of selection at the Mpi allozyme, this locus shows rather low levels of population differentiation at ecological and trans-oceanic scales (F(ST)?～?5%). Pooled population sequencing was performed on samples from Rhode Island (RI), Maine (ME), and Southwold, England (UK). Analysis of more than 650 million reads identified approximately 335,000 high-quality SNPs in 19 million base pairs of the S. balanoides genome. Much variation is shared across the Atlantic, but there are significant examples of strong population differentiation among samples from RI, ME, and UK. An F(ST) outlier screen of more than 22,000 contigs provided a genome-wide context for interpretation of earlier studies on allozymes, mtDNA, and microsatellites. F(ST) values for allozymes, mtDNA and microsatellites are close to the genome-wide average for random SNPs, with the exception of the trans-Atlantic F(ST) for mtDNA. The majority of F(ST) outliers were unique between individual pairs of populations, but some genes show shared patterns of excess differentiation. These data indicate that gene flow is high, that selection is strong on a subset of genes, and that a variety of genes are experiencing diversifying selection at large spatial scales. This survey of polymorphism in S. balanoides provides a number of genomic tools that promise to make this a powerful model for ecological genomics of the rocky intertidal.     

Population structure in Arabidopsis lyrata: evidence for divergent selection on trichome production

Leaf trichomes may serve several biological functions including protection against herbivores, drought, and UV radiation; and their adaptive value can be expected to vary among environments. The perennial, self-incompatible herb Arabidopsis lyrata is polymorphic for trichome production, and occurs in a glabrous and a trichome-producing form. Controlled crosses indicate that the polymorphism is governed by a single gene, with trichome production being dominant. We examined the hypothesis that trichome production is subject to divergent selection (i.e., directional selection favoring different phenotypes in different populations) by comparing patterns of variation at the locus coding for glabrousness and at eight putatively neutral isozyme loci in Swedish populations of A. lyrata. The genetic diversity (He) and allele number at isozyme loci tended to increase with population size and decreased with latitude of origin, whereas genetic diversity at the locus coding for glabrousness did not vary with population size and increased with latitude of origin. The degree of genetic differentiation at the glabrousness locus was much higher than that at isozyme loci. Genetic identity at isozyme loci was negatively related to geographic distance, suggesting isolation by distance. In contrast, there was no significant correlation between genetic identity at the glabrousness locus and at isozyme loci. The results are consistent with the hypothesis that divergent selection contributes to population differentiation in trichome production in A. lyrata.     

Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon

Balancing selection maintains high levels of polymorphism and heterozygosity in genes of the MHC (major histocompatibility complex) of vertebrate organisms, and promotes long evolutionary persistence of individual alleles and strongly differentiated allelic lineages. In this study, genetic variation at the MHC class II DAB-beta1 locus was examined in 31 populations of sockeye salmon (Oncorhynchus nerka) inhabiting the Fraser River drainage of British Columbia, Canada. Twenty-five percent of variation at the locus was partitioned among sockeye populations, as compared with 5% at neutral genetic markers. Geographic heterogeneity of balancing selection was detected among four regions in the Fraser River drainage and among lake systems within regions. High levels of beta1 allelic diversity and heterozygosity, as well as distributions of alleles and allelic lineages that were more even than expected for a neutral locus, indicated the presence of balancing selection in populations throughout much of the interior Fraser drainage. However, proximate populations in the upper Fraser region, and four of six populations from the lower Fraser drainage, exhibited much lower levels of genetic diversity and had beta1 allele frequency distributions in conformance with those expected for a neutral locus, or a locus under directional selection. Pair-wise FST values for beta1 averaged 0.19 and tended to exceed the corresponding values estimated for neutral loci at all levels of population structure, although they were lower among populations experiencing balancing selection than among other populations. The apparent heterogeneity in selection resulted in strong genetic differentiation between geographically proximate populations with and without detectable levels of balancing selection, in stark contrast to observations at neutral loci. The strong partitioning and complex structure of beta1 diversity within and among sockeye populations on a small geographic scale illustrates the value of incorporating adaptive variation into conservation planning for the species.     

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.     

Combining sperm typing and linkage disequilibrium analyses reveals differences in selective pressures or recombination rates across human populations

A previous polymorphism survey of the type 2 diabetes gene CAPN10 identified a segment showing an excess of polymorphism levels in all population samples, coinciding with localized breakdown of linkage disequilibrium (LD) in a sample of Hausa from Cameroon, but not in non-African samples. This raised the possibility that a recombination hotspot is present in all populations and we had insufficient power to detect it in the non-African data. To test this possibility, we estimated the crossover rate by sperm typing in five non-African men; these estimates were consistent with the LD decay in the non-African, but not in the Hausa data. Moreover, resequencing the orthologous region in a sample of Western chimpanzees did not show either an excess of polymorphism level or rapid LD decay, suggesting that the processes underlying the patterns observed in humans operated only on the human lineage. These results suggest that a hotspot of recombination has recently arisen in humans and has reached higher frequency in the Hausa than in non-Africans, or that there is no elevation in crossover rate in any human population, and the observed variation results from long-standing balancing selection.     

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub‐Saharan Africa, has been the subject of long‐term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest–savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.