Retreating or standing: responses of forest species and steppe species to climate change in arid Eastern Central Asia

Comparative evolutionary history of two closely related desert plant, Convolvulus tragacanthoide and Convolvulus gortschakovii (Convolvulaceae) from northwest China

Desert ecosystems are one of the most fragile ecosystems on Earth. The study of the effects of paleoclimatic and geological changes on genetic diversity, genetic structure, and species differentiation of desert plants is not only helpful in understanding the strategies of adaptation of plants to arid habitats, but can also provide reference for the protection and restoration of vegetation in desert ecosystem. Northwest China is an important part of arid regions in the northern hemisphere. Convolvulus tragacanthoides and Convolvulus gortschakovii are closely related and have similar morphology. Through our field investigation, we found that the annual precipitation of the two species distribution areas is significantly different. Thus, C. tragacanthoides and C. gortschakovii provide an ideal comparative template to investigate the evolutionary processes of closely related species, which have adapted to different niches in response to changes in paleogeography and paleoclimate in northwest China. In this study, we employed phylogeographical approaches (two cpDNA spacers: rpl14-rpl36 and trnT-trnY) and species distribution models to trace the demographic history of C. tragacanthoides and C. gortschakovii, two common subshrubs and small shrubs in northwest China. The results showed the following: (1) Populations of C. tragacanthoides in northwest China were divided into three groups: Tianshan Mountains-Ili Valley, west Yin Mountains-Helan Mountains-Qinglian Mountains, and Qinling Mountains-east Yin Mountains. There was a strong correlation between the distribution of haplotypes and the floristic subkingdom. The three groups corresponded to the Eurasian forest subkingdom, Asian desert flora subkingdom, and Sino-Japanese floristic regions, respectively. Thus, environmental differences among different flora may lead to the genetic differentiation of C. tragacanthoides in China. (2) The west Yin Mountains-Helan Mountains-Qinglian Mountains, and Qinling Mountains-east Yin Mountains were thought to form the ancestral distribution range of C. tragacanthoides. (3) C. tragacanthoides and C. gortschakovii adopted different strategies to cope with the Pleistocene glacial cycle. Convolvulus tragacanthoides contracted to the south during the glacial period and expanded to the north during the interglacial period; and there was no obvious north-south expansion or contraction of C. gortschakovii during the glacial cycle. (4) The interspecific variation of C. tragacanthoides and C. gortschakovii was related to the orogeny in northwest China caused by the uplift of the Tibetan Plateau during Miocene. (5) The 200 mm precipitation line formed the dividing line between the niches occupied by C. tragacanthoides and C. gortschakovii, respectively. In this study, from the perspective of precipitation, the impact of the formation of the summer monsoon limit line on species divergence and speciation is reported, which provides a new perspective for studying the response mechanism of species to the formation of the summer monsoon line, and also provides a clue for predicting how desert plants respond to future environmental changes.

Phylogeography of Amygdalus mongolica in relation to Quaternary climatic aridification and oscillations in northwestern China

Quaternary period geological events and climatic oscillations significantly affect the geographic structure and genetic diversity of species distribution in arid northwestern China. Amygdalus mongolica is a relict and endangered shrub that occurs primarily in arid areas of northwestern China. Based on variation patterns present at three cpDNA regions (psbK-psbI, trnL-trnF and trnV) and in one nDNA sequence (ITS1-ITS4) in 174 individuals representing 15 populations, the spatial genetic structure and demographic history of A. mongolica was examined across its entire geographic range. The 17 different haplotypes and 10 ribotypes showed two lineages, distributed across the Western (Mazong Mountains, Hexi Corridor, and Alxa Left Banner) and Eastern regions (Urad Houqi, Yinshan Mountains, Urad Zhongqi, and Daqing Mountains) according to the median-joining network and the BI (Bayesian inference) and ML (Maximum likelihood) trees. AMOVA analysis demonstrated that over 65% of the observed genetic variation was related to this lineage split. The expansions of the Ulanbuhe and Tengger deserts and the eastward extension of the Yinshan Mountains since the Quaternary period likely interrupted gene flow and triggered the observed divergence in the two allopatric regions; arid landscape fragmentation accompanied by local environmental heterogeneity further increased local adaptive differentiation between the Western and Eastern groups. Based on the evidence from phylogeographical patterns and the distribution of genetic variation, A. mongolica distributed in the eastern and western regions are speculated to have experienced eastward migration along the southern slopes of the Lang Mountains and westward migration along the margins of the Ulanbuhe and Tengger deserts to the Hexi Corridor, respectively. For setting a conservation management plan, it is recommended that the south slopes of the Lang Mountains and northern Helan Mountains be identified as the two primary conservation areas, as they have high genetic variation and habitats that are more suitable.

Rear-edge, low-diversity, and haplotypic uniformity in cold-adapted Bupleurum euphorbioides interglacial refugia populations

The high genetic diversity of rear-edge refugia populations is predicted to have resulted from species repeatedly migrating to low latitudes during glacial periods over the course of Quaternary climate change. However, several recent empirical studies of cold-tolerant plants revealed the opposite pattern. We investigated whether current habitats of the cold-adapted and range-restricted Bupleurum euphorbioides in the Baekdudaegan, South Korea, and North Korea could be interglacial refugia, and documented how their rear-edge populations differ genetically from those of typical temperate species. Phylogeographic analysis and ecological niche modeling (ENM) were used. The genetic structure was analyzed using microsatellite markers and chloroplast DNA sequences. The congener B. longiradiatum was included as a typical temperate plant species. Despite having almost identical life history traits, these congeneric species exhibited contrasting patterns of genetic diversity. ENM revealed an apparent maximum range contraction during the last interglacial. In contrast, its range expanded northward to the Russian Far East (Primorsky) during the Last Glacial Maximum. Thus, we hypothesize that B. euphorbioides retreated to its current refugia during interglacial periods. Unlike populations in the central region, the rear-edge populations were genetically impoverished and uniform, both within populations and in pooled regional populations. The rear-edge B. euphorbioides survived at least one past interglacial, contributing to the species' genetic diversity. We believe that such genetic variation in the cold-adapted B. euphorbioides gives the species the necessary adaptations to survive an upcoming favorable environment (the next glacial), unless there is artificial environmental change.

Pleistocene climate change and phylogeographic structure of the Gymnocarpos przewalskii (Caryophyllaceae) in the northwest China: Evidence from plastid DNA, ITS sequences, and Microsatellite

Northwestern China has a wealth of endemic species, which has been hypothesized to be affected by the complex paleoclimatic and paleogeographic history during Quaternary. In this paper, we used Gymnocarpos przewalskii as a model to address the evolutionary history and current population genetic structure of species in northwestern China. We employed two chloroplast DNA fragments (rps16 and psbB-psbI), one nuclear DNA fragment (ITS), and simple sequence repeat (SSRs) to investigate the spatial genetic pattern of G. przewalskii. High genetic diversity (cpDNA: h S = 0.330, h T = 0.866; ITS: h S = 0.458, h T = 0.872) was identified in almost all populations, and most of the population have private haplotypes. Moreover, multimodal mismatch distributions were observed and estimates of Tajima's D and Fu's FS tests did not identify significantly departures from neutrality, indicating that recent expansion of G. przewalskii was rejected. Thus, we inferred that G. przewalskii survived generally in northwestern China during the Pleistocene. All data together support the genotypes of G. przewalskii into three groups, consistent with their respective geographical distributions in the western regions-Tarim Basin, the central regions-Hami Basin and Hexi Corridor, and the eastern regions-Alxa Desert and Wulate Prairie. Divergence among most lineages of G. przewalskii occurred in the Pleistocene, and the range of potential distributions is associated with glacial cycles. We concluded that climate oscillation during Pleistocene significantly affected the distribution of the species.

Contrasting Changes in Vegetation Growth due to Different Climate Forcings over the Last Three Decades in the Selenga-Baikal Basin

The Selenga-Baikal Basin, a transboundary river basin between Mongolia and Russia, warmed at nearly twice the global rate and experienced enhanced human activities in recent decades. To understand the vegetation response to climate change, the dynamic spatial-temporal characteristics of the vegetation and the relationships between the vegetation dynamics and climate variability in the Selenga-Baikal Basin were investigated using the Normalized Difference Vegetation Index (NDVI) and gridded temperature and precipitation data for the period of 1982 to 2015. Our results indicated that precipitation played a key role in vegetation growth across regions that presented multiyear mean annual precipitation lower than 350 mm, although its importance became less apparent over regions with precipitation exceeding 350 mm. Because of the overall temperature-limited conditions, temperature had a more substantial impact on vegetation growth than precipitation. Generally, an increasing trend was observed in the growth of forest vegetation, which is heavily dependent on temperature, whereas a decreasing trend was detected for grassland, for which the predominant growth-limiting factor is precipitation. Additionally, human activities, such as urbanization, mining, increased wildfires, illegal logging, and livestock overgrazing are important factors driving vegetation change.

Phylogeography and conservation genetics of the endangered Tugarinovia mongolica (Asteraceae) from Inner Mongolia, Northwest China

Tugarinovia (Family Asteraceae) is a monotypic genus. It’s sole species, Tugarinovia mongolica Iljin, is distributed in the northern part of Inner Mongolia, with one additional variety, Tugarinovia mongolica var ovatifolia, which is distributed in the southern part of Inner Mongolia. The species has a limited geographical range and declining populations. To understand the phylogeographic structure of T. mongolica, we sequenced two chloroplast DNA regions (psbA-trnH and psbK-psbI) from 219 individuals of 16 populations, and investigated the genetic variation and phylogeographic patterns of T. mongolica. The results identified a total of 17 (H1-H17) chloroplast haplotypes. There were no haplotypes shared between the northern (T. mongolica) and southern groups (T. mongolica var. ovatifolia), and they formed two distinct lineages. The regional split was also supported by AMOVA and BEAST analyses. AMOVA showed the main variation that occurred between the two geographic groups. The time of divergence of the two groups can be dated to the early Pleistocene epoch, when climate fluctuations most likely resulted in the allopatric divergence of T. mongolica. The formation of the desert blocked genetic flow and enhanced the divergence of the northern and southern groups. Our results indicate that the genetic differences between T. mongolica and T. mongolica var. ovatifolia are consistent with previously proposed morphological differences. We speculate that the dry, cold climate and the expansion of the desert during the Quaternary resulted in the currently observed distribution of extant populations of T. mongolica. In the northern group, the populations Chuanjinsumu, Wuliji and Yingen displayed the highest genetic diversity and should be given priority protection. The southern group showed a higher genetic drift (F_ST = 1, G_ST = 1), and the inbreeding load (H_S = 0) required protection for each population. Our results propose that the protection of T. mongolica should be implemented through in situ and ex situ conservation practices to increase the effective population size and genetic diversity.

Decreasing brown bear (Ursus arctos) habitat due to climate change in Central Asia and the Asian Highlands

Around the world, climate change has impacted many species. In this study, we used bioclimatic variables and biophysical layers of Central Asia and the Asian Highlands combined with presence data of brown bear (Ursus arctos) to understand their current distribution and predict their future distribution under the current rate of climate change. Our bioclimatic model showed that the current suitable habitat of brown bear encompasses 3,430,493 km2 in the study area, the majority of which (>65%) located in China. Our analyses demonstrated that suitable habitat will be reduced by 11% (378,861.30 km2) across Central Asia and the Asian Highlands by 2,050 due to climate change, predominantly (>90%) due to the changes in temperature and precipitation. The spatially averaged mean annual temperature of brown bear habitat is currently -1.2°C and predicted to increase to 1.6°C by 2,050. Mean annual precipitation in brown bear habitats is predicted to increase by 13% (from 406 to 459 mm) by 2,050. Such changes in two critical climatic variables may significantly affect the brown bear distribution, ethological repertoires, and physiological processes, which may increase their risk of extirpation in some areas. Approximately 32% (1,124,330 km2) of the total suitable habitat falls within protected areas, which was predicted to reduce to 1,103,912 km2 (1.8% loss) by 2,050. Future loss of suitable habitats inside the protected areas may force brown bears to move outside the protected areas thereby increasing their risk of mortality. Therefore, more protected areas should be established in the suitable brown bear habitats in future to sustain populations in this region. Furthermore, development of corridors is needed to connect habitats between protected areas of different countries in Central Asia. Such practices will facilitate climate migration and connectivity among populations and movement between and within countries.

Phylogeography and genetic effects of habitat fragmentation on endemic Urophysa (Ranunculaceae) in Yungui Plateau and adjacent regions

Urophysa is a Chinese endemic genus with only two species (U. rockii and U. henryi) distributed in Yungui Plateau (Guizhou Province) and adjacent regions (i.e., Provinces of Hunan, Hubei, Chongqing and Sichuan). The aim of this study was to determine the genetic diversity and population differentiation within Urophysa and investigate the effect of the Yungui Plateau uplift and climate oscillations on evolution of Urophysa. In this study, micro-morphological characteristics, nine microsatellite loci (SSR), two nuclear loci (ITS and ETS) and two chloroplast fragments (psbA-trnH and trnL-trnF) were used to analyze the phylogenetic relationships and assess genetic and phylogeographical structure of Urophysa. Isolation by distance (IBD) was performed to research the effects of geographical isolation. We detected high genetic diversity at the species level but low genetic diversity within populations. Striking genetic differentiation (AMOVA) among populations and a significant phylogeographical structure (NST > GST, p < 0.01) were detected among U. henryi populations, along with significant effects of isolation by distance (IBD). Molecular clock estimation using calibration strategy and cpDNA substitution rate indicated that the divergence of U. henryi occurred during late Miocene to early Quaternary, when the orogeny of Yungui Plateau was violent. U. rockii originated at the early Quaternary and further differentiated at early Pleistocene. Our results suggested that habitat fragmentation played an important role in the genetic diversity and population differentiation of U. rockii and U. henryi. Heterogenous geomorphological configuration and complicated environment resulted from rapid uplift of the Yungui Plateau were inferred as important incentives for the modern phylogeograhpical pattern and species divergence of Urophysa. The geographical isolation, limited gene flow, specialized morphologies and the Pleistocene climatic oscillation greatly contributed to the allopatric divergence of U. rockii. Significant genetic drift and inbreeding were detected in these two species, in situ measures should be implemented to protect them.

Spatial Patterns of Species Diversity and Phylogenetic Structure of Plant Communities in the Tianshan Mountains, Arid Central Asia

The Tianshan Mountains, located in arid Central Asia, have a humid climate and are biodiversity hotspots. Here, we aimed to clarify whether the pattern of species diversity and the phylogenetic structure of plant communities is affected by environmental variables and glacial refugia. In this study, plant community assemblies of 17 research sites with a total of 35 sample plots were investigated at the grassland/woodland boundaries on the Tianshan Mountains. Community phylogeny of these plant communities was constructed based on two plant DNA barcode regions. The indices of phylogenetic diversity and phylogenetic community structure were calculated for these sample plots. We first estimated the correlation coefficients between species richness (SR) and environmental variables as well as the presence of glacial refugia. We then mapped the significant values of indices of community phylogeny (PD, RPD, NRI, and NTI) to investigate the correlation between community phylogeny and environmental structure or macrozones in the study area. The results showed that a significantly higher value of SR was obtained for the refugial groups than for the colonizing groups (P < 0.05); presence of refugia and environmental variables were highly correlated to the pattern of variation in SR. Indices of community phylogeny were not significantly different between refugial and colonizing regions. Comparison with the humid western part showed that plant communities in the arid eastern part of the Tianshan Mountains tended to display more significant phylogenetic overdispersion. The variation tendency of the PhyloSor index showed that the increase in macro-geographical and environmental distance did not influence obvious phylogenetic dissimilarities between different sample plots. In conclusion, glacial refugia and environmental factors profoundly influenced the pattern of SR, but community phylogenetic structure was not affected by glacial refugia among different plant communities on the Tianshan Mountains. This pattern of community phylogenetic structure could have resulted from shared ancestry and species pool among these sample plots.

Distribution pattern of poisonous plant species in arid grasslands: a case from Xinjiang, Northwestern China

Poisonous plants threaten the ecosystem health of grasslands and the sustainability of animal husbandry. In arid lands, grassland ecosystems tend to be vulnerable and have been degraded due to the influence of human activities. The total area of the natural grasslands in Xinjiang, a large region in arid north-western China, ranks third in terms of area in China. In the process of grassland degradation, poisonous plants have spread widely and quickly in this region. During recent years, increasing economic losses have been caused by poisonous plants in Xinjiang. Although poisonous plants have been reported at some specific locations, their spatial patterns have rarely been investigated at a large regional scale. To understand the current status of hazards and assess the invasion risks of poisonous plants, we sampled ~150 poisonous plant species from Xinjiang and modelled the present and the future (the 2050s and the 2070s) distribution of 90 species using species distribution modelling. Based on the distribution maps of these poisonous plants, four diversity hotspots of poisonous plants were identified in Xinjiang. The results showed that northern Xinjiang had higher levels of poisonous plant diversity compared with the other part of Xinjiang. The precipitation factors had the most influence on prediction of the poisonous plants distributions in the species distribution modelling. Under the scenarios of future climate change, the results of modelling showed that regions close to the four hotspots of poisonous plants in Xinjiang displayed higher risks of invasion by poisonous plants in the future. In addition, these areas with a high risk of plant invasion will become increasingly large. We propose that policy makers consider implementing monitoring and prevention measures in areas identified as having a high risk of future invasion by poisonous plants.

Long-Distance Dispersal after the Last Glacial Maximum (LGM) Led to the Disjunctive Distribution of Pedicularis kansuensis (Orobanchaceae) between the Qinghai-Tibetan Plateau and Tianshan Region

Quaternary climate fluctuations have profoundly affected the current distribution patterns and genetic structures of many plant and animal species in the Qinghai-Tibetan Plateau (QTP) and adjacent mountain ranges, e.g. Tianshan (TSR), Altay, etc. In this greater area disjunct distributions are prominent but have nevertheless received little attention with respect to the historical processes involved. Here, we focus on Pedicularis kansuensis to test whether the current QTP and TSR disjunction is the result of a recent Holocene range expansion involving dispersal across arid land bridge(s) or a Pleistocene range fragmentation involving persistence in refugia. Two chloroplast DNA spacers were sequenced for 319 individuals from 34 populations covering the entire distribution range of this species in China. We found a total of 17 haplotypes of which all occurred in the QTP, and only five in the TSR. Overall genetic diversity was high (H_T = 0.882, H_S = 0.559) and higher in the QTP than in the TSR. Genetic differentiation among regions and populations was relatively low (G_ST = 0.366) and little evidence for a phylogeographic pattern emerged. The divergence times for the four main lineages could be dated to the early Pleistocene. Surprisingly, the two ubiquitous haplotypes diverged just before or around the Last Glacial Maximum (LGM) and were found in different phylogenetic lineages. The Species Distribution Model suggested a disappearance of P. kansuensis from the TSR during the LGM in contrast to a relatively constant potential distribution in the QTP. We conclude that P. kansuensis colonized the TSR after the LGM. The improbable long-distance dispersal by wind or water across arid land seed flow may well have had birds or men as vector.

Phylogeography of Libanotis buchtormensis (Umbelliferae) in Disjunct Populations along the Deserts in Northwest China

In Northwest China, aridification and desert expansion play significant roles in promoting desert plant diversification and speciation. However, to date, little is known about the effects of the desert barrier on the population structure of montane, non-desert species in the area. In this study, we sequenced chloroplast DNA regions (trnL–trnF and trnS–trnG) and a nuclear gene (rpb2) to investigate the population differentiation and phylogeographical history of Libanotis buchtormensis, a perennial montane species possessing a disjunct distribution at the periphery of the central desert. In total, 23 chloroplast haplotypes and 24 nuclear haplotypes were recovered from the 21 natural populations and six hebarium specimens. Phylogenetic analysis based on the combined plastid and nuclear dataset revealed two distinct lineages of L. buchtormensis, which inhabit the disjunct areas on both sides of the desert zone. The molecular dating analysis indicated that the divergence between the southeastern and the northwestern populations occurred in the middle Pleistocene, concomitantly with the desert expansion. The geographical vicariance likely contributed to the present disjunct distribution of L. buchtormensis across the deserts in Northwest China. Populations in the southeastern region may have migrated from the northwestern region, and seem to be a peripheral distribution of L. buchtormensis.

Phylogeography of the arid shrub Atraphaxis frutescens (Polygonaceae) in northwestern China: evidence from cpDNA sequences

Climatic fluctuations during the Pleistocene are usually considered as a significant factor in shaping intraspecific genetic variation and influencing demographic histories. To well-understand these processes in desert northwest China, we selected arid adapted Atraphaxis frutescens as the study species. Two cpDNA regions (psbK-psbI, psbB-psbH) were sequenced in 272 individuals from 33 natural populations across the range of this shrub, and 10 haplotypes were identified. It was found to contain high levels of total gene diversity (H T = 0.858), and low levels of within-population diversity (H S = 0.092). Analysis of molecular variance (AMOVA) indicates that genetic differentiation primarily occurs among groups of populations. Based on BEAST (Bayesian Evolutionary Analysis Sampling Trees) analysis, we suggest that intraspecific differentiation of the species, resulting from isolated populations, accompanied enhanced desertification during the middle and late Pleistocene. The expansion of the Gurbantunggut and Kumtag deserts in this area appears to have triggered divergence among populations of the western, central, and eastern portions of the region and shaped genetic differentiation among them. Two possible independent glacial refugia were predicted, the Ili Valley and the northern Junggar Basin. Extensive development of arid habitats (desert margin and arid piedmont grassland) coupled with a more equable climate because the early Holocene are factors likely to have generated recent expansion of A. frutescens.

Tertiary montane origin of the Central Asian flora, evidence inferred from cpDNA sequences of Atraphaxis (Polygonaceae)

Atraphaxis has approximately 25 species and a distribution center in Central Asia. It has been previously used to hypothesize an origin from montane forest. We sampled 18 species covering three sections within the genus and sequenced five cpDNA spacers, atpB-rbcL, psbK-psbI, psbA-trnH, rbcL, and trnL-trnF. BEAST was used to reconstruct phylogenetic relationship and time divergences, and S-DIVA and Lagrange were used, based on distribution area and ecotype data, for reconstruction of ancestral areas and events. Our results appear compatible with designation of three taxonomic sections within the genus. The generic stem and crown ages were Eocene, approximately 47 Ma, and Oligocene 27 Ma, respectively. The origin of Atraphaxis is confirmed as montane, with an ancestral area consisting of the Junggar Basin and uplands of the Pamir-Tianshan-Alatau-Altai mountain chains, and ancestral ecotype of montane forest. Two remarkable paleogeographic events, shrinkage of the inland Paratethys Sea at the boundary of the late Oligocene and early Miocene, and the time intervals of cooling and drying of global climate from 24 (22) Ma onward likely facilitated early diversification of Atraphaxis, while rapid uplift of the Tianshan Mountains during the late Miocene may have promoted later diversification.