- Original Paper
- Published:
Genetic variation within and between populations of a desert poplar (Populus euphratica) revealed by SSR markers
Annals of Forest Science volumeΒ 68, ArticleΒ number:Β 1143 (2011)
Abstract
Backgroud
Populus euphratica Oliv., a long-lived woody perennial plant, is the only forest species naturally distributed in desert regions. Severe desertification has resulted in this species becoming endangered, and its populations are seriously fragmented in northwest China.
Objectives
We investigated the within-population and long-distance separated population diversity in northwest China using eight pairs of simple sequence repeat markers.
Results
Our results reveal that this species has high overall genetic diversity with a mean of 12.125 alleles per locus; its expected heterozygosity (H e) ranged from 0.713 to 0.878. However, there was a high level of genetic diversity within the species (when compared within congeneric species) and low genetic differentiation between populations (average F stβ=β0.093). Analyses of molecular variance suggested that 5.21% of the total molecular variance was attributable to between-population diversity (Pβ<β0.001), while the remainder of the variance was associated with differences within populations. There was no distinct correlation between geographical distributions and genetic variation.
Conclusions
We found no evidence to support our initial hypothesis of low genetic diversity within the species and high differentiation between populations separated by long distances. The recent fragmentations of this species due to anthropologic and environmental effects resulted in its endangered status.
1 Introduction
The maintenance of genetic variation provides the potential for tree species to adapt to long-term environmental change (Hedrick 2004). Genetic diversity is also of importance in enabling a population to survive and reproduce under changing environmental conditions. A decrease in genetic variation has been correlated with a reduced adaptive flexibility, which threatens the ability of a species to survive via natural regeneration and has deleterious effects on species fitness. Therefore, assessments of the level and distribution of genetic diversity are important factors informing the management and development of effective conservation strategies for tree species. Most tree species exhibit a high level of genetic variation and a low level of between-population differentiation because they usually have widespread and continuous distributions and outcrossing breeding systems (Petit et al. 2004). However, for those trees species with fragmented distributions, intraspecific diversity tends to be relatively low and the genetic differentiation between populations is very high (Meng et al. 2007).
From an ecological and economic perspective, Populus euphratica is one of the most important poplar species in northwest China and is well known for its phreatophytic habit of growing in deserts (Bruelheide et al. 1999). P. euphratica is the main species responsible for afforestation of saline soil in sandy desert regions in northwest China and occurs at a few sites separated by a distance of 2, 000Β km east to west and 130Β km north to south (Wang et al. 1995). The number of populations and individuals of this species has dropped considerably in recent years to the point where it is now considered to be endangered (Bruelheide et al. 1999). Such limited distribution and the long distances between isolated populations is typically indicative of a low genetic diversity within species and a high between-population differentiation. These characteristic features were observed in part by Saito et al. (2002) based on results from random amplified polymorphic DNA markers (RAPDs) during a study of populations from northern Xinjiang, which are genetically isolated from other populations in southern Xinjiang. One limitation of this study in establishing the level of diversity was the limited distribution of the five populations sampled as these were all collected around the Tianshan mountain range region of Xinjiang. Other representative populations of P. euphratica in northwest China (for example, in Gansu and Mongolia) were not included in the study. Although the RAPD technique can be a valuable tool for obtaining overall estimates of the genetic diversity of a species, it also has limitations, such as poor reproducibility and dominance of molecular markers. In contrast, simple sequence repeats (SSRs) provide more information because they are hypervariable and because of their co-dominance based on the presence of tandem repeated units of nucleotides that are ubiquitous in both protein coding and non-coding regions and are abundant in prokaryotic and eukaryotic genomes. Their uses have included studies of population and conservation genetics, the construction of genetic maps, parenthood analyses (van de Ven and McNicol 1996; Lexer et al. 2010), the assessment of genetic effects of forest management practice, and the development of strategies for conservation and sustainable management of forest genetic resources (Rajora et al. 2000; Rathmacher et al. 2010).
In this study, we investigate genetic variations within and between 10 populations of P. euphratica separated by long distances across northwest China using eight SSR loci (Wu et al. 2008). The two specific aims of this study were: (1) to assess genetic variation within and between populations and (2) to test the hypothesis that the distribution fragmentation has led to a decrease in genetic diversity and an increase in between-population differentiation.
2 Materials and methods
2.1 Plant species and sampled populations
We sampled 10 populations separated by long distances in northwest China. The seeds of P. euphratica are dispersed by both water and wind, and usually travel more than 200Β m. We chose 10β30 adult trees (with diameters at breast height at least 0.5Β m) spaced at least 200Β m apart in each isolated forest stand. Two to three clean leaves were dried and stored in silica gel prior to DNA extraction. Thus, a total of 170 samples were used for DNA analysis, representing populations distributed across the northwest of China (TableΒ 1). Latitude, longitude, and altitude were recorded for each collection site using an Etrex GIS unit (Garmin, Taiwan).
2.2 DNA extraction, PCR amplification, and polyacrylamide gel electrophoresis
Total DNA was extracted from the silica gel-dried leaves using the modified 2Γ CTAB procedure (Doyle and Doyle 1987). Eight SSR markers (Wu et al. 2008) were selected on the basis of their ability to amplify the DNA of all of the plants and the reproducibility of the products. The loci were located on the different genome regions according to the preliminary comparison to the available poplar genome (Tuskan et al. 2006). The SSR reactions were based on the protocol of Wu et al. (2008), with some modifications. Amplification reactions were performed in a volume of 20Β ΞΌl containing 1Γ reaction buffer (Takara), 8β12Β ng genomic DNA, 200Β ΞΌM of each dNTP, 16Β ΞΌM of each primer, and 1Β U of Taq DNA polymerase (Takara). The Cyclerβ’ Thermal Cycler (Bio-Rad, USA) PCR was programmed according to the following profile: 94Β°C for 7Β min; 94Β°C for 40Β s, 50β58Β°C for 30Β s, 72Β°C for 30Β s, 37 cycles; and 72Β°C for 7Β min. After PCR amplification, 2Β ΞΌl of each sample was loaded and electrophoresed on a 2% horizontal agarose gel to check for positive amplification and to determine the approximate amount of the product. The remaining PCR products were separated on 6.5% (w/v) polyacrylamide denaturing gel and visualized by silver staining. The band size was monitored using a 10-bp DNA ladder (Tiangen) as the reference. In order to reduce errors, we maintain all experiments for the same condition.
2.3 Summary statistics
For each locus, genotypes exhibiting one or two bands were scored as being homozygous or heterozygous, respectively. We did not detect null alleles in the sampled loci. We measured intra- and inter-population genetic statistics using the program GENEPOP 3.4 (Raymond and Rousset 1995), including the expected heterozygosity (H e), observed heterozygosity (H o), and the effective number of alleles (N e). We also calculated Shannonβs indices (H pop) to characterize genetic variation within and between populations.
We further calculated the pairwise standard genetic distance (Nei 1972), and an unrooted tree was constructed using a least squares algorithm with 10,000 bootstraps over the microsatellite loci, implemented in the PHYLIP v. 3.6c software package (Felsenstein 1993). We determined the genetic relationship between populations with the STRUCTURE program using the admixture model (Pritchard et al. 2000). This model assumes that the genome of an individual is a mixture of genes originating from K unknown βancestralβ populations that may have undergone introgression events (the unknown K ancestral populations were assumed to have been at HW equilibrium). Eighteen independent runs were performed for each K value, from Kβ=β1 to Kβ=β10. The range of K values was set to include the total number of populations (Kβ=β10). Each run consisted of a burn-in period of 100,000 steps, followed by 106 MCMC (Monte Carlo Markov chain) replicates. The choice of the appropriate K value was based on the recommendations of Pritchard et al. (2000). The runs with the highest lnPr(X|K) values were selected and the proportion of ancestry of each population in each of the clusters was calculated by averaging the estimated membership coefficient of the combined individuals and ancestry estimate. For individuals, the admixture level was analyzed by assigning the individual to a cluster if an arbitrary value of 75% of its genome appeared to belong to that cluster (Matsuoka et al. 2002), while those individuals with membership probabilities of <75% for all clusters were considered to be of βmixed origin.β In order to test for a correlation between Neiβs genetic distance and geographic distance (in meters) between populations, a Mantel test was performed using TFPGA (Miller 1997) for population genetic analysis (999 permutations were calculated). We also derived approximations of the genetic structures using Wrightβs analysis of hierarchical F statistics (Wright 1978). Gene flow (N m) was estimated using the equation N mβ=β0.25 (1βββF st)/F st. To estimate the variance component and to partition the variation within and between populations, we used analysis of molecular variance (AMOVA) implemented in the ARLEQUIN program (Schneider et al. 2000).
3 Results
3.1 Microsatellite polymorphisms
The total number of alleles per locus and the size range of the alleles detected in P. euphratica are given in TableΒ 2. A total of 97 alleles were detected at eight SSR loci within the 10 sampled populations. The number of alleles per locus ranged from 8 to 16, with a mean value of 12.125. Genetic diversity parameters per locus indicated a high number of polymorphisms, which are listed in TableΒ 2. The total H o ranged from 0.761 to 1.000. The average expected H e was slightly lower than the H o, ranging from 0.713 to 0.878. The observed number of alleles per locus (N a) ranged from 8 to 13, with an average of 11, while the effective number of alleles per locus (N e) ranged from 3.45 to 8.00, with an average of 6.56. The maximum value of H o was found in Pop3 and the minimum in Pop9 (TableΒ 3). H e varied from 0.720 (Pop5) to 0.831 (Pop10). Shannonβs information index (H pop) ranged from 1.300 (Pop6) to 1.857 (Pop2), with an average of 1.631 alleles (TableΒ 3).
3.2 Genetic differentiation between populations
To estimate genetic differentiation between populations, the fixation index F st, originally formulated by Wright (1978), was utilized using the program POPGENE 3.4 (Raymond and Rousset 1995), in which modifications were made to avoid errors arising from the small sample size and to enable analyses of multi-alleles at multi-loci. The F st values ranged from 0.025 to 0.162, with an average value of 0.093. The overall gene flow (N m) between populations was estimated to be 2.446, which provided an estimate of the average number of migrants between all studied populations per generation (Table 4).
Long-distance isolation usually results in an increase in genetic differentiation. In order to identify any geographical correlation between the recovered gene pools or genetic groups and sampled populations, we performed a Mantel test with 1,000 permutations. We detected a non-significant correlation between genetic distance and geographic distance (rβ=β0.3217, pβ=β0.058β>β0.05; Fig.Β 1). The most probable division of populations was when Kβ=β4; this received the highest support in terms of log-likelihood values, as suggested by Pritchard et al. (2000). The four clustered groups of populations showed no correspondence to geographic location (Fig.Β 2), nor was there any geographical pattern within the clustered populations. We used AMOVA to examine further the genetic differentiation within and between populations. We conducted two types of group estimate: (1) four genetic groups identified by genetic distances and using STRUCTURE 2.3.2 and (2) two geographical regions isolated by deserts (Xinjiang versus Gansu and Inner Mongolia), with a long distance between them (Fig.Β 1). AMOVA analyses indicated that 5.21% of the total molecular variance was attributable to between-population diversity (Pβ<β0.001), while the remainder (94.79%) was associated with differences within populations (TableΒ 5).
4 Discussion
In this study, SSR markers were used to estimate the level of genetic diversity, an approach that has been shown to be appropriate in many studies (van de Ven and McNicol 1996; Dayanandan et al. 1998; Rajora et al. 2000; Fossati et al. 2003; Rathmacher et al. 2010). In the present study, we examined the genetic variations of the only eight SSR loci and used the polyacrylamide denaturing gel and silver staining to check genetic polymorphisms at each locus. We sampled only 10 individuals for a few populations. The genetic diversity within species and population level will increase with the sampled individuals and populations (e.g., Dayanandan et al. 1998). In addition, our experimental approach may neglect the SSR polymorphisms caused by the small repeats at the same locus between different individuals when the resultant migrating differences are too low. These limitations may underestimate the genetic polymorphisms of this species and the sampled populations, but will not affect the main aim of this study to examine the genetic differentiation between populations. Despite this, we still identified a high level of polymorphism within each population (TableΒ 2) and the species across all loci (TableΒ 3). The observed heterozygosity (H oβ=β0.932) had significantly higher values than the expected heterozygosity (H eβ=β0.787; TablesΒ 2 and 3). All these estimates were much higher than the genetic diversity indexes for this species found using RAPD markers (Saito et al. 2002). Similar values were also reported for congeneric Populus nigra with SSR markers (H oβ=β0.70, H eβ=β0.73; Rathmacher et al. 2010). However, lower levels of polymorphism have been reported for two other poplar species: Populus alba (H o = 0.341, H eβ=β0.368) and Populus tremula (H oβ=β0.483, H eβ=β0.492; Lexer et al. 2005).
The desert poplar studied here is found in northwest China, and the sampled populations were separated by long distances. However, we did not detect high population differentiation (F stβ=β0.093). The low level of differentiation was also confirmed by the AMOVA, which indicated that only 5.21% of the total molecular variance was attributable to between-population diversity (TableΒ 5). The STRUCTURE 2.3.2 analyses and Mantel test also indicated that the sampled populations showed no corresponding geographical groupings and that the genetic variations were not correlated with the geographic isolation (Fig.Β 1 and Table 4). The low differentiation detected in the species is similar to that reported for congeneric species based on the SSR marking. For example, mean G ST values of only 0.15 and 0.017 have been reported for P. nigra and P. tremula, respectively, based on SSR analyses (Fossati et al. 2003; Lexer et al. 2005). These estimates are similar to estimates of corresponding parameters for wide-ranging species including Lychnis flos-cuculi L. (F stβ=β0.022) and Phlebia centrifuga (F stβ=β0.061; Galeuchet et al. 2005; Franzen et al. 2007).
Our genetic analyses indicate that the P. euphratica populations are characterized by a high level of genetic diversity and low population differentiation (TablesΒ 2 and 3). Three possible factors may account for such a genetic pattern. First, P. euphratica is dioecious and is an outcrossing species. Its pollen is dispersed by the wind over long distances, and this is likely to promote gene flow between isolated populations. This is undoubtedly likely to enhance intraspecific diversity and reduce inter-population differentiation (Hedrick 2004). Secondly, the majority of individuals of this species are long-lived (a life span of >200Β years has been estimated) and have extensive clonal growth (Wiehle et al. 2009). The current stands in most populations are, therefore, likely to be well established for a long time through cloning reproductions while these daughter ramets maintained the same genetic diversity as their parents (Vonlanthen et al. 2010). Thirdly, the long isolation between the sampled populations may have occurred recently because of human disturbance and climatic changes. For example, most forests composed of this species were cleared for agricultural purposes during the period from 1950 to 1960 (Chen et al. 2006). In addition, during the past 100Β years, precipitation in the northwest has declined dramatically and most rivers have dried up (Zhang et al. 2003). The shortage of water and expansion of the desert resulted in the death of poplars at a number of sites. However, at the other sites, a few populations persist due to the availability of waters in the microhabitats. The long life cycle and clonal reproduction of this species may have helped these current fragmented populations to retain their historical genetic diversity when these populations were connected. Therefore, it is imaginable that the fragmentations are too short to cause the genetic differentiations between populations. Overall, our findings suggest that the endangered status of the species is mainly due to anthropologic and environmental effects rather than a lack of genetic diversity.
References
Bruelheide H, Jandt U, Gries D, Thomas FM, Foetzki A, Buerkert A, Gang W, Zhang XM, Bussell JD (1999) The distribution of random amplified polymorphic DNA (RAPD) diversity amongst populations of Isotoma petraea (Lobeliaceae). Mol Ecol 8:775β789
Chen YN, Zilliacus H, Li WH, Zhang HF, Chen YP (2006) Ground-water level affects plant species diversity along the lower reaches of the Tarim River, Western China. J Arid Environ 66:231β246
Dayanandan SO, Rajora OP, Bawa KS (1998) Isolation and characterization of microsatellites in trembling aspen (Populus tremuloides). Theor Appl Genet 96:950β956
Doyle JJ, Doyle JL (1987) A rapid DNA isolation procedure from small quantities of fresh leaf tissues. Phyto Bull 19:11β15
Felsenstein J (1993) PHYLIP (Phylogeny Inference Package) version 3.5c. Department of Genetics, University of Washington, Seattle
Fossati T, Grassi F, Sala F, Castiglione S (2003) Molecular analysis of natural populations of Populus nigra L. intermingled with cultivated hybrids. Mol Ecol 12:2033β2043
Franzen I, Vasaitis R, Penttila R, Stenlid J (2007) Population genetics of the wood-decay fungus Phlebia centrifuga P Karst in fragmented and continuous habitats. Mol Ecol 16:3326β3333
Galeuchet DJ, Perret C, Fischer M (2005) Microsatellite variation and structure of 28 populations of the common wetland plant, Lychnis flos-cuculi L., in a fragmented landscape. Mol Ecol 14:991β1000
Hedrick PW (2004) Recent developments in conservation genetics. For Ecol Manag 197:3β19
Lexer C, Fay MF, Joseph A, Nica MS, Heinze B (2005) Barrier to gene flow between two ecologically divergent Populus species, P. alba (white poplar) and P. tremula (European aspen): the role of ecology and life history in gene introgression. Mol Ecol 14:1045β1057
Lexer C, Joseph JA, van Loo M, Barbara T, Heinze B, Bartha D, Castiglione S, Fay MF, Buerkle CA (2010) Genomic admixture analysis in European Populus spp. reveals unexpected patterns of reproductive isolation and mating. Genet 186:699β712
Matsuoka Y, Vigouroux Y, Goodman MM, Sanchez GJ, Buckler E, Doebley J (2002) A single domestication for maize shown by multilocus microsatellite genotyping. Proc Natl Acad Sci USA 99:6080β6084
Meng LH, Yang R, Abbott RJ, Miehe G, Hu TH, Liu JQ (2007) Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the QinghaiβTibetan Plateau and adjacent highlands. Mol Ecol 16:4128β4137
Miller MP (1997) Tools for population genetic analysis. Version 1.3. Department of Biological Sciences, Northern Arizona University, Flagstaff
Nei M (1972) Genetic distance between populations. Am Nat 106:283β292
Petit RJ, Bialozyt R, Garnier-Gereand P, Hampe A (2004) Ecology and genetics of tree invasions: from recent introductions to Quaternary migrations. For Ecol Manag 197:117β137
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945β959
Rajora OP, Rahman MH, Buchert GP, Dancik BP (2000) Microsatellite DNA analysis of genetic effects of harvesting in old-growth eastern white pine (Pinus strobes) Ontario, Canada. Mol Ecol 9:339β348
Rathmacher G, Niggemann M, Kohnen M, Ziegenhagen B, Bialozyt R (2010) Short-distance gene flow in Populus nigra L. accounts for small-scale spatial genetic structures: implications for in situ conservation measures. Conserv Genet 11:1327β1338
Raymond M, Rousset F (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered 86:248β249
Saito Y, Shiraishi S, Tanimoto T, Yin L, Watanabe S, Ide Y (2002) Genetic diversity of Populus euphratica populations in northwestern China determined by RAPD DNA analysis. New For 23:97β103
Schneider S, Roessli D, Excoffier L (2000) Arlequin ver. 2.0: a software for population genetic data analysis. Genetics and Biometry Laboratory, University of Geneva, Geneva
Tuskan GA, DiFazio S, Jansson S et al (2006) The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science 313:1596β1604
Van de Ven WTG, McNicol RJ (1996) Microsatellites as DNA markers in Sitka spruce. Theor Appl Genet 93:613β617
Vonlanthen B, Zhang XM, Bruelheide H (2010) Clonal structure and genetic diversity of three desert phreatophytes. Amer J Bot 97:234β242
Wang SJ, Chen BH, Li HQ (1995) Euphrates poplar forest. China Environmental Science Press, Beijing
Wiehle M, Eusemann P, Thevs N, Schnittler M (2009) Root suckering patterns in Populus euphratica (Euphrates poplar, Salicaceae). Trees 23:991β1001
Wright S (1978) Variability within and among natural populations. The University of Chicago Press, Chicago
Wu YX, Wang J, Liu JQ (2008) Development and characterization of microsatellite markers in Populus euphratica (Populaceae). Mol Ecol Res 8:1142β1144
Zhang H, Wu JW, Zheng QH, Yu YJ (2003) A preliminary study of oasis evolution in the Tarim Basin, Xinjiang, China. J Arid Envir 55:545β553
Acknowledgments
Financial support for this research was provided by NSFC (30800865) and the Program for New Century Excellent Talents in the Ministry of Education in China (NCET-09-0446).
Author information
Authors and Affiliations
Corresponding author
Additional information
Handling Editor: Luc PΓ’ques
Juan Wang and Zhijun Li contributed equally to this work.
Rights and permissions
About this article
Cite this article
Wang, J., Li, Z., Guo, Q. et al. Genetic variation within and between populations of a desert poplar (Populus euphratica) revealed by SSR markers. Annals of Forest Science 68, 1143 (2011). https://0-doi-org.brum.beds.ac.uk/10.1007/s13595-011-0119-6
Received:
Accepted:
Published:
DOI: https://0-doi-org.brum.beds.ac.uk/10.1007/s13595-011-0119-6