Skip to main content

Genetic diversity and population structure of the Tibetan poplar (Populus szechuanica var. tibetica) along an altitude gradient



The Tibetan poplar (Populus szechuanica var. tibetica Schneid), which is distributed at altitudes of 2,000-4,500 m above sea level, is an ecologically important species of the Qinghai-Tibet Plateau and adjacent areas. However, the genetic adaptations responsible for its ability to cope with the harsh environment remain unknown.


In this study, a total of 24 expressed sequence tag microsatellite (EST-SSR) markers were used to evaluate the genetic diversity and population structure of Tibetan poplars along an altitude gradient. The 172 individuals were of genotypes from low-, medium- and high-altitude populations, and 126 alleles were identified. The expected heterozygosity (HE) value ranged from 0.475 to 0.488 with the highest value found in low-altitude populations and the lowest in high-altitude populations. Genetic variation was low among populations, indicating a limited influence of altitude on microsatellite variation. Low genetic differentiation and high levels of gene flow were detected both between and within the populations along the altitude gradient. An analysis of molecular variance (AMOVA) showed that 6.38% of the total molecular variance was attributed to diversity between populations, while 93.62% variance was associated with differences within populations. There was no clear correlation between genetic variation and altitude, and a Mantel test between genetic distance and altitude resulted in a coefficient of association of r = 0.001, indicating virtually no correlation.


Microsatellite genotyping results showing genetic diversity and low differentiation suggest that extensive gene flow may have counteracted local adaptations imposed by differences in altitude. The genetic analyses carried out in this study provide new insight for conservation and optimization of future arboriculture.


Altitude gradients represent one of the most useful natural environments to investigate ecological and evolutionary responses of biota to geophysical influences [1]. For species from habitats which cover different altitudes, differences in their spatial population structure could be due to restricted gene movement, as a result of non-random mating or geographic barriers [2, 3]. Outliers of species found at the boundaries of their distribution zones could be subject to limited gene flow, a small population size and founder effects, all of which lead to a decrease in genetic diversity and an increase in population differentiation [4]. For species living in mountainous areas, altitude changes represent a series of physical factors that can result in the establishment of different populations and species. These factors form barriers, which influence genetic diversity and population structure [57], and include factors such as rainfall [8] and temperature [9]. There is no general rule to summarize the relationship between genetic diversity and altitude; for trees on mountainsides, the pattern of genetic diversity along the altitude gradient is divided into four groups. (1) Populations at an intermediate altitude have greater diversity than populations at lower and higher altitudes, due to local adaptation and milder environmental conditions [10, 11]. (2) Populations at higher altitudes have greater diversity than those at lower altitudes if the higher altitude conditions are similar to their home sites, representing higher fitness [12]. (3) Populations at lower altitudes have greater diversity than those at higher altitudes, as higher altitudes impede growth and the expanding of species countering the bottleneck leaded to decrease of genetic diversity [13]. (4) Populations show no differences in diversity at differing altitudes [14], the pattern may be due to that the sampling area was part of main distribution area, limited number of populations sampled along the gradient may cause the failure to detect altitude-related trends. On the other hand, if the sampled population was large enough, extensive gene flow and other factors also could lead to the similar pattern..

The Qinghai-Tibetan Plateau (QTP) is the highest and largest plateau in the world, with a mean altitude of 4 000 m above sea level, and an area of 2.5 × 106 km2. In recent years, the QTP has become a hotspot for plant phylogeographical studies [15, 16], focusing mainly on the population dynamics that took place during the Quaternary (reviewed in Qiu et al.) [17]. However, genetic variation patterns along altitudinal gradients of the QTP remain unclear.

The Tibetan poplar belongs to Populus sect. Tacamahaca in the genus Populus and is an ecologically important species, mainly distributed in Sichuan and Tibet at altitudes from 2 000 to 4 500 m [18]. Recent studies have focused mainly on the phylogenic and physiological mechanisms responsible for its resistance to the harsh environment where the lowest temperature is -30°C and the annual average temperature is between 4°C to 12°C [19]. However, there is a pressing need to understand the genetic diversity along altitude gradients. In this paper we investigated the genetic variation of the Tibetan poplar along an altitude gradient using microsatellite genotyping. The specific objectives were: (1) to understand the genetic variation and differentiation within and between populations, and (2) to detect any influence of altitude gradients on genetic diversity.

In this study, a total of 24 EST- SSR loci based on Populus euphratica transcriptome [20] were used to analyze the genetic diversity and population structure of Tibetan poplar populations at different altitudes in the Sejila mountain area. The objectives were to provide a complete picture of the genetic diversity of Tibetan poplar populations at different altitudes in the Sejila mountain area, and to identify a relationship between genetic variation and differences in altitude.

Materials and methods

Sampling strategy and DNA extraction

We collected leaves from 64, 34 and 74 individuals from high-, medium- and low-altitude populations, respectively (Figure 1, Table 1). Our sampling scheme was to divide the distribution areas of the Tibetan poplar in the Sejila mountains (in southeastern Tibet) into three altitude-gradient groups (high, medium, and low), even though the trees are distributed continuously throughout the area. We selected individuals at a minimum of 30 m apart to prevent selection of clones. The leaf was rapidly dehydrated using silica gel beads. Total genomic DNA was extracted from approximately 0.5 g of silica-dried leaf using a modified version of the cetyltrimethyl ammonium bromide method [21]. The quality and concentration of the extracted DNA were determined by 1% agarose gel electrophoresis and ultraviolet spectrophotometry. The DNA samples were diluted to 5-10 ng/μL for use as the template for polymerase chain reaction (PCR) amplification.

Figure 1
figure 1

Populus szechuanica population locations. In the Sejila mountain area, we selected three different altitude gradients to select samples.

Table 1 Locations of Populus szechuanica populations.

Primer selection

113 EST-SSR primer pairs based on the Populus euphratica transcriptome [20] were developed and tested for suitability in the Tibetan poplar. DNA extracted from four Tibetan poplar individuals was amplified, and the amplicons were sequenced to confirm the existence of and enumerate repeat motifs. DNA from eight individuals was used to test for polymorphisms of the successfully amplified primers. SSRs were selected if they had at least three alleles and exhibited robust amplification.

SSR amplification

After screening, 24 primer pairs were selected for the PCR analysis. The forward primer of each pair was tagged with a section of the universal M13 sequence (5′-TGTAAAACGACGGCCAGT-3′) during synthesis. Each 10-μL PCR mixture contained 1× Taq buffer, 0.2 µM dNTPs, 10-20 ng template DNA, 1.6 pmol reverse primer, 1.6 pmol fluorescently labeled M13 primer, 0.4 pmol forward primer and 1 U Taq polymerase (BioMed). PCR amplification was performed using a Biometra thermocycler (Biometra, Goettingen, Germany) under the following conditions: 94°C for 5 min; 30 cycles of 94°C for 30 s, annealing at 56°C for 45 s and elongation at 72°C for 45 s; 8 cycles of 94°C for 30 s, annealing at 53°C for 45 s, elongation at 72°C for 45 s; and a final extension at 72°C for 10 min. The PCR products were separated by capillary electrophoresis using an ABI 3730xl DNA Analyzer (Applied Biosystems, Foster City, CA, USA) after confirmation of amplification on a 1.5% agarose gel. Approximately 0.5 μL of the PCR products obtained using each of the four fluorescently labeled primers was then combined. The products were separated using an ABI 3730xl DNA Analyzer with GeneScan-500 LIZ as an internal marker (Applied Biosystems). The amplicon fragments were sized using GeneMarker version 1.75 (Soft Genetics LLC, State College, PA, USA).

Data analysis

The FLEXIBIN software was used for automated binning of the raw molecular data[22], and the Excel Microsatellite Toolkit [23] was used to convert the size data into a format suitable for further analysis. Genetic parameters were estimated under the hypothesis that all the loci were neutral, thereby presenting a true picture of the natural genetic structure affected by neutral forces such as genetic drift and gene flow, etc. There are several methods of investigating whether a particular locus has been under selection pressure. We performed the FST outlier test using LOSITAN [24] to identify candidate SSR loci possibly under selection pressure [25]. After removal of outlier loci, the remaining data were used to estimate the genetic diversity of the population. Genetic diversity parameters used included: number of alleles (Na); observed heterozygosity (HO); expected heterozygosity (HE) within a subpopulation; Wright's fixation indices for within-subpopulation (FIS) and in the total population (FIT); and pair-wise differentiation among subpopulations (FST), according to Weir & Cockerham [26]. F IS measures the deviation from the Hardy-Weinberg equilibrium (HWE) of genotype frequencies in sub-populations, whereas FIT measures the deviation from HWE in the total population. The values of FIT and FIS can be negative, whereas FST is always a positive value. The Shannon's diversity index was conducted using Nei's model, along with the expected heterozygosity [27]. Gene flow (Nm) was calculated to ascertain the conditions of gene communication among populations, and was estimated as follows: Nm = (1- F ST )/4 F ST [28]. Summary statistics were calculated using POPGENE version 1.32 [29]. Inter- and intra-population differentiation was determined by AMOVA analysis using the GenAlEx software version 6.41 [30]. Clustering, based on a Bayesian model which assumed that all the individuals were from K real populations (where K may be unknown), each of which is characterized by a set of allele frequencies at each locus, the method attempts to assign individuals to populations on the basis of their genotypes, while simultaneously estimating population allele frequencies. The model was used to evaluate the genetic structures of the Tibetan poplar populations using STRUCTURE in its extended version 2.3.3 [31, 32]. STRUCTURE is based on a model-based clustering algorithm that applies a Bayesian framework and the Markov chain Monte Carlo (MCMC) algorithm. The optimum number of subpopulations (K) was confirmed after 20 independent runs for each value of K between 1 and 10. The length of the burn-in period and number of MCMC reps after burn-in were set to 25,000 and 100,000, respectively. The K subpopulations identified indicated clusters characterized by a set of allele frequencies at each locus, where individuals were assigned to a subpopulation, or to two or more populations, if the genotype indicated that they are admixed [33]. In this study, the identification of K used the model developed by Evanno et al. [34]. The Bayesian framework was not used to estimate the non-homogeneous original populations, instead we used ΔK, which was based on the rate of change in the log probability of data between successive Ks. STRUCTURE accurately detected the uppermost hierarchical level structure for the scenarios tested.

A Mantel test, performed with GenAlEx version 6.41 [30], was used to calculate the coefficient of association between genetic distance and altitude.


SSR genotyping

SSRs are generally used in genetic diversity studies as evolutionary neutral markers. In this study, 24 SSR primer pairs were developed using the Populus euphratica genome, which were transferable to the Tibetan poplar. Sequencing results were uploaded to GenBank, (Table 2). In total, 114 alleles for 24 loci were amplified (mean = 4.75, SD = 2.71), with locus 7 having 12 alleles and exhibiting the most variation. Locus 18 was detected using LOSITAN based on its FST value(0.16) [24], which showed that it was under positive natural selection ( p = 0 . 01 ) (Figure 2). Subsequently, the sequence was processed using NCBI BLAST [35], and there was high homology with a protein (ID: XM_002311699.1) present in Populus trichocarpa. This implied that this SSR locus could be under selection pressure.

Table 2 Descriptions of and references for the 24 SSR loci analyzed.
Figure 2
figure 2

F ST values of 24 microsatellite loci vs heterozygosity in the Populus szechuanica population. The dot in the red area is the locus U65600 which is positively under nature selection.

Genetic diversity

Previous assessment of genetic diversity among three populations of Tibetan poplar was based on allelic variation observed at 23 neutral microsatellite loci. In this study, a mean of 3.71 alleles per locus were confirmed for all 23 loci in 172 Tibetan poplar individuals. HO and HE are important parameters for assessing genetic diversity of populations, and they ranged from 0.40-0.42 and 0.48-0.49, respectively. In the three populations, as the results were consistent with the Na, it indicated that genetic variation was not significant, and the populations were similar in all parameters. FIS presented a similar pattern in that no significant difference was detected among populations; it ranged from -0.85 to 0.56 (mean = 0.17, SD = 0.18), with only a moderate difference for loci from trees at low altitude. The FST was estimated at each locus for all individuals, and ranged from 0.0004 to 0.16 (table 3), indicating that there was no evidence to support the hypothesis that the populations differed from each othe r.

Table 3 Comparisons of genetic diversity and differentiation among Populus szechuanica populations along altitude gradients.

Genetic structure

The AMOVA indicated different levels of genetic variance among populations and among individuals within populations. Of the total genetic variance, 6.67% was ascribed to population divergence; the remainder was ascribed to the differences between individuals. However, there was a significant difference among populations (p < 0.001). In populations sampled from high, medium, and low altitudes, all genetic diversity parameters were similar, indicating no local adaptation or population differentiation in the study area.

The SSR data was sorted in order of altitude. Population structure analysis was processed according to the known order of individuals, yielding an optimal of K = 2 [34]. Estimated populations of the 172 individuals are shown (Figure 3). The samples plot showed that low- and high-altitude individuals were considered to originate from a single group. However, the medium-altitude group was an admixture of the high- and low-altitude groups, and there was no clear separation between the groups (Figure 4). The FST value showed little differentiation between the populations. Since STRUCTURE could not perform an analysis of K = 1 on populations with no difference, we did not accept the results of K = 2, based on the low FST and the genetic parameters pattern among the three populations. To compare the structure of new clusters (cluster 1 and cluster 2), further STRUCTURE analyses were performed in cluster 1, which contained individuals from a high altitude, and cluster 2 from low altitude. It shows no clear structure despite peaks in K = 3 for cluster 1 and K = 2 for cluster 2 (Figure 5). The ancestry values of all of the individuals revealed that each had an equal probability of being grouped in cluster 1, 2 or 3 for the high- and low-altitude clusters. An analysis based on the Mantel test (Figure 6) showed that genetic distance was not significantly correlated with altitude (r2 = 0.001, p ≤ 0.07), suggesting that altitude was not the principal factor influencing genetic differentiation in the Tibetan poplar.

Figure 3
figure 3

Identification of K. the method of delta K was used to identify the accurate sub-clusters in the population. In this population there is a peak of delta K in K = 2, the population is possibly composed of two sub-clusters.

Figure 4
figure 4

Population structure of Populus szechuanica estimated by STRUCTURE. In the figure, the individuals were sorted as the altitude of sampling distributed area.

Figure 5
figure 5

Structure analysis of two sub-clusters. A,C) the sub-clusters detected at low altitude; B,D) the sub-clusters detected at high altitude

Figure 6
figure 6

Relevance between genetic distance and altitude.


SSR markers and neutrality

In this study, we used GeneMarker version 1.75 to identify the fluorescently labeled PCR products. We selected 24 SSR primer pairs based on the Populus euphratica genome, to analyze genetic diversity within three populations of Tibetan poplar living at different altitudes in Linzhi, Tibet. As the loci were transferable between the two species, this indicates that they may be sited in a conserved region. However, based on the FST using LOSITAN, one locus appeared to be an outlier. SSR loci mutations occasionally occur as a result of the stress of adapting to a change of environment [36] or an external stimulus [37]. Further, studies have shown that some SSR loci are non-neutral [38, 39], and for this reason it is essential that a neutrality test is performed before the SSR data are used in any further analysis. The outlier locus sequence was processed using NCBI BLAST [35], and indicated high homology with a protein in Populus trichocarpa. We conclude that the microsatellite may be linked to expressed genes, and therefore, neutrality should not be assumed, but tested in all of the markers before genetic diversity and structure analysis. This type of marker, however, could be useful for phylogenetic studies of closely related species [40, 41].

Genetic diversity

As expected from perennial and woody species ranging across most areas of the Qinghai-Tibet plateau, the study population contained a high level of genetic diversity, but we did not identify any significant differences among the three populations from different altitudes. The number of alleles per locus in our study was less than in other related Populus species [42]. A mean of 6.1 alleles per locus was identified from the existing literature on Populus genetic diversity [42]. The Na of 3.73 in our study is lower than the Na in P. tremuloides (4.9) as described previously [43]. The difference is most likely due to the limited sampling area. We only collected samples from one mountain area, whereas the Tibetan poplar is distributed throughout southwestern China, of which our samples were from a limited proportion, as we aimed to study adaptation and genetic diversity along an altitude gradient. The samples from high, medium, and low altitudes appeared to be similar in genetic diversity and showed no evidence of local adaptation in the study area. STRUCTURE analysis showed that the population could be divided into two groups (clusters), with individuals from the lower altitude clustered into group 1, and those from the higher altitude clustered into group 2. Altitude appeared to have a direct relationship with the distribution of the groups, but the FST value showed little differentiation between the populations. As STRUCTURE could not provide data for K = 1, we rejected the result showing that the population was divided into two groups. There was no peak in the estimate of the log-likelihood of the cluster number (L(k)) since the lowest likelihood was for K = 1, and L(k) either consistently increased or showed an erratic pattern with increasing variance, with all individuals admixed and the proportion of any individual assigned to each subpopulation remaining roughly similar. The Evanno criterion, ΔK [34], was not relevant as it can only be computed for K ≥ 2 and does not enable comparison of results from K = 1. For K > 2, the value of ΔK remained close to 0 in this study.

The population structure and Mantel test results suggest that the relationship between genetic diversity and altitude is not significant, and hence it is possible to hypothesize that the species has not had sufficient time for evolutionary differentiation to occur along an altitude gradient.

Low FSTand strong gene flow

FST was low for all loci, except for SSR 18. There was no noticeable differentiation among populations at three different altitudes. This may contribute to the local geographic structure and strong gene flow among individuals. The STRUCTURE results showed that the medium-altitude group was an admixture of the low- and high-altitude groups, clearly indicating that the mountain harbored two groups (clusters) of poplars, and that they separated into these clusters at an altitude of ~2700 m. Because the study area altitude ranged from 2 000 to 4 000 m, and the Tibetan poplar is distributed from 2000 to 3096 m, the tree line represented a limiting factor for tree distribution, but it appeared to have had limited impact on gene exchange between individuals and did not hinder pollen or seed dispersal. In this study, gene flow occurred among the populations. Gene flow is a vital element in local adaptation studies, because it can instruct the establishment of the local genetic structure or influence it indirectly. Gene flow among populations can also lead to combining of gene pools, reducing genetic variation among groups [44]. Therefore, gene flow acts strongly against speciation in evolutionary processes [45], by recombining the gene pools of the groups. Gene flow plays a part in evolution through pollen dispersal, seed dispersal, and the establishment of the individual adult. A geographic barrier increases the probability of extinction or local adaptation of a population, as it may push the population to evolve into a different population with a unique genetic structure, or even into a new species [46, 47]. However, gene flow could also be a constraining force of natural evolution by homogenizing populations under a heterogeneous environment, and balancing gene distribution and spread [48]. However, gene flow can also be considered a creative force in evolution, where superior genes or combinations of genes are spread by gene flow [49, 50]. For local adaptation, gene flow and selection are usually considered as the main forces affecting the processes of establishment. This is especially true for high outcrossing trees and perennial species, where there is extensive gene flow [51]. In summary, the factors contributing to the low level of differentiation among populations at different altitudes include: (1) Pollen dispersal and an overlapping flowering period of all three populations (high, medium and low altitude). Generally the flowering phase of Populus is of long duration; for example, flowering in P. × canadensis and P. nigra [52] lasts for 15 and 31 days, respectively. (2) Seed dispersal mechanisms. Most Populus trees live adjacent to rivers and roads, and some in the river channel itself. Therefore, rivers cannot be ignored as an important factor in seed dispersal. Poplar populations are evolutionarily homogeneous. The germplasm and genetic diversity of the Tibetan poplar could be protected by random selection in the future work, which couldprovid all of the genetic diversity to date. An unpublished experiment comparing poplars at two sites showed some differences in the growth rate, leaf characteristics, and branch numbers, etc. of individual clones sampled at different altitudes, indicating that natural selection conserved some fitness types. Genes linked to adaptation mechanisms could contribute to phenotypic variation without genetic structure differentiation which has been proved in this study. Consequently, this makes the population ideal for identifying functional genes and mechanisms of adaptation to high altitudes.


To our knowledge, this is the first genetic analysis of the Tibetan poplar. The results indicate that the Tibetan poplar populations living at different altitudes on the Sejila mountain have a low level of differentiation. They have an excellent ability to adapt to different altitudes; however, local adaptation is not observed due to the lack of a geographic barrier. The high levels of gene flow lead to a low FST, as was observed. We consider the Sejila mountain population to be appropriate for investigation of the mechanisms of adaptation to high altitudes, despite the low level of genetic structure differentiation among populations at different altitudes.


Publication charges for this article came from the Special Fund for Forest Scientific Research in the Public Welfare (201404102), NSF/IOS-0923975, Changjiang Scholars Award and "Thousand-person Plan" Award.

This article has been published as part of BMC Genetics Volume 15 Supplement 1, 2014: Selected articles from the International Symposium on Quantitative Genetics and Genomics of Woody Plants. The full contents of the supplement are available online at


  1. Körner C: The use of 'altitude'in ecological research. Trends in ecology & evolution. 2007, 22 (11): 569-574. 10.1016/j.tree.2007.09.006.

    Article  Google Scholar 

  2. Pickup M, Barrett SC: The influence of demography and local mating environment on sex ratios in a wind-pollinated dioecious plant. Ecology and evolution. 2013, 3 (3): 629-639. 10.1002/ece3.465.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Byars SG, Papst W, Hoffmann AA: Local adaptation and cogradient selection in the alpine plant, Poa hiemata, along a narrow altitudinal gradient. Evolution. 2007, 61 (12): 2925-2941. 10.1111/j.1558-5646.2007.00248.x.

    Article  PubMed  Google Scholar 

  4. Lesica P, Allendorf FW: When are peripheral populations valuable for conservation?. Conservation Biology. 1995, 9 (4): 753-760. 10.1046/j.1523-1739.1995.09040753.x.

    Article  Google Scholar 

  5. Rainey PB, Travisano M: Adaptive radiation in a heterogeneous environment. Nature. 1998, 394 (6688): 69-72. 10.1038/27900.

    Article  CAS  PubMed  Google Scholar 

  6. Ortego J, Riordan EC, Gugger PF, Sork VL: Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak. Molecular Ecology. 2012, 21 (13): 3210-3223. 10.1111/j.1365-294X.2012.05591.x.

    Article  PubMed  Google Scholar 

  7. Mosca E, Eckert A, Di Pierro E, Rocchini D, La Porta N, Belletti P, Neale D: The geographical and environmental determinants of genetic diversity for four alpine conifers of the European Alps. Molecular ecology. 2012, 21 (22): 5530-5545. 10.1111/mec.12043.

    Article  CAS  PubMed  Google Scholar 

  8. Avolio ML, Beaulieu JM, Smith MD: Genetic diversity of a dominant C4 grass is altered with increased precipitation variability. Oecologia. 2013, 171 (2): 571-581. 10.1007/s00442-012-2427-4.

    Article  PubMed  Google Scholar 

  9. Manel S, Gugerli F, Thuiller W, Alvarez N, Legendre P, Holderegger R, Gielly L, Taberlet P: Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation. Molecular Ecology. 2012, 21 (15): 3729-3738. 10.1111/j.1365-294X.2012.05656.x.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Oyama K, Ito M, Yahara T, Ono M: Low genetic differentiation among populations ofArabis serrata (Brassicaceae) along an altitudinal gradient. Journal of plant research. 1993, 106 (2): 143-148. 10.1007/BF02344418.

    Article  Google Scholar 

  11. Taira H, Tsumura Y, Tomaru Y, Ohba K: Regeneration system and genetic diversity of Cryptomeria japonica growing at different altitudes. Canadian Journal of Forest Research. 1997, 27 (4): 447-452. 10.1139/x97-001.

    Article  Google Scholar 

  12. Jump AS, Hunt JM, MARTÍNEZ-IZQUIERDO JA, Penuelas J: Natural selection and climate change: temperature-linked spatial and temporal trends in gene frequency in Fagus sylvatica. Molecular Ecology. 2006, 15 (11): 3469-3480. 10.1111/j.1365-294X.2006.03027.x.

    Article  CAS  PubMed  Google Scholar 

  13. Mathiasen P, Premoli AC: Fine-scale genetic structure of Nothofagus pumilio (lenga) at contrasting elevations of the altitudinal gradient. Genetica. 2013, 141 (1-3): 95-105. 10.1007/s10709-013-9709-6.

    Article  PubMed  Google Scholar 

  14. Korshikov I, Mudrik E: Elevation-dependent genetic variation of plants and seed embryos in the Crimea Mountain population of Pinus pallasiana D. Don. Russian Journal of Ecology. 2006, 37 (2): 79-83. 10.1134/S1067413606020020.

    Article  CAS  Google Scholar 

  15. Zhang Y, Li B, Zheng D: A discussion on the boundary and area of the Tibetan Plateau in China. Geographical Research. 2002, 21 (1): 1-8.

    Google Scholar 

  16. Zheng D: The system of physico-geographical regions of the Qinghai-Tibet (Xizang) Plateau. Science in China (Series D). 1996, 39 (4): 410-417.

    Google Scholar 

  17. Qiu Y-X, Fu C-X, Comes HP: Plant molecular phylogeography in China and adjacent regions: tracing the genetic imprints of Quaternary climate and environmental change in the world's most diverse temperate flora. Molecular Phylogenetics and Evolution. 2011, 59 (1): 225-244. 10.1016/j.ympev.2011.01.012.

    Article  PubMed  Google Scholar 

  18. Chen S-l, Wu Z-y, Raven PH: 1994, Flora of China: Science Press

  19. TangYudan P, CidanZhuoga : Biological Characteristics of Populus szechuanica var. tibetica theRareand Endemic Plant of Qinghai-Tibetan Plateau in the Different Local Environment. Chinese Wild Plant Resources. 2012, 31 (2): 24-32.

    Google Scholar 

  20. Xu F, Feng S, Wu R, Du FK: Two highly validated SSR multiplexes (8-plex) for Euphrates' poplar, Populus euphratica (Salicaceae). Molecular ecology resources. 2013, 13 (1): 144-153. 10.1111/1755-0998.12030.

    Article  CAS  PubMed  Google Scholar 

  21. Doyle JJ: A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987, 19: 11-15.

    Google Scholar 

  22. Amos W, Hoffman J, Frodsham A, Zhang L, Best S, Hill A: Automated binning of microsatellite alleles: problems and solutions. Molecular Ecology Notes. 2007, 7 (1): 10-14. 10.1111/j.1471-8286.2006.01560.x.

    Article  CAS  Google Scholar 

  23. Shaibi T, Lattorff H, Moritz R: A microsatellite DNA toolkit for studying population structure in Apis mellifera. Molecular Ecology Resources. 2008, 8 (5): 1034-1036. 10.1111/j.1755-0998.2008.02146.x.

    Article  CAS  PubMed  Google Scholar 

  24. Antao T, Lopes A, Lopes R, Beja-Pereira A, Luikart G: LOSITAN: a workbench to detect molecular adaptation based on a Fst-outlier method. BMC bioinformatics. 2008, 9 (1): 323-10.1186/1471-2105-9-323.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Beaumont MA: Adaptation and speciation: what can< i> F</i>< sub> st</sub> tell us?. Trends in Ecology & Evolution. 2005, 20 (8): 435-440. 10.1016/j.tree.2005.05.017.

    Article  Google Scholar 

  26. Cockerham CC, Weir B: Covariances of relatives stemming from a population undergoing mixed self and random mating. Biometrics. 1984, 40 (1): 157-164. 10.1007/978-1-4419-1541-2_2.

    Article  CAS  PubMed  Google Scholar 

  27. Nei M: Analysis of gene diversity in subdivided populations. Proceedings of the National Academy of Sciences. 1973, 70 (12): 3321-3323. 10.1073/pnas.70.12.3321.

    Article  CAS  Google Scholar 

  28. Nei M: Molecular evolutionary genetics. 1987, Columbia University Press

    Google Scholar 

  29. Yeh FC, Yang R, Boyle T, Ye Z, Mao JX: POPGENE, version 1.32: the user friendly software for population genetic analysis. Molecular Biology and Biotechnology Centre. 1999, University of Alberta, Edmonton, AB, Canada

    Google Scholar 

  30. Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes. 2006, 6 (1): 288-295. 10.1111/j.1471-8286.2005.01155.x.

    Article  Google Scholar 

  31. Pritchard JK, Donnelly P: Case-control studies of association in structured or admixed populations. Theoretical population biology. 2001, 60 (3): 227-237. 10.1006/tpbi.2001.1543.

    Article  CAS  PubMed  Google Scholar 

  32. Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Molecular ecology resources. 2009, 9 (5): 1322-1332. 10.1111/j.1755-0998.2009.02591.x.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.

    PubMed Central  CAS  PubMed  Google Scholar 

  34. Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Molecular ecology. 2005, 14 (8): 2611-2620. 10.1111/j.1365-294X.2005.02553.x.

    Article  CAS  PubMed  Google Scholar 

  35. Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL: NCBI BLAST: a better web interface. Nucleic acids research. 2008, 36 (suppl 2): W5-W9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Ellegren H: Microsatellites: simple sequences with complex evolution. Nature Reviews Genetics. 2004, 5 (6): 435-445. 10.1038/nrg1348.

    Article  CAS  PubMed  Google Scholar 

  37. Schmidt AL, Mitter V: Microsatellite mutation directed by an external stimulus. Mutation Research/Fundamental and Molecular Mechanisms of Mutagenesis. 2004, 568 (2): 233-243. 10.1016/j.mrfmmm.2004.09.003.

    Article  CAS  Google Scholar 

  38. Li Y-C, Korol AB, Fahima T, Nevo E: Microsatellites within genes: structure, function, and evolution. Molecular biology and evolution. 2004, 21 (6): 991-1007. 10.1093/molbev/msh073.

    Article  CAS  PubMed  Google Scholar 

  39. Ganopoulos I, Aravanopoulos F, Argiriou A, Tsaftaris A: Genome and population dynamics under selection and neutrality: an example of S-allele diversity in wild cherry (Prunus avium L.). Tree Genetics & Genomes. 2012, 8 (6): 1181-1190. 10.1007/s11295-012-0504-9.

    Article  Google Scholar 

  40. Garcia-Lor A, Curk F, Snoussi-Trifa H, Morillon R, Ancillo G, Luro F, Navarro L, Ollitrault P: A nuclear phylogenetic analysis: SNPs, indels and SSRs deliver new insights into the relationships in the 'true citrus fruit trees' group (Citrinae, Rutaceae) and the origin of cultivated species. Annals of botany. 2013, 111 (1): 1-19. 10.1093/aob/mcs227.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Nas MN, Bolek Y, Bardak A: Genetic diversity and phylogenetic relationships of< i> Prunus microcarpa</i> CA Mey.< i> subsp. tortusa</i> analyzed by simple sequence repeats (SSRs). Scientia Horticulturae. 2011, 127 (3): 220-227. 10.1016/j.scienta.2010.09.018.

    Article  CAS  Google Scholar 

  42. Slavov GT, Zhelev P: Salient biological features, systematics, and genetic variation of Populus. Genetics and Genomics of Populus. 2010, Springer, 15-38.

    Chapter  Google Scholar 

  43. Cole CT: Allelic and population variation of microsatellite loci in aspen (Populus tremuloides). New Phytologist. 2005, 167 (1): 155-164. 10.1111/j.1469-8137.2005.01423.x.

    Article  CAS  PubMed  Google Scholar 

  44. Slatkin M: Gene flow in natural populations. Annual review of ecology and systematics. 1985, 16: 393-430. 10.1146/annurev.ecolsys.16.1.393.

    Article  Google Scholar 

  45. Kronforst MR: Gene flow persists millions of years after speciation in Heliconius butterflies. BMC evolutionary biology. 2008, 8 (1): 98-10.1186/1471-2148-8-98.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Vermeij GJ: The dispersal barrier in the tropical Pacific: implications for molluscan speciation and extinction. Evolution. 1987, 1046-1058.

    Google Scholar 

  47. Barnes I, Matheus P, Shapiro B, Jensen D, Cooper A: Dynamics of Pleistocene population extinctions in Beringian brown bears. Science. 2002, 295 (5563): 2267-2270. 10.1126/science.1067814.

    Article  CAS  PubMed  Google Scholar 

  48. Storfer A, Sih A: Gene flow and ineffective antipredator behavior in a stream-breeding salamander. Evolution. 1998, 52 (2): 558-565. 10.2307/2411090.

    Article  Google Scholar 

  49. Hendry AP, Taylor EB, McPhail JD: Adaptive divergence and the balance between selection and gene flow: lake and stream stickleback in the Misty system. Evolution. 2002, 56 (6): 1199-1216. 10.1111/j.0014-3820.2002.tb01432.x.

    Article  PubMed  Google Scholar 

  50. Olson-Manning CF, Wagner MR, Mitchell-Olds T: Adaptive evolution: evaluating empirical support for theoretical predictions. Nature Reviews Genetics. 2012, 13 (12): 867-877. 10.1038/nrg3322.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Muona O, Brown A, Clegg M, Kahler A, Weir B: Population genetics in forest tree improvement. Plant population genetics, breeding, and genetic resources. 1990, 282-298.

    Google Scholar 

  52. Broeck AV, Cox K, Quataert P, Van Bockstaele E, Van Slycken J: Flowering Phenology of Populus nigra L., P. nigra cv. italica and P. × canadensis Moench. and the Potential for Natural Hybridisation in Belgium. Silvae genetica. 2003, 52 (5/6): 280-283.

    Google Scholar 

Download references


The authors thank Dr. Fang Du for comments and for revision of the manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Rongling Wu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

Rongling Wu designed the study. Dengfeng Shen, Wenhao Bo and Fang Xu, contributed extensively to the samples collection. Dengfeng Shen performed PCR experiments and genetic diversity analysis. Dengfeng Shen and Rongling Wu wrote the manuscript. Wenhao Bo and Fang Xu prepared and revised the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shen, D., Bo, W., Xu, F. et al. Genetic diversity and population structure of the Tibetan poplar (Populus szechuanica var. tibetica) along an altitude gradient. BMC Genet 15 (Suppl 1), S11 (2014).

Download citation

  • Published:

  • DOI: