The role of the uplift of the Qinghai‐Tibetan Plateau for the evolution of Tibetan biotas

Decadal shifts in Qingzang Plateau lake carbon dynamics (1970-2020): From predominant carbon sources to emerging sinks

The evasion of carbon dioxide (CO2) from lakes significantly influences the global carbon equilibrium. Amidst global climatic transformations, the role of Qingzang Plateau (QZP) lakes as carbon (C) sources or sinks remains a subject of debate. Furthermore, accurately quantifying their contribution to the global carbon budget presents a formidable challenge. Here, spanning half a century (1970-2020), we utilize a synthesis of literature and empirical field data to assess the CO2 exchange flux of QZP lakes. We find markedly higher CO2 exchange flux in the southeast lakes than that in the northern and western regions from 1970 to 2000. During this time, both freshwater and saltwater lakes served primarily as carbon sources. The annual CO2 exchange flux was estimated at 2.04 ± 0.37 Tg (Tg) C yr-1, mainly influenced by temperature fluctuations. The CO2 exchange flux patterns underwent a geographical inversion between 2000 and 2020, with increased levels in the west and decreased levels in the east. Notably, CO2 emissions from freshwater lakes diminished, and certain saltwater lakes in the QTP transitioned from carbon sources to sinks. From 2000 to 2020, the annual CO2 exchange flux from QZP lakes is estimated at 1.34 ± 0.50 Tg C yr-1, with solar radiation playing a more pronounced role in carbon emissions. Cumulatively, over the past five decades, QZP lakes have generally functioned as carbon sources. Nevertheless, the total annual CO2 emissions have declined since the year 2000, indicating a potential shift trend from being a carbon source to a sink, mirroring broader patterns of global climate change. These findings not only augment our understanding of the carbon cycle in plateau aquatic systems but also provide crucial data for refining China's carbon budget.

Integrative taxonomy, phylogenetics and historical biogeography of subgenus Aeschyntelus Stål, 1872 (Hemiptera: Heteroptera: Rhopalidae)

The subgenus Aeschyntelus includes six species that show variations in body color and shape, thus making it difficult to identify them based on morphological identification alone. To date, no genetic study has evaluated species within this genus. Herein, we collected 171 individuals from 90 localities of Rhopalus and employed an integrative taxonomic approach that incorporated morphological data, mitochondrial genomic data (COI, whole mitochondrial data) and nuclear genomic data (18S + 28S rRNAs, nuclear genome-wide SNPs) to delineate species boundaries. Our analyses confirmed the status of nine described species of Rhopalus and proposed the recognition of one new species known as Rhopalus qinlinganus sp. nov., which is classified within the subgenus Aeschyntelus. Discrepancies arising from nuclear and mitochondrial data suggest the presence of mito-nuclear discordance. Specifically, mitochondrial data indicated admixture within Clade A, comprising R. kerzhneri and R. latus, whereas genome-wide SNPs unambiguously identified two separate species, aligning with morphological classification. Conversely, mitochondrial data clearly distinguished Clade B- consisting of R. sapporensis into two lineages, whereas genome-wide SNPs unequivocally identified a single species. Our study also provides insights into the evolutionary history of Aeschyntelus, thus indicating that it likely originated in East Asia during the middle Miocene. The development of Aeschyntelus biodiversity in the southwestern mountains of China occurred via an uplift-driven diversification process. Our findings highlight the necessity of integrating both morphological and multiple molecular datasets for precise species identification, particularly when delineating closely related species. Additionally, it reveals the important role of mountain orogenesis on speciation within the southwestern mountains of China.

The role of Pleistocene dispersal in shaping species richness of sky island wintergreens from the Himalaya-Hengduan Mountains

In addition to topography and climate, biogeographic dispersal has been considered to influence plant diversity in the Himalaya-Hengduan Mountains (HHM), yet, the mode and tempo of sky island dispersal and its influence on species richness has been little explored. Through phylogenetic analysis of Gaultheria ser. Trichophyllae, a sky island alpine clade within the HHM, we test the hypothesis that dispersal has affected current local species richness. We inferred the dynamics of biogeographic dispersal with correlation tests on direction, distance, occurrence time, and regional species richness. We found that G. ser. Trichophyllae originated at the end of the Miocene and mostly dispersed toward higher longitudes (eastward). In particular, shorter intra-regional eastward dispersals and longer inter-regional westward dispersals were most frequently observed. We detected a prevalence of eastward intra-region dispersals in both glacial periods and interglacials. These dispersals may have been facilitated by the reorganization of paleo-drainages and monsoon intensification through time. We suggest that the timing of dispersal corresponding to glacial periods and the prevalence of intra-region dispersal, rather than dispersal frequency, most influenced the pattern of species richness of G. ser. Trichophyllae. This study facilitates a more comprehensive understanding of biodiversity in the sky islands within the HHM.

Identifying climate refugia for wild yaks (Bos mutus) on the Tibetan Plateau

Climate change is threatening fragile alpine ecosystems and their resident ungulates, particularly the wild yak (Bos mutus) that inhabits alpine areas between the tree line and glaciers on the Tibetan Plateau. Although wild yaks tend to shift habitats in response to changes in climatic factors, the precise impacts of climate change on their habitat distribution and climate refugia remain unclear. Based on over 1000 occurrence records, the maximum entropy (MaxEnt) algorithm was applied to simulate habitat ranges in the last glacial maximum (LGM), Mid-Holocene, current stage, and three greenhouse gas emission scenarios in 2070. Three habitat patches were identified as climate refugia for wild yaks that have persisted from the LGM to the present and are projected to persist until 2070. These stable areas account for approximately 64% of the current wild yak habitat extent and are sufficiently large to support viable populations. The long-term persistence of these climate refugia areas is primarily attributed to the unique alpine environmental features of the Tibetan Plateau, where relatively stable arid or semi-arid climates are maintained, and a wide range of forage resource supplies are available. However, habitat loss by 2070 caused by insufficient protection is predicted to lead to severe fragmentation in the southeastern and northwestern Kunlun, Hengduan, central-western Qilian, and southern Tanggula-northern Himalaya Mountains. Habitat disturbance has also been caused by increasing anthropogenic effects in the southern Tanggula and northern Himalaya Mountains. We suggest that sufficient protection, transboundary cooperation, and community involvement are required to improve wild yak conservation efforts. Our combined modeling method (MaxEnt-Zonation-Linkage Mapper-FRAGSTAT) can be utilized to identify priority areas and linkages between habitat patches while assessing the conservation efficiency of protected areas and analyzing the coupled relationship between climate change and anthropogenic impacts on the habitat distribution of endangered species.

Nocardioides bizhenqiangii sp. nov. and Nocardioides renjunii sp. nov., isolated from soil and faeces of Tibetan antelope (Pantholops hodgsonii) on the Qinghai-Tibet Plateau

Two novel strain pairs (HM61T/HM23 and S-34T/S-58) were isolated from soil and the faeces of Tibetan antelope (Pantholops hodgsonii) collected at the Qinghai-Tibet Plateau of PR China. All four new isolates were aerobic, non-motile, Gram-stain-positive, catalase-positive, oxidase-negative, and short rod-shaped bacteria. The results of phylogenetic analysis based on the full-length 16S rRNA genes and 283 core genomic genes indicated that the four strains were separated into two independent branches belonging to the genus Nocardioides. Strains HM61T and HM23 were most closely related to Nocardioides pelophilus THG T63T (98.58 and 98.65 % 16S rRNA gene sequence similarity). Strains S-34T and S-58 were most closely related to Nocardioides okcheonensis MMS20-HV4-12T (98.89 and 98.89 % 16S rRNA gene sequence similarity). The G+C contents of the genomic DNA of strains HM61T and S-34T were 70.6 and 72.5 mol%, respectively. Strains HM61T, S-34T and the type strains of closely related species in the analysis had average nucleotide identity values of 75.4-90.5 % as well as digital DNA-DNA hybridization values between 20.1 and 40.8 %, which clearly indicated that the four isolates represent two novel species within the genus Nocardioides. The chemotaxonomic characteristics of strains HM61T and S-34T were consistent with the genus Nocardioides. The major fatty acids of all four strains were iso-C16 : 0, C17 : 1  ω8c or C18 : 1  ω9c. For strains HM61T and S-34T, MK-8(H4) was the predominant respiratory quinone, ll-2,6-diaminopimelic acid was the diagnostic diamino acid in the cell-wall peptidoglycan, and the polar lipids profiles were composed of diphosphatidylglycerol and phosphatidylglycerol. Based on phylogenetic, phenotypic, and chemotaxonomic data, we propose that strains HM61T and S-34T represent two novel species of the genus Nocardioides, respectively, with the names Nocardioides bizhenqiangii sp. nov. and Nocardioides renjunii sp. nov. The type strains are HM61T (=GDMCC 4.343T=JCM 36399T) and S-34T (=CGMCC 4.7664T=JCM 33792T).

Complete chloroplast genomes of the hemiparasitic genus Cymbaria: Insights into comparative analysis, development of molecular markers, and phylogenetic relationships

The hemiparasitic tribe Cymbarieae (Orobanchaceae) plays a crucial role in elucidating the initial stage of the transition from autotrophism to heterotrophism. However, the complete chloroplast genome of the type genus Cymbaria has yet to be reported. In addition, the traditional Mongolian medicine Cymbaria daurica is frequently subjected to adulteration or substitution because of the minor morphological differences with Cymbaria mongolica. In this study, the complete chloroplast genomes of the two Cymbaria species were assembled and annotated, and those of other published 52 Orobanchaceae species were retrieved for comparative analyses. We found that the Cymbaria chloroplast genomes are characterized by pseudogenization or loss of stress-relevant genes (ndh) and a unique rbcL-matK inversion. Unlike the high variability observed in holoparasites, Cymbaria and other hemiparasites exhibit high similarity to autotrophs in genome size, guanine-cytosine (GC) content, and intact genes. Notably, four pairs of specific DNA barcodes were developed and validated to distinguish the medicinal herb from its adulterants. Phylogenetic analyses revealed that the genus Cymbaria and the Schwalbea-Siphonostegia clade are grouped into the tribe Cymbarieae, which forms a sister clade to the remaining Orobanchaceae parasitic lineages. Moreover, the diversification of monophyletic Cymbaria occurred during the late Miocene (6.72 Mya) in the Mongol-Chinese steppe region. Our findings provide valuable genetic resources for studying the phylogeny of Orobanchaceae and plant parasitism, and genetic tools to validate the authenticity of the traditional Mongolian medicine "Xinba.".

Morphological observation, molecular identification and evolutionary analysis of Hydatigera kamiyai found in Neodon fuscus from the Qinghai-Tibetan plateau

Hydatigera kamiyai (H. kamiyai) is a new species within Hydatigera that has recently been resurrected. Voles and cats are hosts of H. kamiyai and have a certain impact on its health and economy. Moreover, the Qinghai-Tibetan plateau (QTP) is a research hotspot representing Earth's biodiversity, as its unique geographical environment and climatic conditions support the growth of a variety of mammals and provide favorable conditions for various parasites to complete their life history. The aim of this study was to reveal the phylogenetic relationships and divergence times of H. kamiyai strains isolated from Neodon fuscus on the QTP using morphological and molecular methods. In this study, we morphologically observed H. kamiyai and sequenced the whole mitochondrial genome. Then, we constructed phylogenetic trees with the maximum likelihood (ML) and Bayesian inference (BI) methods. The GTR alternative model was selected for divergence time analysis. These data demonstrated that the results were consistent with the general morphological characteristics of Hydatigera. The whole genome of H. kamiyai was 13,822 bp in size, and the A + T content (73%) was greater than the G + C content (27%). The Ka/Ks values were all <1, indicating that all 13 protein-coding genes (13 PCGs) underwent purifying selection during the process of evolution. The phylogenetic tree generated based on the 13 PCGs, cytochrom oxidase subunit I (COI), 18S rRNA and 28S rRNA revealed close phylogenetic relationships between H. kamiyai and Hydatigera, with high node support for the relationship. The divergence time based on 13 PCGs indicated that H. kamiyai diverged approximately 11.3 million years ago (Mya) in the Miocene. Interestingly, it diverged later than the period of rapid uplift in the QTP. We also speculated that H. kamiyai differentiation was caused by host differentiation due to the favorable living conditions brought about by the uplift of the QTP. As there have been relatively few investigations on the mitochondrial genome of H. kamiyai, our study could provide factual support for further studies of H. kamiyai on the QTP. We also emphasized the importance of further studies of its hosts, Neodon fuscus and cats, which will be important for further understanding the life cycle of H. kamiyai.

Microbial diversity in mountain-dwelling amphibians: The combined effects of host and climatic factors

Comprehending the determinants of host-associated microbiota is pivotal in microbial ecology. Yet, the links between climatic factors and variations in host-associated microbiota necessitate further clarification. Mountain-dwelling amphibians, with limited dispersal abilities, serve as valuable models for addressing these questions. Our study, using 126 amphibian-associated microbial samples (64 gut and 62 skin) and 101 environmental microbial samples (51 soil and 50 water) from the eastern Tibetan Plateau, revealed host factors as primary drivers of the variations in host-associated microbiota. However, climatic factors contributed to additional variations in gut microbial beta-diversity and skin microbial function. Water microbiota were identified as a significant contributor to the amphibian-associated microbiomes, with their climate-driven variations mediating an indirect association between the variations in climatic factors and host-associated microbiota. These findings extend our understanding of the assembly of host-associated microbiota in amphibians, emphasizing the significance of microbiota in evaluating the impact of climate change on animals.

Distinct multitrophic biodiversity composition and community organization in a freshwater lake and a hypersaline lake on the Tibetan Plateau

Alpine lakes play pivotal roles in plateau hydrological processes but are highly sensitive to climate change, yet we lack comprehensive knowledge of their multitrophic biodiversity patterns. Here, we compared the biodiversity characteristics of diverse taxonomic groups across water depths and in surface sediments from a freshwater lake and a hypersaline lake on the northwestern Tibetan Plateau. Using multi-marker environmental DNA metabarcoding, we detected 134 cyanobacteria, 443 diatom, 1,519 invertebrate, and 28 vertebrate taxa. Each group had a substantially different community composition in the two lakes, and differences were also found between water and sediments within each lake. Cooccurrence network analysis revealed higher network complexity, lower modularity, and fewer negative cohesions in the hypersaline lake, suggesting that high salinity may destabilize ecological networks. Our results provide the first holistic view of Tibetan lake biodiversity under contrasting salinity levels and reveal structural differences in the ecological networks that may impact ecosystem resilience.

Landform and lithospheric development contribute to the assembly of mountain floras in China

Although it is well documented that mountains tend to exhibit high biodiversity, how geological processes affect the assemblage of montane floras is a matter of ongoing research. Here, we explore landform-specific differences among montane floras based on a dataset comprising 17,576 angiosperm species representing 140 Chinese mountain floras, which we define as the collection of all angiosperm species growing on a specific mountain. Our results show that igneous bedrock (granitic and karst-granitic landforms) is correlated with higher species richness and phylogenetic overdispersion, while the opposite is true for sedimentary bedrock (karst, Danxia, and desert landforms), which is correlated with phylogenetic clustering. Furthermore, we show that landform type was the primary determinant of the assembly of evolutionarily older species within floras, while climate was a greater determinant for younger species. Our study indicates that landform type not only affects montane species richness, but also contributes to the composition of montane floras. To explain the assembly and differentiation of mountain floras, we propose the 'floristic geo-lithology hypothesis', which highlights the role of bedrock and landform processes in montane floristic assembly and provides insights for future research on speciation, migration, and biodiversity in montane regions.

Historical and dispersal processes drive community assembly of multiple aquatic taxa in glacierized catchments in the Qinghai-Tibet plateau

Understanding the relative role of dispersal dynamics and niche constraints is not only a core task in community ecology, but also becomes an important prerequisite for bioassessment. Despite the recent progress in our knowledge of community assembly in space and time, patterns and processes underlying biotic communities in alpine glacierized catchments remain mostly ignored. To fill this knowledge gap, we combined the recently proposed dispersal-niche continuum index (DNCI) with traditional constrained ordinations and idealized patterns of species distributions to unravel community assembly mechanisms of different key groups of primary producers and consumers (i.e., phytoplankton, epiphytic algae, zooplankton, macroinvertebrates, and fishes) in rivers in the Qinghai-Tibet Plateau, the World's Third Pole. We tested whether organismal groups with contrasting body sizes differed in their assembly processes, and discussed their applicability in bioassessment in alpine zones. We found that community structure of alpine river biotas was always predominantly explained in terms of dispersal dynamics and historical biogeography. These patterns are most likely the result of differences in species-specific functional attributes, the stochastic colonization-extinction dynamics driven by multi-year glacier disturbances and the repeated hydrodynamic separation among alpine catchments after the rising of the Qilian mountains. Additionally, we found that the strength of dispersal dynamics and niche constraints was partially mediated by organismal body sizes, with dispersal processes being more influential for microscopic primary producers. Finding that zooplankton and macroinvertebrate communities followed clumped species replacement structures (i.e., Clementsian gradients) supports the notion that environmental filtering also contributes to the structure of high-altitude animal communities in glacierized catchments. In terms of the applied fields, we argue that freshwater bioassessment in glacierized catchments can benefit from incorporating the metacommunity perspective and applying novel approaches to (i) detect the optimal spatial scale for species sorting and (ii) identify and eliminate the species that are sensitive to dispersal-related processes.

Phylogeographic Structure and Population Dynamics of Baoxing Osmanthus (Osmanthus serrulatus), an Endemic Species from the Southwest Sichuan Basin, China

The mountainous regions of southwest China are recognized as pivotal centers for the origin and evolution of Osmanthus species. Baoxing Osmanthus (Osmanthus serrulatus Rehder), a rare and endemic species known for its spring blooms, is sparsely distributed within the high altitude evergreen broad-leaved forests surrounding the southwestern Sichuan Basin. However, persistent anthropogenic disturbances and habitat fragmentation have precipitated a significant decline in its natural population size, leading to the erosion of genetic resources. To assess the genetic status of O. serrulatus and formulate effective conservation strategies, we conducted sampling across ten wild populations, totaling 148 individuals in their natural habitats. We employed two cpDNA fragments (matK and trnS-trnG) to elucidate the phylogeographic structure and historical population dynamics. The results revealed low species-level genetic diversity, alongside pronounced regional differentiation among populations (FST = 0.812, p < 0.05) and a notable phylogeographic structure (NST = 0.698 > GST = 0.396, p < 0.05). Notably, genetic variation was predominantly observed among populations (81.23%), with no evidence of recent demographic expansion across the O. serrulatus distribution range. Furthermore, divergence dating indicated a timeline of approximately 4.85 Mya, corresponding to the late Miocene to early Pleistocene. This temporal correlation coincided with localized uplift events in the southwestern mountains and heightened Asian monsoons, suggesting pivotal roles for these factors in shaping the current phylogeographic pattern of O. serrulatus. These findings support the effective conservation of O. serrulatus germplasm and offer insights into the impact of Quaternary climate oscillations on companion species within evergreen broad-leaved forests. They also enhance our understanding of the origin and evolution of these forests in the southwestern mountains, aiding biodiversity conservation efforts in the region.

Comparative skin histological and transcriptomic analysis of Rana kukunoris with two different skin colors

This study compares the skin structures of Rana kukunoris with two different skin colors living in the same area of Haibei in the Northeastern Qinghai-Tibet Plateau. The skin thickness of the khaki R. kukunoris was significantly greater than that of the brown R. kukunoris (P < 0.01), and significantly more mucous and granular glands were present on the dorsal skin of the khaki frog (P < 0.05). Meanwhile, the melanocytes on the dorsal skin of the brown frog were significantly larger than those on the khaki one (P < 0.05). Morphological changes in the expansion and aggregation of melanocytes seemed to deepen the skin color of R. kukunoris. Moreover, transcriptome sequencing identified tyrosine metabolism, melanogenesis, and riboflavin metabolism as the main pathways involved in melanin formation and metabolism in brown R. kukunoris. TYR, MC1R was upregulated as the skin color of R. kukunoris was deepened and contributed to melanin production and metabolism. In contrast, the khaki frog had significantly more upregulated genes and metabolic pathways related to autoimmunity. The khaki frog appeared to defend against ultraviolet (UV) radiation-induced damage by secreting mucus and small molecular peptides, whereas the brown frog protected itself by distributing a large amount of melanin. Hence, the different skin colors of R. kukunoris might represent different adaptation strategies for survival in the intense UV radiation environment of the Qinghai-Tibet Plateau.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

The Asian rock-dwelling antlions Gatzara Navás, 1915 and Nepsalus Navás, 1914 (Neuroptera: Myrmeleontidae): new advancements in systematics, biogeography and life history

The antlion genera Gatzara and Nepsalus (Myrmeleontidae: Dendroleontinae) inhabit mountain forests and are characterised by camouflaging larvae. Both genera remain poorly known despite recent findings on systematics and distribution. We report the discovery of new specimens and the previously unknown larvae of the rare species Gatzara jubilaea Navás, 1915, Nepsalus insolitus (Walker, 1860) and N. decorosus (Yang, 1988). These provide new evidence regarding the affinities of these species, and updated knowledge of the distribution, larval morphology and biology. Moreover, a new species of Nepsalus , N. maclachlani Badano, Zheng & Liu, sp. nov. is described from Sri Lanka based on historical museum collections. The discovery of the immature stages of Gatzara shows that the larvae of this genus share the same specialised ecological characteristics and habits as those of Nepsalus but are less morphologically derived. We also reconstruct a molecular phylogeny of this lineage, estimating the divergence time and biogeographical history by adding the new samples. The evolution of the Gatzara + Nepsalus lineage is associated with two major mountain ranges on the southern Tibetan Plateau, i.e. the Himalayas and the Hengduan Mountains. ZooBank: urn:lsid:zoobank.org:pub:68E68211-DFC1-4D98-997B-8A23BA8F9B69.

The Phylogeography of Deciduous Tree Ulmus macrocarpa (Ulmaceae) in Northern China

Disentangling how climate oscillations and geographical events significantly influence plants' genetic architecture and demographic history is a central topic in phytogeography. The deciduous ancient tree species Ulmus macrocarpa is primarily distributed throughout Northern China and has timber and horticultural value. In the current study, we studied the phylogenic architecture and demographical history of U. macrocarpa using chloroplast DNA with ecological niche modeling. The results indicated that the populations' genetic differentiation coefficient (NST) value was significantly greater than the haplotype frequency (GST) (p < 0.05), suggesting that U. macrocarpa had a clear phylogeographical structure. Phylogenetic inference showed that the putative chloroplast haplotypes could be divided into three groups, in which the group Ⅰ was considered to be ancestral. Despite significant genetic differentiation among these groups, gene flow was detected. The common ancestor of all haplotypes was inferred to originate in the middle-late Miocene, followed by the haplotype overwhelming diversification that occurred in the Quaternary. Combined with demography pattern and ecological niche modeling, we speculated that the surrounding areas of Shanxi and Inner Mongolia were potential refugia for U. macrocarpa during the glacial period in Northern China. Our results illuminated the demography pattern of U. macrocarpa and provided clues and references for further population genetics investigations of precious tree species distributed in Northern China.

Across two phylogeographic breaks: Quaternary evolutionary history of a mountain aspen (Populus rotundifolia) in the Hengduan Mountains

Biogeographical barriers to gene flow are central to plant phylogeography. In East Asia, plant distribution is greatly influenced by two phylogeographic breaks, the Mekong-Salween Divide and Tanaka-Kaiyong Line, however, few studies have investigated how these barriers affect the genetic diversity of species that are distributed across both. Here we used 14 microsatellite loci and four chloroplast DNA fragments to examine genetic diversity and distribution patterns of 49 populations of Populus rotundifolia, a species that spans both the Mekong-Salween Divide and the Tanaka-Kaiyong Line in southwestern China. Demographic and migration hypotheses were tested using coalescent-based approaches. Limited historical gene flow was observed between the western and eastern groups of P. rotundifolia, but substantial flow occurred across both the Mekong-Salween Divide and Tanaka-Kaiyong Line, manifesting in clear admixture and high genetic diversity in the central group. Wind-borne pollen and seeds may have facilitated the dispersal of P. rotundifolia following prevalent northwest winds in the spring. We also found that the Hengduan Mountains, where multiple genetic barriers were detected, acted on the whole as a barrier between the western and eastern groups of P. rotundifolia. Ecological niche modeling suggested that P. rotundifolia has undergone range expansion since the last glacial maximum, and demographic reconstruction indicated an earlier population expansion around 600 Ka. The phylogeographic pattern of P. rotundifolia reflects the interplay of biological traits, wind patterns, barriers, niche differentiation, and Quaternary climate history. This study emphasizes the need for multiple lines of evidence in understanding the Quaternary evolution of plants in topographically complex areas.

The complete chloroplast genome and phylogentic results support the species position of Swertia banzragczii and Swertia marginata (Gentianaceae) in Mongolia

BACKGROUND

Swertia banzragczii and S. marginata are important medicinal species in Mongolia. However, their taxonomic positions and genetic backgrounds remain unknown. In this study, we explored the complete chloroplast genomes and DNA barcoding of these species and compared them with those of closely related species within the subgenus to determine their taxonomic positions and phylogenetic relationships.

RESULT

The chloroplast genomes of S. banzragczii and S. marginata encoded 114 genes, including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Among them, 16 genes contained a single intron, and 2 genes had two introns. Closely related species had a conserved genome structure and gene content. Only differences in genome length were noticed, which were caused by the expansion and contraction of the inverted repeat (IR) region and loss of exons in some genes. The trnH-GUG-psbA and trnD-GUC-trnY-GUA intergenic regions had high genetic diversity within Swertia plastomes. Overall, S. banzragczii and S. marginata are true species and belong to the subgenus Swertia.

CONCLUSIONS

These results provide valuable genetic and morphological information on rare and subendemic Swertia species in Mongolia, which can be used for further advanced studies on the Swertia genus.

Successive accumulation of biotic assemblages at a fine spatial scale along glacier-fed waters

Glacier-fed waters create strong environmental filtering for biota, whereby different organisms may assume distinct distribution patterns. By using environmental DNA-based metabarcoding, we investigated the multi-group biodiversity distribution patterns of the Parlung No. 4 Glacier, on the Tibetan Plateau. Altogether, 642 taxa were identified from the meltwater stream and the downstream Ranwu Lake, including 125 cyanobacteria, 316 diatom, 183 invertebrate, and 18 vertebrate taxa. As the distance increased from the glacier terminus, community complexity increased via sequential occurrences of cyanobacteria, diatoms, invertebrates, and vertebrates, as well as increasing taxa numbers. The stream and lake showed different community compositions and distinct taxa. Furthermore, the correlations with environmental factors and community assembly mechanisms showed group- and habitat-specific patterns. Our results reveal the rapid spatial succession and increasing community complexity along glacial flowpaths and highlight the varying adaptivity of different organisms, while also providing insight into the ecosystem responses to global change.

Population Genetic Assessment Model Reveals Conservation Priorities for Gymnocypris Species Resources on the Qinghai-Tibetan Plateau

The Qinghai-Tibetan Plateau (QTP) has nurtured a rich diversity of species because of its unique geographical and environmental conditions. Gymnocypris species (subfamily Schizopygopsinae) are primitive fishes that live in the special environment of the plateau, and their evolution and distribution are inseparable from the historical changes of the QTP. Recently, the resources of Gymnocypris species have been decreasing due to habit deterioration and the intensification of human activities. Therefore, the scientific conservation of the genetic resources of Gymnocypris species is urgently required. In this study, we established two models for the priority conservation assessment of germplasm resources of Gymnocypris species on the basis of the genetic diversity and phylogenetic relationships of 674 individuals from eight Gymnocypris species populations. The results show that the Gymnocypris potanini (GPO), Gymnocypris eckloni (GE), and Gymnocypris przewalskii (GPR) populations are the most genetically diverse in terms of combined genetic diversity values and should be prioritized for conservation. In terms of genetic contribution, the GPO, GE, and GPR populations have a positive impact on maintaining the distinctiveness and diversity of the entire Gymnocypris species population and should be prioritized for conservation. However, in terms of different evolutionary clades, the Gymnocypris namensis, Gymnocypris waddellii, Gymnocypris dobula, and GE populations in clade A should be given priority for protection, the GE population in clade B should be given priority, and the GPR population in clade C should be given priority. In conclusion, the two models and assessment of conservation priorities will provide a scientific basis for the conservation of Gymnocypris species.

Phylogenomic analyses sheds new light on the phylogeny and diversification of Corydalis DC. in Himalaya-Hengduan Mountains and adjacent regions

The Himalaya-Hengduan Mountains (HHM), a renowned biodiversity hotspot of the world, harbors the most extensive habitats for alpine plants with extraordinary high levels of endemism. Although the general evolution pattern has been elucidated, the underlying processes driving spectacular radiations in many species-rich groups remain elusive. Corydalis DC. is widely distributed throughout the Northern Hemisphere containing more than 500 species, with high diversity in HHM and adjacent regions. Using 95 plastid genes, 3,258,640 nuclear single nucleotide polymorphisms (SNPs) and eight single-copy nuclear genes (SCNs) generated from genome skimming data, we reconstructed a robust time-calibrated phylogeny of Corydalis comprising more than 100 species that represented all subgenera and most sections. Molecular dating indicated that all main clades of Corydalis began to diverge in the Eocene, with the majority of extant species in HHM emerged from a diversification burst after the middle Miocene. Global pattern of mean divergence times indicated that species distributed in HHM were considerably younger than those in other regions, particularly for the two most species-rich clades (V and VI) of Corydalis. The early divergence and the recent diversification of Corydalis were most likely promoted by the continuous orogenesis and climate change associated with the uplift of the Qinghai-Tibetan Plateau (QTP). Our study demonstrates the effectivity of phylogenomic analyses with genome skimming data on the phylogeny of species-rich taxa, and sheds lights on how the uplift of QTP has triggered the evolutionary radiations of large plant genera in HHM and adjacent regions.

Genetic differentiation and phylogeography of Erythroneurini (Hemiptera, Cicadellidae, Typhlocybinae) in the southwestern karst area of China

Erythroneurini is the largest tribe of the microleafhopper subfamily Typhlocybinae. Most prior research on this tribe has focused on traditional classification, phylogeny, and control of agricultural pests, and the phylogeography of the group remains poorly understood. In this study, the mitochondrial genomes of 10 erythroneurine species were sequenced, and sequences of four genes were obtained for 12 geographical populations of Seriana bacilla. The new sequence data were combined with previously available mitochondrial DNA sequence data and analyzed using Bayesian and Maximum-Likelihood-based phylogenetic methods to elucidate relationships among genera and species and estimate divergence times. Seriana was shown to be derived from within Empoascanara. Phylogeographic and population genetic analysis of the endemic Chinese species Seriana bacilla suggest that the species diverged about 54.85 Mya (95% HPD: 20.76-66.23 million years) in the Paleogene period and that population divergence occurred within the last 14 million years. Ancestral area reconstruction indicates that Seriana bacilla may have originated in the central region of Guizhou, and geographical barriers are the main factors affecting gene flow among populations. Ecological niche modeling using the MaxEnt model suggests that the distribution of the species was more restricted in the past but is likely to expand in the future years 2050 and 2070.

Cryptic diversity and rampant hybridization in annual gentians on the Qinghai-Tibet Plateau revealed by population genomic analysis

Understanding the evolutionary and ecological processes involved in population differentiation and speciation provides critical insights into biodiversity formation. In this study, we employed 29,865 single nucleotide polymorphisms (SNPs) and complete plastomes to examine genomic divergence and hybridization in Gentiana aristata, which is endemic to the Qinghai-Tibet Plateau (QTP) region. Genetic clustering revealed that G. aristata is characterized by geographic genetic structures with five clusters (West, East, Central, South and North). The West cluster has a specific morphological character (i.e., blue corolla) and higher values of FST compared to the remaining clusters, likely the result of the geological barrier formed by the Yangtze River. The West cluster diverged from the other clusters in the Early Pliocene; these remaining clusters diverged from one another in the Early Quaternary. Phylogenetic reconstructions based on SNPs and plastid data revealed substantial cyto-nuclear conflicts. Genetic clustering and D-statistics demonstrated rampant hybridization between the Central and North clusters, along the Bayankala Mountains, which form the geological barrier between the Central and North clusters. Species distribution modeling demonstrated the range of G. aristata expanded since the Last Interglacial period. Our findings provide genetic and morphological evidence of cryptic diversity in G. aristata, and identified rampant hybridization between genetic clusters along a geological barrier. These findings suggest that geological barriers and climatic fluctuations have an important role in triggering diversification as well as hybridization, indicating that cryptic diversity and hybridization are essential factors in biodiversity formation within the QTP region.

Historical biogeography and evolutionary diversification of Lilium (Liliaceae): New insights from plastome phylogenomics

Here, we infer the historical biogeography and evolutionary diversification of the genus Lilium. For this purpose, we used the complete plastomes of 64 currently accepted species in the genus Lilium (14 plastomes were newly sequenced) to recover the phylogenetic backbone of the genus and a time-calibrated phylogenetic framework to estimate biogeographical history scenarios and evolutionary diversification rates of Lilium. Our results suggest that ancient climatic changes and geological tectonic activities jointly shaped the distribution range and drove evolutionary radiation of Lilium, including the Middle Miocene Climate Optimum (MMCO), the late Miocene global cooling, as well as the successive uplift of the Qinghai-Tibet Plateau (QTP) and the strengthening of the monsoon climate in East Asia during the late Miocene and the Pliocene. This case study suggests that the unique geological and climatic events in the Neogene of East Asia, in particular the uplift of QTP and the enhancement of monsoonal climate, may have played an essential role in formation of uneven distribution of plant diversity in the Northern Hemisphere.

Landscape Heterogeneity Explains the Genetic Differentiation of a Forest Bird across the Sino-Himalayan Mountains

Mountains are the world's most important centers of biodiversity. The Sino-Himalayan Mountains are global biodiversity hotspot due to their extremely high species richness and endemicity. Ample research investigated the impact of the Qinghai-Tibet Plateau uplift and Quaternary glaciations in driving species diversification in plants and animals across the Sino-Himalayan Mountains. However, little is known about the role of landscape heterogeneity and other environmental features in driving diversification in this region. We utilized whole genomes and phenotypic data in combination with landscape genetic approaches to investigate population structure, demography, and genetic diversity in a forest songbird species native to the Sino-Himalayan Mountains, the red-billed leiothrix (Leiothrix lutea). We identified 5 phylogeographic clades, including 1 in the East of China, 1 in Yunnan, and 3 in Tibet, roughly consistent with differences in song and plumage coloration but incongruent with traditional subspecies boundaries. Isolation-by-resistance model best explained population differentiation within L. lutea, with extensive secondary contact after allopatric isolation leading to admixture among clades. Ecological niche modeling indicated relative stability in the extent of suitable distribution areas of the species across Quaternary glacial cycles. Our results underscore the importance of mountains in the diversification of this species, given that most of the distinct genetic clades are concentrated in a relatively small area in the Sino-Himalayan Mountain region, while a single shallow clade populates vast lower-lying areas to the east. This study highlights the crucial role of landscape heterogeneity in promoting differentiation and provides a deep genomic perspective on the mechanisms through which diversity hotspots form.

Disentangling the effects of geographic distance, environment and history on beta diversity of freshwater fish at a biogeographical crossroads

Examining assemblage turnover and variation along geographic and environmental distances is a useful approach to evaluate beta diversity patterns and associated driving mechanisms. However, such studies are relatively limited in freshwater systems. Here, we compared the relationships between freshwater fish beta diversity and geographic distances among 165 hydrological units (HUs) in four zoogeographical regions (PA, Palearctic Region; CA, High Central Asia; EA, East Asia, SA, South Asia) across China and adjacent areas. This area can be considered a biogeographical crossroads, where faunal composition shares elements with different biogeographic and evolutionary origins. We found a considerably high level of between-HU overall dissimilarity (βsor, range from ca. 0.60 to 0.85) in all four regions, mainly due to the turnover component (the relative contribution of βsim to βsor ranged from 60% to 90%). In general, βsor and βsim both significantly increased with geographic distance (except in PA), whereas the nestedness-resultant component (βsne) decreased with geographic distance. The intercepts and slopes of the relationships between dissimilarities and distance (RDDs) both varied significantly among the four regions. The intercepts of βsor and βsim were both highest in SA, followed by CA, PA and EA, implying different levels of fish faunal heterogeneity at short distances. In contrast, the slopes of these two dissimilarities followed the decreasing trend from EA > CA > SA > PA, suggesting different environmental suitability and dispersal ability of fish species among regions. Variation partitioning in distance-based redundancy analysis showed that the spatial and historical factors were more important than area-heterogeneity and energy factors across all HUs and within three individual ecoregions (EA, SA and CA), but spatial factors were non-significant in PA. Our study highlighted the usefulness of RDDs in understanding biogeographical patterns and enhancing the biodiversity conservation of freshwater fishes.

Comparative Analysis of Chloroplast Genome of Meconopsis (Papaveraceae) Provides Insights into Their Genomic Evolution and Adaptation to High Elevation

The Meconopsis species are widely distributed in the Qinghai-Tibet Plateau, Himalayas, and Hengduan Mountains in China, and have high medicinal and ornamental value. The high diversity of plant morphology in this genus poses significant challenges for species identification, given their propensity for highland dwelling, which makes it a question worth exploring how they cope with the harsh surroundings. In this study, we recently generated chloroplast (cp) genomes of two Meconopsis species, Meconopsis paniculata (M. paniculata) and M. pinnatifolia, and compared them with those of ten Meconopsis cp genomes to comprehend cp genomic features, their phylogenetic relationships, and what part they might play in plateau adaptation. These cp genomes shared a great deal of similarities in terms of genome size, structure, gene content, GC content, and codon usage patterns. The cp genomes were between 151,864 bp and 154,997 bp in length, and contain 133 predictive genes. Through sequence divergence analysis, we identified three highly variable regions (trnD-psbD, ccsA-ndhD, and ycf1 genes), which could be used as potential markers or DNA barcodes for phylogenetic analysis. Between 22 and 38 SSRs and some long repeat sequences were identified from 12 Meconopsis species. Our phylogenetic analysis confirmed that 12 species of Meconopsis clustered into a monophyletic clade in Papaveraceae, which corroborated their intrageneric relationships. The results indicated that M. pinnatifolia and M. paniculata are sister species in the phylogenetic tree. In addition, the atpA and ycf2 genes were positively selected in high-altitude species. The functions of these two genes might be involved in adaptation to the extreme environment in the cold and low CO2 concentration conditions at the plateau.

Predicting the potential distribution of 12 threatened medicinal plants on the Qinghai-Tibet Plateau, with a maximum entropy model

Climate change is a vital driver of biodiversity patterns and species distributions, understanding how organisms respond to climate change will shed light on the conservation of endangered species. In this study, the MaxEnt model was used to predict the potential suitable area of 12 threatened medicinal plants in the QTP (Qinghai-Tibet Plateau) under the current and future (2050s, 2070s) three climate scenarios (RCP2.6, RCP4.5, RCP8.5). The results showed that the climatically suitable habitats for the threatened medicinal plants were primarily found in the eastern, southeast, southern, and some parts of the central regions on the QTP. Moreover, 25% of the threatened medicinal plants would have reduced suitable habitat areas within the next 30-50 years in the different future global warming scenarios. Among these medicinal plants, RT (Rheum tanguticum) would miss the most habitat (98.97%), while the RAN (Rhododendron anthopogonoides) would miss the least habitat (10.15%). Nevertheless, 33.3% of the threatened medicinal plants showed an increase in their future habitat area because of their physiological characteristics which are more adaptable to a wide range of climates. The climatic suitable habitat for 50% of the threatened medicinal plants would migrate to higher altitudes or higher latitudes regions. This study provides a data foundation for the conservation of biodiversity and wild medicinal plants on the QTP.

Drivers of desert plant beta diversity on the Qinghai-Tibet plateau

The desert ecosystem of the Qinghai-Tibet Plateau (QTP) is an important component of China's desert ecosystem. Studying the mechanisms shaping the taxonomic, phylogenetic, and functional beta diversity of plant communities in the QTP desert will help us to promote scientific conservation and management of the region's biodiversity. This study investigated the effects of environmental (including altitude, climate factors, and soil factors) and geographic distances on three facets of beta diversity as well as their turnover and nestedness components based on field survey data. The results showed that turnover components dominate the three facets of beta diversity. However, the turnover contributions to phylogenetic and functional beta diversity were lower than for taxonomic beta diversity. Environmental distance had a greater influence than geographic distance, with the former uniquely explaining 15.2%-22.8% of beta diversity and the latter explaining only 1.7%-2.4%. Additionally, the explanatory power of different factors for beta diversity differed between herbs and shrubs, with environmental distance being more important for the latter. Distance-based redundancy analysis suggested that soil total potassium content had a substantial impact on the beta diversity of three dimensions, with mean temperature of the coldest month and soil total phosphorus content having a substantial impact on taxonomic and functional beta diversity as well. Our results support that environmental sorting plays a predominant role in shaping plant community composition across QTP desert ecosystems. To maintain the plant diversity of this region, it is crucial to prioritize the conservation of its diverse environmental conditions and actively mitigate its degradation by anthropogenic pressures.

The influence of climate oscillations and geological events on population differentiation of Camponotus japonicus in the Chinese mainland

Camponotus japonicus (Hymenoptera: Formicidae) is an omnivorous social insect which builds sizable colonies in sparse woodlands or cropland and spreads across multiple climatic zones in the Chinese mainland. This study aims to reveal the role of climate changes and geological events in driving the genetic structure of social insect populations by investigating the phylogenetics and historical demography of C. japonicus in the Chinese mainland. Phylogenetic analyses were conducted based on the mitochondria DNA dataset using MrBayes and IQ-TREE. We constructed a haplotype network, calculated analyses of molecular variance, estimated the divergence time, and reconstructed the maximum clade credibility tree. Mismatch distribution and Bayesian skyline plots were used to infer historical population fluctuations. Additionally, ecological niche modeling was employed to predict the potential distribution of the species during the present, mid-holocene, and last glacial maximum periods in the Chinese mainland. The phylogenetic tree and median-joining network analyses support the presence of four distinct lineages in C. japonicus. These lineages exhibit significant genetic differentiation and limited gene flow. The divergence among the four lineages began in the early Pleistocene, approximately 1.41 million years ago (Ma). Subsequently, the central lineage diverged from both the northern and southern lineages around 1.16 Ma, while the northern and southern lineages diverged from each other at approximately 1.07 Ma. Population expansion was observed in the southern, central, and northern lineages prior to the last glacial maximum, while the Yunnan-Sichuan lineage experienced a slight increase in population size in more recent times. The predicted distribution of the species corresponds well with the actual distribution. Furthermore, the current suitable habitat areas in northern Xinjiang, southern Tibet, and the southeast coastal regions have significantly decreased compared to the last glacial maximum and the mid-holocene periods. Our results suggest that climate oscillations and geological events play an important role in driving genetic patterns and differentiation of C. japonicus. Mountain barriers isolate populations from each other, hinder the flow of genes, and effectively prevent the spread of this species. But at the same time, it also formed refugia at low altitudes areas such as Qinling-Bashan Mountains and Yanshan-Taihang Mountains and provide suitable habitats during glaciation. This study provides a good model for understanding how complex climate changes and geological events affect population genetic differentiation of social insects in the Chinese mainland.

Diversification of freshwater crabs on the sky islands in the Hengduan Mountains Region, China

The numerous naturally-fragmented sky islands (SIs) in the Hengduan Mountains Region (HMR) of southwestern China constitute discontinuous landscapes where montane habitats are isolated by dry-hot valleys which have fostered exceptional species diversification and endemicity. However, studies documenting the crucial role of SI on the speciation dynamics of native freshwater organisms are scarce. Here we used a novel set of comprehensive genetic markers (24 nuclear DNA sequences and complete mitogenomes), morphological characters, and biogeographical information to reveal the evolutionary history and speciation mechanisms of a group of small-bodied montane potamids in the genus Tenuipotamon. Our results provide a robustly supported phylogeny, and suggest that the vicariance events of these montane crabs correlate well with the emergence of SIs due to the uplift of the HMR during the Late Oligocene. Furthermore, ancestrally, mountain ridges provided corridors for the dispersal of these montane crabs that led to the colonization of moist montane-specific habitats, aided by past climatic conditions that were the crucial determinants of their evolutionary history. The present results illustrated that the mechanisms isolating SIs are reinforced by the harsh-dry isolating climatic features of dry-hot valleys separating SIs and continue to affect local diversification. This offers insights into the causes of the high biodiversity and endemism shown by the freshwater crabs of the HMR-SIs in southwestern China.

Genome-wide phylogenetic and genetic evolutionary analyses of mitochondria in Hypoderma bovis and H. sinense on the Qinghai-Tibetan Plateau

Hypoderma bovis (H. bovis) and Hypoderma sinense (H. sinense) are insects that cause hypodermosis in yaks and Bos taurus. Hypodermosis is a severe skin condition that not only impairs the development of local animal husbandry but also poses threats to human health as a zoonosis. The Qinghai-Tibetan Plateau (QTP) is known as the "Roof of the World." Its unique geographical environment and climate conditions have supported the growth of a wide range of mammals, providing favorable conditions for Hypoderma spp. to complete their life cycles. In this study, the whole mitochondrial genomes of H. bovis and H. sinense collected from the QTP were sequenced and phylogenetically analyzed. We found that the whole genomes of H. bovis and H. sinense are 16,283 bp and 16,300 bp in length, respectively. Both the H. bovis and H. sinense genomes have 37 mitochondrial genes, which include two rRNA genes (16S rRNA and 12S rRNA), 22 tRNA genes, the control region (D-loop region), the light chain replication initiation region, and 13 protein-coding genes (PCGs). The phylogenetic tree generated based on the 13 PCGs revealed close phylogenetic relationships between H. sinense, H. bovis, and Hypoderma lineatum. A similar result was also found in our phylogenetic analysis based on 18S rRNA and 28S rRNA. However, analysis of cytochrome oxidase subunit I (COI) showed cluster of H. bovis, H. sinense, and Cuterebra spp. on the same branch, all belonging to Oestridae. The differentiation time generated based on 13 PCGs indicates that H. bovis and H. sinense differentiated and formed ~4.69 million years ago (Mya) and ~4.06 Mya, respectively. This timing coincides with the differentiation and appearance of yak and Bos taurus in the Pliocene (~4.7 Mya), indicating that the parasites and mammals diverged in close temporal proximity. Of note, this period also witnessed a rapid uplift of the QTP, causing significant climate and environmental changes. Thus, we conjecture that the differentiation of Hypoderma spp. is potentially related to the differentiation of their host species, as well as climate changes caused by the uplift of the QTP. Overall, our study can provide valuable data to support further studies on the phylogeny and differentiation of Hypoderma spp. on the QTP.

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

The emergence of modern zoogeographic regions in Asia examined through climate-dental trait association patterns

The complex and contrasted distribution of terrestrial biota in Asia has been linked to active tectonics and dramatic climatic changes during the Neogene. However, the timings of the emergence of these distributional patterns and the underlying climatic and tectonic mechanisms remain disputed. Here, we apply a computational data analysis technique, called redescription mining, to track these spatiotemporal phenomena by studying the associations between the prevailing herbivore dental traits of mammalian communities and climatic conditions during the Neogene. Our results indicate that the modern latitudinal zoogeographic division emerged after the Middle Miocene climatic transition, and that the modern monsoonal zoogeographic pattern emerged during the late Late Miocene. Furthermore, the presence of a montane forest biodiversity hotspot in the Hengduan Mountains alongside Alpine fauna on the Tibetan Plateau suggests that the modern distribution patterns may have already existed since the Pliocene.

Past climate cooling and orogenesis of the Hengduan Mountains have influenced the evolution of Impatiens sect. Impatiens (Balsaminaceae) in the Northern Hemisphere

BACKGROUND

Impatiens sect. Impatiens is distributed across the Northern Hemisphere and has diversified considerably, particularly within the Hengduan Mountains (HDM) in southwest China. Yet, the infra-sectional phylogenetic relationships are not well resolved, largely due to limited taxon sampling and an insufficient number of molecular markers. The evolutionary history of its diversification is also poorly understood. In this study, plastome data and the most complete sampling to date were used to reconstruct a robust phylogenetic framework for this section. The phylogeny was then used to investigate its biogeographical history and diversification patterns, specifically with the aim of understanding the role played by the HDM and past climatic changes in its diversification.

RESULTS

A stable phylogeny was reconstructed that strongly supported both the monophyly of the section and its division into seven major clades (Clades I-VII). Molecular dating and ancestral area reconstruction suggest that sect. Impatiens originated in the HDM and Southeast China around 11.76 Ma, after which different lineages dispersed to Northwest China, temperate Eurasia, and North America, mainly during the Pliocene and Pleistocene. An intercontinental dispersal event from East Asia to western North America may have occurred via the Bering Land Bridge or Aleutian Islands. The diversification rate was high during its early history, especially with the HDM, but gradually decreased over time both within and outside the HDM. Multiple linear regression analysis showed that the distribution pattern of species richness was strongly associated with elevation range, elevation, and mean annual temperature. Finally, ancestral niche analysis indicated that sect. Impatiens originated in a relatively cool, middle-elevation area.

CONCLUSIONS

We inferred the evolutionary history of sect. Impatiens based on a solid phylogenetic framework. The HDM was the primary source or pump of its diversity in the Northern Hemisphere. Orogeny and climate change may have also shaped its diversification rates, as a steady decrease in the diversification rate coincided with the uplift of the HDM and climate cooling. These findings provide insights into the distribution pattern of sect. Impatiens and other plants in the Northern Hemisphere.

Elucidation of the evolutionary history of Stipa in China using comparative transcriptomic analysis

Phylogenetic analysis provides crucial insights into the evolutionary relationships and diversification patterns within specific taxonomic groups. In this study, we aimed to identify the phylogenetic relationships and explore the evolutionary history of Stipa using transcriptomic data. Samples of 12 Stipa species were collected from the Qinghai-Tibet Plateau and Mongolian Plateau, where they are widely distributed, and transcriptome sequencing was performed using their fresh spikelet tissues. Using bidirectional best BLAST analysis, we identified two sets of one-to-one orthologous genes shared between Brachypodium distachyon and the 12 Stipa species (9397 and 2300 sequences, respectively), as well as 62 single-copy orthologous genes. Concatenation methods were used to construct a robust phylogenetic tree for Stipa, and molecular dating was used to estimate divergence times. Our results indicated that Stipa originated during the Pliocene. In approximately 0.8 million years, it diverged into two major clades each consisting of native species from the Mongolian Plateau and the Qinghai-Tibet Plateau, respectively. The evolution of Stipa was closely associated with the development of northern grassland landscapes. Important external factors such as global cooling during the Pleistocene, changes in monsoonal circulation, and tectonic movements contributed to the diversification of Stipa. This study provided a highly supported phylogenetic framework for understanding the evolution of the Stipa genus in China and insights into its diversification patterns.

Ancient allopatry and ecological divergence act together to promote plant diversity in mountainous regions: evidence from comparative phylogeography of two genera in the Sino-Himalayan region

BACKGROUND

How geographical isolation and ecological divergence act together to promote plant diversity in mountainous regions remains largely unknown. In this study, we chose two genera comprising a small number of species distributed in the Sino-Himalayan region, Megacodon (Gentianaceae) and Beesia (Ranunculaceae), which both exhibit a fragmented distribution pattern and are found across a wide range of elevations. By summarizing their common patterns of speciation and/or divergence processes, we aim to understand how environmental changes accelerated lineage diversification in the Sino-Himalayan region through ancient allopatry and ecological divergence.

RESULTS

Using ddRAD-seq, chloroplast genome sequences, and specific molecular markers, we studied the phylogenetic relationships, population structure, and historical biogeography of Beesia and Megacodon. Both genera began to diverge from the late Miocene onwards, with ancient allopatry at lower elevations formed narrow-range species or relict populations. Mantel tests between genetic distance and climatic, elevational, or geographic distance revealed an isolation-by-distance pattern in Beesia and Megacodon stylophorus. Megacodon showed two clades occupying entirely different altitudinal ranges, whereas Beesia calthifolia exhibited a genetic divergence pattern along an elevation gradient. Furthermore, we conducted morphological measurements on Beesia calthifolia and found that different elevational groups had distinct leaf shapes.

CONCLUSIONS

The regional disjunctions of plant groups in the Sino-Himalayan region are drastic and closely related to several biogeographic boundaries. As a consequence of major geological and climate change, ecological divergence when different elevations are colonized often happens simultaneously within plant groups. Although habitat fragmentation and parapatric ecological divergence each spur speciation to different extents, a combined effect of these two factors is a common phenomenon in the Sino-Himalayan region.

Comparative metabolomics analysis reveals high-altitude adaptations in a toad-headed viviparous lizard, Phrynocephalus vlangalii

Extreme environmental conditions at high altitude, such as hypobaric hypoxia, low temperature, and strong UV radiation, pose a great challenge to the survival of animals. Although the mechanisms of adaptation to high-altitude environments have attracted much attention for native plateau species, the underlying metabolic regulation remains unclear. Here, we used a multi-platform metabolomic analysis to compare metabolic profiles of liver between high- and low-altitude populations of toad-headed lizards, Phrynocephalus vlangalii, from the Qinghai-Tibet Plateau. A total of 191 differential metabolites were identified, consisting of 108 up-regulated and 83 down-regulated metabolites in high-altitude lizards as compared with values for low-altitude lizards. Pathway analysis revealed that the significantly different metabolites were associated with carbohydrate metabolism, amino acid metabolism, purine metabolism, and glycerolipid metabolism. Most intermediary metabolites of glycolysis and the tricarboxylic acid cycle were not significantly altered between the two altitudes, but most free fatty acids as well as β-hydroxybutyric acid were significantly lower in the high-altitude population. This may suggest that high-altitude lizards rely more on carbohydrates as their main energy fuel rather than lipids. Higher levels of phospholipids occurred in the liver of high-altitude populations, suggesting that membrane lipids may undergo adaptive remodeling in response to low-temperature stress at high altitude. In summary, this study demonstrates that metabolic profiles differ substantially between high- and low-altitude lizard populations, and that these differential metabolites and metabolic pathways can provide new insights to reveal mechanisms of adaptation to extreme environments at high altitude.

Phytodiversity is associated with habitat heterogeneity from Eurasia to the Hengduan Mountains

The geographic distribution of plant diversity matches the gradient of habitat heterogeneity from lowlands to mountain regions. However, little is known about how much this relationship is conserved across scales. Using the World Checklist of Vascular Plants and high-resolution biodiversity maps developed by species distribution models, we investigated the associations between species richness and habitat heterogeneity at the scales of Eurasia and the Hengduan Mountains (HDM) in China. Habitat heterogeneity explains seed plant species richness across Eurasia, but the plant species richness of 41/97 HDM families is even higher than expected from fitted statistical relationships. A habitat heterogeneity index combining growing degree days, site water balance, and bedrock type performs better than heterogeneity based on single variables in explaining species richness. In the HDM, the association between heterogeneity and species richness is stronger at larger scales. Our findings suggest that high environmental heterogeneity provides suitable conditions for the diversification of lineages in the HDM. Nevertheless, habitat heterogeneity alone cannot fully explain the distribution of species richness in the HDM, especially in the western HDM, and complementary mechanisms, such as the complex geological history of the region, may have contributed to shaping this exceptional biodiversity hotspot.

The Sino-Himalayan flora evolved from lowland biomes dominated by tropical floristic elements

BACKGROUND

The Sino-Himalayan flora harbors highly diverse high-elevation biotas, but our understanding of its evolutionary history in temporal and spatial dimensions is limited. In this study, we integrated a dated phylogenetic tree with comprehensive species distribution data to investigate changes over time and space in floristic elements, including the tropical, Tethys, northern temperate, and East Asian floristic elements, across the entire Sino-Himalaya and its three floristic regions: the Yunnan Plateau, Hengduan Mountains, and East Himalaya regions.

RESULTS

Our results revealed that the Sino-Himalayan flora developed from lowland biomes and was predominantly characterized by tropical floristic elements before the collision between the Indian subcontinent and Eurasia during the Early Cenozoic. Subsequently, from the Late Eocene onwards, the uplifts of the Himalaya and Hengduan Mountains transformed the Sino-Himalayan region into a wet and cold plateau, on which harsh and diverse ecological conditions forced the rapid evolution of local angiosperms, giving birth to characteristic taxa adapted to the high altitudes and cold habitat. The percentage of temperate floristic elements increased and exceeded that of tropical floristic elements by the Late Miocene.

CONCLUSIONS

The Sino-Himalayan flora underwent four significant formation periods and experienced a considerable increase in endemic genera and species in the Miocene, which remain crucial to the present-day patterns of plant diversity. Our findings support the view that the Sino-Himalayan flora is relatively young but has ancient origins. The three major shifts in the divergence of genera and species during the four formation periods were primarily influenced by the uplifts of the Himalaya and Hengduan Mountains and the onset and intensification of the Asian monsoon system. Additionally, the temporal patterns of floristic elements differed among the three floristic regions of the Sino-Himalaya, indicating that the uplift of the Himalaya and surrounding areas was asynchronous. Compared to the Yunnan Plateau region, the East Himalaya and Hengduan Mountains experienced more recent and drastic uplifts, resulting in highly intricate topography with diverse habitats that promoted the rapid radiation of endemic genera and species in these regions.

Evolution of Cherries (Prunus Subgenus Cerasus) Based on Chloroplast Genomes

Cherries (Prunus Subgenus Cerasus) have economic value and ecological significance, yet their phylogeny, geographic origin, timing, and dispersal patterns remain challenging to understand. To fill this gap, we conducted a comprehensive analysis of the complete chloroplast genomes of 54 subg. Cerasus individuals, along with 36 additional genomes from the NCBI database, resulting in a total of 90 genomes for comparative analysis. The chloroplast genomes of subg. Cerasus exhibited varying sizes and consisted of 129 genes, including protein-coding, transfer RNA, and ribosomIal RNA genes. Genomic variation was investigated through InDels and SNPs, showcasing distribution patterns and impact levels. A comparative analysis of chloroplast genome boundaries highlighted variations in inverted repeat (IR) regions among Cerasus and other Prunus species. Phylogeny based on whole-chloroplast genome sequences supported the division of Prunus into three subgenera, I subg. Padus, II subg. Prunus and III subg. Cerasus. The subg. Cerasus was subdivided into seven lineages (IIIa to IIIg), which matched roughly to taxonomic sections. The subg. Padus first diverged 51.42 Mya, followed by the separation of subg. Cerasus from subg. Prunus 39.27 Mya. The subg. Cerasus started diversification at 15.01 Mya, coinciding with geological and climatic changes, including the uplift of the Qinghai-Tibet Plateau and global cooling. The Himalayans were the refuge of cherries, from which a few species reached Europe through westward migration and another species reached North America through northeastward migration. The mainstage of cherry evolution was on the Qing-Tibet Plateau and later East China and Japan as well. These findings strengthen our understanding of the evolution of cherry and provide valuable insights into the conservation and sustainable utilization of cherry's genetic resources.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants

Introduction

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification.

Methods

In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized.

Results

exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago.

Discussion

The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

The Himalayan uplift and evolution of aquatic biodiversity across Asia: Snowtrout (Cyprininae: Schizothorax) as a test case

Global biodiversity hotspots are often remote, tectonically active areas undergoing climatic fluctuations, such as the Himalaya Mountains and neighboring Qinghai-Tibetan Plateau (QTP). They provide biogeographic templates upon which endemic biodiversity can be mapped to infer diversification scenarios. Yet, this process can be somewhat opaque for the Himalaya, given substantial data gaps separating eastern and western regions. To help clarify, we evaluated phylogeographic and phylogenetic hypotheses for a widespread fish (Snowtrout: Cyprininae; Schizothorax) by sequencing 1,140 base pair of mtDNA cytochrome-b (cytb) from Central Himalaya samples (Nepal: N = 53; Bhutan: N = 19), augmented with 68 GenBank sequences (N = 60 Schizothorax/N = 8 outgroups). Genealogical relationships (N = 132) were analyzed via maximum likelihood (ML), Bayesian (BA), and haplotype network clustering, with clade divergence estimated via TimeTree. Snowtrout seemingly originated in Central Asia, dispersed across the QTP, then into Bhutan via southward-flowing tributaries of the east-flowing Yarlung-Tsangpo River (YLTR). Headwaters of five large Asian rivers provided dispersal corridors from Central into eastern/southeastern Asia. South of the Himalaya, the YLTR transitions into the Brahmaputra River, facilitating successive westward colonization of Himalayan drainages first in Bhutan, then Nepal, followed by far-western drainages subsequently captured by the (now) westward-flowing Indus River. Two distinct Bhutanese phylogenetic groups were recovered: Bhutan-1 (with three subclades) seemingly represents southward dispersal from the QTP; Bhutan-2 apparently illustrates northward colonization from the Lower Brahmaputra. The close phylogenetic/phylogeographic relationships between the Indus River (Pakistan) and western tributaries of the Upper Ganges (India/Nepal) potentially implicate an historic, now disjunct connection. Greater species-divergences occurred across rather than within-basins, suggesting vicariance as a driver. The Himalaya is a component of the Earth's largest glacial reservoir (i.e., the "third-pole") separate from the Arctic/Antarctic. Its unique aquatic biodiversity must be defined and conserved through broad, trans-national collaborations. Our study provides an initial baseline for this process.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

Himalayan upliftment and Shiwalik succession act as a cradle for divergence in Bengal monitor lizard Varanus bengalensis (Reptilia: Varanidae) in India

The Himalayan foothills and associated environment are well-known for driving the rapid diversification of many species and the formation of biodiversity hotspots. The effects of environmental change since the Miocene have accelerated species diversification, and hence are useful for studying population genetic structure, and evolutionary relationships via genetic approaches. To date, the effects of climatic fluctuations on the biogeography of large-bodied lizards have not been assessed comprehensively. Herein, we examine the diversification of Varanus bengalensis, focusing on its genetic structure to provide insights into how landscape structure and climatic fluctuations have shaped species differentiation. We confirm the existence of two distinct lineages within V. bengalensis distributed across the Himalayan foothills and the remainder of mainland India. Divergence analyses revealed the split between the Himalayan foothills and the remainder of the mainland lineages of V. bengalensis in the mid-Pliocene ~3.06 Ma, potentially as a consequence of the Siwalik broadening and climatic fluctuations across the Himalayan foothills. The results suggest recognition of a new lineage of V. bengalensis from the Himalayan foothills as a distinctive evolutionarily significant unit.

Himalayan Marmot (Marmota himalayana) Redistribution to High Latitudes under Climate Change

Climate warming and human activities impact the expansion and contraction of species distribution. The Himalayan marmot (Marmota himalayana) is a unique mammal and an ecosystem engineer in the Qinghai-Tibet Plateau (QTP). This pest aggravates grassland degradation and is a carrier and transmitter of plagues. Therefore, exploring the future distribution of Himalayan marmots based on climate change and human activities is crucial for ecosystem management, biodiversity conservation, and public health safety. Here, a maximum entropy model was explored to forecast changes in the distribution and centroid migration of the Himalayan marmot in the 2050s and 2070s. The results implied that the human footprint index (72.80%) and altitude (16.40%) were the crucial environmental factors affecting the potential distribution of Himalayan marmots, with moderately covered grassland being the preferred habitat of the Himalayan marmot. Over the next 30-50 years, the area of suitable habitat for the Himalayan marmot will increase slightly and the distribution center will shift towards higher latitudes in the northeastern part of the plateau. These results demonstrate the influence of climate change on Himalayan marmots and provide a theoretical reference for ecological management and plague monitoring.

Geographical distribution and conservation strategy of national key protected wild plants of China

National key protected wild plants (NKPWPs) are considered flagship species for plant diversity conservation in China. Using data for 1101 species, we characterized NKPWPs distribution patterns in China and assessed conservation effectiveness and conservation gaps. In total, 4880 grid cells at a 20 × 20 km resolution were filled with occurrence records for NKPWPs. We identified 444 hotspot grid cells and 27 diversity hotspot regions, containing 92.37% of NKPWPs. However, 43.24% of these hotspot grid cells were fully or partially covered by national nature reserves (NNRs), where 70.21% of species were distributed. Approximately 61.49% of the NKPWPs species were protected by NNRs, but the populations or habitats of 963 species were partially or fully outside of NNRs. With global warming, the overall change in the extent of suitable habitats for NKPWPs is expected to be small, however, habitat quality in some areas with a high habitat suitability index will decrease.

Molecular phylogeny of mega-diverse Carabus attests late Miocene evolution of alpine environments in the Himalayan-Tibetan Orogen

The timing, sequence, and scale of uplift of the Himalayan-Tibetan Orogen (HTO) are controversially debated. Many geoscientific studies assume paleoelevations close to present-day elevations and the existence of alpine environments across the HTO already in the late Paleogene, contradicting fossil data. Using molecular genetic data of ground beetles, we aim to reconstruct the paleoenvironmental history of the HTO, focusing on its southern margin (Himalayas, South Tibet). Based on a comprehensive sampling of extratropical Carabus, and ~ 10,000 bp of mitochondrial and nuclear DNA we applied Bayesian and Maximum likelihood methods to infer the phylogenetic relationships. We show that Carabus arrived in the HTO at the Oligocene-Miocene boundary. During the early Miocene, five lineages diversified in different parts of the HTO, initially in its southern center and on its eastern margin. Evolution of alpine taxa occurred during the late Miocene. There were apparently no habitats for Carabus before the late Oligocene. Until the Late Oligocene elevations must have been low throughout the HTO. Temperate forests emerged in South Tibet in the late Oligocene at the earliest. Alpine environments developed in the HTO from the late Miocene and, in large scale, during the Pliocene-Quaternary. Findings are consistent with fossil records but contrast with uplift models recovered from stable isotope paleoaltimetry.

Impacts of Quaternary glaciation, geological history and geography on animal species history in continental East Asia: A phylogeographic review

Continental East Asia has a mild Pleistocene climate and a complex recent geological history. Phylogeographic studies of animals over the last 30 years have produced several distinctive patterns. Glaciation refugia are numerous and are not restricted to any particular regions. Most of them are localized and species-specific, although several large refugia, for example the mountains of SW China, are shared by multiple species and have refugia-within-refugia. Furthermore, postglaciation range expansion events vary greatly in time, scale and direction. Large-scale south-to-north post-LGM expansions are few and mostly occurred in the northern regions. Additionally, several unique geographic features, including the three-step terrain of China and the northern arid belt, have significant impacts on many species histories. Overall, the impacts of Pleistocene glaciations, particularly the LGM, on species history vary drastically from nondetectable to significant. The impacts are the least for species from the southwestern region and are most dominant for species from the north. Geological events play a more significant role in shaping species history than Pleistocene climatic changes. Phylogeographic patterns among animals species are highly consistent with those of plants. Future phylogeographic endeavour in East Asia should be hypothesis-driven and seek processes that underlie common patterns. The wide use of genomic data allow accurate estimates of historical population processes and exploration of older history beyond the Pleistocene.

Phylogenomics and historical biogeography of Hydrangeeae (Hydrangeaceae) elucidate the effects of geologic and climatic dynamics on diversification

Demonstrating the process of transregional biogeography and mechanisms underlying evolutionary radiations is crucial to understanding biological evolution. Here, we use Hydrangeeae (Hydrangeaceae), a tribe with a unique disjunct distribution and complex trait variations, using a solid phylogenetic framework, to investigate how geographical and climatic factors interact with functional traits to trigger plant evolutionary radiations. We constructed the first highly supported and dated phylogenetic framework using 79 protein-coding genes obtained from 81 plastomes, representing 63 species and all major clades, and found that most extant species originated from asynchronous diversification of two lineages undergoing repeated expansion and retraction, at middle and high latitudes of the Northern Hemisphere between East Asia and North America, during the Eocene to Pleistocene (driven by geologic and climatic dynamics). In accordance with these drivers, interactions of flora between central-eastern China and Japan occurred frequently after the Late Tertiary. We found that resource limitation and range fragmentation probably accelerated the diversification of Hydrangeeae, which supports the resource-use hypothesis. Our study sheds light on the evolutionary radiation and assembly of flora within East Asia, and the East Asian-North American disjunction, through integration of phylogenomic and biogeographic data with functional trait and ecological data.