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

Phylogenomic analyses shed new light on the spatiotemporal evolution of global larches: Implications for the dynamics of boreal forests

As the Earth warms, understanding the long-term dynamics of forest ecosystems is essential for guiding forest management and biodiversity conservation. Insights from past dynamics may provide valuable lessons for managing today's forests. Here, we investigated the spatiotemporal evolution of global larches to gain further insights into how boreal forests change over time. We first reconstructed a highly resolved and robust phylogeny of Larix covering all widely recognized species, using both transcriptome-based 1,301 orthologous genes (OGs) and plastid genomes. In sharp contrast to previous studies, an unexpected deep split between the circumboreal and Qinghai-Tibetan Plateau (QTP) larches was revealed in our study. Within each lineage, two geographically distinct clades were further resolved. Biogeographical analyses suggest that Larix might have an origin of Eocene in high-latitude uplands, and during the Miocene, all extant species have appeared. Cenozoic climate- and orogeny-triggered vicariance likely played a major role in the divergence of global larches. Our results also demonstrate that the proto-boreal forest biome may have a relatively old origin back to the early Miocene, and significant winnowing and species alteration would have occurred as the climate shifted to much colder and drier conditions during the Neogene. Ecological niche analyses show various responses of the circumboreal and QTP larches under different climate scenarios, but both lineages are negatively impacted by warming climates. These findings have important conservation implications given the sensitivity of boreal forests in the face of global warming. Our work further emphasizes the importance of a solid phylogenetic framework for evolutionary and biogeographical inferences.

Genome sequencing and transcriptome analysis provide an insight into the hypoxia resistance of Channa asiatica

Channa asiatica is an economically valuable fish species and excellent model for studying hypoxic tolerance. However, the underlying genetic and molecular mechanisms are poorly understood. In this study, we assembled a high-quality C. asiatica genome (23 chromosomes, totaling 722 Mb) using a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing. Repetitive elements accounted for 28.39%of the C. asiatica genome, and 23,949 protein-coding genes were predicted, with 96.63 % of these functionally annotated. Moreover, a comparative genomic analysis of 12 fish genomes showed that gene families associated with oxygen binding and transport were expanded in C. asiatica. In addition, transcriptome analysis revealed that multiple oxidative stress pathways were activated when C. asiatica was exposed to air. In conclusion, this study provided high-quality genome assembly and transcriptome data, both serving as critical resources for researching the genetic basis of hypoxic tolerance in C. asiatica.

Recent and Recurrent Autopolyploidization Fueled Diversification of Snow Carp on the Tibetan Plateau

Whole-genome duplication (WGD), or polyploidization, is a major contributor to biodiversity. However, the establishment and survival of WGDs are often considered to be stochastic, since elucidating the processes of WGD establishment remains challenging. In the current study, we explored the processes leading to polyploidy establishment in snow carp (Cyprinidae: Schizothoracinae), a predominant component of the ichthyofauna of the Tibetan Plateau and its surrounding areas. Using large-scale genomic data from isoform sequencing, we analyzed ohnolog genealogies and divergence in hundreds to thousands of gene families across major snow carp lineages. Our findings demonstrated that independent autopolyploidization subsequent to speciation was prevalent, while autopolyploidization followed by speciation also occurred in the diversification of snow carp. This was further supported by matrilineal divergence and drainage evolution evidence. Contrary to the long-standing hypothesis that ancient polyploidization preceded the diversification of snow carp, we determined that polyploidy in extant snow carp was established by recurrent autopolyploidization events during the Pleistocene. These findings indicate that the diversification of extant snow carp resembles a coordinated duet: first, the uplift of the Tibetan Plateau orchestrated the biogeography and diversification of their diploid progenitors; then, the extensive Pliocene-Pleistocene climate changes acted as relay runners, further fueling diversification through recurrent autopolyploidization. Overall, this study not only reveals a hitherto unrecognized recent WGD lineage in vertebrates but also advances current understanding of WGD processes, emphasizing that WGD establishment is a nonstochastic event, emerging from numerous adaptations to environmental challenges and recurring throughout evolutionary history rather than merely in plants.

Contrasting demographic history and mutational load in three threatened whitebark pines (Pinus subsect. Gerardianae): implications for conservation

Demographic history and mutational load are of paramount importance for the adaptation of the endangered species. However, the effects of population evolutionary history and genetic load on the adaptive potential in endangered conifers remain unclear. Here, using population transcriptome sequencing, whole chloroplast genomes and mitochondrial DNA markers, combined with niche analysis, we determined the demographic history and mutational load for three threatened whitebark pines having different endangered statuses, Pinus bungeana, P. gerardiana and P. squamata. Demographic inference indicated that severe bottlenecks occurred in all three pines at different times, coinciding with periods of major climate and geological changes; in contrast, while P. bungeana experienced a recent population expansion, P. gerardiana and P. squamata maintained small population sizes after bottlenecking. Abundant homozygous-derived variants accumulated in the three pines, particularly in P. squamata, while the species with most heterozygous variants was P. gerardiana. Abundant moderately and few highly deleterious variants accumulated in the pine species that have experienced the most severe demographic bottlenecks (P. gerardiana and P. squamata), most likely because of purging effects. Finally, niche modeling showed that the distribution of P. bungeana might experience a significant expansion in the future, and the species' identified genetic clusters are also supported by differences in the ecological niche. The integration of genomic, demographic and niche data has allowed us to prove that the three threatened pines have contrasting patterns of demographic history and mutational load, which may have important implications in their adaptive potential and thus are also key for informing conservation planning.

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.

Plant-Soil Moisture Positive Feedback Maintaining Alternative Stable States in the Alpine Marsh Ecosystem

A self-reinforcing positive feedback is regarded as a critical process for maintaining alternative stable states (ASS); however, identification of ASS and quantification of positive feedbacks remain elusive in natural ecosystems. Here, we used large-scale field surveys to search for ASS and a positive feedback mechanism under a wide range of habitats on the Tibetan Plateau. Using multiple methods, we proved that three stable states exist that accompany alpine marsh degradation. Positive feedbacks between changing soil moisture and plant community composition forced the ecosystem into another stable state, and the alteration of water use efficiency (WUE) of the component species contributed to this shift. This study provides the first empirical evidence that positive feedback loops maintain ASS in the alpine marsh ecosystem on the Tibetan Plateau. Our research revealed the powerful driving role of plants in transitions between states, which may support the conservation and restoration of global alpine marsh ecosystems.

Mitochondrial phylogeography of grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) endemic to the Qinghai-Tibetan plateau

Grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) are the most damaging pests to alpine meadows in the Qinghai-Tibetan Plateau (QTP). Here, we conducted extensive sampling from 39 geographic populations covering almost the entire distribution of the eight QTP Gynaephora (Hübner) species to investigate phylogeographic patterns and speciation based on two mitochondrial genes (COI and ND5). A total of 40 haplotypes were detected in the 39 populations, with >70% of all haplotypes not shared between populations. The monophyletic QTP Gynaephora migrated from non-QTP regions during the Pliocene, corresponding to the uplift of the QTP, suggesting a mode of transport into the QTP. Among the eight QTP Gynaephora species described by morphological characteristics, two species (G. alpherakii and G. menyuanensis) were recovered as monophyletic groups (Clades B and C), while the remaining six formed two monophyletic clades: Clade A (G. qinghaiensis, G. jiuzhiensis, and G. qumalaiensis) and Clade D (G. aureata, G. ruoergensis, and G. minora). These results suggested that the number of the QTP Gynaephora species may be overestimated and further studies based on both morphological and nuclear gene data are needed. Genetic differentiation and speciation of the QTP Gynaephora were likely driven by the QTP uplifts and associated climate fluctuations during the Pleistocene, indicated by divergence time estimation, suggesting that isolation and subsequent divergence was the dominant mode of speciation. The Sanjiangyuan region (i.e., Clade A, characterized by high genetic diversity) may have been a glacial refugium of the QTP Gynaephora, as supported by analyses of gene flow and biogeography. High levels of genetic diversity were found in QTP Gynaephora, without population expansion, which may explain the high-altitude adaptation and outbreaks of grassland caterpillars in alpine meadows of the QTP. This study provides the largest phylogeographic analysis of QTP Gynaephora and improves our understanding of the diversity and speciation of QTP insects.

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.

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.

Amphibians rise to flourishing under climate change on the Qinghai-Tibetan Plateau

Amphibian populations are declining globally due to climate change. However, the impacts on the geographic distribution of amphibians on the Qinghai-Tibetan Plateau (QTP), a global biodiversity hotspot with 112 species of amphibians that is sensitive to global climate change, remains unclear. In this study, MaxEnt and barycentre shift analyses were performed to reveal the impact of climate change on the potential future habitats of amphibians on the QTP using the BCC-CSM2-MR global climate model of the Coupled Model Intercomparison Projects Phase 6 (CMIP6) climate pattern with three shared socioeconomic pathways (SSP). In contrast to the widespread decline in the amphibian population, the future scenarios projected an increase in most amphibian habitats on the QTP, accompanied by migration to higher elevations or latitudes under three climatic projections (SSP 1-2.6, 3-7.0, and 5-8.5). Average annual precipitation was the most crucial environmental variable impacting the future distribution of amphibians. The findings indicate that amphibians would flourish under climate change on the QTP, which is of great significance for the protection of amphibians and biodiversity on the QTP.

Incorporating New Datatypes to Enhance Species Delimitation: A Case Study in Rice Paddy Snakes (Homalopsidae: Hypsiscopus)

Homalopsids (Old World Mud Snakes) include 59 semiaquatic species in Asia and Australasia that display an array of morphological adaptations, behaviors, and microhabitat preferences. These attributes make homalopsids an ideal model system for broader questions in evolutionary biology, but the diversity of this understudied group of snakes is still being described. Recognized species diversity in rice paddy snakes (Hypsiscopus) has recently doubled after nearly 200 years of taxonomic stability. However, the evolutionary distinctiveness of some populations remains in question. In this study, we compare mainland Southeast Asian populations of Hypsiscopus east and west of the Red River Basin in Vietnam, a known biogeographic barrier in Asia, using an iterative approach with molecular phylogenetic reconstruction, machine-learning morphological quantitative statistics, and ecological niche modeling. Our analyses show that populations west of the Red River Basin represent an independent evolutionary lineage that is distinct in genetics, morphospace, and habitat suitability, and so warrants species recognition. The holotype of H. wettsteini, a species originally described in error from Costa Rica, grouped morphometrically with the population at the Red River Basin and eastward, and those west of the Red River Basin are referred to the recently described H. murphyi. The two species may have diversified due to a variety of geological and environmental factors, and their recognition exemplifies the importance of multifaceted approaches in taxonomy for downstream biogeographic studies on speciation scenarios.

Comparative and phylogenetic analysis of chloroplast genomes from ten species in Quercus section Cyclobalanopsis

The genus Quercus L. is widely acknowledged as a significant assemblage within East Asia tropical and subtropical broadleaf evergreen forests, possessing considerable economic importance. Nevertheless, the differentiation of Quercus species is deemed arduous, and the interrelations among these species remain enigmatic. Leveraging Illumina sequencing, we undertook the sequencing and assembly of the chloroplast (cp) genomes of seven species belonging to Quercus section Cyclobalanopsis (Quercus argyrotricha, Q. augustinii, Q. bambusifolia, Q. bella, Q. edithiae, Q. jenseniana, and Q. poilanei). Furthermore, we collated three previously published cp genome sequences of Cyclobalanopsis species (Q. litseoides, Q. obovatifolia, and Q. saravanensis). Our primary objective was to conduct comparative genomics and phylogenetic analyses of the complete cp genomes of ten species from Quercus section Cyclobalanopsis. This investigation unveiled that Quercus species feature a characteristic circular tetrad structure, with genome sizes ranging from 160,707 to 160,999 base pairs. The genomic configuration, GC content, and boundaries of inverted repeats/single copy regions exhibited marked conservation. Notably, four highly variable hotspots were identified in the comparative analysis, namely trnK-rps16, psbC-trnS, rbcL-accD, and ycf1. Furthermore, three genes (atpF, rpoC1, and ycf2) displayed signals of positive selection pressure. Phylogenetic scrutiny revealed that the four sections of Cyclobalanopsis clustered together as sister taxa. The branch support values ranged from moderate to high, with most nodes garnering 100% support, underscoring the utility of cp genomic data in elucidating the relationships within the genus. Divergence time analysis revealed that Section Cyclobalanopsis represents the earliest type of Quercus genus. The outcomes of this investigation establish a foundation for forthcoming research endeavors in taxonomy and phylogenetics.

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.".

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.

On eight species of Chrysillini from Xizang, China (Araneae: Salticidae: Salticinae)

Six new species of Chrysillini are described from Xizang, China: Chrysilla yarlungzangbo sp. nov. (♂), Icius faker sp. nov. (♂♀), I. han sp. nov. (♂♀), I. kui sp. nov. (♂♀), I. kulakangri sp. nov. (♂♀), I. zhengi sp. nov. (♂♀). In addition, a new combination is proposed: Nandicius gyirongensis (Hu, 2001) comb. nov. (♂♀) (transferred from Psenuc Prószyński, 2016), and the males of this species are described for the first time. Siler niser Caleb, Parag & Datta-Roy, 2023 (♂♀) is newly recorded in 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.