Landscape genomics reveal that ecological character determines adaptation: a case study in smoke tree (Cotinus coggygria Scop.)
© The Author(s). 2017
Received: 27 June 2017
Accepted: 18 August 2017
Published: 23 August 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.
Details of population locations, sample size, genetic diversity of 15 population for C. coggygria
Population no. and code
Wudang Mt., Hubei
Song Mt., Henan
Jiulian Mt., Henan
Yuan Mt., Shandong
Lingkong Mt., Shanxi
Laojun Mt., Henan
Taibai Mt., Shaanxi
Tiantai Mt., Shaanxi
Laojun Mt., Shaanxi
Yuntai Mt., Henan
Hua Mt., Shaanxi
Wuzhi Mt., Hebei
Tianlong Mt., Shanxi
Hierarchical AMOVAs for SCOT variation surveyed in C. coggygria
Source of variation
F ST = 0.115
P < 0.001
F ST = 0.053
P < 0.001
F ST = 0.117
P < 0.001
Among two groups
F CT = 0.053
P < 0.05
F SC = 0.091
P < 0.001
F ST = 0.139
P < 0.001
Nineteen environmental variables used in this study
Bio1: Annual mean temperature (°C × 10)
Bio2: Mean diurnal range (Mean of monthly (max temp - min temp))
Bio3: Isothermality (Bio2/Bio7) (×100)
Bio4: Temperature seasonality (standard deviation ×100)
Bio5: Max temperature of warmest month (°C × 10)
Bio6: Min temperature of coldest month (°C × 10)
Bio7: Temperature annual range (E5-E6)
Bio8: Mean temperature of wettest quarter (°C × 10)
Bio9: Mean temperature of driest quarter (°C × 10)
Bio10: Mean temperature of warmest quarter (°C × 10)
Bio11: Mean temperature of coldest quarter (°C × 10)
Bio12: Annual precipitation (mm)
Bio13: Precipitation of wettest month (mm)
Bio14: Precipitation of driest month (mm)
Bio15: Precipitation seasonality (coefficient of variation)
Bio16: Precipitation of wettest quarter (mm)
Bio17: Precipitation of driest quarter (mm)
Bio18: Precipitation of warmest quarter (mm)
Bio19: Precipitation of coldest quarter (mm)
Correlations between environmental variables and the ordination axes
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 EAL as indicated by Pearson’s correlation coefficients
Pearson’s correlation coefficients
The EAL as indicated by |z|-score
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–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 126.96.36.199  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 188.8.131.52 . 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.
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.
YL conceived and designed the experiments; CYM and JY performed the experiments; CYM, JY and RLM analyzed the data; YL wrote the paper. All authors read and approved the final manuscript.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hoffmann AA, Sgrò CM. Climate change and evolutionary adaptation. Nature. 2011;470(7335):479–85.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Davis MB, Shaw RG. Range shifts and adaptive responses to quaternary climate change. Science. 2001;292(5517):673–9.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Korte A, Farlow A. The advantages and limitations of trait analysis with GWAS: a review. Plant Methods. 2013;9:29.View ArticlePubMedPubMed CentralGoogle Scholar
- Stinchcombe JR, Hoekstra HE. Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits. Heredity. 2008;100(2):158–70.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ć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.View ArticleGoogle Scholar
- Kültz D. Physiological mechanisms used by fish to cope with salinity stress. J Exp Biol. 2015;218(Pt 12):1907–14.View ArticlePubMedGoogle Scholar
- Givnish TJ: Convergent evolution, adaptive radiation, and species diversification in plants. Encyclopedia Evol Biol. 2016;362–373.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- Hall LA, Beissinger SR. A practical toolbox for design and analysis of landscape genetics studies. Landsc Ecol. 2014;29(9):1487–504.View ArticleGoogle Scholar
- Hewitt GM. Genetic consequences of climatic oscillations in the quaternary. Philos Trans Roy Soc Lond B Biol Sci. 2004;359(1442):183–95.View ArticleGoogle Scholar
- Nybom H. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants. Mol Ecol. 2004;13(5):1143–55.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Nei M. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci U S A. 1973;70(12):3321–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Jensen JL, Bohonak AJ, Kelley ST. Isolation by distance, web service. BMC Genet. 2005;6:13.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar