Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai–Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.     

Different responses of ecosystem carbon exchange to warming in three types of alpine grassland on the central Qinghai–Tibetan Plateau

Climate is a driver of terrestrial ecosystem carbon exchange, which is an important product of ecosystem function. The Qinghai–Tibetan Plateau has recently been subjected to a marked increase in temperature as a consequence of global warming. To explore the effects of warming on carbon exchange in grassland ecosystems, we conducted a whole‐year warming experiment between 2012 and 2014 using open‐top chambers placed in an alpine meadow, an alpine steppe, and a cultivated grassland on the central Qinghai–Tibetan Plateau. We measured the gross primary productivity, net ecosystem CO₂ exchange (NEE), ecosystem respiration, and soil respiration using a chamber‐based method during the growing season. The results show that after 3 years of warming, there was significant stimulation of carbon assimilation and emission in the alpine meadow, but both these processes declined in the alpine steppe and the cultivated grassland. Under warming conditions, the soil water content was more important in stimulating ecosystem carbon exchange in the meadow and cultivated grassland than was soil temperature. In the steppe, the soil temperature was negatively correlated with ecosystem carbon exchange. We found that the ambient soil water content was significantly correlated with the magnitude of warming‐induced change in NEE. Under high soil moisture condition, warming has a significant positive effect on NEE, while it has a negative effect under low soil moisture condition. Our results highlight that the NEE in steppe and cultivated grassland have negative responses to warming; after reclamation, the natural meadow would subject to loose more C in warmer condition. Therefore, under future warmer condition, the overextension of cultivated grassland should be avoided and scientific planning of cultivated grassland should be achieved.     

Occurrence of Christella (Thelypteridaceae) in Southwest China and its indications of the paleoenvironment of the Qinghai–Tibetan Plateau and adjacent areas

The uplift of the Qinghai–Tibetan Plateau dramatically changed the regional topography and climate, profoundly impacting the distribution of many plant lineages. Plant responses to environmental changes are particularly prominent in lineages that require ecological factors differentiated from those present before the uplift of the QTP. Two fossil occurrences of Christella H. Lév., Fl. Kouy–Tchéou (Thelypteridaceae), a fern genus now distributed mainly at low elevations of the pantropics with warm and moist habitats, are described based on fossilized Cenozoic leaf fronds recovered from SW China: late Paleocene Christella nervosa (J. R. Tao) C. L. Xu, T. Su & Z. K. Zhou comb. nov. found in Liuqu, southern Tibet and middle Miocene Christella sp. recovered from the Jinggu Basin in western Yunnan. The frond fossils from both sites share key morphological characteristics that diagnose these fossils as Christella. After detailed comparisons, we further clarified Christella papilio (C. Hope) Holttum, a species distributed in warm, humid habitats at altitudes no more than 1300?m, as the nearest living relative of C. nervosa. This finding suggested that southern Tibet had not reached its present elevation during the late Paleocene (ca. 56 Ma). We propose that the uplift, accompanied by severe cooling and aridification after the late Paleocene, caused the disappearance of Christella in southern Tibet, whereas paleoenvironmental conditions enabled the genus to survive in Yunnan. Our study provides the first example of distributional constraints of ferns in SW China in response to paleoenvironmental changes in the Qinghai–Tibetan Plateau and nearby areas.     

Exploring the community phylogenetic structure along the slope aspect of subalpine meadows in the eastern Qinghai–Tibetan Plateau,China

Exploring the community assembly has been important for explaining the maintenance mechanisms of biodiversity and species coexistence, in that it is a central issue in community ecology. Here, we examined patterns of the community phylogenetic structure of the subalpine meadow plant community along the slope gradient in the Qinghai–Tibetan Plateau of China. We surveyed all species and constructed the phylogenetic tree of the plant community based on data from the Angiosperm Phylogeny Group III. We selected the net relative index (NRI) and evaluated the community phylogenetic structure along the five slope plants communities. We found that the phylogenetic structure varied from phylogenetic clustering to phylogenetic overdispersion with the slope aspect from north to south. In the north slope, the community phylogenetically cluster indicated that the limiting similarity played a leading role in the community assembly and the maintenance of biodiversity. Community phylogenetic overdispersion in the east, southeast, and south slopes indicated that habitat filtration was the driving force for community assembly. The NRI index of the northeast slope was close to zero, implying random dispersion. But it may be driven by the neutral process or limiting similarity, in that the community assembly process was the result of a combination of several ecological factors and thus required further study.     

Leaf and infructescence fossils of Alnus (Betulaceae) from the late Eocene of the southeastern Qinghai–Tibetan Plateau

Plant fossils from the Qinghai–Tibetan Plateau (QTP), China are critical to understand not only the diversification history of plants there, but also the paleoenvironmental conditions. Alnus are deciduous trees, mainly distributed in temperate and subtropical regions of Eurasia and North America, and they are well known in the fossil records throughout the Cenozoic in the Northern Hemisphere. We collected numerous well‐preserved Alnus leaf and infructescence fossils from the Lawula Formation (～34.6 Ma with ⁴⁰Ar/³⁹Ar dating) at the present elevation of 3910 m a.s.l. in the southeastern QTP. Based on detailed morphological comparisons with existing and fossil species, these fossils show closest affinity to Alnus ferdinandi‐coburgii C. K. Schneid., and we refer to these fossils as A. cf. ferdinandi‐coburgii. These specimens comprise the oldest megafossil record of Alnus in the QTP, and provide solid evidence for the distribution of Alnus there as early as the late Eocene. Extant A. ferdinandi‐coburgii is distributed in areas with mean annual temperature values between 9.7 °C and 16.9 °C, and mean annual precipitation values ranging from 896.2 mm to 1161.2 mm; therefore, fossils of A. cf. ferdinandi‐coburgii suggest a much warmer and wetter climate during the late Eocene than today in the southeastern QTP. This finding is consistent with other evidence for continued uplift of the southeastern QTP after the late Eocene that might be due to the eastward extension of the QTP.     

Contrasting responses to Pleistocene climate changes: a case study of two sister species Allium cyathophorum and A. spicata (Amaryllidaceae) distributed in the eastern and western Qinghai–Tibet Plateau

It has been hypothesized that species occurring in the eastern and the western Qinghai–Tibet Plateau (QTP) responded differently to climate changes during the Pleistocene. Here, we test this hypothesis by phylogeographic analysis of two sister species, Allium cyathophorum and A. spicata. We sequenced two chloroplast DNA (cpDNA) fragments (accD‐psaI and the rpl16 intron) of 150 individuals, and the nuclear (ITS) region of 114 individuals, from 19 populations throughout the distributional ranges of these species. The divergence between the two species was dated at 779 ‐ 714 thousand years before the present and was likely initiated by the most major glaciation in the QTP. Analysis of chlorotype diversity showed that A. spicata, the species occurring in the western QTP, contains much lower genetic diversity (0.25) than A. cyathophorum (0.93), which is distributed in the eastern QTP. Moreover, multiple independent tests suggested that the A. spicata population had expanded recently, while no such expansion was detected in A. cyathophorum, indicating a contrasting pattern of responses to Pleistocene climate changes. These findings highlight the importance of geographical topography in determining how species responded to the climate changes that took place in the QTP during the Pleistocene.     

Diploid hybrid origin of Ostryopsis intermedia (Betulaceae) in the Qinghai‐Tibet Plateau triggered by Quaternary climate change

Despite the well‐known effects that Quaternary climate oscillations had on shaping intraspecific diversity, their role in driving homoploid hybrid speciation is less clear. Here, we examine their importance in the putative homoploid hybrid origin and evolution of Ostryopsis intermedia, a diploid species occurring in the Qinghai‐Tibet Plateau (QTP), a biodiversity hotspot. We investigated interspecific relationships between this species and its only other congeners, O. davidiana and O. nobilis, based on four sets of nuclear and chloroplast population genetic data and tested alternative speciation hypotheses. All nuclear data distinguished the three species clearly and supported a close relationship between O. intermedia and the disjunctly distributed O. davidiana. Chloroplast DNA sequence variation identified two tentative lineages, which distinguished O. intermedia from O. davidiana; however, both were present in O. nobilis. Admixture analyses of genetic polymorphisms at 20 SSR loci and sequence variation at 11 nuclear loci and approximate Bayesian computation (ABC) tests supported the hypothesis that O. intermedia originated by homoploid hybrid speciation from O. davidiana and O. nobilis. We further estimated that O. davidiana and O. nobilis diverged 6–11 Ma, while O. intermedia originated 0.5–1.2 Ma when O. davidiana is believed to have migrated southward, contacted and hybridized with O. nobilis possibly during the largest Quaternary glaciation that occurred in this region. Our findings highlight the importance of Quaternary climate change in the QTP in causing hybrid speciation in this important biodiversity hotspot.     

Responses of biomass allocation across two vegetation types to climate fluctuations in the northern Qinghai–Tibet Plateau

The Qinghai–Tibet Plateau (QTP) is particularly sensitive to global climate change, especially to elevated temperatures, when compared with other ecosystems. However, few studies use long‐term field measurements to explore the interannual variations in plant biomass under climate fluctuations. Here, we examine the interannual variations of plant biomass within two vegetation types (alpine meadow and alpine shrub) during 2008–2017 and their relationships with climate variables. The following results were obtained. The aboveground biomass (AGB) and belowground biomass (BGB) response differently to climate fluctuations, the AGB in KPM was dominated by mean annual precipitation (MAP), whereas the AGB in PFS was controlled by mean annual air temperature (MAT). However, the BGB of both KPM and PFS was only weakly affected by climate variables, suggesting that the BGB in alpine ecosystems may remain as a stable carbon stock even under future global climate change. Furthermore, the AGB in PFS was significantly higher than KPM, while the BGB and R/S in KPM were significantly higher than PFS, reflecting the KPM be more likely to allocate more photosynthates to roots. Interestingly, the proportion of 0–10 cm root biomass increased in KPM and PFS, whereas the other proportions both decreased, reflecting a shift in biomass toward the surface layer. Our results could provide a new sight for the prediction how alpine ecosystem response to future climate change.     

Gradual genome size evolution and polyploidy in Allium from the Qinghai–Tibetan Plateau

Background and AimsGenome size is an important plant trait, with substantial interspecies variation. The mechanisms and selective pressures underlying genome size evolution are important topics in evolutionary biology. There is considerable diversity in Allium from the Qinghai–Tibetan Plateau, where genome size variation and related evolutionary mechanisms are poorly understood.MethodsWe reconstructed the Allium phylogeny using DNA sequences from 71 species. We also estimated genome sizes of 62 species, and determined chromosome numbers in 65 species. We examined the phylogenetic signal associated with genome size variation, and tested how well the data fit different evolutionary models. Correlations between genome size variations and seed mass, altitude and 19 bioclimatic factors were determined.Key Results Allium genome sizes differed substantially between species and within diploids, triploids, tetraploids, hexaploids and octaploids. Size per monoploid genome (1Cx) tended to decrease with increasing ploidy levels. Allium polyploids tended to grow at a higher altitude than diploids. The phylogenetic tree was divided into three evolutionary branches. The genomes in Clade I were mostly close to the ancestral genome (18.781 pg) while those in Clades II and III tended to expand and contract, respectively. A weak phylogenetic signal was detected for Allium genome size. Furthermore, significant positive correlations were detected between genome size and seed mass, as well as between genome size and altitude. However, genome size was not correlated with 19 bioclimatic variables.Conclusions Allium genome size shows gradual evolution, followed by subsequent adaptive radiation. The three well-supported Allium clades are consistent with previous studies. The evolutionary patterns in different Allium clades revealed genome contraction, expansion and relative stasis. The Allium species in Clade II may follow adaptive radiation. The genome contraction in Clade III may be due to DNA loss after polyploidization. Allium genome size might be influenced by selective pressure due to the conditions on the Qinghai–Tibetan Plateau (low temperature, high UV irradiation and abundant phosphate in the soil).     

Genetic Variation Within and Among Populations of Rhodiola alsia (Crassulaceae) Native to the Tibetan Plateau as Detected by ISSR Markers

Genetic variation of 10 Rhodiola alsia (Crassulaceae) populations from the Qinghai–Tibet Plateau of China was investigated using intersimple sequence repeat (ISSR) markers. R. alsia is an endemic species of the Qinghai–Tibet Plateau. Of the 100 primers screened, 13 were highly polymorphic. Using these primers, 140 discernible DNA fragments were generated with 112 (80%) being polymorphic, indicating pronounced genetic variation at the species level. Also there were high levels of polymorphism at the population level with the percentage of polymorphic bands (PPB) ranging from 63.4 to 88.6%. Analysis of molecular variance (AMOVA) showed that the genetic variation was mainly found among populations (70.3%) and variance within populations was 29.7%. The main factors responsible for the high level of differentiation among populations are probably the isolation from other populations and clonal propagation of this species. Occasional sexual reproduction might occur in order to maintain high levels of variation within populations. Environmental conditions could also influence population genetic structure as they occur in severe habitats. The strong genetic differentiation among populations in our study indicates that the conservation of genetic variability in R. alsia requires maintenance of as many populations as possible.     

Phylogeography and evolution of a fungal–insect association on the Tibetan Plateau

Parasitoidism refers to a major form of interspecies interactions where parasitoids sterilize and/or kill their hosts typically before hosts reach reproductive age. However, relatively little is known about the evolutionary dynamics of parasitoidism. Here, we investigate the spatial patterns of genetic variation of Chinese cordyceps, including both the parasitoidal fungus Ophiocordyceps sinensis and its host insects. We sampled broadly from alpine regions on the Tibetan Plateau and obtained sequences on seven fungal and three insect DNA fragments from each of the 125 samples. Seven and five divergent lineages/cryptic species were identified within the fungus and host insects, respectively. Our analyses suggested that O. sinensis and host insects originated at similar geographic regions in southern Tibet/Yunnan, followed by range expansion to their current distributions. Cophylogenetic analyses revealed a complex evolutionary relationship between O. sinensis and its host insects. Significant congruence was found between host and parasite phylogenies and the time estimates of divergence were similar, raising the possibility of the occurrence of cospeciation events, but the incongruences suggested that host shifts were also prevalent. Interestingly, one fungal genotype was broadly distributed, consistent with recent gene flow. In contrast, the high‐frequency insect genotypes showed limited geographic distributions. The dominant genotypes from both the fungus and the insect hosts may represent ideal materials from which to develop artificial cultivation of this important Chinese traditional medicine. Our results demonstrate that both historical and contemporary events have played important roles in the phylogeography and evolution of the O. sinensis–ghost moth parasitoidism on the Tibetan Plateau.     

Concentrations of the Active Constituents of the Tibetan Folk Medicine Qinjiao (Gentiana sect. Cruciata) within and between Taxonomic Species across the Qinghai‐Tibetan Plateau

The Tibetan folk medicine Qinjiao is traditionally used to treat various conditions, and its main active constituents comprise four iridoid glycosides, i.e., loganic acid, swertiamarin, gentiopicroside, and sweroside. The traditional crude medicine Qinjiao is derived from the dried roots of three species belonging to Gentiana sect. Cruciata (Gentianaceae) growing in the Qinghai‐Tibetan Plateau (QTP). In this study, we determined by HPLC the contents of the four main active constituents in the dried roots collected from 83 localities at different altitudes across the QTP. The material was classified under the seven taxonomic species G. straminea, G. dahurica, G. crassicaulis, G. waltonii, G. officinalis, G. ihassica, and G. macrophylla. Our results suggested that the four constituents were present in the roots of all seven species for all localities, but their concentrations varied greatly within and between species. The level of gentiopicroside revealed to be the most dominant for all examined localities (2.1–12.4 mg/g), and G. macrophylla Pall . contained the highest concentration of all the four constituents at the species level. Except for loganic acid in G. officinalis, there was no significant correlation between the contents of these constituents and the altitude of the sampling localities. These results suggest that all species of all origins can be used as reliable resource for the crude medicine Qinjiao. However, a few species contain higher concentrations of the main active constituents, irrespective of their origin.     

Mountains too high and valleys too deep drive population structuring and demographics in a Qinghai–Tibetan Plateau frog Nanorana pleskei (Dicroglossidae)

Pleistocene glacial–interglacial climatic oscillations greatly shaped the current genetic structure of many species. However, geographic features may influence the impact of climatic cycling. Distinct geographic and environmental characters between northern and southern parts of the eastern Qinghai–Tibetan Plateau (EQTP) facilitate explorations into the impacts of geographic features on species. The northern parts of EQTP contain large areas of marsh, and the environment is rather homogeneous. In contrast, the southern EQTP harbors complex alpine valleys and a much more heterogeneous setting. We evaluate DNA sequence variation from both the mitochondrial and nuclear genomes in Nanorana pleskei, a species endemic to the EQTP. Hypothesis testing on the evolutionary history of N. pleskei indicates that northern populations can disperse freely, but alpine valleys isolate southern populations. Demographic histories between northern and southern populations also differ. Northern populations appear to have experienced population expansions, while southern frogs exhibit a far more stable demographic history. By combining climatic analyses and species' distribution models, our study suggests that geographic and environmental features drive the differences between the northern and southern EQTP.     

Grazing intensifies degradation of a Tibetan Plateau alpine meadow through plant–pest interaction

Understanding the plant–pest interaction under warming with grazing conditions is critical to predict the response of alpine meadow to future climate change. We investigated the effects of experimental warming and grazing on the interaction between plants and the grassland caterpillar Gynaephora menyuanensis in an alpine meadow on the Tibetan Plateau in 2010 and 2011. Our results showed that grazing significantly increased nitrogen concentration in graminoids and sward openness with a lower sward height, sward coverage, and plant litter mass in the community. Grazing significantly increased G. menyuanensis body size and potential fecundity in 2010. The increases in female body size were about twofold greater than in males. In addition, grazing significantly increased G. menyuanensis density and its negative effects on aboveground biomass and graminoid coverage in 2011. We found that G. menyuanensis body size was significantly positively correlated with nitrogen concentration in graminoids but negatively correlated with plant litter mass. Even though warming did not significantly increased G. menyuanensis performance and the negative effects of G. menyuanensis on alpine meadow, the increases in G. menyuanensis growth rate and its negative effect on aboveground biomass under the warming with grazing treatment were significantly higher than those under the no warming with grazing treatment. The positive effects of grazing on G. menyuanensis performance and its damage were exacerbated by the warming treatment. Our results suggest that the fitness of G. menyuanensis would increase under future warming with grazing conditions, thereby posing a greater risk to alpine meadow and livestock production.