History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau:Aconitum gymnandrum(Ranunculaceae)

Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens)

BACKGROUND

There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity.

RESULTS

A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks.

CONCLUSIONS

Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

Cryptic divergences and repeated hybridizations within the endangered "living fossil" dove tree (Davidia involucrata) revealed by whole genome resequencing

The identification and understanding of cryptic intraspecific evolutionary units (lineages) are crucial for planning effective conservation strategies aimed at preserving genetic diversity in endangered species. However, the factors driving the evolution and maintenance of these intraspecific lineages in most endangered species remain poorly understood. In this study, we conducted resequencing of 77 individuals from 22 natural populations of Davidia involucrata, a "living fossil" dove tree endemic to central and southwest China. Our analysis revealed the presence of three distinct local lineages within this endangered species, which emerged approximately 3.09 and 0.32 million years ago. These divergence events align well with the geographic and climatic oscillations that occurred across the distributional range. Additionally, we observed frequent hybridization events between the three lineages, resulting in the formation of hybrid populations in their adjacent as well as disjunct regions. These hybridizations likely arose from climate-driven population expansion and/or long-distance gene flow. Furthermore, we identified numerous environment-correlated gene variants across the total and many other genes that exhibited signals of positive evolution during the maintenance of two major local lineages. Our findings shed light on the highly dynamic evolution underlying the remarkably similar phenotype of this endangered species. Importantly, these results not only provide guidance for the development of conservation plans but also enhance our understanding of evolutionary past for this and other endangered species with similar histories.

Phylogeography and ecological niche modeling implicate multiple microrefugia of Swertia tetraptera during quaternary glaciations

BACKGROUND

Climate fluctuations during the Pleistocene and mountain uplift are vital driving forces affecting geographic distribution. Here, we ask how an annual plant responded to the Pleistocene glacial cycles.

METHODS

In this study, we analyzed the population demographic history of the annual herb Swertia tetraptera Maxim (Gentianaceae) endemic to Qinghai-Tibetan Plateau (QTP). A total of 301 individuals from 35 populations of S. tetraptera were analyzed based on two maternally inherited chloroplast fragments (trnL-trnF and trnS-trnG). Phylogeographic analysis was combined with species distribution modeling to detect the genetic variations in S. tetraptera.

RESULTS

The genetic diversity of S. tetraptera was high, likely due to its wide natural range, high proportion of endemic haplotypes and evolutionary history. Fifty-four haplotypes were identified in S. tetraptera. Only a few haplotypes were widespread (Hap_4, Hap_1, Hap_3), which were dispersed throughout the present geographical range of S. tetraptera, while many haplotypes were confined to single populations. The cpDNA dataset showed that phylogeographic structuring was lacking across the distribution range of S. tetraptera. Analyses of molecular variance showed that most genetic variation was found within populations (70.51%). In addition, the relationships of the haplotypes were almost completely unresolved by phylogenetic reconstruction. Both mismatch distribution analysis and neutrality tests showed a recent expansion across the distribution range of S. tetraptera. The MaxEnt analysis showed that S. tetraptera had a narrow distribution range during the Last Glacial Maximum (LGM) and a wide distribution range during the current time, with predictions into the future showing the distribution range of S. tetraptera expanding.

CONCLUSION

Our study implies that the current geographic and genetic distribution of S. tetraptera is likely to have been shaped by Quaternary periods. Multiple microrefugia of S. tetraptera existed during Quaternary glaciations. Rapid intraspecific diversification and hybridization and/or introgression may have played a vital role in shaping the current distribution patterns of S. tetraptera. The distribution range of S. tetraptera appeared to have experienced contraction during the LGM; in the future, when the global climate becomes warmer with rising carbon dioxide levels, the distribution of S. tetraptera will expand.

The known, unknown, and the intriguing about members of a critically endangered traditional medicinal plant genus Aconitum

Humanity will always be indebted to plants. In the ongoing scientific era, the 'Herbal Revolution' has helped discover several valuable medicinal plants and associated novel secondary metabolites from the diverse unexplored ecosystems, treating several diseases via phytotherapy. The Aconitum genus comprises several economically-important poisonous mountainous medicinal plant species whose unique biodiversity is on the verge of extinction due to illegal human intervention triggered habitat loss, over-harvesting, and unrestricted trading. Owing to its vast diversity of diterpene alkaloids, most species are extensively used to treat several ailments in rural parts of the world. Irrespective of this, many unexplored and intriguing prospects exist to understand and utilize this critical plant for human benefit. This systematic review tries to fill this gap by compiling information from the sporadically available literature known for ~300 Aconitum spp. regarding its nomenclature and classification, endangerment, plant morphology, ploidy, secondary metabolites, drug pharmacokinetics, conservation, and omics-based computational studies. We also depicted the disparity in the studied model organisms for this diverse genus. The absence of genomic/metagenomic data is becoming a limiting factor in understanding its plant physiology, metabolic pathways, and plant-microbes interactions, and therefore must be promoted. Additionally, government support and public participation are crucial in establishing conservation protocols to save this plant from endangerment.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

Neutral and adaptive genetic diversity in plants: An overview

Genetic diversity is a prerequisite for evolutionary change in all kinds of organisms. It is generally acknowledged that populations lacking genetic variation are unable to evolve in response to new environmental conditions (e.g., climate change) and thus may face an increased risk of extinction. Although the importance of incorporating genetic diversity into the design of conservation measures is now well understood, less attention has been paid to the distinction between neutral (NGV) and adaptive (AGV) genetic variation. In this review, we first focus on the utility of NGV by examining the ways to quantify it, reviewing applications of NGV to infer ecological and evolutionary processes, and by exploring its utility in designing conservation measures for plant populations and species. Against this background, we then summarize the ways to identify and estimate AGV and discuss its potential use in plant conservation. After comparing NGV and AGV and considering their pros and cons in a conservation context, we conclude that there is an urgent need for a better understanding of AGV and its role in climate change adaptation. To date, however, there are only a few AGV studies on non-model plant species aimed at deciphering the genetic and genomic basis of complex trait variation. Therefore, conservation researchers and practitioners should keep utilizing NGV to develop relevant strategies for rare and endangered plant species until more estimates of AGV are available.

Genus Aconitum (Ranunculaceae) in the Ukrainian Carpathians and adjacent territories

Background

The dataset represents a comprehensive collection of occurrence records concerning the genus Aconitum (Ranunculaceae) in the Ukrainian Carpathians and adjacent territories. It is based primarily on the results of critical revision of the main herbarium collections of the Carpathian region (i.e. LW, LWS, LWKS, KRA, KRAM, CHER, KW, UU and KWHU). Besides this, the dataset contains the data parsed (and taxonomically revised) from the published materials and other available sources (e.g. Karel Domin's Card Index).

New information

In total, 2,280 occurrence records of the genus Aconitum representatives distributed in the Ukrainian Carpathians were published.

Out of the Qinghai-Tibetan plateau: Origin, evolution and historical biogeography of Morchella (both Elata and Esculenta clades)

Introduction

Morchella has become a research hotspot because of its wide distribution, delicious taste, and phenotypic plasticity. The Qinghai-Tibet Plateau subkingdoms (QTPs) are known as the cradle of Ice age biodiversity. However, the diversity of Morchella in the QTPs has been poorly investigated, especially in phylogenetic diversity, origin, and biogeography.

Methods

The genealogical concordance phylogenetic species recognition (GCPSR, based on Bayesian evolutionary analysis using sequences from the internal transcribed spacer (ITS), nuclear large subunit rDNA (nrLSU), translation elongation factor 1-α (EF1-α), and the largest and second largest subunits of RNA polymerase II (RPB1 and RPB2)), differentiation time estimation, and ancestral region reconstruction were used to infer Morchella's phylogenetic relationships and historical biogeography in the QTPs.

Results

Firstly, a total of 18 Morchella phylogenetic species are recognized in the QTPs, including 10 Elata clades and 8 Esculenta clades of 216 individuals Secondly, the divergences of the 18 phylogenetic species were 50.24-4.20 Mya (Eocene-Pliocene), which was closely related to the geological activities in the QTPs. Furthermore, the ancestor of Morchella probably originated in the Northern regions (Qilian Shan, Elata cade) and southwestern regions (Shangri-La, Esculenta clade) of QTPs and might have migrated from North America (Rufobrunnea clade) via Beringian Land Bridge (BLB) and Long-Distance Dispersal (LDD) expansions during the Late Cretaceous. Moreover, as the cradle of species origin and diversity, the fungi species in the QTPs have spread out and diffused to Eurasia and South Africa starting in the Paleogene Period.

Conclusion

This is the first report that Esculenta and Elata clade of Morchella originated from the QTPs because of orogenic, and rapid differentiation of fungi is strongly linked to geological uplift movement and refuge in marginal areas of the QTPs. Our findings contribute to increasing the diversity of Morchella and offer more evidence for the origin theory of the QTPs.

Integrating Metabolomics and Transcriptomics to Unveil Atisine Biosynthesis in Aconitum gymnandrum Maxim

Diterpene alkaloids (DAs) are characteristic compounds in Aconitum, which are classified into four skeletal types: C18, C19, C20, and bisditerpenoid alkaloids. C20-DAs are thought to be the precursor of the other types. Their biosynthetic pathway, however, is largely unclear. Herein, we combine metabolomics and transcriptomics to unveil the methyl jasmonate (MJ) inducible biosynthesis of DAs in the sterile seedling of A. gymnandrum, the only species in the Subgenus Gymnaconitum (Stapf) Rapaics. Target metabolomics based on root and aerial portions identified 51 C19-DAs and 15 C20-DAs, with 40 inducible compounds. The highest content of C20-DA atisine was selected for further network analysis. PacBio Isoform sequencing integrated with RNA sequencing not only provided the full-length transcriptome but also their response to induction, revealing 1994 genes that exhibited up-regulated expression. Further, 38 genes involved in terpenoid biosynthesis were identified, including 7 diterpene synthases. In addition to the expected function of the four diterpene synthases, AgCPS5 was identified to be a new ent-8,13-CPP synthase in Aconitum and could also combine with AgKSL1 to form the C20-DAs precursor ent-atiserene. Combined with multiple network analyses, six CYP450 and seven 2-ODD genes predicted to be involved in the biosynthesis of atisine were also identified. This study not only sheds light on diterpene synthase evolution in Aconitum but also provides a rich dataset of full-length transcriptomes, systemic metabolomes, and gene expression profiles, setting the groundwork for further investigation of the C20-DAs biosynthesis pathway.

Genetic diversity and population dynamic of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow in Central China

Phylogeographic research concerning Central China has been rarely conducted. Population genetic and phylogeography of Ziziphus jujuba var. spinosa (also called sour jujube) were investigated to improve our understanding of plant phylogeographic patterns in Central China. Single-copy nuclear gene markers and complete chloroplast genome data were applied to 328 individuals collected from 21 natural populations of sour jujube in China. Nucleotide variation of sour jujube was relatively high (π = 0.00720, θ w = 0.00925), which resulted from the mating system and complex population dynamics. Analysis of molecular variation analysis revealed that most of the total variation was attributed to variation within populations, and a high level of genetic differentiation among populations was detected (F st = 0.197). Relatively low long-distance dispersal capability and vitality of pollen contributed to high genetic differentiation among populations. Differences in the environmental conditions and long distance among populations further restricted gene flow. Structure clustering analysis uncovered intraspecific divergence between central and marginal populations. Migrate analysis found a high level of gene flow between these two intraspecific groups. Bayesian skyline plot detected population expansion of these two intraspecific groups. Network and phylogeny analysis of chloroplast haplotypes also found intraspecific divergence, and the divergence time was estimated to occur at about 55.86 Ma. Haplotype native to the Loess Plateau was more ancient, and multiple glacial refugia of sour jujube were found to locate at the Loess Plateau, areas adjacent to the Qinling Mountains and Tianmu Mountains. Species distribution model analysis found a typical contraction-expansion model corresponding to the Quaternary climatic oscillations. In the future, the distribution of sour jujube may shift to high-latitude areas. This study provides new insights for phylogeographic research of temperate plant species distributed in Central China and sets a solid foundation for the application of the scientific management strategy of Z. jujuba var. spinosa.

Dispersal into the Qinghai-Tibet plateau: evidence from the genetic structure and demography of the alpine plant Triosteum pinnatifidum

Triosteum pinnatifidum Maxim., an alpine plant, is traditionally used for several medicinal purposes. Here, both chloroplast DNA sequences and nuclear low copy sequence markers were used to investigate the genetic diversity and population structure of T. pinnatifidum. Materials were collected from thirteen localities in the northeast Qinghai-Tibet Plateau (QTP) and adjacent highlands and advanced analytical toolkits were used to access their origin and range shifts. The results revealed a higher level of population differentiation based on chloroplast DNA (cpDNA) concatenated sequences compared with the nuclear DNA sequences (F _ST = 0.654 for cpDNA, F _ST = 0.398 for AT103), indicating that pollen flow was still extensive in T. pinnatifidum. A decline in haplotype variation was observed from the plateau edge and adjoining highlands toward the platform of the QTP. The hypothesis "dispersal into the QTP," proposing that T. pinnatifidum experienced migration from the plateau edge and adjacent highlands to the platform, was supported. These results were in line with the hypothesis that multiple refugia exist on the plateau edge and adjacent highlands rather than on the plateau platform. Our unimodal mismatch distribution, star-like network supported a recent expansion in T. pinnatifidum.

Phylogeography of Iris loczyi (Iridaceae) in Qinghai-Tibet Plateau revealed by chloroplast DNA and microsatellite markers

Quaternary climate oscillations and complex topography have tremendous effects on current distribution and genetic structure of species, and hence the Qinghai-Tibet Plateau (QTP), the largest plateau in the world, has become a hotspot for many phylogeographic studies. However, little is known about the phylogeographic pattern of herbaceous plants in QTP. Here, we investigate the genetic diversity, population structure and historical dynamics of Iris loczyi, using five chloroplast DNA (cpDNA) fragments and seven microsatellite markers. A total of 15 populations, and 149 individuals were sampled throughout the QTP. High genetic diversity was detected both in cpDNA (H _d = 0.820) and SSR (H _o = 0.689, H _e = 0.699). Ten cpDNA haplotypes and 163 alleles were identified. AMOVA and clustering analyses revealed obvious differentiation between regions. The N _st, G _st and Mantel test showed significant phylogeographic structure of I. loczyi. The neutrality test and mismatch distribution analyses indicated that I. loczyi could not have undergone a historical population expansion, but population XS from the Qilian Mountain area could have experienced a local expansion. Bottleneck analyses indicated that I. loczyi had not experienced bottleneck recently. Based on cpDNA and SSR results, the Qilian Mountain area was inferred as a potential glacial refuge, and the southern Tibet valley was considered as a 'microrefugia' for I. loczyi. These findings provided new insights into the location of glacial refuges for the species distributed in QTP, and supplemented more plant species data for the response of QTP species to the Quaternary climate.

Evolutionary history of two rare endemic conifer species from the eastern Qinghai-Tibet Plateau

BACKGROUND AND AIMS

Understanding the population genetics and evolutionary history of endangered species is urgently needed in an era of accelerated biodiversity loss. This knowledge is most important for regions with high endemism that are ecologically vulnerable, such as the Qinghai-Tibet Plateau (QTP).

METHODS

The genetic variation of 84 juniper trees from six populations of Juniperus microsperma and one population of Juniperus erectopatens, two narrow-endemic junipers from the QTP that are sister to each other, was surveyed using RNA-sequencing data. Coalescent-based analyses were used to test speciation, migration and demographic scenarios. Furthermore, positively selected and climate-associated genes were identified, and the genetic load was assessed for both species.

KEY RESULTS

Analyses of 149 052 single nucleotide polymorphisms showed that the two species are well differentiated and monophyletic. They diverged around the late Pliocene, but interspecific gene flow continued until the Last Glacial Maximum. Demographic reconstruction by Stairway Plot detected two severe bottlenecks for J. microsperma but only one for J. erectopatens. The identified positively selected genes and climate-associated genes revealed habitat adaptation of the two species. Furthermore, although J. microsperma had a much wider geographical distribution than J. erectopatens, the former possesses lower genetic diversity and a higher genetic load than the latter.

CONCLUSIONS

This study sheds light on the evolution of two endemic juniper species from the QTP and their responses to Quaternary climate fluctuations. Our findings emphasize the importance of speciation and demographic history reconstructions in understanding the current distribution pattern and genetic diversity of threatened species in mountainous regions.

Plastome sequencing reveals phylogenetic relationships among Comastoma and related taxa (Gentianaceae) from the Qinghai-Tibetan Plateau

Genus Comastoma (subt. Swertiinae, Gentianaceae) contains species, such as "Zangyinchen," that are important herbs in Tibetan medicine. The phylogenetic relationship of this within Gentianaceae and the circumscriptions of its species have long been controversial with conflicting morphological and molecular data reported. Here, we used whole chloroplast genome sequences for Comastoma species and related taxa to reconstruct their phylogeny and clarify their taxonomic relationships. The results revealed that the length of all plastome sequenced varied from 149 to 151 kb and have high similarity in structure and gene content. Phylogenomic analysis showed that Comastoma is a monophyletic group, closely related to the genus Lomatogonium. The divergence time estimation showed that Gentianaceae diverged at about 21.81 Ma, while the split of Comastoma occurred at 7.70 Ma. However, the results suggested the crown age of species formation in this genus is after 4.19 Ma. Our results suggest that QTP uplift, the alternation of Quaternary glaciation and interglaciation, and monsoon changes might have acted as drivers of speciation in Comastoma.

The function of floral traits and phenotypic selection in Aconitum gymnandrum (Ranunculaceae)

Floral evolution in angiosperms is thought to be driven by pollinator-mediated selection. Understanding flower integration and adaptation requires resolving the additive and nonadditive contributions of floral pollinator attraction and pollination efficiency traits to fitness components. In this study, a flower manipulation experiment with a factorial design was used to study the adaptive significance of galea height (a putative attraction trait) and entrance width (a putative efficiency trait) in Aconitum gymnandrum Maxim. flowers. Simultaneously, phenotypic selection analysis was conducted to examine selection by pollinators on galea height, entrance width, nectar production and plant height. Increased galea height increased the pollinator visitation rate, which confirmed its attractiveness function. Increasing floral entrance width by spreading the lower sepals increased the seed number per fruit without affecting pollinator visitation. This suggests a pollination efficiency role for the entrance width. The phenotypic selection analysis, however, did not provide evidence of pollinator-mediated selection for either of these traist, but it did for plant height. According to the manipulation treatment and correlational selection results, the combined variation in galea height and entrance width of A. gymnandrum flowers did not have nonadditive effects on female reproductive success. This study demonstrated the adaptive value of A. gymnandrum flowers through manipulation of an attractiveness trait and an efficiency trait. However, neither trait was associated with pollinator-mediated selection. A combination of manipulating traits and determining current phenotypic selection could help to elucidate the mechanism of selection on floral traits involved in different functions and flower integration.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

The fate of Meconopsis species in the Tibeto-Himalayan region under future climate change

High-mountain areas such as the Tibeto-Himalayan region (THR) host cold-adapted biota expected to be sensitive to anthropogenic climate change. Meconopsis is a representative endangered genus confined to alpine meadow or subnival habitats in the THR. We used climate-niche factor analysis to study the vulnerability of ten Meconopsis species to climate change, comparing current climate (representative of 1960-1990) to future climate scenarios (2070: average 2061-2080). For these ten Meconopsis species, we then identified potential future climate refugia and determined optimal routes for each species to disperse to the proposed refugia. Our results indicate that for the ten Meconopsis species, the regions with low vulnerability to climate change in the THR are the central Qinghai-Tibet Plateau, the Hengduan Mountains (HDM), the eastern Himalayas, and the West Qinling Mountain (WQL), and can be considered potential future climate refugia. Under future climate change, we found for the ten Meconopsis species potential dispersal routes to three of the four identified refugia: the HDM, the eastern Himalayas, and the WQL. Our results suggest that past refugia on the THR will also be the future climate refugia for the ten Meconopsis species, and these species may potentially persist in multiple future climate refugia, likely reducing risks from climate change. Furthermore, climate change may affect the threat ranking of Red Listed Species for Meconopsis species, as Least Concern species were estimated to become more vulnerable to climate change than the only Near Threatened species.

Phylogeographic analyses of poplar revealed potential glacial refugia and allopatric divergence in southwest China

The distribution pattern and genetic structure of plant species have been profoundly influenced by climate oscillations. Phylogeographic analyses have been numerously conducted in biodiversity hotspot regions and some general conclusions have been documented. However, other regions have received less attentions and these places may serve as potential glacial refugia for plant species to survive the Quaternary glaciation. Here, we used six nuclear and three cpDNA markers to estimate the phylogeographic pattern of Populus davidiana, a widespread species distributed in southwest China. As a widely distributed species in southwest China, the nucleotide diversity of P. davidiana was relatively high (N _a = 6.28, H _O = 0.534, and H _E = 0.658). Genetic differentiation (F _ST) between the two main distribution regions, Yunnan and Guizhou provinces, was 0.21221. According to the composition of chloroplast haplotypes and the result of structure in these populations, we clearly distinguished two distantly sublineages corresponding to two distribution regions. Results of the Mantel test showed that there was a significant correlation between genetic distance and geographical distance (R ² = 0.8252, p<.05). The topographically heterogeneous regions and the low dispersal ability of seed and pollen may lead to high genetic differentiation between these two regions. A potential glacial refugia for P. davidiana located in adjacent regions to the Hengduan range was revealed and allopatric divergence in separated glacial refugia may directly lead to the present phylogeographic pattern of this species.

Deciphering the archaeal communities in tree rhizosphere of the Qinghai-Tibetan plateau

BACKGROUND

The Qinghai-Tibetan Plateau represents one of the most important component of the terrestrial ecosystem and a particularly vulnerable region, which harbouring complex and diverse microbiota. The knowledge about their underground microorganisms have largely been studied, but the characteristics of rhizosphere microbiota, particularly archaeal communities remains unclear.

RESULTS

High-throughput Illumina sequencing was used to investigate the rhizosphere archaeal communities of two native alpine trees (Picea crassifolia and Populus szechuanica) living on the Qinghai-Tibetan Plateau. The archaeal community structure in rhizospheres significantly differed from that in bulk soil. Thaumarchaeota was the dominant archaeal phylum in all soils tested (92.46-98.01%), while its relative abundance in rhizospheres were significantly higher than that in bulk soil. Ammonium nitrogen, soil organic matter, available phosphorus and pH were significantly correlated with the archaeal community structure, and the deterministic processes dominated the assembly of archaeal communities across all soils. In addition, the network structures of the archaeal community in the rhizosphere were less complex than they were in the bulk soil, and an unclassified archaeal group (Unclassified_k_norank) was identified as the keystone species in all archaeal networks.

CONCLUSIONS

Overall, the structure, assembly and co-occurrence patterns of archaeal communities are significantly affected by the presence of roots of alpine trees living on the Qinghai-Tibetan Plateau. This study provides new insights into our understanding of archaeal communities in vulnerable ecosystems.

The complete plastome of Nannoglottis ravida, an extremely endangered species in the Qinghai-Tibet Plateau

Nannoglottis ravida is an extremely endangered species in the Qinghai-Tibet Plateau. Based on the second-generation high-throughput genome sequencing, we assembled the plastome of this species. The length of the total plastome is 152,324 bp with a typical quadripartite structure including a large single-copy region of 83,708 bp, a small single-copy region of 29,882 bp and two reverse repeat regions of 19,367 bp respectively. A total of 131 genes were annotated including 85 protein-coding genes (PCG), 36 tRNA genes, 8 rRNA genes and 2 pseudogenes. The constructed phylogenetic tree with other species of two tribes Senecioneae and Astereae based on plastomes suggests that N. ravida has a close relationship with the Astereae, but diverged early from this tribe.

Spatial genetic structure and historical demography of East Asian wild boar

Pleistocene climatic fluctuations may have had a profound impact on the evolutionary history of many species. The geographical pattern of European wild boar (Sus scrofa) is clearly studied, and it was greatly influenced by ancient climatic events, especially the Last Glacial Maximum. Previous research on genetic variation has mainly focused on the origin and distribution histories of domestic pigs. However, some questions have not been answered, including those concerning the genetic diversity, geographical pattern and possible historic influence of climate on East Asian wild boar (EAWB). Employing the control region of mtDNA (511 bp), we investigated the contributions of historic climate, which possibly shaped the genetic pattern of wild boar. Given that the level of genetic diversity of wild boars is higher in East Asia than in Europe, 172 haplotypes were detected from 680 individuals. Phylogenetic analysis demonstrated the complex phylogeographic structure of EAWB. Mismatch analysis, neutrality tests and the Bayesian Skyline Plot results all retrieved signals of a rapid population expansion, which might have played an important role in driving the formation of complex spatial genetic structure. Genetic data and species distribution modelling showed that the Last Glacial Maximum had weak effect on the distribution of the EAWB. We suggest that, in shaping spatial genetic structure in East Asian, long-term gene flow and population history played more important roles than Pleistocene climate fluctuations.

A comparative plastomics approach reveals available molecular markers for the phylogeographic study of Dendrobium huoshanense, an endangered orchid with extremely small populations

Comparative plastomics approaches have been used to identify available molecular markers for different taxonomic level studies of orchid species. However, the adoption of such methods has been largely limited in phylogeographic studies. Therefore, in this study, Dendrobium huoshanense, an endangered species with extremely small populations, was used as a model system to test whether the comparative plastomic approaches could screen available molecular markers for the phylogeographic study. We sequenced two more plastomes of D. huoshanense and compared them with our previously published one. A total of 27 mutational hotspot regions and six polymorphic cpSSRs have been screened for the phylogeographic studies of D. huoshanense. The cpDNA haplotype data revealed that the existence of haplotype distribution center was located in Dabieshan Mts. (Huoshan). The genetic diversity and phylogenetic analyses showed that the populations of D. huoshanense have been isolated and evolved independently for long period. On the contrary, based on cpSSR data, the genetic structure analysis revealed a mixed structure among the populations in Anhui and Jiangxi province, which suggested that the hybridization or introgression events have occurred among the populations of D. huoshanense. These results indicated that human activities have played key roles in shaping the genetic diversity and distributional patterns of D. huoshanense. According to our results, both two markers showed a high resolution for the phylogeographic studies of D. huoshanense. Therefore, we put forth that comparative plastomic approaches could revealed available molecular markers for phylogeographic study, especially for the species with extremely small populations.

Revealing Fungal Communities in Alpine Wetlands Through Species Diversity, Functional Diversity and Ecological Network Diversity

The biodiversity of fungi, which are extremely important in maintaining the ecosystem balance in alpine lakeside wetlands, has not been fully studied. In this study, we investigated the fungal communities of three lakeside wetlands from different altitudes in the Qinghai-Tibet Plateau and its edge. The results showed that the fungi of the alpine lakeside wetland had higher species diversity. Functional annotation of fungi by FUNGild software showed that saprophytic fungi were the most abundant type in all three wetlands. Further analysis of the microbial phylogenetic molecular ecological network (pMEN) showed that saprophytic fungi are important species in the three wetland fungal networks, while symbiotic fungi and pathotrophic fungi have different roles in the fungal networks in different wetlands. Community diversity was high in all three lakeside wetlands, but there were significant differences in the composition, function and network structure of the fungal communities. Contemporary environmental conditions (soil properties) and historical contingencies (geographic sampling location) jointly determine fungi community diversity in this study. These results expand our knowledge of fungal biodiversity in the alpine lakeside wetlands.

Molecular Phylogeography and Evolutionary History of the Endemic Species Corydalis hendersonii (Papaveraceae) on the Tibetan Plateau Inferred From Chloroplast DNA and ITS Sequence Variation

In response to past climatic changes, the species with different habits or adaptive traits likely have experienced very different evolutionary histories, especially for species that restricted to high mountain areas. In order to trace how Quaternary climatic oscillations affected range distributions and intraspecific divergence of such alpine plants on the Tibetan Plateau, here, we investigated maternally inherited chloroplast DNA (cpDNA) markers and biparentally inherited nuclear ribosomal internal transcribed spacer (ITS) DNA variations and aimed to explore the phylogeographic history of the endemic alpine species Corydalis hendersonii Hemsl. (Papaveraceae). We sequenced four cpDNA fragments (trnS-trnG, trnT-trnL, atpH-atpI, and psbE-petL) and also the nuclear (ITS) region in 368 individuals from 30 populations across the species' range. The network and phylogenetic analysis based on cpDNA variations identified 15 chlorotypes that cluster into three distinct clades. However, our nuclear DNA results demonstrated that there were four genetic/geographical groups within C. hendersonii. Some common and highly divergent cpDNA and ITS haplotypes were distributed in the populations of central and northeastern Tibetan Plateau, and the highest nucleotide diversity and genetic differentiation were detected in the central region. Demographic tests further indicated that the populations of southwestern and western Tibet may have experienced recent range expansion, which most likely occurred during the last glacial maximum (LGM) and continued its expansion after the beginning of the Holocene. These two different groups of this species may have derived from potential refugia that existed in the central and/or northeastern regions of Tibet during recent interglacial periods. In addition, our AMOVA analyses detected high genetic differentiation along with the whole sampling range. Also, distinct phylogeographic structures were detected among populations of C. hendersonii based on both of cpDNA and ITS variation. These findings shed new light on the importance of climatic oscillations during Quaternary and complex local topography as causes of intraspecific diversification and demographic changes within cold-tolerant herbs in the Tibetan Plateau biodiversity hotspot.

Evolutionary history of a relict conifer, Pseudotaxus chienii (Taxaceae), in south-east China during the late Neogene: old lineage, young populations

BACKGROUND AND AIMS

Many monotypic gymnosperm lineages in south-east China paradoxically remain in relict status despite long evolutionary histories and ample opportunities for allopatric speciation, but this paradox has received little attention and has yet to be resolved. Here, we address this issue by investigating the evolutionary history of a relict conifer, Pseudotaxus chienii (Taxaceae).

METHODS

DNA sequences from two chloroplast regions and 14 nuclear loci were obtained for 134 samples. The demographic history was inferred and the contribution of isolation by environment (IBE) in patterning genetic divergence was compared with that of isolation by distance (IBD).

KEY RESULTS

Three genetic clusters were identified. Approximate Bayesian computation analyses showed that the three clusters diverged in the late Pliocene (~3.68 Ma) and two admixture events were detected. Asymmetric gene flow and similar population divergence times (~ 3.74 Ma) were characterized using the isolation with migration model. Neither IBD nor IBE contributed significantly to genetic divergence, and the contribution of IBE was much smaller than that of IBD.

CONCLUSIONS

These results suggest that several monotypic relict gymnosperm lineages like P. chienii in south-east China did not remain in situ and undiversified for millions of years. On the contrary, they have been evolving and the extant populations have become established more recently, having insufficient time to speciate. Our findings provide a new perspective for understanding the formation and evolution of the relict gymnosperm flora of China as well as of the Sino-Japanese Flora.

Floral Organogenesis in Three Members of the Tribe Delphinieae (Ranunculaceae)

Three species (Aconitum taipeicum, Delphinium giraldii, and Consolida ajacis) of the tribe Delphinieae (Ranunculaceae) were examined using scanning electron microscopy and histological methods. The results showed that members of Delphinieae differ from their polysymmetrical relatives by four unique features: (1) a spiral phyllotaxis of their perianth and stamens, and a series of carpels, which initiated superficially in a whorl-liked arrangement; (2) sepal 2 being the largest one among the five sepals and becoming helmet-shaped or having a spur; (3) petals 2 and 5 initiated adaxially of sepal 2 and also becoming spurred; and (4) the monosymmetry of the first flower becoming established when sepal 2 becomes the largest. Major differences among the species include the timing of development of the second series; the fusion of two petals into a single one in C. ajacis; and, during early developmental stages, the two young spurred petals giving rise to a stalk and two bulges in A. taipeicum, a single bulge in D. giraldii, or an arch blade in C. ajacis. The unequal growth of the perianth, together with the reduction and the rearrangement of the carpels, are critical in inducing the symmetrical transformation of the flowers.

Development of genomic microsatellite markers for Aconitum gymnandrum (Ranunculaceae) by next generation sequencing (NGS)

Mating plays key roles in the demographic and genetic dynamics of populations. Estimates of mating portfolios and system based on progeny array (PA) method required highly polymorphic genetic markers, of which microsatellite is a good choice. In this study, we reported 19 polymorphic microsatellite loci for Aconitum gymnandrum. The number of alleles per locus ranged from 2 to 12. Observed and expected heterozygosity ranged from 0.000 to 1.000 and from 0.219 to 0.842, respectively. Seven loci showed significant deviation from Hardy-Weinberg equilibrium. These markers will provide a useful tool for pollination ecology and population genetic studies of A. gymnandrum in Qinghai-Tibet plateau.

Spatial Genetic Patterns and Distribution Dynamics of the Rare Oak Quercus chungii: Implications for Biodiversity Conservation in Southeast China

A rapidly changing climate and frequent human activity influences the distribution and community structure of forests. Increasing our knowledge about the genetic diversity and distribution patterns of trees is helpful for forest conservation and management. In this study, nSSRs (nuclear simple sequence repeats) were integrated with a species distribution model (SDM) to investigate the spatial genetic patterns and distribution dynamics of Quercus chungii F.P.Metcalf, a rare oak in the subtropics of southeast China. A total of 188 individuals from 11 populations distributed across the natural range of Q. chungii were genotyped using nine nSSRs. The STRUCTURE analysis indicated that genetic admixture was present in all populations, but the population genetic variation and genetic differentiation were related to their geographical distributions. The SDM result indicated that Q. chungii retreated to the Nanling Mountains and adjacent areas during the Last Glacial Maximum (LGM) period, which corresponds to higher genetic diversity for populations in this region. Landscape genetic analysis showed that the Nanling Mountains served as a corridor for organism dispersal at the glacial and interglacial periods within the Quaternary. Based on these results, we propose that establishing nature reserves to protect the ecological corridor across the Nanling Mountains is necessary for the conservation of regional species genetic diversity, as well as the ecosystem of evergreen broadleaved forests in southern China. The study combines species distribution models and genetic diversity to provide new insight into biodiversity conservation and forest management under future climate change.

Complex population evolutionary history of four cold-tolerant Notopterygium herb species in the Qinghai-Tibetan Plateau and adjacent areas

Historical geological and climatic events are the most important drivers of population expansions/contractions, range shifts, and interspecific divergence in plants. However, the species divergence and spatiotemporal population dynamics of alpine cold-tolerant herbal plants in the high-altitude Qinghai-Tibetan Plateau (QTP) and adjacent areas remain poorly understood. In this study, we investigated population evolutionary history of four endangered Notopterygium herb species in the QTP and adjacent regions. We sequenced 10 nuclear loci, 2 mitochondrial DNA regions, and 4 chloroplast DNA regions in a total of 72 natural populations from the 4 species, and tested the hypothesis that the population history of these alpine herbs was markedly affected by the Miocene-Pliocene QTP uplifts and Quaternary climatic oscillations. We found that the four Notopterygium species had generally low levels of nucleotide variability within populations. Molecular dating and isolation-with-migration analyses suggested that Notopterygium species diverged ~1.74-7.82 million years ago and their differentiation was significantly associated with recent uplifts of the eastern margin of the QTP. In addition, ecological niche modeling and population history analysis showed that N. incisum and N. franchetii underwent considerable demographic expansions during the last glacial period of the Pleistocene, whereas a demographic contraction and a expansion occurred for N. forrestii and N. oviforme during the antepenultimate interglacial period and penultimate glacial period, respectively. These findings highlight the importance of geological and climatic changes during the Miocene-Pliocene and Pleistocene as causes of species divergence and changes in population structure within cold-tolerant herbs in the QTP biodiversity hotspot.

Quaternary climate change and habitat preference shaped the genetic differentiation and phylogeography of Rhodiola sect. Prainia in the southern Qinghai-Tibetan Plateau

There are two long-standing biogeographic hypotheses regarding the glacial survival of plant species in the Qinghai-Tibetan Plateau (QTP): the in situ survival hypothesis and the tabula rasa hypothesis. We tested these two hypotheses in a phylogeographic study of Rhodiola sect. Prainia, a monophyletic section with ecologically divergent lineages. Molecular data from the nuclear internal transcribed spacer, six plastid markers and 13 nuclear microsatellite loci were analyzed for 240 individuals from 19 populations of this section. Environmental data were used to analyze the niches of major phylogenetic lineages within this section and to model changes in their distributions since the Last Glacial Maximum (LGM). We found that Rhodiola sect. Prainia consists of three evolutionary lineages: all populations of R. stapfii, R. prainii populations at the southern edge of the QTP, and R. prainii populations in the interior part of the QTP. During the LGM, the survival of R. prainii in the interior part of the QTP corresponded with the in situ survival hypothesis, while R. stapfii most probably survived the LGM in a manner corresponding with the tabula rasa hypothesis. The evolutionary history of different lineages of this section was shaped by topography, climate change, and lineage-specific habitat preferences.

Recent Trends in Research on the Genetic Diversity of Plants: Implications for Conservation

Genetic diversity and its distribution, both within and between populations, may be determined by micro-evolutionary processes, such as the demographic history of populations, natural selection, and gene flow. In plants, indices of genetic diversity (e.g., k, h and π) and structure (e.g., FST) are typically inferred from sequences of chloroplast markers. Given the recent advances and popularization of molecular techniques for research in population genetics, phylogenetics, phylogeography, and ecology, we adopted a scientometric approach to compile evidence on the recent trends in the use of cpDNA sequences as markers for the analysis of genetic diversity in botanical studies, over the years. We also used phylogenetic modeling to assess the relative contribution of relatedness or ecological and reproductive characters to the genetic diversity of plants. We postulated that genetic diversity could be defined not only by microevolutionary factors and life history traits, but also by relatedness, so that species more closely related phylogenetically would have similar genetic diversities. We found a clear tendency for an increase in the number of studies over time, confirming the hypothesis that the advances in the area of molecular genetics have supported the accumulation of data on the genetic diversity of plants. However, we found that the vast majority of these data have been produced by Chinese authors, and refer specifically to populations of Chinese plants. Most of the data on genetic diversity have been obtained for species in the International Union for Conservation of Nature (IUCN) category NE (Not Evaluated), which indicates a relative lack of attention on threatened species. In general, we observed very high FST values in the groups analyzed and, as we focused primarily on species that have not been evaluated by the IUCN, the number of plant species that are threatened with extinction may be much greater than that indicated by the listing of this organization. We also found that the number of haplotypes (k) was influenced by the type of geographic distribution of the plant, while haplotype diversity (h) was affected by the type of flower, and the fixation index (FST), by the type of habitat. The plant species most closely-related phylogenetically have similar levels of genetic diversity. Overall, then, it will important to consider phylogenetic dependence in future studies that evaluate the effects of life-history traits on plant genetic diversity.

Evolutionary History of Rhus chinensis (Anacardiaceae) From the Temperate and Subtropical Zones of China Based on cpDNA and Nuclear DNA Sequences and Ecological Niche Model

To explore the origin and evolution of local flora and vegetation, we examined the evolutionary history of Rhus chinensis, which is widely distributed in China’s temperate and subtropical zones, by sequencing three maternally inherited chloroplast DNAs (cpDNA: trnL-trnF, psbA-trnH, and rbcL) and the biparentally inherited nuclear DNA (nuDNA: LEAFY) from 19 natural populations of R. chinensis as well as the ecological niche modeling. In all, 23 chloroplast haplotypes (M1–M23) and 15 nuclear alleles (N1–N15) were detected. The estimation of divergence time showed that the most recent common ancestor dated at 4.2 ± 2.5 million years ago (Mya) from cpDNA, and the initial divergence of genotypes occurred at 4.8 ± 3.6 Mya for the nuDNA. Meanwhile, the multimodality mismatch distribution curves and positive Tajima’s D values indicated that R. chinensis did not experience population expansion after the last glacial maximum. Besides, our study was also consistent with the hypothesis that most refugia in the temperate and subtropical zones of China were in situ during the glaciation.

Mountains as Evolutionary Arenas: Patterns, Emerging Approaches, Paradigm Shifts, and Their Implications for Plant Phylogeographic Research in the Tibeto-Himalayan Region

Recently, the "mountain-geobiodiversity hypothesis" (MGH) was proposed as a key concept for explaining the high levels of biodiversity found in mountain systems of the Tibeto-Himalayan region (THR), which comprises the Qinghai-Tibetan Plateau, the Himalayas, and the biodiversity hotspot known as the "Mountains of Southwest China" (Hengduan Mountains region). In addition to the MGH, which covers the entire life span of a mountain system, a complementary concept, the so-called "flickering connectivity system" (FCS), was recently proposed for the period of the Quaternary. The FCS focuses on connectivity dynamics in alpine ecosystems caused by the drastic climatic changes during the past ca. 2.6 million years, emphasizing that range fragmentation and allopatric speciation are not the sole factors for accelerated evolution of species richness and endemism in mountains. I here provide a review of the current state of knowledge concerning geological uplift, Quaternary glaciation, and the main phylogeographic patterns ("contraction/recolonization," "platform refugia/local expansion," and "microrefugia") of seed plant species in the THR. In addition, I make specific suggestions as to which factors future avenues of phylogeographic research should take into account based on the fundamentals presented by the MGH and FCS, and associated complementary paradigm shifts.

Molecular Phylogeography and Ecological Niche Modeling of Sibbaldia procumbens s.l. (Rosaceae)

The phylogeographical analysis and ecological niche modeling (ENM) of the widely distributed Northern Hemisphere Sibbaldia procumbens s.l. can help evaluate how tectonic motion and climate change helped shape the current distribution patterns of this species. Three chloroplast regions (the atpI-atpH and trnL-trnF intergenic spacers and the trnL intron) were obtained from 332 (156 from present study and 176 from the previous study) individuals of S. procumbens s.l. An unrooted haplotype network was constructed using the software NETWORK, while BEAST was used to estimate the divergence times among haplotypes. ENM was performed by MAXENT to explore the historical dynamic distribution of S. procumbens s.l. The haplotype distribution demonstrates significant phylogeographical structure (N_ST > G_ST; P < 0.01). The best partitioning of genetic diversity by SAMOVA produced three groups, while the time to the most recent common ancestor of all haplotypes was estimated to originate during the Miocene, with most of the haplotype diversity having occurred during the Quaternary. The MAXENT analysis showed S. procumbens s.l. had a wider distribution range during the last glacial maximum and a narrower distribution range during the last interglacial, with predictions into the future showing the distribution range of S. procumbens s.l. shrinking.

A Phytogeographic Divide Along the 500 mm Isohyet in the Qinghai-Tibet Plateau: Insights From the Phylogeographic Evidence of Chinese Alliums (Amaryllidaceae)

The Qinghai-Tibet Plateau (QTP) has been biogeographically divided into the eastern monsoonal and the western continental climatic zones along the 500 mm isohyet. However, this biogeographic hypothesis has been rarely tested using a phylogeographic approach. The members of the genus Allium subgenus Cyathophora coincidentally distribute across this biogeographical divide. Intriguingly, Allium fasciculatum of subgenus Amerallium co-occurs in the distribution range of subgenus Cyathophora. To illuminate the role of this biogeographic divide on the genetic divergence, we genotyped 466 individuals of 52 populations of subgenus Cyathophora and 110 individuals of 19 populations of A. fasciculatum using three chloroplast DNA fragments, whole nrITS and nine nuclear microsatellite loci, supplemented with the present environmental space and paleo-distribution modeling. Our phylogeographical evidence recovered the concordant east-west genetic breaks both for subgenus Cyathophora and A. fasciculatum along the 500 mm isohyet. The divergence time estimations and environmental niche differentiations suggested this east-west genetic breaks could have been triggered by the climatic-induced vicariance during the early Pleistocene. Noticeably, this split within subgenus Cyathophora could have been deepened by the morphological vicariance from the eastern umbel to the western spicate, while that within A. fasciculatum could have been obscured due to the pollen flows from the east to west caused by the postglacial expansion. The genetic structures and ecological niche modelings (ENMs) recovered the distinct responses to the Quaternary climatic oscillations for species constricted to different climatic zones, further highlighting the profound effect of the climatic differences and tectonic uplifts on the genetic diversification. Overall, our findings offer strong evidence for the existence of a biogeographic divide between the eastern monsoonal and the west continental climatic zones of the QTP nearly along the 500 mm isohyet.

Geophysical upheavals and evolutionary diversification of plant species in the Himalaya

The Himalaya is one of the youngest and the loftiest mountain chains of the world; it is also referred to as the water tower of Asia. The Himalayan region harbors nearly 10,000 plant species constituting approximately 2.5% of the global angiosperm diversity of which over 4,000 are endemics. The present-day Himalayan flora consists of an admixture of immigrant taxa and diversified species over the last 40 million years. The interesting questions about the Himalayan flora discussed here are: how did the Himalaya achieve high endemic plant diversity starting with immigrant taxa and what were the main drivers of this diversity? This contribution aims to answer these questions and raise some more. We review and analyze existing information from diverse areas of earth and climate sciences, palaeobiology and phytogeography to evolve a bio-chronological record of plant species divergence and evolution in the Himalaya. From the analysis we infer the effects of major environmental upheavals on plant diversity in the region. The understanding developed in the following discussion is based on the idea that Himalaya experienced at least five phases of major geophysical upheavals, namely: (i) mega-collision between India and Eurasian plates, (ii) tectonic uplift in phases and progressive landform elevation, (iii) onset of southwest (SW) Indian monsoon, (iv) spurring of arid conditions in Central Asia, and (v) cyclic phases of cooling and warming in the Quaternary. The geophysical upheavals that were potentially disrupting for the ecosystem stability had a key role in providing impetus for biological diversification. The upheavals produced new geophysical environments, new ecological niches, imposed physical and physiological isolation barriers, acted as natural selection sieves and led to the formation of new species. This contribution aims to develop a comprehensive understanding of the plant biodiversity profile of the Himalaya in the context of complex, interconnected and dynamic relationship between earth system processes, climate and plant diversity.

Deep Intraspecific Divergence in the Endemic Herb Lancea tibetica (Mazaceae) Distributed Over the Qinghai-Tibetan Plateau

Qinghai-Tibetan Plateau (QTP) is an important biodiversity hub, which is very sensitive to climate change. Here in this study, we investigated genetic diversity and past population dynamics of Lancea tibetica (Mazaceae), an endemic herb to QTP and adjacent highlands. We sequenced chloroplast and nuclear ribosomal DNA fragments for 429 individuals, collected from 29 localities, covering their major distribution range at the QTP. A total of 19 chloroplast haplotypes and 13 nuclear genotypes in two well-differentiated lineages, corresponding to populations into two groups isolated by Tanggula and Bayangela Mountains. Meanwhile, significant phylogeographical structure was detected among sampling range of L. tibetica, and 61.50% of genetic variations was partitioned between groups. Gene flow across the whole region appears to be restricted by high mountains, suggesting a significant role of geography in the genetic differences between the two groups. Divergence time between the two lineages dated to 8.63 million years ago, which corresponded to the uplifting of QTP during the late Miocene and Pliocene. Ecological differences were found between both the lineages represent species-specific characteristics, sufficient to keep the lineages separated to a high degree. The simulated distribution from the last interglacial period to the current period showed that the distribution of L. tibetica experienced shrinkage and expansion. Climate changes during the Pleistocene glacial-interglacial cycles had a dramatic effect on L. tibetica distribution ranges. Multiple refugia of L. tibetica might have remained during the species history, to south of the Tanggula and north of Bayangela Mountains, both appeared as topological barrier and contributed to restricting gene flow between the two lineages. Together, geographic isolation and climatic factors have played a fundamental role in promoting diversification and evolution of L. tibetica.

Nunataks or massif de refuge? A phylogeographic study of Rhodiola crenulata (Crassulaceae) on the world's highest sky islands

Background

Quaternary climatic oscillations had tremendous effects on the current distribution of species. Here, we aim to elucidate the glacial history of Rhodiola crenulata, a perennial herb almost exclusively restricted to rock crevices on mountain peaks, and to test whether the nunatak or massif de refuge hypotheses could explain its distribution pattern.

Results

Six haplotypes and six ribotypes were detected in the cpDNA data set and the ITS data set, respectively. The divergence of R. crenulata and its closest relatives was dated have occurred ca. 0.65 Mya, during the Naynayxungla glaciation on the QTP. Mismatch distribution analysis suggested that the species experienced a range expansion around 0.31 Mya. Populations with high genetic and haplotype diversity were found on the QTP platform as well in the Hengduan Mountains. The ecological niche modeling results showed that there were suitable habitats on both the QTP platform and in the Hengduan Mountains during the LGM.

Conclusion

Our results support a scenario that both nunataks and the massif de refuge hypotheses could explain the distribution of R. crenulata. We also confirmed that Quaternary climatic oscillations could promote plant speciation in some circumstances. This study adds to a growing body of evidence suggesting that the QTP plant lineages exhibited diverse reactions to the Quaternary climatic oscillations.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1270-6) contains supplementary material, which is available to authorized users.

Rapid Intraspecific Diversification of the Alpine Species Saxifraga sinomontana (Saxifragaceae) in the Qinghai-Tibetan Plateau and Himalayas

An increasing number of phylogeographic studies have been conducted for plant species in the Qinghai-Tibetan Plateau (QTP) and its flanking mountains. However, these studies have mainly focused on the determination of glacial refugia and routes of inter-/post-glacial expansions. Rapid intraspecific diversification of plants in this region have not been thoroughly discussed. Herein, we investigate the effects of the Quaternary climate changes on population genetic structure and diversifications of a herbaceous alpine species, Saxifraga sinomontana, which may have an evolutionary time scale <5 million years in the QTP and Himalayan regions. Using a total of 350 individuals from 29 populations, we studied the evolutionary history of S. sinomontana by analyzing cpDNA trnL-trnF, rpl16 and nrDNA ITS sequences. A total of 89 haplotypes and 158 genotypes were detected for cpDNA and ITS sequences, respectively. Only a few haplotypes/genotypes were widespread, while an extremely large number of haplotypes/genotypes were restricted to single populations, which were scattered throughout the current geographical range of S. sinomontana. This suggests the existence of microrefugia of this species during the Quaternary glaciations. In addition, the relationships of the haplotypes/genotypes were almost completely not resolved by phylogenetic reconstruction. Combining characteristics in terms of high haplotype richness, large proportion of private haplotypes, and shallow haplotype divergence, we speculate that recent intraspecific diversification has occurred in S. sinomontana. Molecular clock analysis estimated that the onset diversification within S. sinomontana to be 1.09 Ma (95% HPD = 0.80–1.45), coinciding with the extensive Quaternary glaciations on the QTP which started ca. 1.17 Ma. The Quaternary climatic oscillations may have triggered rapid intraspecific diversification in this QTP-Himalayan species. However, large niche breadth, as well as introgression/hybridization between the studied species and its closely related sympatric saxifrages, may also played a role to some extent on the current genetic structure of S. sinomontana, which need to be further studied.

The effect of Hengduan Mountains Region (HMR) uplift to environmental changes in the HMR and its eastern adjacent area: Tracing the evolutionary history of Allium section Sikkimensia (Amaryllidaceae)

Exploring the effects of orographic events and climatic shifts on geographic distribution of organism in the Hengduan Mountains Region (HMR) and its eastern adjacent area is crucial to the understanding of the environmental changes to organismal evolution. To gain further insight into these processes, we reconstruct evolutionary history of ten species in Allium section Sikkimensia, distributed across regions abovementioned. Using chloroplast and nuclear sequence variation of 79 populations of these ten Allium species with known morphological preferences, we elucidate the phylogenetic relationship, divergence time, ancestral area and genetic structures. Climatic variables analysis, Isolation by distance (IBD) and environment (IBE) and Species distribution modeling (SDM) were analyzed along different genetic clades. These analyses indicated that the initial split of Sikkimensia was triggered by climate changes following Qinghai-Tibet Plateau sensu lato (QTPsl) uplift during the late Miocene. Subsequently, divergences within lineage (lineage A)/among lineages (lineage C and D) in Sikkimensia may be induced by the intense uplift of the HMR around 3-4 Ma and abrupt intensifying of the Asian monsoon regimes. Furthermore, Sikkimensia populations exhibited lopsided demographic history in the Last Glacial Maximum (LGM), as was indicated by the expansion of their range in the QDM and contraction in the HMR. Our findings appear to suggest that the HMR uplift could have strengthened the orographic difference between the HMR and its eastern adjacent area and led to a colder climate in the HMR, while geological topography also played an important role for taxa to respond the climate change that had taken place in the HMR and its eastern adjacent area during the Pleistocene.

Evolutionary history of enigmatic bears in the Tibetan Plateau-Himalaya region and the identity of the yeti

Although anecdotally associated with local bears (Ursus arctos and U. thibetanus), the exact identity of ‘hominid’-like creatures important to folklore and mythology in the Tibetan Plateau–Himalaya region is still surrounded by mystery. Recently, two purported yeti samples from the Himalayas showed genetic affinity with an ancient polar bear, suggesting they may be from previously unrecognized, possibly hybrid, bear species, but this preliminary finding has been under question. We conducted a comprehensive genetic survey of field-collected and museum specimens to explore their identity and ultimately infer the evolutionary history of bears in the region. Phylogenetic analyses of mitochondrial DNA sequences determined clade affinities of the purported yeti samples in this study, strongly supporting the biological basis of the yeti legend to be local, extant bears. Complete mitochondrial genomes were assembled for Himalayan brown bear (U. a. isabellinus) and black bear (U. t. laniger) for the first time. Our results demonstrate that the Himalayan brown bear is one of the first-branching clades within the brown bear lineage, while Tibetan brown bears diverged much later. The estimated times of divergence of the Tibetan Plateau and Himalayan bear lineages overlap with Middle to Late Pleistocene glaciation events, suggesting that extant bears in the region are likely descendants of populations that survived in local refugia during the Pleistocene glaciations.