- Research article
- Open Access
Landscape genomics reveal that ecological character determines adaptation: a case study in smoke tree (Cotinus coggygria Scop.)
BMC Evolutionary Biology volume 17, Article number: 202 (2017)
The adaptive evolution of species response to environment are the key issues in molecular ecology and evolutionary biology. The direction of adaptive differentiation of species in regions lacking strong selection pressure is usually diverse. However, the driving mechanism of the diverse adaptive differentiation for regional species is still undetermined to date. In this study, we used landscape genomics modelling to infer the adaptive evolution of Cotinus coggygria in China’s warm-temperate zone.
Using fifteen natural populations and nine start codon targeted (SCoT) markers, a total of 1131 unambiguous loci were yielded. Our results showed two genetic groups existed in the fifteen natural populations of C. coggygria, which is due to the divergent selection driven by six environmental factors. Environmental association analyses revealed the environmental variables related to precipitation were associated with high numbers of environment-associated loci.
Our results indicated that the ecological characters of C. coggygria, i.e. avoiding wetness and tolerating drought, determine its adaptive evolution. This study provides a reference that ecological character determines the adaptive evolution of species in regions lacking strong selection pressure.
The adaptive evolution of species response to environment are the fundamental issues in molecular ecology and evolutionary biology. The global climate is changing rapidly and affecting the global ecosystem and biodiversity. Under the pressure of climate change, species either adapt or become extinct . Species adaptation response to rapid climate change can be divided into two. On the one hand, species adapt through migration to adjust the distribution . On the other hand, for species that cannot adapt to climate change through migration, they resort to local adaptation. Thus, these species require adaptive changes, especially in terms of phenology, reproductive behavior, and phenotypic characteristics . Adaptive phenotypic changes for fitting to the changing climate are usually based on the adaptive evolution of species genome . Therefore, the identification of these adaptive genes is the key in understanding species adaptive evolution.
To date, two strategies can be used to identify the adaptive genes. One is the “top-down” strategy, which measures the adaptive phenotype and phenological data using common garden experiment or reciprocal transplant experiment; then, this strategy links these data to genetic variation via genome-wide association studies  or quantitative trait locus . The other is the “bottom-up” strategy, which searches the selected signals of adaptive genetic evolution using genomics scanning and then associates these signals with climatic data to determine the adaptive genes . An example of the latter strategy is landscape genomics, which is less cost and more time efficient than the methods belonging to the former strategy [7, 8].
In recent years, landscape genomics studies have provided information on the interactions between environmental variations and adaptive genetic variations in natural populations . The universal pattern of adaptive evolution is especially popular for regional landscape genomics studies. However, the universal pattern on adaptation needs extreme selection pressure, such as drought of desert and high salinity of ocean, for most species in the same region [10, 11]. However, the strong selection pressure produced by extreme environmental conditions does not exist in most regions. Thus, the direction of adaptive differentiation of regional species is usually diverse [9, 12]. However, the determinant of the diverse adaptive differentiation for regional species is still unclear to date. Determining the reason why these environmental variables drive the comprehensive adaptive differentiation of species genomes is interesting.
In this study, we sampled Cotinus coggygria Scop. (smoke tree) in China’s warm-temperate zone to infer the relationship between environmental variables and adaptive genetic variations in the plant’s genome. This deciduous tree species is widely distributed in China’s warm-temperate zone. This species prefers light; tolerates semi-shade, cold, and drought; and avoids wetness. Although previous population genetics study on C. coggygria had been conducted , the used neutral markers of chloroplast DNA (cpDNA) reflected more demography history events rather than adaptive evolution driving by environmental variations. Here, novel molecular markers, start codon targeted (SCoT) polymorphisms, were used for genome scanning. Start codon targeted polymorphism (SCoT) is a kind of gene targeted marker that was developed based on the conserved region flanking the ATG start codon . SCoT markers have certain advantages, such as simplicity, reproducibility, abundant polymorphism, high throughput and not require priori genomic information . Therefore, this marker is suitable for landscape genomics research.
We employed SCoT markers to detect environment-associated loci (EAL) in response to environmental variations in natural populations of C. coggygria. The present study aimed (i) to reveal the population genetic structure of C. coggygria, (ii) identify the EAL in the genome of C. coggygria, and (iii) detect the key environmental factors that drive the adaptive differentiation of C. coggygria.
Nine SCoT primers were selected in investigating the population genetic structure in C. coggygria. A total of 1131 unambiguous loci were identified with sizes ranging from 60 bp to 1000 bp. The numbers of loci of the nine primers ranged from 93 (SCoT31) to 163 (SCoT2). The number of polymorphic alleles (N A) of each population ranged from 134 (P9) to 325 (P14). The percentage of polymorphic alleles (PPA) of each population ranged from 11.8 (P9) to 28.7 (P14). The level of genetic diversity (H E) of each population ranged from 0.038 (P9) to 0.098 (P4). Summary statistics for the genetic diversity analyses for each population are shown in Table 1.
The Bayesian analysis with all loci of population structure (Fig. 1a) clearly demonstrated that the highest ΔK value (Fig. 2) was obtained when populations were clustered into two groups. One is the East group (P1 to P5), the other is the West group (P6 to P15) (Fig. 3). The non-hierarchical AMOVA (Table 2) revealed that these populations were significantly structured at the species-range scale (F ST = 0.115, P < 0.001). Although two groups were divided, only 5.30% genetic variation occurred among groups (F CT = 0.053, P < 0.05) and most genetic variation occurred within populations (86.06%, F ST = 0.139, P < 0.001). In addition, significant patterns of isolation by distance were detected by comparing F ST values with geographical distances at the species-range scale (r = 0.2452, P < 0.05).
A total of 15 natural populations of C. coggygria and environmental variables were used as subjects and explanatory variables to perform redundancy analysis (RDA). To avoid overestimation of the contribution of environmental variables to population structure, the strong correlated environmental variables were excluded from 19 environmental variables (Table 3). After removing the strong correlated environmental variables, nine remaining environmental variables (Bio3, Bio5, Bio6, Bio7, Bio11, Bio12, Bio14, Bio15 and Bio18) were selected for RDA and environmental association analyses. Figure 4 shows the results of RDA performed using 1131 SCoT allele frequencies as response variables. Correlations of genetic variables with environmental variables in axes 1 and 2 were 0.934 and 0.965, respectively. The ratios of the total eigenvalues of axes 1 and 2 were 37.4% and 15.8%, respectively. RDA showed that six environmental variables were significantly associated with RDA axes 1 and 2 (Table 4), which suggested that the two axes represented more of the changes in the six environmental variables. Among these six environmental variables, isothermality (Bio3), max temperature of warmest month (Bio5), and temperature annual range (Bio7) were related to temperature. Meanwhile, annual precipitation (Bio12), precipitation seasonality (Bio15), and precipitation of warmest quarter (Bio18) were related to precipitation. Bio18 was the highest contributor among the nine environmental variables because of the high ratios of the total eigenvalues of axis 1.
Using the FDIST2 method, 100 outlier loci with a P-value under 0.05 were identified (Additional file 1: Figure S5A). The BayeScan method detected 74 loci as outliers with a log10PO above 0.5; this value is considered a substantial evidence for selection and corresponds to a posterior probability above 0.76 (Additional file 1: Figure S5A). Among the 74 loci detected using the second approach, 27 common loci were detected using the first approach (Additional file 1: Figure S5A). The detected outlier loci all showed a positive alpha value, thereby indicating positive or directional selection. To reduce the false discovery rate, we used the 27 common loci for further environmental association analyses.
The Pearson’s correlation analysis detected thirteen EAL among the 27 outlier loci, which associated with at least one environmental variable (Table 5). This number accounted to 1.14% of the total number of SCoT loci. Among the 13 detected loci, three were significantly related to temperature and precipitation, one was significantly related to temperature, and nine were significantly related to precipitation (Table 5). Among these associated environmental variables, Bio14 and Bio18 were associated with the highest numbers of EAL. LFMM identified seventeen EAL among the 27 outlier loci, which associated with at least one environmental variable (Table 6). This number accounted to 1.50% of the total number of SCoT loci. Among the 17 detected loci, seven were significantly associated with temperature and precipitation, and ten were significantly associated with precipitation (Table 6). Similar to the Pearson’s correlation analysis, LFMM showed that Bio14 and Bio18 were associated with the highest numbers of EAL. However, Bio15 was also identified as environmental variable with high numbers of EAL. A total of 27 EAL with 12 common EAL were detected using the combined two detection methods (Fig. 5b). To further test the contribution of environmental variables to the spatial genetic structure, the Bayesian analysis with all loci except EAL of population structure was performed. After removing the EAL, the two genetic groups also dissolved (Fig. 1b).
This study analyzed the adaptive evolution of C. coggygria to environmental factors through SCoT markers. The overall outlier detection rates in the present study were 8.84% (100 out of 1131) in FDIST and 6.54% in BayeScan. The detection rates are consistent with previously reported rates on landscape genomics studies (2.85% to 10%), such as 2.85% in Alnus glutinosa  and 4.5% in Picea mariana  from SNPs, 9% in Arabis alpina  and 10% in 13 alpine species  from AFLPs, and 4.22% in Cephalotaxus oliveri  from ISSRs. Although SCoT markers are non-neutral biased nature, their detection rate in this study is not significantly higher than that of other molecular markers. Due to lacking DNA sequence information, these loci could not be validated and might be suspected as false-positive loci. In order to minimize false positive rates, we used the loci that were jointly identified by two methods at the same time.
Spatial population genetic structure
Landscape genomics studies have focused mostly on the spatial population genetic structure of species [20, 21]. The influence of environmental variables on population genetic structure has been increasingly revealed. However, addressing spatial population genetic structure contributed by environmental variables is a serious challenge because of the complex reciprocal interactions of multifactors (i.e., gene flow, natural selection, and historical events) [22,23,24]. Our survey on SCoT data demonstrated significant hierarchical population genetic structure across all the studied populations. The two genetic groups, the East group (P1 to P5) and the West group (P6 to P15) are geographically separated. The spatial pattern is quite different from previous findings in Wang et al. . This is mainly due to the difference in markers. The neutral cpDNA markers mainly reflect seed mediated gene flow and population dynamics. SCOT markers can simultaneously reflect seed and pollen mediated gene flow, as well as population dynamics and species adaptive evolution. Our results showing significant IBD pattern in C. coggygria indicated restricted gene flow. However, the restricted gene flow can only increase the degree of genetic differentiation among populations and cannot explain the spatial separation of the two genetic groups. Three hypotheses can be used to explain this large-scale intraspecific genetic disjunction. First, C. coggygria was compressed to two refugia (i.e., gene pools) during climate fluctuations in the past, and the current distribution pattern resulted from the redistribution of the two gene pools. Second, a geographical barrier existed between the two groups, and long-term blocking of gene flow by geographical barrier led to the genetic divergence of the two groups. Third and last, significant environmental differences occurred between the two regions, and these heterogeneous environmental conditions resulted in divergent selection and eventually led to the genetic divergence of the two groups. Previous phylogeographical study on C. coggygria suggested that this species survives in situ and occupies multiple localized glacial refugia rather than compressing to two refugia during the Pleistocene glaciations . Therefore, the first hypothesis failed to explain the genetic divergence of C. coggygria. Assuming that C. coggygria agreed with the second hypothesis, significant genetic differentiation could be expected between the two groups with sufficiently strong geographic barrier. However, our AMOVA analysis based on all loci and STRUCTURE analysis with all loci except EAL (i.e., most of the neutral loci) both showed weak genetic differentiation between the two groups. Thus the second hypothesis was also not appropriate for C. coggygria. Considering that the third hypothesis was appropriate for C. coggygria, a shallow genetic divergence might occur between the two groups because of the interaction between natural selection and gene flow. This expectation is consistent with our detection results (F CT = 0.053). To verify that the differences between the groups were caused by environmental factors, we compared the STRUCTURE analyses with all loci and all loci except EAL. When the EAL driven by environmental variables were excluded, the two genetic groups were also dissolved (Fig. 1a and b). Therefore, the genetic divergence of the two groups is due to the divergent selection driven by environmental factors. To further infer the relative contribution of environmental variables in driving population genetic structure, RDA was performed. The RDA results suggested that six environmental variables related to temperature and precipitation remarkably influenced the spatial population genetic structure. Among these environmental variables, Bio18 was the most important environmental factor in driving population genetic structure. Furthermore, our results showed that these environmental variables could significantly subdivide the populations into two groups (Fig. 4). In general, the third hypothesis holds true for C. coggygria.
EAL driven by environmental factors
Previous phylogeographical studies have confirmed that the warm temperate vegetation adapts to climate change through migration or local adaptation . Plants with long-term local adaptation often face the divergent selection driven by environmental variables, thereby leading to the adaptive evolution of species genome . For some regions under extreme selection pressures, such as desert areas, species undergo convergent evolution at the genomic or phenotypic scale . However, for most regions, such as the distribution area of C. coggygria, extreme selection pressures that drive species convergent evolution do not exist. In this case, species follows a variety of evolution directions. However, the environmental variables that significantly affect the genome and play a decisive role on the direction of species evolution are still unknown.
In this study, we selected C. coggygria as a model in addressing the abovementioned issues. We hypothesized that species ecological characters would drive adaptive evolution and produce large number of EAL. Here, the species ecological characters refer to the sensitivity, adaptability and resistance of species to environmental factors during local adaptation. Thus, examining the ecological characters of C. coggygria is urgently necessary. This species tolerates cold and drought, and avoids wetness. According to the adopted nine environmental variables, min temperature of coldest month (Bio6) and mean temperature of coldest quarter (Bio11) were associated with the ecological character of tolerating cold; precipitation of driest month (Bio14) was associated with tolerating drought, precipitation of warmest quarter (Bio18) were associated with avoiding wetness. In China’s warm-temperate zone, rainy and hot seasons overlap most of the time. Thus, precipitation of warmest quarter (Bio18) is tantamount to precipitation of wettest quarter (Bio16). The results of auto correlation analysis of environmental variables also confirmed this climate characteristic in this region. Therefore, we expected that the large numbers of EAL were associated with the environmental variables associated with the ecological characters, i.e. Bio6, Bio11, Bio14 and Bio18. As expected, the results of the Pearson’s correlation analysis and LFMM both showed that Bio14 and Bio18 were associated with the highest number of EAL. However, the environmental variables related to tolerating cold, Bio6 and Bio11, were not associated with high number of EAL. Our results suggested the environmental variables associated with the ecological characters of tolerating drought and avoiding wetness played more important roles in adaptive evolution in C. coggygria. Thus, most aspects of the characterized EAL of C. coggygria agreed with the hypothesis that ecological characters determine adaptation.
Whether the deduction that ecological character determines adaptation is universal must also be discussed. To date, landscape genomics studies in China’s warm-temperate zone are rare. Therefore, comparison with other species cannot help in confirming the universality of our results in this region. Thus, we reviewed some published works on landscape genomics in other regions in recent years. Prunier et al.  argued that the adaptive SNPs in Picea mariana are related to temperature and correspond to the nature of cold resistance of budset. De Kort et al.  mentioned that the detection of several temperature-dependent SNPs in Alnus glutinosa is related to its resistance to drought. Wang et al.  stated that the ecological characters of Cephalotaxus oliveri in response to temperature and precipitation sensitivity determine its adaptive evolution. Roschanski et al.  reported that the detected adaptive SNPs in Abies alba are associated with winter and drought and correspond to the ecological characters related to drought and cold tolerance. On the basis of this considerable evidence, we suggested that the ecological characters of species might be related to species adaptive evolution. Under high intensity of selection pressure, species usually loses ecological character and evolves toward extreme environments for survival. Without the strong selection pressure, species ecological character appears and evolution develops toward variety to ensure a better survival of species. This study provides a reference that ecological character, i.e. species sensitivity, adaptability and resistance to environmental factors, determines the adaptive evolution of species in regions lacking strong selection pressure.
The differences in the intensity of selection pressure affect the direction of species evolution. In this study, we sampled C. coggygria from China’s warm-temperate zone, a region lacking strong selection pressure. Nine SCoT markers were used to investigate the adaptive genetic variation in C. coggygria. Our results showed that significant hierarchical population genetic structure of C. coggygria is due to the divergent selection driven by environmental factors. The ecological characters of C. coggygria, tolerating drought and avoiding wetness, determine its adaptive evolution. Therefore, species ecological character determines the adaptive evolution of species in regions lacking strong selection pressure.
A total of 156 individuals from fifteen natural populations of C. coggygria were collected from the entire distribution range in China (Fig. 3). Population samples included five to 12 individuals, and each sample was collected at least 10 meters apart. All individuals were collected when the population size was less than ten. Young, healthy leaves were collected and stored in silica gel at room temperature until DNA extraction and genotyping. The geographical coordinates for each sampled population are presented in Table 1.
Genomic DNA was isolated from approximately 30 mg of dried leaves using Plant DNA Extraction Kit DP305 (Tiangen, Beijing, China) following the protocols of the manufacturer. DNA concentration was measured using Microcolume Spectrophotometer ND5000 (BioTeke, Beijing, China). After preliminary screening, nine SCoT primers (SCoT2, SCoT3, SCoT6, SCoT14, SCoT16, SCoT22, SCoT30, SCoT31, and SCoT33) from Collard and Mackill  were selected for polymerase chain reaction (PCR). SCoT2, SCoT16, and SCoT33 were 5′ fluorescent primers labeled with FAM; SCoT3, SCoT22, and SCoT30 were labeled with HEX; SCoT6, SCoT14, and SCoT31 were labeled with TAMRA. PCR was conducted in a 20 μL-reaction mixture consisting of 20 ng template DNA, 1 × reaction buffer (pH 8.3), 0.2 mM dNTPs, 0.3 μM primer, 1 unit of Taq polymerase, and DNA-free water. In an iCycler gene amplification system (Bioteke, Beijing, China), PCR was started with an initial denaturation at 94 °C for 5 min followed by 35 cycles at 94 °C for 40 s, primer-specific annealing temperature (50 °C for SCoT2; 52 °C for SCoT3, SCoT6, SCoT22, SCoT30, and SCoT31; 56 °C for SCoT16 and SCoT33; and 60 °C for SCoT14) for 45 s and 72 °C for 1 min, a final extension at 72 °C for 5 min, and termination by a final hold at 4 °C. PCR products were mixed with 10 μL of HiDi formamide and 0.1 μL of ROX1000 size standard (Applied Biosystems, Foster City, USA). These products were then separated on an ABI 3730 DNA Analyzer at BGI (Beijing, China).
Electropherograms were viewed with GeneMarker 2.2.0 (SoftGenetics, State College, Pennsylvania, USA). To minimize scoring false alleles, peaks between 60 and 1000 bp and heights above 300 relative fluorescent units were scored as a presence (1) or absence (0) matrix for each sample. Subsequent statistical analyses were performed on the basis of this matrix.
AFLPSURV 1.0  was used to calculate the genetic parameters for each population. The estimates included the number of polymorphic alleles (N A), percentage of polymorphic alleles (PPA), gene diversity of Nei (H E) , pairwise Fst between populations, and gene frequencies per allele.
The analysis of hierarchical population structure was implemented using the Bayesian-based program STRUCTURE 2.3.4 . For the analysis, a no-admixture model with independent allele frequencies was selected. K values were tested from 1 to 10, and 10 replicates were performed for each K. Burn-in periods of 1 × 105 and 2 × 104 Monte Carlo and Markov chains were specified. The choice of the optimal value of K was based on the method introduced by Evanno et al. , and this method was implemented in STRUCTURE HARVESTER . Hierarchical and non-hierarchical AMOVA were calculated in ARLEQUIN 22.214.171.124  in inferring the distribution of genetic differentiation at various levels. Mantel tests of isolation-by-distance (IBD) were performed in IBD 3.23  in determining the relation of geographical distance (km) to genetic differentiation (F ST). RDA was conducted using CANOCO 4.5  in disentangling the relative contribution of environmental variables in driving population genetic structure. In RDA, allele frequencies per population (Additional file 2) were used as the response variable and environmental variables (Additional file 3) were used as explanatory variables. Environmental data from 1950 to 2000 at 2.5 arcmin resolution were downloaded from the world climate database (http://www.diva-gis.org/climate). Additional data for each population were extracted using the DIVA-GIS 7.5.0 . In China’s warm-temperate zone, rainy and hot are over the same period, cold and drought are over the same period. Thus, some environmental variables of temperature and precipitation might be significantly acossiated. To avoid overestimation of the contribution of environmental variables to population structure, the strong correlated environmental variables (r > 0.95) are excluded. Auto correlation analysis of environmental variables was performed using Pearson’s regression in SPSS 19 (SPSS Inc., Chicago, IL, USA).
Two approaches were used to identify the outlier loci. The first approach was based on the FDIST2 approach proposed by Beaumont and Nichols  and was implemented in the program Arlequin 126.96.36.199 . The hierarchical island model in Arlequin was selected. The running parameters were set as follows: 100 simulated demes and 20,000 coalescent simulations. The loci outside the 95% confidence interval were regarded as outlier loci. To reduce the false discovery rate, loci with minor allele frequency < 5% were excluded. The second approach was based on the Bayesian approach and implemented in BayeScan 2.01 . The running parameters were set as follows: sample size of 5000, thinning interval of 10, 20 pilot runs with a run length of 5000, and additional burn-in of 50,000 iterations. The loci with posterior probability >0.76 were regarded as outlier loci.
In further detecting the EAL potentially driven by environmental variations, two methods of environmental association analyses were performed. The first method is Pearson’s correlation analysis, which was implemented using SPSS 19 (SPSS Inc., Chicago, IL, USA). This regression analysis ignored the population structure, which might produce a relaxed result of EAL. Similar to RDA, allele frequencies per population were used as the response variable and environmental variables were utilized as explanatory variables. The loci with |r| > 0.50 and P < 0.05 were regarded as EAL. The second method is LFMM, which was implemented in LFMM 1.2 . This Bayesian mixed model considered the population structure, thereby avoiding the bias caused by population history and isolation by distance and producing a robust result of EAL. The running parameters were set as follows: 10,000 sweeps, 1000 burn-in sweeps, and number of latent factors as suggested by STRUCTURE. The loci with |z| > 3 and P < 0.001 were regarded as EAL.
Hoffmann AA, Sgrò CM. Climate change and evolutionary adaptation. Nature. 2011;470(7335):479–85.
Aitken SN, Yeaman S, Holliday JA, Wang TL, Curtis-McLane S. Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol Appl. 2008;1(1):95–111.
Davis MB, Shaw RG. Range shifts and adaptive responses to quaternary climate change. Science. 2001;292(5517):673–9.
Rellstab C, Gugerli F, Eckert AJ, Hancock AM, Holderegger R. A practical guide to environmental association analysis in landscape genomics. Mol Ecol. 2015;24(17):4348–70.
Korte A, Farlow A. The advantages and limitations of trait analysis with GWAS: a review. Plant Methods. 2013;9:29.
Stinchcombe JR, Hoekstra HE. Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits. Heredity. 2008;100(2):158–70.
Lepais O, Bacles CF. Two are better than one: combining landscape genomics and common gardens for detecting local adaptation in forest trees. Mol Ecol. 2014;23(19):4671–3.
Shryock DF, Havrilla CA, Defalco LA, Esque TC, Custer NA, Wood TE. Landscape genomics of Sphaeralcea ambigua in the Mojave Desert: a multivariate, spatially-explicit approach to guide ecological restoration. Conserv Genet. 2015;16(6):1303–17.
Ćalić I, Bussotti F, Martínez-García PJ, Neale DB. Recent landscape genomics studies in forest trees—what can we believe? Tree Genet Genomes. 2016;12(1):3.
Kültz D. Physiological mechanisms used by fish to cope with salinity stress. J Exp Biol. 2015;218(Pt 12):1907–14.
Givnish TJ: Convergent evolution, adaptive radiation, and species diversification in plants. Encyclopedia Evol Biol. 2016;362–373.
Manel S, Gugerli F, Thuiller W, Alvarez N, Legendre P, Holderegger R, Gielly L, Taberlet P, IntraBioDiv C. Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation. Mol Ecol. 2012;21(15):3729–38.
Wang W, Tian CY, Li YH, Li Y. Molecular data and ecological niche modelling reveal the phylogeographic pattern of Cotinus coggygria (Anacardiaceae) in China's warm-temperate zone. Plant Biol. 2014;16(6):1114–20.
Collard BCY, Mackill DJ. Start Codon targeted (SCoT) polymorphism: a simple, novel DNA marker technique for generating gene-targeted markers in plants. Plant Mol Biol Rep. 2009;27(1):86–93.
Deng L, Liang Q, He X, Luo C, Chen H, Qin Z. Investigation and analysis of genetic diversity of diospyros germplasms using SCoT molecular markers in Guangxi. PLoS One. 2015;10(8):e0136510.
De Kort H, Vandepitte K, Bruun HH, Closset-Kopp D, Honnay O, Mergeay J. Landscape genomics and a common garden trial reveal adaptive differentiation to temperature across Europe in the tree species Alnus glutinosa. Mol Ecol. 2014;23(19):4709–21.
Prunier J, Laroche J, Beaulieu J, Bousquet J. Scanning the genome for gene SNPs related to climate adaptation and estimating selection at the molecular level in boreal black spruce. Mol Ecol. 2011;20(8):1702–16.
Poncet BN, Herrmann D, Gugerli F, Taberlet P, Holderegger R, Gielly L, Rioux D, Thuiller W, Aubert S, Manel S. Tracking genes of ecological relevance using a genome scan in two independent regional population samples of Arabis alpina. Mol Ecol. 2010;19(14):2896–907.
Wang T, Wang Z, Xia F, Su YJ. Local adaptation to temperature and precipitation in naturally fragmented populations of Cephalotaxus oliveri, an endangered conifer endemic to China. Sci Rep. 2016;6:25031.
Schoville S, Bonin A, Francois O, Lobreaux S, Melodelima C, Manel S. Adaptive genetic variation on the landscape: methods and cases. Annu Rev Ecol Syst. 2012;43:23–43.
Hall LA, Beissinger SR. A practical toolbox for design and analysis of landscape genetics studies. Landsc Ecol. 2014;29(9):1487–504.
Hewitt GM. Genetic consequences of climatic oscillations in the quaternary. Philos Trans Roy Soc Lond B Biol Sci. 2004;359(1442):183–95.
Nybom H. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants. Mol Ecol. 2004;13(5):1143–55.
Ohsawa T, Ide Y. Global patterns of genetic variation in plant species along vertical and horizontal gradients on mountains. Glob Ecol Biogeogr. 2008;17(2):152–63.
Fu ZZ, Li YH, Zhang KM, Li Y. Molecular data and ecological niche modeling reveal population dynamics of widespread shrub Forsythia suspensa (Oleaceae) in China’s warm-temperate zone in response to climate change during the Pleistocene. BMC Evol Biol. 2014;14:114.
Azua-Bustos A, González-Silva C, Arenas-Fajardo C, Vicuña R. Extreme environments as potential drivers of convergent evolution by exaptation: the Atacama Desert coastal range case. Front Microbiol. 2012;3(2):426.
Roschanski AM, Csilléry K, Liepelt S, Oddou-Muratorio S, Ziegenhagen B, Huard F, Ullrich KK, Postolache D, Vendramin GG, Fady B. Evidence of divergent selection for drought and cold tolerance at landscape and local scales in Abies alba mill. In the French Mediterranean alps. Mol Ecol. 2016;25(3):776–94.
Vekemans X, Beauwens T, Lemaire M, Roldán-Ruiz I. Data from amplified fragment length polymorphism (AFLP) markers show indication of size homoplasy and of a relationship between degree of homoplasy and fragment size. Mol Ecol. 2002;11(1):139–51.
Nei M. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci U S A. 1973;70(12):3321–3.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005;14(8):2611–20.
Earl DA, BM vH. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol Ecol Resour. 2010;10(3):564–7.
Jensen JL, Bohonak AJ, Kelley ST. Isolation by distance, web service. BMC Genet. 2005;6:13.
Ter Braak CJF, Smilauer P. CANOCO reference manual and CanoDraw for windows user's guide: software for canonical community ordination (version 4.5). New York: Microcomputer Power; 2002.
Hijmans RJ, Guarino L, Cruz M, Rojas E. Computer tools for spatial analysis of plant genetic resources data: 1. DIVA-GIS. Plant Genet Resour Newsl. 2001;127:15–9.
Beaumont MA, Nichols RA. Evaluating loci for use in the genetic analysis of population structure. P Roy Soc Lond B-Biol Sci. 1996;263(1377):1619–26.
Foll M, Gaggiotti OE. A genome scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics. 2008;180(2):977–93.
Frichot E, Schoville SD, Bouchard G, Francois O. Testing for associations between loci and environmental gradients using latent factor mixed models. Mol Biol Evol. 2013;30(7):1687–99.
We thank professor Xue-Jun Ge (South China Botanical Garden, Chinese Academy of Sciences) for suggestions and two anonymous reviewers for their input.
This work was supported by the National Natural Science Foundation of China (31770225), the Henan Agricultural University Science & Technology Innovation Fund (KJCX2016A2), the Funding Scheme of Young Backbone Teachers of Higher Education Institutions in Henan Province (2015GGJS-081), and the Key Scientific Research Projects of Henan Higher School (16A220002).
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
Ethics approval and consent to participate
Consent for publication
The authors declare that have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.