Climate change profoundly influences species distributions. These effects are evident in poleward latitudinal range shifts for many taxa, and upward altitudinal range shifts for alpine species, that resulted from increased annual global temperatures since the Last Glacial Maximum (LGM, ca. 22,000 BP). For the latter, the ultimate consequence of upward shifts may be extinction as species in the highest alpine ecosystems can migrate no further, a phenomenon often characterized as “nowhere to go”. To predict responses to climate change of the alpine plants on the Qinghai-Tibetan Plateau (QTP), we used ecological niche modelling (ENM) to estimate the range shifts of 14 Rhodiola species, beginning with the Last Interglacial (ca. 120,000–140,000 BP) through to 2050. Distributions of Rhodiola species appear to be shaped by temperature-related variables. The southeastern QTP, and especially the Hengduan Mountains, were the origin and center of distribution for Rhodiola, and also served as refugia during the LGM. Under future climate scenario in 2050, Rhodiola species might have to migrate upward and northward, but many species would expand their ranges contra the prediction of the “nowhere to go” hypothesis, caused by the appearance of additional potential habitat concomitant with the reduction of permafrost with climate warming.

Estimating the potential of species to cope with rapid environmental climatic modifications is of vital importance for determining their future viability and conservation. The variation between existing populations along a climatic gradient may predict how a species will respond to future climate change. Stipa purpurea is a dominant grass species in the alpine steppe and meadow of the Qinghai-Tibetan Plateau (QTP). Ecological niche modelling was applied to S. purpurea, and its distribution was found to be most strongly correlated with the annual precipitation and the mean temperature of the warmest quarter. We established a north-to-south transect over 2000 km long on the QTP reflecting the gradients of temperature and precipitation, and then we estimated the morphological by sampling fruited tussocks and genetic divergence by using 11 microsatellite markers between 20 populations along the transect. Reproductive traits (the number of seeds and reproductive shoots), the reproductive-vegetative growth ratio and the length of roots in the S. purpurea populations varied significantly with climate variables. S. purpurea has high genetic diversity (He = 0.585), a large effective population size (Ne >1,000), and a considerable level of gene flow between populations. The S. purpurea populations have a mosaic genetic structure: some distant populations (over 1000 km apart) clustered genetically, whereas closer populations (< 100 km apart) had diverged significantly, suggesting local adaptation. Asymmetrical long-distance inter-population gene flow occurs along the sampling transect and might be mediated by seed dispersal via migratory herbivores, such as the chiru (Pantholops hodgsonii). These findings suggest that population performance variation and gene flow both facilitate the response of S. purpurea to climate change.

Rhodiola crenulata, a well-known medicinal Tibetan herb, is mainly grown in high-altitude regions of the Tibet, Yunnan, and Sichuan provinces in China. In the past few years, increasing numbers of studies have been published on the potential pharmacological activities of R. crenulata, strengthening our understanding into its putitive active ingredient composition, pharmacological activity, and mechanism of action. These findings also provide strong evidence supporting the important medicinal and economical value of R. crenulata. Consequently, some Rhodiola species are becoming endangered because of overexploitation and environmental destruction. However, little is known about the genetic and genomic information of any Rhodiola species. Here we report the first draft assembly ofthe R. crenulata genome, which was 344.5 Mb (25.7 Mb Ns), accounting for 82% of the estimated genome size, with a scaffold N50 length of 144.7 kb and a contig N50 length of 25.4 kb. The R. crenulata genome is not only highly heterozygous but also highly repetitive, with ratios of 1.12% and 66.15%, respectively, based on the k-mer analysis. Furthermore, 226.6 Mb of transposable elements were detected, of which 77.03% were long terminal repeats. In total, 31 517 protein-coding genes were identified, capturing 86.72% of expected plant genes in BUSCO. Additionally, 79.73% of protein-coding genes were functionally annotated. R. crenulata is an important medicinal plant and also a potentially interesting model species for studying the adaptability of Rhodiola species to extreme environments. The genomic sequences of R. crenulata will be useful for understanding the evolutionary mechanism of the stress resistance gene and the biosynthesis pathways of the different medicinal ingredients, for example, salidroside in R. crenulata.