Localized environmental heterogeneity drives the population differentiation of two endangered and endemic Opisthopappus Shih species
BMC Ecology and Evolution volume 21, Article number: 56 (2021)
Climate heterogeneity not only indirectly shapes the genetic structures of plant populations, but also drives adaptive divergence by impacting demographic dynamics. The variable localized climates and topographic complexity of the Taihang Mountains make them a major natural boundary in Northern China that influences the divergence of organisms distributed across this region. Opisthopappus is an endemic genus of the Taihang Mountains that includes only two spatially partitioned species Opisthopappus longilobus and Opisthopappus taihangensis. For this study, the mechanisms behind the genetic variations in Opisthopappus populations were investigated.
Using SNP and InDel data coupled with geographic and climatic information, significant genetic differentiation was found to exist either between Opisthopappus populations or two species. All studied populations were divided into two genetic groups with the differentiation of haplotypes between the groups. At approximately 17.44 Ma of the early Miocene, O. taihangensis differentiated from O. longilobus under differing precipitation regimes due to the intensification of the Asian monsoon. Subsequently, intraspecific divergence might be induced by the dramatic climatic transformation from the mid- to late Miocene. During the Pleistocene period, the rapid uplift of the Taihang Mountains coupled with violent climatic oscillations would further promote the diversity of the two species. Following the development of the Taihang Mountains, its complex topography created geographical and ecological heterogeneity, which could lead to spatiotemporal isolation between the Opisthopappus populations. Thus the adaptive divergence might occur within these intraspecific populations in the localized heterogeneous environment of the Taihang Mountains.
The localized environmental events through the integration of small-scale spatial effects impacted the demographic history and differentiation mechanism of Opisthopappus species in the Taihang Mountains. The results provide useful information for us to understand the ecology and evolution of organisms in the mountainous environment from population and species perspective.
Understanding the processes that drive differentiation between populations and elucidating the mechanisms that underlie the origins and maintenance of genetic variations are major aims and fundamental tasks in evolutionary biology [1,2,3,4,5], which are also core issues in conservation biology [6, 7]. Myriad factors may impact the evolution and genetic differentiation of plant populations, where geological events and climate oscillations have been suggested as critical drivers [8,9,10]. In terms of geological events, mountain uplifts lead to complex topographies that can segregate large plant populations into multiple smaller sub-populations and enhance differentiation between species or populations through geographic isolation. Further, climatic oscillations can shift the ranges of species, resulting in novel environments with increased variability [11, 12]. To adapt to different environments, organisms evolve corresponding phenotypic variations and genetic differentiation [13, 14].
During this process, intensifying climate change has left an indelible imprint on the composition and divergence of populations or species [15, 16], which further greatly influenced the distribution patterns and shaped the genetic structures of populations [5, 17]. In general, the geographic processes of mountainous regions may influence the genetic makeup of plant populations over large spatial scales, whereas ecological processes from climate change may impact the genetic structures of plant populations at small spatial scales [5, 18,19,20].
The Taihang Mountains, with a north–south orientation (36–40 ºN, 112–115 ºE), are a prominent natural boundary in Northern China [21, 22], which have a geological developmental history of more than 2.5 billion years with a typical platform type crustal structure spanning the Mesoproterozoic to Paleozoic Eras. Their distinct geotectonic positioning has produced a unique geological and geomorphic landscape. The Southern Taihang Mountains have existed as a major boundary of neotectonic deformation, represented by the Yuntai Landform [21,22,23]. The Northern Taihang Mountains, represented by Zhangshiyan, Cnagyan, and Linlv Landforms, are higher than their southern counterparts with an average elevation of 1500 m [21, 24, 25].
Climatically, the southern region is home to a warm temperate semi-humid climate with a mean annual temperature of 12.7 °C and precipitation of 606.4 mm, while the northern region has a temperate continental monsoon climate with a mean annual temperature of ~ 10 °C and precipitation of 700 mm. The topographic complexity of the Taihang Mountains coupled with increasing climate variability have significantly impacted many organisms [21, 24]. Being an important germplasm resource, Pyrus betulaefolia exhibits abundant genetic diversity and variation, which might have been derived from the diversified environments of different populations located in the Taihang Mountains . In contrast Episyrphus balteatus presented a non-obvious phylogeographical structure, which resulted from invalid geographical barriers of the Taihang and Yashan Mountains, where its population division was driven by the climatic changes following the uplift of the Taihang and Yashan Mountains . The Taihang Mountains have been regarded as a distribution and diversity center for numerous genera [21, 24].
As an important perennial herbal germplasm resource of Asteraceae (Opisthopappus), which has been listed as a second-class protected plant in China , grows only on the steep slopes and cliffs of the Taihang Mountains that span Shanxi, Hebei, and Henan Provinces [28,29,30]. Being a diploid species (2n = 18) , comprised of Opisthopappus taihangensis and Opisthopappus longilobus [25, 31, 32], this genus is endemic in China and possesses significant ornamental and medicinal value [25, 32, 33]. Between O. taihangensis and O. longilobus, the morphological distinctions are primarily manifested in the leaves and bracts. For the former, there is appressed puberulent on both surfaces of the leaves, two pinnatisect stem leaves, and no bracteal leaves. For the latter, which has hairless leaves, there is one pinnatisect, except the basal stem leaves, and a pair of bracteal leaves beneath the involucres. Furthermore, several other categories of morphological differentiation have been observed, such as the leaf pinnatisect, sparsely pubescent or glabrous surfaces, stoma size and density, pollen colpus depth, and ostiole density [32, 34, 35].
Meanwhile, genetic variations emerged between these two species, whether in the forms of nuclear molecular markers or chloroplast gene sequences [36,37,38,39,40]. Remarkably, phylogenetic analyses based on chloroplast microsatellites (cpSSR) and sequence related amplified polymorphisms (SRAP) revealed that neither O. longilobus nor O. taihangensis formed a monophyletic clade [37, 38]. In particular, some populations of O. longilobus always integrated with O. taihangensis populations. Nevertheless, the interspecific hybridization of these two species have never been reported.
In previous researches, it was revealed that geographical distance had significant correlations with genetic differentiation among populations across Opisthopappus species [37, 38]. However, precisely how the climatic heterogeneity of the Taihang Mountains influenced genetic differentiation in populations of this genus had not been addressed as yet. Consequently, we hypothesized that the differentiation between these two species would be a hierarchically comprehensive process that might be initially impacted by climate shifts, subsequently by the geographical topography of China, and finally by the environmental heterogeneity of the Taihang Mountains.
For the present study (according to the hypothesis above), the roles and influences of environmental factors on species and population differentiation were investigated through the combination of geographic and climatic data, using single nucleotide polymorphisms (SNP) and insertion-deletion (InDel) markers of nuclear genes developed by Chai et al.  based on transcriptome data. The aims of this study were to: (i) analyze the genetic variations between species and between all studied populations; (ii) investigate the evolutionary processes and histories of O. taihangensis and O. longilobus species; (iii) estimate the effects of geographical and climatic variables; (iv) identify the possible key environmental factors that drive this differentiation. These results evaluated the demographic dynamics of Opisthopappus during the evolutionary process, explored the underlying mechanisms of inter-/intra-species differentiation, and provided some clues for the investigation of additional plant species in the Taihang Mountains.
Genetic variation of Opisthopappus populations
Among the SNP primers, eight pairs produced repeatable, clear, and stable bands. The total length of eight SNP combination segments was 1921 bp, which contained 1870 conservative sites and 51 polymorphic sites. Based on one hundred twenty sequences, seventy-five haplotypes were identified (Table 1). Therein, 47 haplotypes (H1–H47) were detected in O. longilobus and 28 haplotypes (H48–H75) in O. taihangensis. No shared haplotypes were detected between O. longilobus and O. taihangensis. For O. longilobus, the H5 haplotype was the most widely distributed, which was shared by three populations. Five haplotypes (H2, H3, H8, H22, and H30) were shared by two populations, whereas the other 41 haplotypes were distributed only among a single population (Table 1, Fig. 1). For O. taihangensis, the H50 and H52 haplotypes were the most widely distributed (both shared by six populations), followed by H56 (shared by five populations), H53 (shared by four populations), and H51, 54, 55, 60 (shared by two populations). The remaining 20 haplotypes were detected only in a single population.
The value of the genetic differentiation coefficient NST was 0.743 across species for the haplotypes with SNP data, which was significantly larger than the value of GST = 0.105 (P < 0.05). This indicated that there were significant phylogeographic structures in the Opisthopappus genus.
Both the SNP sequences and InDel data revealed a high genetic diversity in Opisthopappus (Table 1). For O. longilobus, the haplotype diversity (Hd) was 0.99327 and the total nucleotide polymorphism (π) was 0.00308 according to the SNP data. As relates to the InDel data, the Nei’s gene diversity index (H), polymorphic loci ratio (PPL), and Shannon’s polymorphism information index (I) were 0.1695, 76.92%, and 0.2777, respectively. For O. taihangensis, the genetic indices, Hd, π, H, PPL and I, were 0.95399, 0.00178, 0.1592, 73.85% and 0.2585 respectively (Table 1).
Significant genetic variations occurred either between Opisthopappus populations or two species (Table 2). Based on the SNP sequences, 80% of the mutations was found between O. longilobus and O. taihangensis (FCT = 0.8003, P < 0.01), 15% of the molecular variations within the populations (FST = 0.8460, P < 0.01), and only 5% of the molecular variations between populations within species (FSC = 0.2287, P < 0.01). For the InDel data, the genetic variation distribution trend was similar to the SNP sequences. The results verified that molecular variations existed primarily between the two Opisthopappus species (Table 2).
A maximum likelihood (ML) phylogenetic tree was constructed based on eight SNP combination fragments (Additional file 1: Fig. S1A), which revealed that all individuals were clearly divided into two groups corresponding to two species. Further, UPGMA cluster analysis (Additional file 1: Fig. S1B) performed based on Nei's genetic distance revealed that twenty-four populations were separated into two (O. longilobus and O. taihangensis) clusters.
For the structural analysis, when △K (mean (|L’(K)|/sd(L(K))) attained a maximum value, K = 2 was taken on both SNP and InDel data. The most significant possibilities were gathered into two groups (Fig. 2). When K ranged from 3 to 6, the genetic structural pattern was similar to that when K = 2, only more mixed individuals presented within O. longilobus or O. taihangensis populations (Additional file 2: Fig. S2).
By the DAPC analysis, the conserved first ten principal components (PCs) represented a 91.1% variation of the total genetic components. The first three linear discriminant functions (LDs) explained 56.9%, 7.5%, and 7.5% of the eigenvalues of the remaining PCs. All O. longilobus and O. taihangensis populations could be thoroughly separated based on LD1 & LD2. From the plotting of LD2 & LD3, the two species were seemingly not well separated, which was attributed to the limited eigenvalues of variation proportion contained by LD2 and LD3 (Fig. 3). The results were confirmed by the K-W test (Additional file 3: Fig. S3).
Historical dynamics of Opisthopappus populations
A Bayesian inference tree was developed for seventy-five haplotypes (Fig. 4), which were segregated within two distinct branches. One branch contained the H1-H47 haplotypes from O. longilobus, whereas the other contained the H48-H75 haplotypes from O. taihangensis. Further, the haplotype network (Additional file 4: Fig. S4) presented a similar haplotype distribution pattern within the haplotype phylogenetic tree (Fig. 4).
At the Miocene-Oligocene boundary (~ 25.08 Ma, 95% HPD: 22.14–27.92 Ma), the Opisthopappus genus diverged from the outgroup (Fig. 4). In early Miocene, approximately 17.44 Ma (95% HPD: 15.58–19.43 Ma), Opisthopappus began to differentiate into two major lineages (O. longilobus and O. taihangensis). Within O. taihangensis, the differentiation time was at 12.82 Ma (95% HPD: 10.96–14.79 Ma) between haplotypes, while the haplotypes of O. longilobus were differentiated at ~ 13.67 Ma (95% HPD: 11.82–15.58 Ma). From the later Miocene to Pliocene, the intraspecies divergence continuously occurred, where the approximate more recent differentiation time of the intraspecific haplotypes for both species was during the Quaternary Era (e.g. H1 and H17, H53 and H57).
The neutral test and mismatch distribution analysis (MDA) (Additional file 5: Table S1) suggested that both Opisthopappus genus and two species had experienced a recent expansion based on significantly negative Fu's Fs (− 25.2963, − 18.5566, and − 24.1000 for O. longilobus, O. taihangensis, and Opisthopappus, respectively, P < 0.05) values and a non-significant sum of square deviation (SSD) and raggedness index (Rag) values (P > 0.05).
According to ABC analysis, there was a significant difference between the observed and simulated data based on the posterior distributions of all scenarios. Scenario 3, namely, O. longilobus was ancestral, and O. taihangensis differentiated from O. longilobus, was considered to be the most unambiguously supported evolutionary model with the highest posterior probability under direct estimate (0.3900, 95% CI 0.0000–0.8175) and logistic regression tests (0.7187, 95% CI 0.7093–0.7282). Furthermore, a 95% CI of the PP for this scenario did not overlap with other scenarios under logistic estimation.
The PCA of the posterior distributions of model checking analysis revealed that the summary statistics from the observed data produced eigenvectors that were within, or at the margins of 1000 simulated pseudo-observed data sets (PODs), which indicated that scenario 3 was generally suitable for the observed data (Fig. 5). For Scenario 3, the low values of type I (direct estimate: 0.182; logistic estimate: 0.210) and type II errors (direct estimate: 0.094; logistic estimate: 0.093) were obtained based on 500 PODs. Under this best supported model, the estimated effective population sizes of O. longilobus and O. taihangensis were 2.21E + 04 (95% CI 4.22E + 03–1.24e + 05) and 7.40E + 04 (95% CI 1.36E + 04–3.51E + 05).
The historical gene flow generated using MIGRATE were low between the two species, Nm O. longilobus → O. taihangensis = 0.3813 (95% CIs 0–1.8790) and Nm O. taihangensis → O. longilobus = 0.7860 (95% CIs 0–3.1748). The mean contemporary gene flow (m) between the two species was also low. Migration occurred from O. longilobus to O. taihangensis and in turn was 0.0059 (95% CIs 0–0.0173) and 0.0052 (95% CIs 0–0.0156), respectively. However, the migration rates within each species were relatively high (O. longilobus: 0.9941 (95% Cis 0.9827–1.0055), O. taihangensis: 0.9948 (95% CIs 0.9844–1.0052)) as estimated by BAYESASS.
Influences of environmental heterogenicity on Opisthopappus populations
One-way ANOVA following the extraction of bioclimatic variables found that most of them, distributed along the two species were significantly different, including Mean diurnal range (bio2), Isothermality (bio3), Temperature seasonality (bio4), Min temperature of coldest month (bio6), Temperature annual range (bio7), Mean temperature of driest quarter (bio9), Mean temperature of coldest quarter (bio11), Precipitation of wettest month (bio13), Precipitation of driest month (bio14), Precipitation seasonality (bio15), Precipitation of wettest quarter (bio16), Precipitation of driest quarter (bio17), and Precipitation of coldest quarter (bio19) (Additional file 6: Table S2). The partial correlation of bioclimatic variables via PCA revealed that the explanatory direction was different for these variables. The PCA plot drawn on the first two axes explained 55.25% and 31.85% of the variations in the climate variables, respectively (Fig. 6).
To assess whether geographic or environmental differences might drive genetic differentiation, Barrier analysis, Mantel and partial Mantel tests were conducted. Among the Opisthopappus populations, geographical barriers existed (Fig. 1). In particular, more significant barriers were found between the boundary populations of O. longilobus and O. taihangensis (Fig. 1), which reflected a pattern of geographical isolation.
Significant associations between geographical and environmental distances appeared (r = 0.3588, P = 0.002) across all populations of Opisthopappus. Moreover, significant correlations were discovered in both the genetic and geographic distances (r = 0.5039, P = 0.001), and genetic and environmental distances (r = 0.3132, P = 0.001). The partial Mantel tests also detected significant correlations between genetic and environmental distances conditioned on geographic effects (r = 0.1631, P = 0.011).
Similar results, which revealed that geographic and environmental distances influenced genetic distances, were obtained by multiple matrix regression with randomization (MMRR) analyses. The effects were that both the geographic (coefficient = 0.2189, r2 = 0.2944, P = 0.001) and environmental (coefficient = 0.2060, r2 = 0.0944, P = 0.001) distances significantly related to the genetic distance. The joint effects of both the geographical and environmental distances also significantly impacted the genetic distance (r2 = 0.3108, coefficient geo = 0.1955, P = 0.001, coefficient env = 0.0751, P = 0.011).
The above results indicated that the genetic differentiation of populations across the two species was significantly and linearly correlated with geographic and/or climatic differentiation. Scatterplots were subsequently constructed to show further details of the relationships between the genetic, geographical, and environmental distances (Fig. 7).
The full RDA model, including geographic distribution and climatic factors, explained 72.42% (conditioned: 12.22%, constrained: 60.20%, Table 3) of the variations between the genetic components. The partial RDA, while conditioned on the geographic distribution (coordinates) of sites, found a significant climate variable effect following the removal of the isolation by distance effects (Proportion = 60.20%, adj R2 = 0.5497, P = 0.03). The ANOVA indicated that the PC2 of bioclimatic variables significantly explained the genetic components with the highest explanatory proportions (Proportion = 56.75%, adj R2 = 0.5451, P = 0.01).
The distribution of the PC1-3 of climatic variables along the ordination axis was further examined by GLM (Table 4). The PC2 (adj R2 = 0.1687, F = 5.6682, P = 0.0264) and PC3 (adj R2 = 0.1171, F = 4.0500, P = 0.0456), had significant F statistics. PC2 correlated significantly with the ordination axis1 of dbRDA, while PC3 was significantly correlated with axis2. The high adjusted R2 indicated that PC2 and PC3 were sufficiently explained by the two dbRDA axes. Consistent influence estimates were obtained even when the order of predictors was altered; thus, the two variables were highly relevant for explaining the genetic differentiation of the populations across the two species. However, the adjusted R2 of PC1 might be too small to be meaningful in either dbRDA or GLM.
Owing to its high explanatory proportions, the scatter and ordisurf plots (Fig. 8) for PC2 were further developed. The results congruently revealed broader ranges of environmental contours in O. taihangensis populations in contrast to O. longilobus populations. The populations of the two species demonstrated a significantly different distribution of the dbRDA space along axis 1, but not axis 2, which were similar to the DAPC clustering patterns. The distribution of populations within species was too close to be distinguished. The ordisurf plots clearly illustrated the climatic differentiation between O. longilobus and O. taihangensis.
Genetic differentiation in localized heterogeneous environments of the Taihang Mountains
Abundant genetic variations were observed in O. longilobus or O. taihangensis (Table 1), which was consistent with previous researches [36,37,38]. As an ancestor of O. taihangensis, O. longilobus had a relatively higher genetic diversity than O. taihangensis (Table 1, Fig. 5) . The perenniality and insect-pollination of these two species likely led to high genetic diversity, enabling them to persist across a range of environmental conditions on the Taihang Mountains.
The studied Opisthopappus populations could be well split into O. longilobus and O. taihangensis (Fig. 2–3), which was slightly inconsistent with the previous results [36,37,38]. This might be from different molecular markers used. The markers used in the previous researches might reveal less genetic information than the markers used in this study that obtained across the whole genome of Opisthopappus. The EST-SSR markers developed from the transcriptome data also divided Opisthopappus populations into two clusters .
Overall genetic variations existed primarily among populations within Opisthopappus species (Table 2, Fig. 2–3, Additional file 3: Fig. S1-4). The genetic differentiation among populations or between species might be an outcome contributed to by both geography and environment (Figs. 1 and 7) [44,45,46,47].
With peaks reaching ~ 2800 m in elevation, the Taihang Mountains have many deep gullies and immense valleys. Its complex topography serves as a significant geographical barrier in Northern China , which can cause distributional gaps between species, leading not only to disjointed populations but also to cascading effects associated with interrupted gene flow and habitat heterogeneity . Following limited gene exchange, O. taihangensis and O. longilobus, which are distributed only on the Taihang Mountains, were segmented into spatially isolated subpopulations (Figs. 1 and 7). Contemporary gene flow between the two species was found to be fairly low, and migration occurred primarily within the species. With the decline in gene flow, population differentiation would increase due to genetic drift or/and local adaptation effects under the conditions of heterogeneous habitats [47, 49].
From north to south, the weather of the Taihang Mountains changes from a temperate continental monsoon climate to a warm temperate semi-humid climate [21, 25]. The habitat for each population of Opisthopappus was diversified due to variable localized climates, such as temperature and precipitation.
Commonly, temperature and precipitation were found to play prominent roles as selective drivers for the variations in various plant species [17, 50,51,52]. Significant temperature fluctuations could contribute to the physiological states, metabolic levels, and genetic alterations of plants, to further drive the genetic makeup of populations. Meanwhile, precipitation during seed germination and growth would impact the demographic size as well as influencing successful seed colonization.
According our field investigations, various sample sites had different climatic conditions . These differences in localized climates lead to the environmental isolation among Opisthopappus populations (Fig. 7). A driver of this pattern may be attributed to a neutral process of temporally disrupted gene flow between individuals growing in environmentally distinct habitats [53, 54]. In turn, these environmentally distinct habitats can serve as a barrier to gene flow, causing genetic differentiation between spatially close populations . Consequently, these would promote isolated subpopulations to eventually evolve into genetically distinct lineages, while adapting to local environments [16, 17, 55].
In the Taihang Mountains, all extracted climatic variables might well be employed to explore weather differentiations between O. taihangensis and O. longilobus (Fig. 6, Additional file 6: Table S2). The PC2 of bioclimatic variables possessed an overwhelming explanatory ability toward genetic differentiation (Table 3, Fig. 8). Thereinto, bio6 and bio13 contained the first two longest projections with the rotation values in PC2 (Fig. 6). Bio6, the minimum temperature of the coldest month belonged to the temperature dimension, while bio13, the precipitation of the wettest month, belonged to the precipitation dimension.
Single or simple environmental variations could typically initiate the adaptive divergence between populations, followed by the expansion of accumulating genetic differentiation [5, 50, 56]. The significantly distinct distribution conditions and patterns of the two climate variables in the Taihang Mountains, corresponding to O. longilobus and O. taihangensis, might be regarded as original primordial indicators with driving forces toward causing and promoting the genetic differentiation and diversification of the two species.
Evolutionary demographic dynamics from climatic transformation and topographic events
The most recent common ancestor of Asteraceae originated 76–66 Ma ago . According to the fossil record, the Asteraceae family occurred at the Eocene–Oligocene boundary (42–49 Ma) [58, 59]. Zhao  pointed out that the Anthemideae tribe of Asteraceae originated in Eurasia, which then gradually dispersed eastward toward Asia. Some initial taxa of the Artemisiinae subtribe in central Asia evolved separately into Dendranthema and Artemisia groups. Within the Artemisia group, the oldest fossil pollen of the Artemisia genus was recorded at Eocene–Oligocene boundary [61, 62]. In the Dendranthema group, Opisthopappus is a relatively close taxon with Ajania, being the primary ancestral genus . In our study, the temporal divergence between Opisthopappus and the outgroup Chrysanthemum indicum occurred at the Oligocene–Miocene boundary (25.08 Ma, 95% HPD: 22.14–27.92 Ma, Fig. 4). This period was coincidental with the evolution of the Dendranthema group.
Due to the collision of the India Plate with Eurasia, the Qinghai-Tibetan Plateau (QTP) began to lift during the Eocene (50–55 Ma) and then experienced different stages of growth to attain its current elevation [63,64,65,66,67,68,69,70]. The formation of the QTP dramatically modified the global and East Asian climate [70, 71] and triggered the Asian monsoon that had a significant impact on the weather in China . The climatic systems of China were transformed from the original planetary to the monsoon system during the early Miocene, which created profound ecological heterogeneity . It may be that this ecological heterogeneity drove the divergence of Opisthopappus from Asteraceae.
In the Burdigalian of Miocene (17.44 Ma, 95% HPD: 15.58–19.43 Ma), the Opisthopappus genus began to differentiate from O. longilobus to O. taihangensis (Fig. 4), which coincided with the second progressive and heterogeneous uplift of the QTP (15–13 Ma) [72,73,74]. During this period, the Asian monsoon was intensified owing to the extensive uplifting of the QTP. The monsoon characteristics of different areas resulted in the segregation of species populations .
During the early mid-Miocene period, the Asian monsoon was further enhanced, while global cooling events occurred in the ambient ocean and atmosphere . However, following cooling events, the global temperature rapidly returned and the Middle Miocene Climate Optimum period (MMCO) emerged (14.7–16.9 Ma) . Subsequently, as the Antarctic ice sheet expanded, the global weather changed from MMCO to a colder period. These climatic shifts served as a stimulus to promote the divergence of different plant populations . Under this scenario, O. longilobus or O. taihangensis began to undergo intraspecific differentiation (Fig. 4), from 14.3 to 11.0 Ma during the Serravallian of Miocene periods.
Later in the Miocene (5.30–11.0 Ma) period, the continuously cooled global climate and the progressive extended QTP uplift brought the topographic and vegetational changes in China [79, 80]. In Eastern China, different weather systems (e.g., tropical humid, subtropical humid, warm temperate humid, and temperate humid) were gradually established from south to north . The uplift process coupled with climatic changes initiated habitat diversification as well as that of the two species of Opisthopappus. This was verified by not only the divergence period of different intraspecific haplotypes (8.42–11.22 Ma), but the diverse haplotypes in each species (Fig. 4).
Toward the emergence of the Taihang Mountains the stage was being set for a neotectonics movement, from the late Miocene to early Pliocene . The intermittent activity of the QTP, from 3.5 to 1.6 Ma (late Miocene to Pleistocene periods), drove the rapid uplift of Taihang Mountains during the Pleistocene. And the Taihang Mountain regions were within the ranges of monsoonal system at that time [82, 83]. The monsoon system interacted with the interglacial cycle to produce a more variable monsoon climate during the Pleistocene period [70, 83, 84]. The habitats of O. longilobus and O. taihangensis became ever more fragmented, where originally large and continuous populations may be separated into multiple smaller subpopulations. Thus, the most recent diversification of haplotypes within species occurred during the Pleistocene of the Quaternary , which also implied divergent selection between environments and localized adaptations to their respective habitats (Figs. 1 and 4).
Additionally, both O. taihangensis and O. longilobus exhibited signs of recent expansion during the evolutionary process (Additional file 5: Table S1). This was confirmed by the distribution of the Opisthopappus haplotypes (Fig. 1) and network (Additional file 4: Fig. S4). During the Quaternary period, the paleovegetation of Taihang Mountains repeatedly appeared as replacement species during grassland and forest cycles . Emerging grasslands might have served as a transitional corridor that provided opportunities for populations to expand and colonize.
Implication of the outcome from this research
Spatial environmental heterogeneity is typically proposed as a critical driver that leads to population differentiation, and even the acceleration of speciation. Here, we provided comprehensive evidence, including genetics, geographical conditions, climate variables, and evolutionary processes to interpret the differentiation of two Opisthopappus species.
Based on the above results, the divergence and intraspecies variations of Opisthopappus primarily resulted from climate fluctuations, the intensification of Asian monsoon, and the topographic complexity of China with the extensive uplift of the QTP. Subsequently, the ecological stratification and environmental heterogeneity of different climatic systems and the rapid rise of the Taihang Mountains shaped the contemporary geographical distribution pattern of the two Opisthopappus species.
Our results indicated that ecological factors play important role in shaping the physiological states, metabolic levels, and genetic alterations of species and populations and might drive the genetic makeup of populations and species. The results provide useful information for us to understand the ecology and evolution of organisms in the mountainous environment from population and species perspective.
In summary, when genetics, geographical conditions, climate variables, and evolutionary processes were all considered, O. taihangensis and O. longilobus were clearly distinct. At ~ 17.44 Ma during the early Miocene, the establishment of differing monsoon regimes due to the enhanced Asian monsoon from the QTP uplift triggered the derivation of O. taihangensis from O. longilobus. During the mid- late Miocene period, dramatic climatic shifts coupled with the progressive and heterogeneous uplift of the QTP initiated the intraspecific differentiation of these two species. Up until the Pleistocene, the rapid uplift of the Taihang Mountains coupled with violent climatic oscillations further promoted the diversity of the two species. With the formation of the Taihang Mountains, this complex topography led to localized environments and ecological heterogeneity, which established spatiotemporal isolation between populations. Under this scenario, O. taihangensis and O. longilobus underwent adaptive divergence, which gradually shaped current genetic structures and distribution patterns. The results of this study explored the differentiation mechanisms of these two species of the Opisthopappus genus, revealing the impacts of environmental events by taking small-scale spatial niches into consideration, while providing clues for the further investigation of other germplasm resources of the Taihang Mountains.
Our study was conducted in accordance with the laws of the People’s Republic of China, and field collection was approved by the Chinese Government. All researchers received permission letters from the College of Life Science, Shanxi Normal University, to collect the samples, which were taxonomically identified based on their phenotype by Junxia Su (Associate Professor of systematic botany) at Shanxi Normal University. The voucher specimens were deposited in the herbarium of College of Life Science, Shanxi Normal University (No:20170105030–20170105050).
Eleven populations of O. longilobus and thirteen populations of O. taihangensis were sampled, which covered the Opisthopappus distribution ranges (Table 1, Fig. 1). Individuals growing at a common site were regarded as a single "population". Fresh young leaves devoid of disease or insect pests were selected for each of the sample sites, where 10–15 individuals from each population were collected. These samples were placed into sealed bags filled with silica gel, dehydrated/quickly dried, and stored at 20 °C for later use. A global positioning system (GPS) was employed to demarcate each sample site and record the longitude, latitude, and elevation of each population (Table 1).
PCR amplification, sequencing, and genotyping
The total genomic DNA was extracted using the modified 2 × CTAB method . The quality of DNA was measured using an ultraviolet spectrophotometer and 0.8% agarose gel electrophoresis, and stored at − 20 °C for further use.
The SNP and InDel primers (Additional file 8: Table S3) of nuclear genes of Opisthopappus were obtained from a pervious study . For the SNP primers, the 20 µL PCR reaction contained 10 µL 2 × MasterMix, 2 µL template DNA (30 ng/µL), 1 µL primer S (10 µM), 1 µL primer A (10 µM), and 6 µL ddH2O. The PCR procedure proceeded as follows: pre-denaturation at 94 °C for 5 min., denaturation at 94 °C for 1 min, annealing temperature based on each primer setting for 1 min, elongation at 72 °C for 1.5 min., repeated for 35 cycles, last elongation at 72 °C for 10 min, and preservation at 4 °C. The PCR products detected using 2% agarose gel electrophoresis were confirmed via an automatic analysis electrophoresis gel imaging system, which were then sent to Sangon Biotech (Shanghai) for sequencing.
For the InDel primers, the PCR reaction was 20 µL, which contained 10 µL 2 × MasterMix, 3 µL template DNA (30 ng/µL), 1 µL primer S (10 µM), 1 µL primer A (10 µM), and 5 µL ddH2O. The PCR procedure was as follows: pre-denaturation at 94 °C for 1 min, denaturation at 94 °C for 1 min, annealing temperature based on each primer setting for 1 min, elongation at 72 °C for 1 min, repeated for 35 cycles, last elongation at 72 °C for 10 min, preservation at 4 °C. The PCR products were detected using 8% polyacrylamide gel electrophoresis. The presence or absence of each InDel fragment were coded as ‘1′and ‘0′ respectively. The details for the numbers of individuals for SNP sequencing and InDel genotyping are shown in Table 1.
Population genetic differentiation analyses
Prior to population genetic analysis, the partition homogeneity test (PHT) were initially conducted by PAUP  to identify whether the SNP sequences were suitable to be combined. The non-significant (P > 0.05) of the results revealed that the combined SNP sequences were suitable.
The haplotypes, haplotype frequencies, haplotype diversity (Hd), and nucleotide diversity (π) were calculated using DNASP 5.10 . The genetic GST and NST differentiation parameters were examined by PERMUT 2.0  based on the haplotype frequency.
For the InDel data, the genetic characteristics, Nei's gene diversity index (H), Shannon’s information index (I), and the percentage of polymorphic loci (PPL), were calculated by POPGENE 1.31 . An analysis of molecular variance (AMOVA) was implemented by ARLEQUIN 3.5  and GENALEX 6.5  to detect the distribution of genetic variations within and between populations or species. Subsequently, the FST, FCT, and FSC values  were calculated based on hierarchical AMOVA, and the permutation test was set to 1000.
Cluster analysis based on the maximum likelihood (ML) method and Nei’s genetic distance, respectively, was performed using MEGA 7.0 . Bayesian clustering analysis (BCA) was employed to examine the similarity and divergence of genetic components between populations and performed using STRUCTURE 2.2  for both the SNP sequencing and InDel data. The posterior probability of grouping number (K = 2–24) was estimated through 10 independent runs using 500,000 step Markov chain Monte Carlo (MCMC) replicates, following a 1,000,000-step burn-in for each run to evaluate consistency. The best grouping number was evaluated by ΔK  in STRUCTURE HARVESTER 0.6.94 . These 10 runs were aligned and summarized using CLUMPP 1.1.2  and the visualization of the results was plotted using DISTRUCT 1.1 .
To test the genetic differentiation between populations or species, a discriminant analysis of principal components (DAPC) was implemented by the function dapc in the R package ‘adegenet’ , which initially transformed the genetic data using principal component analysis (PCA) results, and subsequently performed discriminant analysis on the retained principal components. The properties of the “without a priori”, using partial synthetic variables to minimize variations within groups , might assist with objectively evaluating the artificial classification of O. taihangensis and O. longilobus. Kruskal–Wallis tests for the first two principal components (PCs), and the first two linear discriminants (LDs) of DAPC, were conducted to examine the genetic differentiation between the populations and species.
Inference of population demographic history
A network relationship was generated through the median-joining method in POPART 1.7 , to investigate the evolutionary relationships between the Opisthopappus haplotypes. BEAST 1.84  was employed to estimate the differentiation and diversification time between haplotypes. Chrysanthemum indicum, belonging to the same subtribe of Chrysantheminae with Opisthopappus (holding identified genomic information) was selected as the outgroup in BEAST analysis. The haplotype sequence of each primer was aligned to the NT (Nucleotide Sequence) database followed by manual splicing. Owing to the absence of the record of the Opisthopappus fossil data, the divergence time of Chrysanthemum and Opisthopappus (25.40 Ma) referred to the Time Tree website (http://www.timetree.org/) was adopted as a prerequisite for calibrating the age of most recent common ancestor (tMRCA).
The Akaike Information Criterion (AIC) with a “greedy” algorithm in PartitionFinder 2.1.1  was employed to select the best-fit partitioning schemes and evolutionary models. Based on the AIC results, the dataset was partitioned into four groups (group1: SNP2 + SNP29, group2: SNP4 + SNP26, group3: SNP13 + SNP32, and group4: SNP19 + SNP23), and the phylogenetic relationships were inferred based on four optimal evolutionary models, namely HKY + I + G + X, HKY + I + G, SYM + I + G and GTR + I + X, corresponding to group1 to group4, respectively. The generic average mutation rate of 6.1 × 10–9 (5.1 and 7.1 × 10–9) for the nuclear DNA of the Asteraceae species was employed according to the present study . Markov chain Monte Carlo (MCMC) was repeated 8 × 107 times by sampling every 80,000 generations. TRACER 1.5  was used to check the convergence of the framework, which ensured that every tested parameter was greater than 200.
To assess whether the species had experienced a significant expansion, we utilized ARLEQUIN 3.5  to calculate the Tajima’s D  and Fu’s FS  values. Moreover, the sum of square deviation (SSD) and raggedness index (Rag) in the mismatch distribution analysis (MDA) was performed in ARLEQUIN 3.5. The process employed a 1000 step permutation test.
Approximate Bayesian computation (ABC) analysis, provided by DIY-ABC 2.1.0 , enabled the estimation of complex evolutionary population histories. Based on the estimated genetic variations, genetic structures, and current geographic distributions, three evolutionary scenarios were proposed. Scenario 1: O. longilobus and O. taihangensis were differentiated from a common ancestral population during the same period. Scenario 2: O. taihangensis was an ancestral population, and O. longilobus was differentiated from O. taihangensis. Scenario 3: O. longilobus was the ancestral population, and O. taihangensis was differentiated from O. longilobus.
Each scenario was performed with 1,000,000 simulations and six summary statistics (number of haplotypes, number of segregating sites, mean of pairwise differences, Tajima’s D and private segregating sites) were selected. The substitution rates of nuclear genes were the same as those used in the BEAST analysis. To identify the best-supported scenario under direct and logistic approaches, we selected 1% of the simulated datasets closest to the observed data to evaluate model accuracy and estimate the relative posterior probability (PP) with 95% confidence intervals (95% CI) for each scenario. Further, the parameters including effective population size and divergence generation was estimated under the optimal scenario. The goodness of fit of the best supported scenario was evaluated by the option ‘model checking’ with principal component analysis (PCA). To estimate type I and II errors on the power of model selection, we assessed confidence in scenario choice with 500 simulated pseudo-observed data sets (PODs) for the seven plausible scenarios.
Additionally, the historical and contemporary gene flow were estimated within the two Opisthopappus species by MIGRATE-N 3.6  and BAYESASS 3.0 , respectively. In MIGRATE-N 3.6, maximum-likelihood analyses were performed using 10 short chains (104 trees) and three long chains (105 trees) with 104 trees discarded as an initial burn-in’ and astatic heating scheme at four temperatures (1, 1.5, 3, and 1000,000). To ensure the consistency of estimates, we repeated this procedure five times and reported average maximum-likelihood estimates with 95% confidence intervals. The parameters θ and M were estimated using a Bayesian method, which could be employed to estimate the number of migrants per generation (Nm) into each population using the Eq. 4Nm = θ*M.
When estimating the contemporary gene flow using BAYESASS 3.0, the parameters were examined including migration rates (m), allele frequencies (a) and inbreeding coefficients (f) to ensure that the optimal acceptance rates of the three parameters fell within the 20–60% range. Ten independent runs were executed to minimize the convergence problem. The result with the lowest deviance was adopted according to the method of Meirmans , where the 95% credible interval was estimated as m ± 1.96 × standard deviation (SD).
Environmental variables influence analyses
Nineteen bioclimatic variables (Bioclim) representing Grinnellian niches [110, 111], which are defined as the scenopoetic environmental variables of a species required to survive, were downloaded from the WorldClim database (http://www.worldclim.org/) with a resolution of 30 arc-sec (~ 1 × 1 km) and extracted using the R package ‘raster’ . Subsequently, the significance test of the distribution of climate factors along the two species was tested by one-way ANOVA. A principal component analysis (PCA) of independent climatic variables to reduce the dimensionality that defined the niche space, allowed for the comparison of the integrity of climate variables between O. longilobus and O. taihangensis, after which the PC1–PC3 were reserved for further analysis.
To test how the geographical and environmental differences impacted genetic differentiation, the Mantel test, partial Mantel test, and Barrier analysis were applied in this study. Further, a multiple matrix regression with randomization (MMRR) was performed to explore whether the genetic distance responded to variations in geographic and/or environmental distances.
Pairwise FST distance calculated in ARLEQUIN 3.5 was used as the genetic distance. The geographic distance was estimated using the GENALEX 6.5 according to three-dimensional factors (latitude, longitude, and elevation). The environmental distance was calculated using the Euclidean distance with PASSAGE 2.0 based on the first three PCs .
The Mantel test was performed in the R package ‘vegan’ , whereas the MMRR analysis was performed using the R package ‘PopGenReport’ [115, 116]. Logarithmic transformation of the distance matrices was conducted to ensure that they are in the same or similar order of magnitude. Regression coefficients of the Mantel test (r) and MMRR (r2) and their significance were determined based on 9,999 random permutations. Scatterplots to reveal the relationships between genetic, environmental, and geographic distances were conducted using GraphPad Prism 8 .
The biogeographic boundaries between population pairs were calculated by the Monmonier’s maximum-difference algorithm in BARRIER 2.2  based on the multiple distance matrix. Permutation and bootstrap tests were conducted with 1000 replicates for each case (Fig. 1).
In addition, distance based redundancy analyses (dbRDA) were performed to elucidate whether the climatic variables conditioned on the geographic distribution explained the genetic differentiation of the populations using the R package ‘vegan’. Firstly, a distance-based principal coordinate analysis (PCoA) of the genetic data at the species level was performed to generate several principal coordinates (PCs) using the R package ‘ape’ . Next, the PC1-3 of climatic variables were employed as explanatory variables conditioned on geographic factors, and significance tests were performed using the “anova. cca”  function in the R package ‘vegan’ with 999 permutations. The distribution pattern of the PC1-3 of climate variables along the ordination axes1-2 was further analyzed using a generalized linear model (GLM). Finally, the first two RDA axes and the explanatory variables were employed to construct the ordination and ordisurf plots of the dbRDA.
Availability of data and materials
The datasets generated and/or analyzed during the current study was available in the Dryad repository, https://doi.org/10.5061/dryad.p5hqbzkpd. The datasets used and analyzed during the current study was also available from the corresponding author on reasonable request.
Analysis of molecular variance
Bayesian clustering analysis
Principal components analysis
- K-W test:
Mismatch distribution analysis
Approximate Bayesian computation
Multiple matrix regression with randomization
Distance based redundancy analysis
Generalized linear model
Akaike information criterion
Mayr E: Animal species and evolution. 1963.
Knowles LL. Tests of Pleistocene speciation in Montane grasshoppers (genus Melanoplus) from the sky islands of western north Ametica. Evolution. 2000;54(4):1337–48.
Callahan CM, Rowe CA, Ryel RJ, Shaw JD, Madritch MD, Mock KE. Continental-scale assessment of genetic diversity and population structure in quaking aspen (Populus tremuloides). J Biogeogr. 2013;40(9):1780–91.
Carvalho CDS, Ballesteros-Mejia L, Ribeiro MC, Côrtes MC, Santos AS, Collevatti RG. Climatic stability and contemporary human impacts affect the genetic diversity and conservation status of a tropical palm in the Atlantic forest of Brazil. Conservation Genetics. 2017;18(2):467–78.
Ben-Menni Schuler S, López-Pujol J, Blanca G, Vilatersana R, Garcia-Jacas N, Suárez-Santiago VN. Influence of the quaternary glacial cycles and the mountains on the reticulations in the subsection Willkommia of the genus Centaurea. Front Plant Sci. 2019. https://doi.org/10.3389/fpls.2019.00303.
Lavergne S, Mouquet N, Thuiller W, Ronce O. Biodiversity and climate change: integrating evolutionary and ecological responses of species and communities. Annu Rev Ecol Evol Syst. 2010;41(1):321–50.
Mouquet N, Devictor V, Meynard CN, Munoz F, Bersier LF, Chave J, Couteron P, Dalecky A, Fontaine C, Gravel D, et al. Ecophylogenetics: advances and perspectives. Biol Rev. 2012;87(4):769–85.
Young A, Boyle T, Brown T. The population genetic consequences of habitat fragmentation for plants. Trends Ecol Evol. 1996;11(10):413–8.
Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405(6789):907–13.
Xie DF, Li MJ, Tan JB, Price M, Xiao QY, Zhou SD, Yu Y, He XJ. Phylogeography and genetic effects of habitat fragmentation on endemic Urophysa (Ranunculaceae) in Yungui Plateau and adjacent regions. PLoS ONE. 2017;12(10):e0186378.
Drake JM. Population effects of increased climate variation. Proc R Society B Biol Sci. 2005;272(1574):1823–7.
Stojanova B, Šurinová M, Klápště J, Koláříková V, Hadincová V, Münzbergová Z. Adaptive differentiation of Festuca rubra along a climate gradient revealed by molecular markers and quantitative traits. PLoS ONE. 2018;13(4):e0194670.
Sexton JP, McIntyre PJ, Angert AL, Rice KJ. Evolution and ecology of species range limits. Annu Rev Ecol Evol Syst. 2009;40(1):415–36.
Bernatchez L. On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes. J Fish Biol. 2016;89(6):2519–56.
Hewitt GM. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Lin Soc. 2008;58(3):247–76.
Deli T, Kiel C, Schubart CD. Phylogeographic and evolutionary history analyses of the warty crab Eriphia verrucosa (Decapoda, Brachyura, Eriphiidae) unveil genetic imprints of a late Pleistocene vicariant event across the Gibraltar Strait, erased by postglacial expansion and admixture among refugial lineages. BMC Evol Biol. 2019;19(1):105.
Contreras-Moreira B, Serrano-Notivoli R, Mohammed NE, Cantalapiedra CP, Beguería S, Casas AM, Igartua E. Genetic association with high-resolution climate data reveals selection footprints in the genomes of barley landraces across the Iberian Peninsula. Mol Ecol. 2019;28(8):1994–2012.
Sacks BN, Brown SK, Ernest HB. Population structure of California coyotes corresponds to habitat-specific breaks and illuminates species history. Mol Ecol. 2004;13(5):1265–75.
He SL, Wang YS, Li DZ, Yi TS. Environmental and historical determinants of patterns of genetic differentiation in wild soybean (Glycine soja Sieb. et Zucc). Sci Rep. 2016;6(1):22795.
Ye JW, Zhang ZK, Wang HF, Bao L, Ge JP. Phylogeography of Schisandra chinensis (Magnoliaceae) reveal multiple refugia with ample gene flow in Northeast China. Front Plant Sci. 2019;10:199.
Gong M. Uplifting process of southern Taihang Mountain in Cenozoic. Chinese Academy of Geological Science Thesis for Doctor Degree. 2010.
Zhang Y, Ma Y, Yang N, Shi W, Dong S. Cenozoic extensional stress evolution in North China. J Geodyn. 2003;36(5):591–613.
Zhang M, Li P. Discussion on the main uplift period of the Southern segment of Taihang Mountains. Territory Nat Res Study. 2014;4:20.
Zhu L. Spider community structure in fragmented habitats of Taihang Mountain area, China. Master of Dissertation. Hebei University; 2008.
Wu C, Zhang X, Ma Y. The Taihang and Yan mountains rose mainly in Quaternary. Norht China Earthquake Sciences. 1999;17(3):1–7.
Yan S: The investigation and collection of Pyrus betulaefolia in Taihang Mountains and evaluation of genetic diversity. Master of Dissertation. Agricultural University of Hebei; 2015.
Bai QQ, Pan Z, Ren GD. Phylogeographical analysis of Episyrphus balteatus (Diptera: Syrphidae) in Yanshan-Taihang Mountains Area. Chin J Ecol. 2018;37(1):157–63.
Zhao HB, Chen FD, Chen SM, Wu G-S, Guo WM. Molecular phylogeny of Chrysanthemum, Ajania and its allies (Anthemideae, Asteraceae) as inferred from nuclear ribosomal ITS and chloroplast trnL-F IGS sequences. Plant Syst Evol. 2010;284(3):153–69.
Yang D, Hu X, Liu Z, Zhao H. Intergeneric hybridizations between Opisthopappus taihangensis and Chrysanthemum lavandulifolium. Sci Hortic. 2010;125(4):718–23.
Tang F, Wang H, Chen S, Chen F, Teng N, Liu Z. First intergeneric hybrids within the tribe Anthemideae Cass. III. Chrysanthemum indicum L. Des Moul. × Opisthopappus taihangensis (Ling) Shih. Biochem Systematics Ecol. 2012;43:87–92.
Hu X: Preliminary studies on inter-generic hybridization within Chrysanthemum in broad sense (III). Thesis for Master’ Degree, Beijing Forestry University; 2008.
Ding BZ, Wang SY. Flora of Henan. Zhengzhou: Henan Science & Technology Press; 1998.
Wei DW, Xu MM, Sun WY, Jia CY, Zhang XW. Antioxidant activity of aqueous extracts from different organs of Opisthopappus Shih. J Chin Institute Food Sci Technol. 2015;15(2):56–63.
Wu ZY. Compositae. Flora of China. Beijing: Science Press; 1993.
Jia R, Wang Y. Leaves micromorphological characteristics of Opisthopappus taihangensis and Opisthopappus longilobus from Taihang Mountain, China. Vegetos. 2015;28(2):82–9.
Guo R, Zhou L, Zhao H, Chen F. High genetic diversity and insignificant interspecific differentiation in Opisthopappus Shih, an endangered cliff genus endemic to the Taihang Mountains of China. Sci World J. 2013;2013:275753.
Wang Y. Chloroplast microsatellite diversity of Opisthopappus Shih (Asteraceae) endemic to China. Plant Syst Evol. 2013;299(10):1849–58.
Wang Y, Yan G. Genetic diversity and population structure of Opisthopappus longilobus and Opisthopappus taihangensis (Asteraceae) in China determined using sequence related amplified polymorphism markers. Biochem Syst Ecol. 2013;49:115–24.
Wang Y, Yan G. Molecular phylogeography and population genetic structure of O. longilobus and O. taihangensis (Opisthopappus) on the Taihang Mountains. Plos One. 2014;9(8):e104773.
Wang Y, Zhang C, Lin L, Yuan L. ITS sequence analysis of Opisthopappus taihangensis and O. longilobus. Acta Horticulturae Sinica. 2015;42(1):86–94.
Chai M, Wang S, He J, Chen W, Fan Z, Li J, Wang Y. De novo assembly and transcriptome characterization of Opisthopappus (Asteraceae) for population differentiation and adaption. Front Genet. 2018;9:371.
Geng Q, Sun L, Zhang P, Wang Z, Qiu Y, Liu H, Lian C. Understanding population structure and historical demography of Litsea auriculata (Lauraceae), an endangered species in east China. Sci Rep. 2017;7(1):17343.
Chai M, Ye H, Wang Z, Zhou YC, Wu JH, Gao Y, Han W, Zang E, Zhang H, Ru WM, Sun GL, Wang YL. Genetic divergence and relationship among Opisthopappus species identified by development of EST-SSR markers. Front Genetics. 2020;11:177.
Lenormand T. Gene flow and the limits to natural selection. Trends Ecol Evol. 2002;17(4):183–9.
Shih KM, Chang CT, Chung JD, Chiang YC, Hwang SY. Adaptive genetic divergence despite significant isolation-by-distance in populations of Taiwan cowtail fir (Keteleeria davidiana Var. formosana). Front Plant Sci. 2018;9:92.
Endler JA. Gene flow and population differentiation. Science. 1973;179(4070):243–50.
Sexton JP, Hangartner SB, Hoffmann AA. Genetic isolation by environment of distance: which pattern of gene flow is most common? Evolution. 2014;68(1):1–15.
Liu W, Zhao Y, Qi D, You J, Zhou Y, Song Z. The Tanggula Mountains enhance population divergence in Carex moorcroftii: a dominant sedge on the Qinghai-Tibetan Plateau. Sci Rep. 2018;8(1):2741.
Star B, Spencer HG. Effects of genetic drift and gene flow on the selective maintenance of genetic variation. Genetics. 2013;194(1):235–44.
Huang BH, Huang CW, Huang CL, Liao PC. Continuation of the genetic divergence of ecological speciation by spatial environmental heterogeneity in island endemic plants. Sci Rep. 2017;7(1):5465.
Yang J, Vázquez L, Feng L, Liu Z, Zhao G. Climatic and soil factors shape the demographical history and genetic diversity of a deciduous oak (Quercus liaotungensis) in Northern China. Front Plant Sci. 2018. https://doi.org/10.3389/fpls.2018.01534.
Manel S, Poncet BN, Legendre P, Gugerli F, Holderegger R. Common factors drive adaptive genetic variation at different spatial scales in Arabis alpina. Mol Ecol. 2010;19(17):3824–35.
Wang IJ, Bradburd GS. Isolation by environment. Mol Ecol. 2014;23(23):5649–62.
Mosca E, González-Martínez SC, Neale DB. Environmental versus geographical determinants of genetic structure in two subalpine conifers. New Phytol. 2014;201(1):180–92.
Meng L, Chen G, Li Z, Yang Y, Wang Z, Wang L. Refugial isolation and range expansions drive the genetic structure of Oxyria sinensis (Polygonaceae) in the Himalaya-Hengduan Mountains. Sci Rep. 2015;5(1):10396.
Wu CI. The genic view of the process of speciation. J Evol Biol. 2001;14(6):851–65.
Barreda VD, Palazzesi L, Tellería MC, Olivero EB, Raine JI, Forest F. Early evolution of the angiosperm clade Asteraceae in the Cretaceous of Antarctica. Proc Natl Acad Sci. 2015;112(35):10989–94.
Graham A. A contribution to the geologic history of the Compositae. In: Compositae: systematics Proceedings of the international Compositae conference, Kew. Royal Botanic Gardens Kew, 1994. pp. 123–140.
Zhu Z, Wu L, Xi P, Song Z, Zhang Y. A research on Tertiary palynology from the Qaidam Basin, Qinghai Province. The Petroleum Industry Press; 1985. pp. 1–297.
Zhao HB. Phylogeny of tribe Anthemideae (Asteraceae) from east Asia and intergeneric cross between Dendranthema × Grandiflorum (Ramat.) Kitam. and Ajania pacifica (Nakai) K. Bremer & Humphries. Doctor dissertation. Nanjing Agricultural University; 2007.
Wang WM. On the origin and development of Artemisia (Asteraceae) in the geological past. Bot J Linn Soc. 2004;145(3):331–6.
Hobbs CR, Baldwin BG. Asian origin and upslope migration of Hawaiian Artemisia (Compositae–Anthemideae). J Biogeogr. 2013;40(3):442–54.
Li J, Fang X. Uplift of the Tibetan Plateau and environmental changes. Chin Sci Bull. 1999;44(23):2117–24.
An ZS, Zhang PZ, Wang EQ, Wang SM, Qaing XK, Li L, Song YG, Chang H, Liu XD, Zhou WJ, et al. Changes of the monsoon-arid environment in China and growth of the Tibetan Plateau since the Miocene. Quaternary Sci. 2006;26(5):678–93.
Harrison TM, Copeland P, Kidd WSF, Lovera OM. Activation of the nyainqentanghla shear zone: implications for uplift of the southern Tibetan Plateau. Tectonics. 1995;14(3):658–76.
Shi Y, Li J, Li B. Uplift and environmental changes of Qinghai-Tibetan Plateau in the Late Cenozoic. Guangzhou: Guangdong Science and Technology Press; 1998.
Spicer RA, Harris NBW, Widdowson M, Herman AB, Guo S, Valdes PJ, Wolfe JA, Kelley SP. Constant elevation of southern Tibet over the past 15 million years. Nature. 2003;421(6923):622–4.
Molnar P, England P, Martinod J. Mantle dynamics, uplift of the Tibetan Plateau, and the Indian Monsoon. Rev Geophys. 1993;31(4):357–96.
Wan S, Li A, Clift PD, Stuut JBW. Development of the East Asian monsoon: mineralogical and sedimentologic records in the northern South China Sea since 20 Ma. Palaeogeogr Palaeoclimatol Palaeoecol. 2007;254(3):561–82.
Zhisheng A, Kutzbach JE, Prell WL, Porter SC. Evolution of Asian monsoons and phased uplift of the Himalaya-Tibetan plateau since Late Miocene times. Nature. 2001;411(6833):62–6.
Bloemendal J, Demenocal P. Evidence for a change in the periodicity of tropical climate cycles at 2.4 Myr from whole-core magnetic susceptibility measurements. Nature. 1989;342(6252):897–900.
Favre A, Päckert M, Pauls SU, Jähnig SC, Uhl D, Michalak I, Muellner-Riehl AN. The role of the uplift of the Qinghai-Tibetan Plateau for the evolution of Tibetan biotas. Biol Rev. 2015;90(1):236–53.
Mulch A, Chamberlain CP. The rise and growth of Tibet. Nature. 2006;439(7077):670–1.
Wang C, Zhao X, Liu Z, Lippert PC, Graham SA, Coe RS, Yi H, Zhu L, Liu S, Li Y. Constraints on the early uplift history of the Tibetan Plateau. Proc Natl Acad Sci. 2008;105(13):4987–92.
Zhao YJ, Gong X. Genetic divergence and phylogeographic history of two closely related species (Leucomeris decora and Nouelia insignis) across the “Tanaka Line” in Southwest China. BMC Evol Biol. 2015;15(1):134.
Zachos J, Pagani M, Sloan L, Thomas E, Billups K. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science. 2001;292(5517):686–93.
Li G, Pettke T, Chen J. Increasing Nd isotopic ratio of Asian dust indicates progressive uplift of the north Tibetan Plateau since the middle Miocene. Geology. 2011;39(3):199–202.
Zhao Y, Yin G, Pan Y, Gong X. Ecological and genetic divergences with gene flow of two sister species (Leucomeris decora and Nouelia insignis) driving by climatic transition in Southwest China. Front Plant Sci. 2018. https://doi.org/10.3389/fpls.2018.00031.
Ge J, Guo Z, Zhan T, Yao Z, Deng C, Oldfield F. Magnetostratigraphy of the Xihe loess-soil sequence and implication for late Neogene deformation of the West Qinling Mountains. Geophys J Int. 2012;189(3):1399–408.
Rost KT. Paleoclimatic field studies in and along the Qinling Shan (Central China). Geo J. 1994;34(1):107–20.
Lei Q, Liu DK, Li S, Ji W. Geomorphological characteristics and cause of Taihang Mountains. Technol Innovation Appl. 2019;09:78–9.
Zhang Z, Zhang JL. Discussion on the uplift of the south section of Taihang Mountain in Quaternary period. J Arid Land Resources Environ. 2020;34(10):87–92.
Yao YF, Bruch AA, Mosbrugger V, Li CS. Quantitative reconstruction of Miocene climate patterns and evolution in Southern China based on plant fossils. Palaeogeogr Palaeoclimatol Palaeoecol. 2011;304(3):291–307.
Su T, Jacques FMB, Spicer RA, Liu YS, Huang YJ, Xing YW, Zhou ZK. Post-Pliocene establishment of the present monsoonal climate in SW China: evidence from the late Pliocene Longmen megaflora. Clim Past. 2013;9(4):1911–20.
Yang XL, Xu QH, Zhao HP, Liang WD, Sun LM. Vegetation changes of the Taihang Mountains since the last glacial. Chin Geogra Sci. 2000;10(3):261–9.
Swofford DL. Paup*: phylogenetic analysis using parsimony (and other methods) 4.0. B5. 2001.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.
Pons O, Petit RJ. Measwring and testing genetic differentiation with ordered versus unordered alleles. Genetics. 1996;144(3):1237–45.
Yeh FC, Y RC, Boyle T. POPGENE Version 1.31. Microsoft windows-based freeware for population genetic analysis. University of Alberta and Centre for International Forestry Research 1998:11–23.
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 Resources. 2010;10(3):564–7.
Peakall R, Smouse PE. Genalex 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Hubisz MJ, Falush D, Stephens M, Pritchard JK. Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009;9(5):1322–32.
Evanno G, Reganut 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, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.
Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23(14):1801–6.
Rosenberg NA. Distruct: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4(1):137–8.
Jombart T. Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24(11):1403–5.
Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11(1):94.
Leigh JW, Bryant D, Nakagawa S. Popart: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6(9):1110–6.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.
Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. Partition finder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2016;34(3):772–3.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123(3):585–95.
Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147(2):915–25.
Cornuet JM, Pudlo P, Veyssier J, DehneGarcia A, Gautier M, Leblois R, Marin JM, Estoup A. DIYABC v20: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism DNA sequence and microsatellite data. Bioinformatics. 2014;30(8):1187–9.
Beerli P. Comparison of Bayesian and maximum-likelihood inference of population genetic parameters. Bioinformatics. 2005;22(3):341–5.
Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003;163(3):1177–91.
Meirmans PG. Nonconvergence in Bayesian estimation of migration rates. Mol Ecol Resour. 2014;14(4):726–33.
Grinnell J. The niche-relationships of the California thrasher. Auk. 1917;34(4):427–33.
Stephen EF, Robert JH. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int J Climatol. 2017. https://doi.org/10.1002/joc.5086.
Hijmans RJ: raster: Geographic data analysis and modeling. R package version 3.0–2. https://CRAN.R-project.org/package=raster. 2019.
Rosenberg MS, Anderson CD. PASSaGE: pattern analysis, spatial statistics and geographic exegesis. version 2. Methods Ecol Evol. 2011;2(3):229–32.
Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14(6):927–30.
Wang IJ. Examining the full effects of landscape heterogeneity on spatial genetic variation: a multiple matrix regression approach for quantifying grographic and ecological isolation. Evolution. 2013;67(12):3403–11.
Adamack AT, Gruber B. PopGenReport: simplifying basic population genetic analyses in R. Methods Ecol Evol. 2014;5(4):384–7.
Swift ML. GraphPad prism, data analysis, and scientific graphing. J Chem Inf Comput Sci. 1997;37(2):411–2.
Manni F, Gu XC, Rard E, Heyer E. Geographic patterns of (genetic, morphologic, linguistic) variation: how barriers can be detected by using Monmonier’s Algorithm. Hum Biol. 2004;76(2):173–90.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20(2):289–90.
Legendre P, Oksanen J, ter Braak CJF. Testing the significance of canonical axes in redundancy analysis. Methods Ecol Evol. 2011;2(3):269–77.
The authors thank Junxia Su for her assistance with the sample identification.
This study was funded by the National Natural Science Foundation of China (31970358 to Y.W.) and Research Project Supported by Shanxi Scholarship Council of China (2020-090 to Y.W.). The funding agency did not play a role in the experimental design, results analysis, or writing of the manuscript, but did provide financial support for the manuscript.
Ethics approval and consent to participate
Our study was conducted in accordance with the laws of People’s Republic of China, and the field collection was approved by Chinese Government. All researchers received permission letters from College of Life Science, Shanxi Normal University to collect the samples, which were taxonomically identified based on their morphological characteristics by Junxia Su (Associate professor of the systematic botany) at Shanxi Normal University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Phylogenetic relationships between O. longilobus and O. taihangensis. (A): Individual ML clustering of Opisthopappus. Blue branches presented individuals of O. longilobus and red branches presented individuals of O. taihangensis. (B): UPGMA clustering for 24 populations of Opisthopappus based on Nei’s genetic distance. Blue for populations of O. longilobus and red for populations of O. taihangensis.
Structure analysis from K = 3 to K = 6 for SNP and InDel, respectively.
The Kruskal–Wallis test of the first two principal components and the first two linear discriminants of the genetic variation revealed significant genetic divergence between species but no or little population differentiation within species. (A–D): Comparisons between species. (E–H): Comparisons among populations.
Haplotypes network of Opisthopappus. 47 haplotypes (H1-H47) were detected in O. longilobus and 28 haplotypes (H48-H75) in O. taihangensis. No shared haplotypes were detected between O. longilobus and O. taihangensis. The color of each haplotype corresponded to Fig. 1 The size of the circles corresponds to the frequency of each haplotype and each solid line represents one mutational step.
The results of neutrality tests (Tajima’s D and Fu’s FS tests) and mismatch distribution analyses.
ANOVA analysis for the nineteen bioclimatic variables grouped by two different species.
The average temperature of every month (A) and the average precipitation of every month (B) of the studied populations of the distribution of Opisthopappus.
Information of primer pairs.
About this article
Cite this article
Ye, H., Wang, Z., Hou, H. et al. Localized environmental heterogeneity drives the population differentiation of two endangered and endemic Opisthopappus Shih species. BMC Ecol Evo 21, 56 (2021). https://doi.org/10.1186/s12862-021-01790-0
- Genetic differentiation
- Environmental heterogeneity
- Taihang Mountains