Does fragmentation of wetlands affect gene flow in sympatric Acrocephalus warblers with different migration strategies?

Wetlands are naturally patchy habitats, but patchiness has been accentuated by the extensive wetlands loss due to human activities. In such a fragmented habitat, dispersal ability is especially important to maintain gene flow between populations. Here we studied population structure, genetic diversity and demographic history of Iberian and North African populations of two wetland passerines, the Eurasian reed warbler Acrocephalus scirpaceus and the moustached warbler Acrocephalus melanopogon. These species are closely related and sympatric in our study sites, but the reed warbler is a widespread long‐distance migrant while the moustached warbler's breeding range is patchier and it is resident or migrates over short distances. Using microsatellite and mtDNA data, we found higher differentiation in moustached than in reed warblers, indicating higher dispersal capability of the latter species. Our results also suggest that the sea limits dispersal in the moustached warbler. However, we found evidence of gene flow between the study sites in both species, indicating a capability to compensate for habitat fragmentation. In most cases, the gene flow was restricted, possibly because of the large distances between study sites (from ca 290 to 960 km) or breeding site fidelity. The reed warbler had higher haplotype diversity, likely due to dispersal from different populations, past admixture event and a larger population size. We found also signs of postglacial population growth for both species and evidence of a recent colonization or re‐colonization of the Mallorca Island by the moustached warbler.     

Molecular evidence of a reed warbler × great reed warbler hybrid ( Acrocephalus scirpaceus × A. arundinaceus) in Belgium

We report on the second case of a reed warbler × great reed warbler hybrid (Acrocephalus scirpaceus and A. arundinaceus). The bird was captured during a standardised ringing session in Belgium in autumn 1999, and fell between the parental species in all measurements. Molecular analyses of two microsatellite loci verified that the bird was a female that had a reed warbler father and a great reed warbler mother.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

Morphometric differentiation across House Sparrow Passer domesticus populations in Finland in comparison with the neutral expectation for divergence

To understand the biology of organisms it is important to take into account the evolutionary forces that have acted on their constituent populations. Neutral genetic variation is often assumed to reflect variation in quantitative traits under selection, though with even low neutral divergence there can be substantial differentiation in quantitative genetic variation associated with locally adapted phenotypes. To study the relative roles of natural selection and genetic drift in shaping phenotypic variation, the levels of quantitative divergence based on phenotypes (P_ST) and neutral genetic divergence (F_ST) can be compared. Such a comparison was made between 10 populations of Finnish House Sparrows (n = 238 individuals) collected in 2009 across the whole country. Phenotypic variation in tarsus‐length, wing‐length, bill‐depth, bill‐length and body mass were considered and 13 polymorphic microsatellite loci were analysed to quantify neutral genetic variation. Calculations of P_ST were based on Markov‐Chain Monte Carlo Bayesian estimates of phenotypic variances across and within populations. The robustness of the conclusions of the P_ST–F_ST comparison was evaluated by varying the proportion of variation due to additive genetic effects within and across populations. Our results suggest that body mass is under directional selection, whereas the divergence in other traits does not differ from neutral expectations. These findings suggest candidate traits for considering gene‐based studies of local adaptation. The recognition of locally adapted populations may be of value in the conservation of this declining species.     

Genetic differentiation between collections of hatchery and wild masu salmon (<Emphasis Type="Italic">Oncorhynchus masou</Emphasis>) inferred from mitochondrial and microsatellite DNA analyses

There has been very little effort to understand genetic divergence between wild and hatchery populations of masu salmon (Oncorhynchus masou). In this study, we used mitochondrial (mt) NADH dehydrogenase subunit 5 gene (ND5) and six polymorphic nuclear microsatellite DNA loci to compare the genetic variability in three hatchery broodstocks of masu salmon with the variability in eight putative wild masu populations sampled in five rivers including one known source river for the hatchery broodstocks. Both ND5 and microsatellites showed no significant genetic divergence (based on F_ST estimates) between four annual collections from the source river population, suggesting no change in genetic diversity over this time period. The F_ST estimates, an analysis of molecular variance (AMOVA), and a neighbor-joining tree using both DNA markers suggested significant differentiation between the three hatchery and all eight putative wild populations. We conclude that genetic diversity of hatchery populations are low relative to putative wild populations of masu salmon, and we discuss the implications for conservation and fisheries management in Hokkaido.     

Spatial Genetic Structure in Natural Populations of Phragmites australis in a Mosaic of Saline Habitats in the Yellow River Delta, China

Determination of spatial genetic structure (SGS) in natural populations is important for both theoretical aspects of evolutionary genetics and their application in species conservation and ecological restoration. In this study, we examined genetic diversity within and among the natural populations of a cosmopolitan grass Phragmites australis (common reed) in the Yellow River Delta (YRD), China, where a mosaic of habitat patches varying in soil salinity was detected. We demonstrated that, despite their close geographic proximity, the common reed populations in the YRD significantly diverged at six microsatellite loci, exhibiting a strong association of genetic variation with habitat heterogeneity. Genetic distances among populations were best explained as a function of environmental difference, rather than geographical distance. Although the level of genetic divergence among populations was relatively low (F’_ST = 0.073), weak but significant genetic differentiation, as well as the concordance between ecological and genetic landscapes, suggests spatial structuring of genotypes in relation to patchy habitats. These findings not only provided insights into the population dynamics of common reed in changing environments, but also demonstrated the feasibility of using habitat patches in a mosaic landscape as test systems to identify appropriate genetic sources for ecological restoration.     

Calculations of population differentiation based on GST and D: forget GST but not all of statistics!

G _ST‐values and its relatives (F_ST) belong to the most used parameters to define genetic differences between populations. Originally, they were developed for allozymes with very low number of alleles. Using highly polymorphic microsatellite markers it was often puzzling that G_ST‐values were very low but statistically significant. In their papers, Jost (2008) and Hedrick (2005) explained that G_ST‐values do not show genetic differentiation, and Jost suggested calculating D‐values instead. Theoretical mathematical considerations are often difficult to follow; therefore, we chose an applied approach comparing two artificial populations with different number of alleles at equal frequencies and known genetic divergence. Our results show that even for more than one allele per population G_ST‐values do not calculate population differentiation correctly; in contrast, D‐values do reflect the genetic differentiation indicating that data based on G_ST‐values need to be re‐evaluated. In our approach, statistical evaluations remained similar. We provide information about the impact of different sample sizes on D‐values in relation to number of alleles and genetic divergence.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Microgeographic structure of Anopheles gambiae in western Kenya based on mtDNA and microsatellite loci

The population genetic structure of the Anopheles gambiae in western Kenya was studied using length variation at five microsatellite loci and sequence variation in a 648-nt mtDNA fragment. Mosquitoes were collected from houses in villages spanning up to 50 km distance, The following questions were answered, (i) Are mosquitoes in a house more related genetically to each other than mosquitoes between houses? (ii) What degree of genetic differentiation occurs on these geographical scales? (iii) How consistent are the results obtained with both types of genetic markers? At the house level, no differentiation was detected by F_ST and R_ST, and the band sharing index test revealed no significant associations of alleles across loci. Likewise, indices of kinship based on mtDNA haplotypes in houses were even lower than in the pooled sample. Therefore, the hypothesis that mosquitoes in a house are more related genetically was rejected. At increasing geographical scales, microsatellite allele distributions were similar among all population samples and no subdivision of the gene pool was detected by F_ST or R_ST. Likewise, estimates of haplotype divergence of mtDNA between populations were not higher than the within population estimates, and mtDNA-based F_ST values were not significantly different from zero. That sequence variation in mtDNA provided matching results with microsatellite loci (while high genetic variation was observed in all loci), suggested that this pattern represents the whole genome. The minimum area associated with a deme of A. gambiae in western Kenya is therefore larger than 50 km in diameter.     

Extensive gene flow along the urban–rural gradient in a migratory colonial bird

Urban colonization by wildlife involves a combination of several different mechanisms, including phenotype or genotype sorting, phenotypic plasticity and microevolutionary adaptation. Combination of these processes can produce a rapid phenotypic, but also genetic divergence of urban versus rural populations. Here, we examined the pattern of genetic differentiation between urban and rural populations of a colonial migratory bird, the black‐headed gull Chroicocephalus ridibundus. To this end, we sampled ca 170 individuals from six (two urban and four rural) colonies in northern Poland, and genotyped them at ten microsatellite loci. Our analysis provided evidence for negligible genetic divergence of urban and rural colonies, as assessed with fixation index F_ST and Nei's unbiased genetic distance D (mean pairwise urban‐rural comparisons: F_ST = 0.003 ± 0.001 [SE] and D = 0.012 ± 0.006 [SE]). Bayesian clustering methods provided support for homogeneous genetic structure across all urban and rural populations. Also, we found no support for reduced allelic diversity in urban versus rural colonies. These results stand in a stark contrast to the previous findings on the genetic consequences of urbanization in birds. We hypothesize that this pattern could possibly be attributed to the important life‐history characters of the black‐headed gull, including coloniality, migratoriness, and high dispersal propensity. Our study provides a novel insight into the urban landscape genetics, underlining large variation in the mechanisms of urban colonization and its genetic consequences in wild animal populations.     

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of F_ST and to evaluate whether highly variable microsatellite loci resulted in deflated F_ST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (F_CT = 0.010) and within (F_ST = 0.014–0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (P_CT > F_CT). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (P_ST > F_ST), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated F_ST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for P_ST–F_ST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

Morphological and microsatellite diversity associated with ecological factors in natural populations of Medicago laciniata Mill. (Fabaceae)

Genetic variability in 10 natural Tunisian populations of Medicago laciniata were analysed using 19 quantitative traits and 12 polymorphic microsatellite loci. A large degree of genetic variability within-populations and among-populations was detected for both quantitative characters and molecular markers. High genetic differentiation among populations for quantitative traits was seen, with Q _ST = 0.47, and F _ST = 0.47 for microsatellite markers. Several quantitative traits displayed no statistical difference in the levels of Q _ST and F _ST . Further, significant correlations between quantitative traits and eco-geographical factors suggest that divergence in the traits among populations may track environmental differences. There was no significant correlation between genetic variability at quantitative traits and microsatellite markers within populations. The site-of-origin of eco-geographical factors explain between 18.13% and 23.40% of genetic variance among populations at quantitative traits and microsatellite markers, respectively. The environmental factors that most influence variation in measured traits among populations are assimilated phosphorus (P₂0₅) and mean annual rainfall, followed by climate and soil texture, altitude and organic matter. Significant associations between eco-geographical factors and gene diversity, H _e , were established in five-microsatellite loci suggesting that these simple sequence repeats (SSRs) are not necessarily biologically neutral.     

Comparison of genetic differentiation at marker loci and quantitative traits

The comparison of the degree of differentiation in neutral marker loci and genes coding quantitative traits with standardized and equivalent measures of genetic differentiation (F_ST and Q_ST, respectively) can provide insights into two important but seldom explored questions in evolutionary genetics: (i) what is the relative importance of random genetic drift and directional natural selection as causes of population differentiation in quantitative traits, and (ii) does the degree of divergence in neutral marker loci predict the degree of divergence in genes coding quantitative traits? Examination of data from 18 independent studies of plants and animals using both standard statistical and meta‐analytical methods revealed a number of interesting points. First, the degree of differentiation in quantitative traits (Q_ST) typically exceeds that observed in neutral marker genes (F_ST), suggesting a prominent role for natural selection in accounting for patterns of quantitative trait differentiation among contemporary populations. Second, the F_ST – Q_ST difference is more pronounced for allozyme markers and morphological traits, than for other kinds of molecular markers and life‐history traits. Third, very few studies reveal situations were Q_ST < F_ST, suggesting that selection pressures, and hence optimal phenotypes, in different populations of the same species are unlikely to be often similar. Fourth, there is a strong correlation between Q_ST and F_ST indices across the different studies for allozyme (r=0.81), microsatellite (r=0.87) and combined (r=0.75) marker data, suggesting that the degree of genetic differentiation in neutral marker loci is closely predictive of the degree of differentiation in loci coding quantitative traits. However, these interpretations are subject to a number of assumptions about the data and methods used to derive the estimates of population differentiation in the two sets of traits.     

Historical vicariance and male-mediated gene flow in the toad-headed lizards Phrynocephalus przewalskii

Using mitochondrial and microsatellite DNA data and a population genetic approach, we tested male‐mediated gene flow in the toad‐headed lizards Phrynocephalus przewalskii. The mitochondrial DNA (ND2 gene), on the one hand, revealed two major lineages and a strong population genetic structure (F_ST = 0.692; F_ST′ = 0.995). The pairwise differences between the two lineages ranged from 2.1% to 6.4% and the geographical division of the two lineages coincided with a mountain chain consisting of the Helan and Yin Mountains, suggesting a historical vicariant pattern. On the other hand, the nuclear microsatellite DNA revealed a significant but small population genetic structure (F_ST = 0.017; F_ST′ = 0.372). The pairwise F_ST among the nine populations examined with seven microsatellite DNA loci ranged from 0.0062 to 0.0266; the assignment test failed to detect any naturally occurring population clusters. Furthermore, the populations demonstrated a weak isolation by distance and a northeast to southwest clinal variation, rather than a vicariant pattern. A historical vicariant event followed by male‐mediated gene flow appears to be the best explanation for the data. Approximately 2–5 Ma, climatic change may have created an uninhabitable zone along the Helan‐Yin mountain chain and initiated the divergence between the two mitochondrial lineages. With further climatic changes, males were able to disperse across the mountain chain, causing sufficient gene flow that eventually erased the vicariant pattern and drastically reduced the population genetic structure, while females remained philopatric and maintained the mitochondrial DNA (mtDNA) divergence. Although polygyny mating system and female philopatry may partially contribute to the reduced movement of females, other hypotheses, such as female intrasexual aggression, should also be explored.     

Genetic divergence in nuclear genomes between populations of <Emphasis Type="Italic">Fagus crenata</Emphasis> along the Japan Sea and Pacific sides of Japan

Genetic diversity and structure in Fagus crenata were studied by analyzing 14 nuclear microsatellite loci in 23 populations distributed throughout the species’ range. Although population differentiation was very low (F _ST = 0.027; R _ST = 0.041), both neighbor-joining tree and Bayesian clustering analyses provided clear evidence of genetic divergence between populations along the Japan Sea (Japan Sea lineage) and Pacific (Pacific lineage) sides of Japan, indicating that physical barriers to migration and gene flow, notably the mountain ranges separating the populations along the Japan Sea and Pacific sides, have promoted genetic divergence between these populations. The two lineages of the nuclear genome are generally consistent with those of the chloroplast genome detected in a previous study, with several discrepancies between the two genomes. Within-population genetic diversity was generally very high (average H _E = 0.839), but decreased in a clinal fashion from southwest to northeast, largely among populations of the Japan Sea lineage. This geographical gradient may have resulted from the late-glacial and postglacial recolonization to the northeast, which led to a loss of within-population genetic diversity due to cumulative founder effects.     

Genetic structure of the Painted Bunting and its implications for conservation of migratory populations

The Painted Bunting Passerina ciris is a Neotropical songbird which breeds primarily in the USA during the summer and migrates to Mexico, Central America, southern Florida and the Caribbean over the winter. Male Painted Buntings are brightly coloured, which makes them highly sought after as pets, particularly in Mexico, Central America and Europe. We used short sequence repeats (microsatellite DNA) to investigate the population genetic structure of the Painted Bunting and its implications in conservation management of migratory populations. We found a detectable level of population differentiation as revealed by pairwise F_ST and R_ST comparisons and Bayesian clustering analyses, with strong support for differentiation between eastern and western Painted Buntings (e.g. Oklahoma and Georgia F_ST = 0.1; P = 0.005; R_ST = 0.18; P = 0.04) in accordance with previous mitochondrial DNA analysis. We recovered additional support for two sub‐groups within the western clade. While linking migrant songbirds captured outside of the USA to their breeding populations remains a challenge, we show that natural levels of population genetic differentiation can be detected via microsatellite DNA markers and exploited in migratory connectivity studies. We also demonstrate the potential utility of our low‐cost markers for population identification of birds recovered from the pet trade by screening a small subset of samples (n = 5) collected as part of wildlife tracking. We discuss the implications of our results for future efforts to understand patterns of population decline in Painted Buntings more generally, as well as how we might expand this methodology to combat illegal pet‐trade activity in this and other songbird species.     

Genetic structure of Polytrichum formosum in relation to the breeding system as revealed by microsatellites

Microsatellite variation was determined for three Danish and three Dutch populations of the haploid moss species Polytrichum formosum to gain insight into the relative importance of sexual vs. asexual reproduction for the amount and structure of genetic variation. In general, low levels of microsatellite variation were observed within this species. Even when estimated for polymorphic loci only, the levels of microsatellite variability (P=90.6, A=4.3 and H_S=0.468) within populations were on average lower than those reported for most other plant species. In contrast, genotypic diversity was high within each of the examined populations, indicating that sexual reproduction is a very important determinant of the genetic structure of P. formosum within populations. In agreement with previous findings for allozyme data, no significant genetic differentiation (F_ST=0.028, R_ST=0.015) was observed neither between populations nor between regions approximately 450 km apart (Denmark vs. the Netherlands). These low levels of population differentiation observed for both types of genetic markers are probably best explained by a high level of effective spore dispersal (gene flow) between populations. Therefore, also on a large geographical scale sexual reproduction is the most important determinant of the genetic structure of P. formosum, despite the high potential to reproduce clonally.     

Population structure of bottlenose dolphins (Tursiops aduncus) impacted by bycatch along the east coast of South Africa

Shark nets placed along the coast of KwaZulu-Natal to protect bathers from shark attacks result in an incidental by-catch of bottlenose dolphins (Tursiops aduncus) at twice the level suggested by the IWC as the maximum sustainable capture rate for a cetacean population. Observational data reported the presence of at least two putative populations, one a seasonal migratory population moving into the KwaZulu-Natal area following the movement of sardines, and the other a coastal, apparently resident population. Some observations indicated possible further subdivision of the coastal population. We analysed nine microsatellite loci and 599 bps of the mitochondrial control region and found small but significant differentiation within the putative coastal population (F_ST = 0.022) separated north and south of Ifafa. There was weaker evidence for differentiation between the southern coastal and migratory populations (F_ST = 0.012). Genetic diversity was low for both types of markers in all populations, and the pattern of mtDNA variation was consistent with a founder event. Considering the high capture rate, the low level of genetic diversity, and the evidence of possible population differentiation along the coast, we suggest that particular care should be taken in managing the take from shark nets especially in the north coastal area of KZN where the population showed the higher degree of differentiation.     

Low levels of quantitative and molecular genetic differentiation among natural populations of Medicago ciliaris Kroch. (Fabaceae) of different Tunisian eco-geographical origin

In this paper, we analyze the genetic variability in four Tunisian natural populations of Medicago ciliaris using 19 quantitative traits and six polymorphic microsatellite loci. We investigated the amplification transferability of 30 microsatellites developed in the model legume M. truncatula to M. ciliaris. Results revealed that about 56.66% of analyzed markers are valuable genetic markers for M. ciliaris. The most genetic diversity at quantitative traits and microsatellite loci was found to occur within populations (>80%). Low differentiations among populations at quantitative traits Q _ST  = 0.146 and molecular markers F _ST  = 0.18 were found. The majority of measured traits exhibited no significant difference in the level of Q _ST and F _ST . Furthermore, significant correlations established between these traits and eco-geographical factors suggested that natural selection should be invoked to explain the level of phenotypic divergence among populations rather than drift. There was no significant correlation between population differentiation at quantitative traits and molecular markers. Significant spatial genetic structure consistent with models of isolation by distance was detected within all studied populations. The site-of-origin environmental factors explain about 9.07% of total phenotypic genetic variation among populations. The eco-geographical factors that influence more the variation of measured traits among populations are the soil texture and altitude. Nevertheless, there were no consistent pattern of associations between gene diversity (He) and environmental factors.