Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau

Pleistocene glaciation advances the cryptic speciation of Stellera chamaejasme L. in a major biodiversity hotspot

The mountains of Southwest China comprise a significant large mountain range and biodiversity hotspot imperiled by global climate change. The high species diversity in this mountain system has long been attributed to a complex set of factors, and recent large-scale macroevolutionary investigations have placed a broad timeline on plant diversification that stretches from 10 million years ago (Mya) to the present. Despite our increasing understanding of the temporal mode of speciation, finer-scale population-level investigations are lacking to better refine these temporal trends and illuminate the abiotic and biotic influences of cryptic speciation. This is largely due to the dearth of organismal sampling among closely related species and populations, spanning the incredible size and topological heterogeneity of this region. Our study dives into these evolutionary dynamics of speciation using genomic and eco-morphological data of Stellera chamaejasme L. We identified four previously unrecognized cryptic species having indistinct morphological traits and large metapopulation of evolving lineages, suggesting a more recent diversification (~2.67-0.90 Mya), largely influenced by Pleistocene glaciation and biotic factors. These factors likely influenced allopatric speciation and advocated cyclical warming-cooling episodes along elevational gradients during the Pleistocene. The study refines the evolutionary timeline to be much younger than previously implicated and raises the concern that projected future warming may influence the alpine species diversity, necessitating increased conservation efforts.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

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.

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.

Evolutionary history of an Irano-Turanian cushion-forming legume (Onobrychis cornuta)

The Irano-Turanian region is one of the largest floristic regions in the world and harbors a high percentage of endemics, including cushion-like and dwarf-shrubby taxa. Onobrychis cornuta is an important cushion-forming element of the subalpine/alpine flora of the Irano-Turanian floristic region. To specify the genetic diversity among the populations of this species (including individuals of O. elymaitica), we employed nrDNA ITS and two noncoding regions of plastid DNA (rpl32-trnL(UAG) and trnT(UGU)-trnL(UAA)). The most striking feature of O. cornuta assemblages was the unexpectedly high nucleotide diversity in both the nDNA and cpDNA dataset. In the analyses of nuclear and plastid regions, 25 ribotypes and 42 haplotypes were found among 77 and 59 accessions, respectively, from Iran, Turkey, and Afghanistan. Network analysis of the datasets demonstrated geographic differentiation within the species. Phylogenetic analyses of all dataset retrieved O. cornuta as a non-monophyletic species due to the inclusion of O. elymaitica, comprising four distinct lineages. In addition, our analyses showed cytonuclear discordance between both nuclear and plastid topologies regarding the position of some O. cornuta individuals. The underlying causes of this inconsistency remain unclear. However, we speculate that chloroplast capture, incomplete lineage sorting, and introgression were the main reasons for this event. Furthermore, molecular dating analysis indicated that O. cornuta originated in the early Pliocene (around 4.8 Mya) and started to diversify throughout the Pliocene and in particular the Pleistocene. Moreover, O. elymaitica was reduced to a subspecific rank within the species.

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.

Repeated upslope biome shifts in Saxifraga during late-Cenozoic climate cooling

Mountains are among the most biodiverse places on Earth, and plant lineages that inhabit them have some of the highest speciation rates ever recorded. Plant diversity within the alpine zone - the elevation above which trees cannot grow-contributes significantly to overall diversity within mountain systems, but the origins of alpine plant diversity are poorly understood. Here, we quantify the processes that generate alpine plant diversity and their changing dynamics through time in Saxifraga (Saxifragaceae), an angiosperm genus that occurs predominantly in mountain systems. We present a time-calibrated molecular phylogenetic tree for the genus that is inferred from 329 low-copy nuclear loci and incorporates 73% (407) of known species. We show that upslope biome shifts into the alpine zone are considerably more prevalent than dispersal of alpine specialists between regions, and that the rate of upslope biome shifts increased markedly in the last 5 Myr, a timeframe concordant with a cooling and fluctuating climate that is likely to have increased the extent of the alpine zone. Furthermore, alpine zone specialists have lower speciation rates than generalists that occur inside and outside the alpine zone, and major speciation rate increases within Saxifraga significantly pre-date increased rates of upslope biome shifts. Specialisation to the alpine zone is not therefore associated with speciation rate increases. Taken together, this study presents a quantified and broad scale perspective of processes underpinning alpine plant diversity.

Comparatively analyzing of chloroplast genome and new insights into phylogenetic relationships regarding the genus Stephania

The medicinal plant of the genus Stephania holds significant economic importance in the pharmaceutical industry. However, accurately classifying and subdividing this genus remains a challenge. Herein, the chloroplast (cp) genomes of Stephania and Cyclea were sequenced, and the primary characteristics, repeat sequences, inverted repeats regions, simple sequence repeats, and codon usage bias of 17 species were comparatively analyzed. Twelve markers were identified through genome alignment and sliding window analysis. Moreover, a molecular clock analysis revealed the divergence between subgenus (subg.) Botryodiscia and the combined Cyclea, subg. Stephania and Tuberiphania during the early Oligocene epoch. Notably, the raceme-type inflorescence represents the ancestral state of the Stephania and Cyclea. The genetic relationships inferred from the cp genome and protein-coding genes exhibited similar topologies. Additionally, the paraphyletic relationship between the genera Cyclea and Stephania was confirmed. Bayesian inference, maximum likelihood, and neighbor-joining trees consistently showed that section Tuberiphania and Transcostula were non-monophyletic. In conclusion, this research provides valuable insights for further investigations into species identification, evolution, and phylogenetics within the Stephania genus.

Hedysarumqilianshanense sp. nov. (Fabaceae, Hedysareae), a new species from the Qilianshan Mountains in Gansu, China

Hedysarumqilianshanensesp. nov. (Fabaceae, Hedysareae) is described and illustrated from the Qilianshan Mountains in Gansu, China. This new species is similar to H.przewalskii, but can be distinguished by its corolla being light purple to purple, standard 15-19 mm long, wings 14-16 mm long, keels 16-19 mm long, and the ovary and legume being glabrous. The new species can be easily distinguished from H.neglectum Ledeb. by its bract being shorter than the pedicel, and the ovary and legume being glabrous. Phylogenetic tree based on the nuclear ITS and ETS sequences shows that H.qilianshanense is sister to H.przewalskii, while the tree based on the plastid psbA-trnH, trnC-petN, trnL-F, trnS-G and petN-psbM sequences shows H.qilianshanense as sister to a clade consisting of H.hedysaroides, H.inundatum, H.americanum and H.neglectum. The new species is a diploid with the chromosome number 2n = 14. Based on morphological, phylogenetic and karyotypic evidence, the new species may originate from an ancient homoploid hybrid speciation event.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae)

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

Reticulate evolution of the tertiary relict Osmanthus

When interspecific gene flow is common, species relationships are more accurately represented by a phylogenetic network than by a bifurcating tree. This study aimed to uncover the role of introgression in the evolution of Osmanthus, the only genus of the subtribe Oleinae (Oleaceae) with its distribution center in East Asia. We built species trees, detected introgression, and constructed networks using multiple kinds of sequencing data (whole genome resequencing, transcriptome sequencing, and Sanger sequencing of nrDNA) combined with concatenation and coalescence approaches. Then, based on well-understood species relationships, historical biogeographic analyses and diversification rate estimates were employed to reveal the history of Osmanthus. Osmanthus originated in mid-Miocene Europe and dispersed to the eastern Tibetan Plateau in the late Miocene. Thereafter, it continued to spread eastwards. Phylogenetic conflict is common within the 'Core Osmanthus' clade and is seen at both early and late stages of diversification, leading to hypotheses of net-like species relationships. Incomplete lineage sorting proved ineffective in explaining phylogenetic conflicts and thus supported introgression as the main cause of conflicts. This study elucidates the diversification history of a relict genus in the subtropical regions of eastern Asia and reveals that introgression had profound effects on its evolutionary history.

Plastid phylogenomics of the tribe potentilleae (Rosaceae)

The tribe Potentilleae comprises approximately 1700 species in 13 genera, making it one of the largest of the 16 tribes in Rosaceae. Our understanding of the composition and relationships among members of Potentilleae has advanced dramatically with the application of molecular markers in the last two decades. Yet there is still much work remaining toward a robust phylogenetic framework for the entire Potentilleae and a comprehensive genus-level dating framework for the tribe. The goals of the present study were to establish a phylogenetic framework for Potentilleae, infer the origin and diversification of the tribe using a temporal framework, and explore the taxonomic implications in light of the updated phylogenetic framework. We used the plastome sequences from 158 accessions representing 139 taxa covering all 13 recognized genera of the tribe to reconstruct the Potentilleae phylogeny. High phylogenetic resolution was recovered along the Potentilleae backbone. Two major clades were recovered within Potentilleae, corresponding to the two subtribes Fragariinae and Potentillinae. Within Fragariinae, two subclades were recovered. In one subclade, Sibbaldia sensu stricto is sister to a clade containing Sibbaldianthe, Comarum, Farinopsis, and Alchemilla sensu lato. In the other subclade, Fragaria is sister to a clade comprising Chamaerhodos, Chamaecallis, Drymocallis, Dasiphora, and Potaninia. Within Potentillinae, Argentina is sister to Potentilla sensu stricto. Within Potentilla sensu stricto, clade Himalaya is sister to Alba, and the Himalaya-Alba clade together is sister to a clade comprising Reptans, Potentilla ancistrifolia Bunge, Fragarioides, Ivesioid, and Argentea. Divergence time estimates indicated that tribe Potentilleae originated during the middle Eocene, and subtribes Fragariinae and Potentillinae diverged around the Eocene-Oligocene transition, and divergence times dated for Potentilleae genera ranged from the early Miocene to the late Pleistocene.

Hybridization and divergent climatic preferences drive divergence of two allopatric Gentiana species on the Qinghai-Tibet Plateau

BACKGROUND AND AIMS

Exploring how species diverge is vital for understanding the drivers of speciation. Factors such as geographical separation and ecological selection, hybridization, polyploidization and shifts in mating system are all major mechanisms of plant speciation, but their contributions to divergence are rarely well understood. Here we test these mechanisms in two plant species, Gentiana lhassica and G. hoae, with the goal of understanding recent allopatric species divergence on the Qinghai-Tibet Plateau (QTP).

METHODS

We performed Bayesian clustering, phylogenetic analysis and estimates of hybridization using 561 302 nuclear genomic single nucleotide polymorphisms (SNPs). We performed redundancy analysis, and identified and annotated species-specific SNPs (ssSNPs) to explore the association between climatic preference and genetic divergence. We also estimated genome sizes using flow cytometry to test for overlooked polyploidy.

KEY RESULTS

Genomic evidence confirms that G. lhassica and G. hoae are closely related but distinct species, while genome size estimates show divergence occurred without polyploidy. Gentiana hoae has significantly higher average FIS values than G. lhassica. Population clustering based on genomic SNPs shows no signature of recent hybridization, but each species is characterized by a distinct history of hybridization with congeners that has shaped genome-wide variation. Gentiana lhassica has captured the chloroplast and experienced introgression with a divergent gentian species, while G. hoae has experienced recurrent hybridization with related taxa. Species distribution modelling suggested range overlap in the Last Interglacial Period, while redundancy analysis showed that precipitation and temperature are the major climatic differences explaining the separation of the species. The species differ by 2993 ssSNPs, with genome annotation showing missense variants in genes involved in stress resistance.

CONCLUSIONS

This study suggests that the distinctiveness of these species on the QTP is driven by a combination of hybridization, geographical isolation, mating system differences and evolution of divergent climatic preferences.

Peripatric speciation within Torreya fargesii (Taxaceae) in the Hengduan Mountains inferred from multi-loci phylogeography

BACKGROUND

The Hengduan Mountains (HDM) are one of the major global biodiversity hotspots in the world. Several evolutionary scenarios, especially in-situ diversification, have been proposed to account for the high species richness of temperate plants. However, peripatric speciation, an important mode of allopatric speciation, has seldom been reported in this region.

RESULTS

Here, two chloroplast DNA regions and 14 nuclear loci were sequenced for 112 individuals from 10 populations of Torreya fargesii var. fargesii and 63 individuals from 6 populations of T. fargesii var. yunnanensis. Population genetic analyses revealed that the two varieties are well differentiated genetically (FST, 0.5765) and have uneven genetic diversity (π, 0.00221 vs. 0.00073 on an average of nuclear loci). The gene genealogical relationship showed that T. fargesii var. yunnanensis is inferred as derived from T. fargesii var. fargesii, which was further supported by the coalescent simulations (DIYABC, fastsimcoal2 and IMa2). By the coalescent simulations, the divergence time (~ 2.50-3.65 Ma) and the weak gene flow between the two varieties were detected. The gene flow was asymmetrical and only occurred in later stages of divergence, which is caused by second contact due to the population expansion (~ 0.61 Ma) in T. fargesii var. fargesii. In addition, niche modeling indicated that the two varieties are differentiated geographically and ecologically and have unbalanced distribution range.

CONCLUSIONS

Overall, T. fargesii var. fargesii is always parapatric with respect to T. fargesii var. yunnanensis, and the latter derived from the former in peripatry of the HDM following a colonization from central China during the late Pliocene. Our findings demonstrate that peripatric speciation following dispersal events may be an important evolutionary scenario for the formation of biodiversity hotspot of the HDM.

Impact of the climatic changes in the Pliocene-Pleistocene transition on Irano-Turanian species. The radiation of genus Jurinea (Compositae)

The Irano-Turanian region is one of the world's richest floristic regions and the centre of diversity for numerous xerophytic plant lineages. However, we still have limited knowledge on the timing of evolution and biogeographic history of its flora, and potential drivers of diversification remain underexplored. To fill this knowledge gap, we focus on the Eurasian genus Jurinea (ca. 200 species), one of the largest plant radiations that diversified in the region. We applied a macroevolutionary integrative approach to explicitly test diversification hypotheses and investigate the relative roles of geography vs. ecology and niche conservatism vs. niche lability in speciation processes. To do so, we gathered a sample comprising 77% of total genus richness and obtained data about (1) its phylogenetic history, recovering 502 nuclear loci sequences; (2) growth forms; (3) ecological niche, compiling data of 21 variables for more than 2500 occurrences; and (4) paleoclimatic conditions, to estimate climatic stability. Our results revealed that climate was a key factor in the evolutionary dynamics of Jurinea. The main diversification and biogeographic events that occurred during past climate changes, which led to colder and drier conditions, are the following: (1) the origin of the genus (10.7 Ma); (2) long-distance dispersals from the Iranian Plateau to adjacent regions (∼7-4 Ma); and (3) the diversification shift during Pliocene-Pleistocene Transition (ca. 3 Ma), when net diversification rate almost doubled. Our results supported the pre-adaptation hypothesis, i.e., the evolutionary success of Jurinea was linked to the retention of the ancestral niche adapted to aridity. Interestingly, the paleoclimatic analyses revealed that in the Iranian Plateau long-term climatic stability favoured old-lineage persistence, resulting in current high species richness of semi-arid and cold adapted clades; whereas moderate climate oscillations stimulated allopatric diversification in the lineages distributed in the Circumboreal region. In contrast, growth form lability and high niche disparity among closely related species in the Central Asian clade suggest adaptive radiation to mountain habitats. In sum, the radiation of Jurinea is the result of both adaptive and non-adaptive processes influenced by climatic, orogenic and ecological factors.

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.

Genetic Diversity and Phylogenetic Analysis of Zygophyllum loczyi in Northwest China's Deserts Based on the Resequencing of the Genome

In order to study the genetics of local adaptation in all main deserts of northwest China, whole genomes of 169 individuals were resequenced, which covers 20 populations of Zygophyllum loczyi (Zygophyllales: Zygophylaceae). We describe more than 15 million single nucleotide polymorphisms and numerous InDels. The expected heterozygosity and PIC values associated with local adaptation varied significantly across biogeographic regions. Variation in environmental factors contributes largely to the population genetic structure of Z. loczyi. Bayesian analysis performed with STRUCTURE defined four genetic clusters, while the results of principle component analysis were similar. Our results shows that the Qaidam Desert group appears to be diverging into two branches characterized by significant geographic separation and gene flow with two neighboring deserts. Geological data assume that it is possible that the Taklamakan Desert was the original distribution site, and Z. loczyi could have migrated later on and expanded within other desert areas. The above findings provide insights into the processes involved in biogeography, phylogeny, and differentiation within the northwest deserts of China.

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.

Unveiling the genetic landscape of high-altitude adaptive ethnic groups with polymorphic markers: Implications of comprehensive forensic appraisals and population genetic investigations

Background

Deletion/insertion polymorphisms (DIPs), a novel class of biomarker, have been widely utilized in forensic areas for individual identification, paternity tests, and ancestral origin inference due to its applicability to degraded samples and low mutation rates. Despite the availability of a well-established commercial kit, the Investigator® DIPplex kit (Qiagen), certain loci exhibit limited levels of polymorphisms in East Asian populations, particularly in Chinese populations.

Objective

This dissertation seeks to undertake a comprehensive evaluation about the forensic efficiency of a self-developed multiplex amplification system in high-altitude adaptive ethnic groups of China. Healthy unrelated Tibetan individuals residing in Tibet Autonomous Region and Qinghai Province were genotyped using previously reported 43 deletion/insertion polymorphism loci. Forensic statistical analyses including allele frequencies and forensic parameters were conducted in the two Tibetan groups, and the genetic relatedness of the studied groups with reference populations from the 1000 Genomes Project Phase 3 were investigated.

Results

Forensic statistical results showed that the polymorphism information content values of the 43 deletion/insertion polymorphism loci in the two Tibetan groups exceeded 0.35. Moreover, the combined power of discrimination using the 43 deletion/insertion polymorphism loci was calculated to be 0.9999999999999999984 in the Qinghai Tibetan group and 0.9999999999999999921 in the Tibet Tibetan group. The cumulative power of exclusion using the 43 deletion/insertion polymorphism loci was calculated to be 0.999782512 in the Qinghai Tibetan group and 0.999886205 in the Tibet Tibetan group. Analysis of population genetics demonstrated that the two studied Tibetan groups shared close genetic relationships with East Asia populations.

Conclusion

The set of 43 deletion/insertion polymorphism loci exhibited remarkable forensic efficacy, rendering it a promising tool for forensic practice. Population genetic analyses indicated that the two Tibetan groups had closer genetic affinities to East Asian populations.

Investigation of phylogenetic relationships within Saxifraga diversifolia complex (Saxifragaceae) based on restriction-site associated DNA sequence markers

Subsect. Hirculoideae Engl. & Irmsch., belonging to Saxifraga sect. Ciliatae Haw., has high species richness. It can be divided into S. diversifolia, S. pseudohirculus, and S. sinomontana complexes based on morphological characteristics. The species with prominent leaf veins on the posterior leaf edge were placed in the S. diversifolia complex, which is mainly distributed on the eastern and southern margins of the Qinghai-Tibetan Plateau. In this study, 53 samples, representing 15 of the 33 described species in the S. diversifolia complex, were sequenced using the Restriction-site Associated DNA Sequence (RAD-seq) technique. A total of 111,938 high-quality SNP loci were screened to investigate the phylogenetic relationships within the S. diversifolia complex. The result of the neighbor-joining (NJ) tree shows that the S. diversifolia complex is a paraphyletic group. Despite of some inconsistencies as revealed by genetic structural analysis, clustering results of representative species reconstructed by both NJ and principal component analysis analyses support previous biogeographic and morphological evidences. In addition, long-distance gene flow events for 11 taxa were detected in the S. diversifolia complex, respectively from S. implicans 1 to S. implicans 2, S. diversifolia and S. maxionggouensis, and from S. maxionggouensis to S. nigroglandulifera. These findings may improve our comprehension of the phylogeny, classification, and evolution of the S. diversifolia complex.

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.

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.

Parallel evolution, atavism, and extensive introgression explain the radiation of Epimedium sect. Diphyllon (Berberidaceae) in southern East Asia

East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.

Exploring the effects of palm kernel meal feeding on the meat quality and rumen microorganisms of Qinghai Tibetan sheep

Palm kernel meal (PKM) has been shown to be a high-quality protein source in ruminant feeds. This study focused on the effects of feed, supplemented with different amounts of PKM (ZL-0 as blank group, and ZL-15, ZL-18, and ZL-21 as treatment group), on the quality and flavor profile of Tibetan sheep meat. Furthermore, the deposition of beneficial metabolites in Tibetan sheep and the composition of rumen microorganisms on underlying regulatory mechanisms of meat quality were studied based on ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry as well as 16S rDNA sequencing. The results of the study showed that Tibetan sheep in the ZL-18 group exhibited superior eating quality and flavor profile while depositing more protein and fat relative to the other groups. The ZL-18 group also changed significantly in terms of the concentration and metabolic pathways of meat metabolites, as revealed by metabolomics. Metabolomics and correlation analyses finally showed that PKM feed mainly affected carbohydrate metabolism in muscle, which in turn affects meat pH, tenderness, and flavor. In addition, 18% of PKM increased the abundance of Christensenellaceae R-7 group, Ruminococcaceae UCG-013, Lachnospiraceae UCG-002, and Family XIII AD3011 group in the rumen but decreased the abundance of Prevotella 1; the above bacteria groups regulate meat quality by regulating rumen metabolites (succinic acid, DL-glutamic acid, etc.). Overall, the addition of PKM may improve the quality and flavor of the meat by affecting muscle metabolism and microorganisms in the rumen.

Target sequence capture data shed light on the deeper evolutionary relationships of subgenus Chamaecerasus in Lonicera (Caprifoliaceae)

The genus Lonicera L. is widely distributed in the north temperate zone and is well-known for its high species richness and morphological diversity. Previous studies have suggested that many sections of Lonicera are not monophyletic and phylogenetic relationships within the genus are still poorly resolved. In this study, we sampled 37 accessions of Lonicera, covering four sections of subgenus Chamaecerasus plus six outgroup taxa, to recover the main clades of Lonicera based on sequences of nuclear loci generated by target enrichment and cpDNA from genome skimming. We found extensive cytonuclear discordance across the subgenus. Both nuclear and plastid phylogenetic analyses supported subgenus Chamaecerasus sister to subgenus Lonicera. Within subgenus Chamaecerasus, sections Isika and Nintooa were each polyphyletic. Based on the nuclear and chloroplast phylogenies, we propose to merge Lonicera korolkowii into section Coeloxylosteum and Lonicera caerulea into section Nintooa. In addition, Lonicera is estimated to have originated in the mid Oligocene (26.45 Ma). The stem age of section Nintooa was estimated to be 17.09 Ma (95% HPD: 13.30-24.45). The stem age of subgenus Lonicera was estimated to be 16.35 Ma (95% HPD: 14.12-23.66). Ancestral area reconstruction analyses indicate that subgenus Chamaecerasus originated in East Asia and Central Asia. In addition, sections Coeloxylosteum and Nintooa originated in East Asia, with subsequent dispersals into other areas. The aridification of the Asian interior likely promoted the rapid radiation of sections Coeloxylosteum and Nintooa within this region. Moreover, our biogeographic analysis fully supports the Bering and the North Atlantic Land Bridge hypotheses for the intercontinental migrations in the Northern Hemisphere. Overall, this study provides new insights into the taxonomically complex lineages of subgenus Chamaecerasus and the process of speciation.

Phylogenomics, phylogeography and germplasms authentication of the Rheum palmatum complex based on complete chloroplast genomes

As a traditional Chinese medicine, rhubarb is used to treat several diseases such as severe acute pancreatitis, sepsis and chronic renal failure. However, few studies focused on the authentication of germplasm for the Rheum palmatum complex, and no studies have been conducted to elucidate the evolutionary history of the R. palmatum complex using plastome datasets. Hence, we aim to develop the potential molecular markers to identify the elite germplasms of rhubarb and explore the divergence and biogeographic history of the R. palmatum complex based on the newly sequenced chloroplast genome datasets. Chloroplast genomes of thirty-five the R. palmatum complex germplasms were sequenced, and the length ranged from 160,858 to 161,204 bp. The structure, gene content and gene order were highly conserved across all genomes. Eight InDels and sixty-one SNPs loci could be used to authenticate the high-quality germplasms of rhubarb in specific areas. Phylogenetic analysis revealed that all rhubarb germplasms were clustered in the same clade with high bootstrap support values and Bayesian posterior probabilities. According to the molecular dating result, the intraspecific divergence of the complex occurred in the Quaternary, which might be affected by climatic fluctuation. The biogeography reconstruction indicated that the ancestor of the R. palmatum complex might originate from the Himalaya-Hengduan Mountains or/and Bashan-Qinling Mountains, and then spread to surrounding areas. Several useful molecular markers were developed to identify rhubarb germplasms, and our study will provide further understanding on speciation, divergence and biogeography of the R. palmatum complex.

Origin and diversification of a Himalayan orchid genus Pleione

Pleione is an orchid endemically distributed in high mountain areas across the Hengduan Mountains (HDM), Himalayas, Southeast Asia and South of China. The unique flower shapes, rich colors and immense medicinal importance of Pleione are valuable ornamental and economic resources. However, the phylogenetic relationships and evolutionary history of the genus have not yet been comprehensively resolved. Here, the evolutionary history of Pleione was investigated using single-copy gene single nucleotide polymorphisms and chloroplast genome datasets. The data revealed that Pleione could be divided into five clades. Discordance in topology between the two phylogenetic trees and network and D-statistic analyses indicated the occurrence of reticulate evolution in the genus. The evolution could be attributed to introgression and incomplete lineage sorting. Ancestral area reconstruction suggested that Pleione was originated from the HDM. Uplifting of the HDM drove rapid diversification by creating conditions favoring rapid speciation. This coincided with two periods of consolidation of the Asian monsoon climate, which caused the first rapid diversification of Pleione from 8.87 to 7.83 Mya, and a second rapid diversification started at around 4.05 Mya to Pleistocene. The interaction between Pleione and climate changes, especially the monsoons, led to the current distribution pattern and shaped the dormancy characteristic of the different clades. In addition to revealing the evolutionary relationship of Pleione with orogeny and climate changes, the findings of this study provide insights into the speciation and diversification mechanisms of plants in the East Asian flora.

Simplified Genomic Data Revealing the Decline of Aleuritopteris grevilleoides Population Accompanied by the Uplift of Dry-Hot Valley in Yunnan, China

Understanding the evolutionary history of endangered species is crucial for identifying the main reasons for species endangerment in the past and predicting the changing trends and evolutionary directions of their future distribution. In order to study the impact of environmental changes caused by deep valley incision after the uplift of the Qinghai-Tibet Plateau on endangered species, we collected 23 samples belonging to four populations of Aleuritopteris grevilleoides, an endangered fern endemic to the dry-hot valleys (DHV) of Yunnan. Single-nucleotide variation sites (SNPs) were obtained by the genotyping-by-sequencing (GBS) method, and approximately 8085 SNP loci were identified. Through the reconstruction and analysis of genetic diversity, population structure, population dynamics, evolution time, and ancestral geographical distribution, combined with geological historical events such as the formation of dry-hot valleys, this study explores the formation history, current situation, reasons for endangerment and scientifically sound measures for the protection of A. grevilleoides. In our study, A. grevilleoides had low genetic diversity (Obs_Het = 0.16, Exp_Het = 0.32, Pi = 0.33) and a high inbreeding coefficient (Fis = 0.45). The differentiation events were 0.18 Mya, 0.16 Mya, and 0.11 Mya in the A. grevilleoides and may have been related to the formation of terraces within the dry-hot valleys. The history of population dynamics results shows that the diversion of the river resulted in a small amount of gene flow between the two clades, accompanied by a rapid increase in the population at 0.8 Mya. After that, the effective population sizes of A. grevilleoides began to contract continuously due to topographic changes resulting from the continuous expansion of dry-hot valleys. In conclusion, we found that the environmental changes caused by geological events might be the main reason for the changing population size of A. grevilleoides.

Integrated phylogenomic analyses unveil reticulate evolution in Parthenocissus (Vitaceae), highlighting speciation dynamics in the Himalayan-Hengduan Mountains

Hybridization caused by frequent environmental changes can lead both to species diversification (speciation) and to speciation reversal (despeciation), but the latter has rarely been demonstrated. Parthenocissus, a genus with its trifoliolate lineage in the Himalayan-Hengduan Mountains (HHM) region showing perplexing phylogenetic relationships, provides an opportunity for investigating speciation dynamics based on integrated evidence. We investigated phylogenetic discordance and reticulate evolution in Parthenocissus based on rigorous analyses of plastome and transcriptome data. We focused on reticulations in the trifoliolate lineage in the HHM region using a population-level genome resequencing dataset, incorporating evidence from morphology, distribution, and elevation. Comprehensive analyses confirmed multiple introgressions within Parthenocissus in a robust temporal-spatial framework. Around the HHM region, at least three hybridization hot spots were identified, one of which showed evidence of ongoing speciation reversal. We present a solid case study using an integrative methodological approach to investigate reticulate evolutionary history and its underlying mechanisms in plants. It demonstrates an example of speciation reversal through frequent hybridizations in the HHM region, which provides new perspectives on speciation dynamics in mountainous areas with strong topographic and environmental heterogeneity.

The evolutionary history of rice azaleas (Rhododendron tschonoskii alliance) involved niche evolution to a montane environment

PREMISE

The formation of isolated montane geography on islands promotes evolution, speciation, and then radiation if there are ecological changes. Thus, investigating evolutionary histories of montane species and associated ecological changes may help efforts to understand how endemism formed in islands' montane floras. To explore this process, we investigated the evolutionary history of the Rhododendron tschonoskii alliance, which grows in montane environments of the Japanese archipelago and the Korean Peninsula.

METHODS

We studied the five species in the R. tschonoskii alliance and 30 outgroup species, using genome-wide single-nucleotide polymorphisms and cpDNA sequences, in association with environmental analyses.

RESULTS

The monophyletic R. tschonoskii alliance diverged since the late Miocene. Species in the alliance currently inhabit a cold climatic niche that is largely different from that of the outgroup species. We observed clear genetic and niche differentiations between the taxa of the alliance.

CONCLUSIONS

The association of the alliance's evolution with the formation of cooler climates on mountains indicates that it was driven by global cooling since the mid-Miocene and by rapid uplift of mountains since the Pliocene. The combination of geographic and climatic isolation promoted high genetic differentiation between taxa, which has been maintained by climatic oscillations since the Quaternary.

Phylogenomics of Aralia sect. Aralia (Araliaceae): Signals of hybridization and insights into its species delimitations and intercontinental biogeography

Genome-scale data have significantly increased the number of informative characters for phylogenetic analyses and recent studies have also revealed widespread phylogenomic discordance in many plant lineages. Aralia sect. Aralia is a small plant lineage (14 spp.) of the ginseng family Araliaceae with a disjunct distribution between eastern Asia (11 spp.) and North America (3 spp.). We herein employ sequences of hundreds of nuclear loci and the complete plastomes using targeted sequence capture and genome skimming to reconstruct the phylogenetic and biogeographic history of this section. We detected substantial conflicts among nuclear genes, yet different analytical strategies generated largely congruent topologies from the nuclear data. Significant cytonuclear discordance was detected, especially concerning the positions of the three North American species. The phylogenomic results support two intercontinental disjunctions: (1) Aralia californica of western North America is sister to the eastern Asian clade consisting of A. cordata and A. continentalis in the nuclear tree, and (2) the eastern North American A. racemosa forms a clade with A. bicrenata from southwestern North America, and the North American A. racemosa - A. bicrenata clade is then sister to the eastern Asian clade consisting of A. glabra (Japan), A. fargesii (C China), and A. apioides and A. atropurpurea (the Hengduan Mountains). Aralia cordata is supported to be disjunctly distributed in Japan, Taiwan, the Ulleung island of Korea, and in Central, Southwest and South China, and Aralia continentalis is redefined with a narrower distribution in Northeast China, eastern Russia and peninsular Korea.

Insights into the differentiation and adaptation within Circaeasteraceae from Circaeaster agrestis genome sequencing and resequencing

Circaeaster agrestis and Kingdonia uniflora are sister species that reproduce sexually and mainly asexually respectively, providing a good system for comparative genome evolution between taxa with different reproductive models. Comparative genome analyses revealed the two species have similar genome size, but C. agrestis encodes many more genes. The gene families specific to C. agrestis show significant enrichment of genes associated with defense response, while those gene families specific to K. uniflora are enriched in genes regulating root system development. Collinearity analyses revealed C. agrestis experienced two rounds of whole-genome duplication. Fst outlier test across 25 C. agrestis populations uncovered a close inter-relationship between abiotic stress and genetic variability. Genetic feature comparisons showed K. uniflora presents much higher genome heterozygosity, transposable element load, linkage disequilibrium degree, and π_N/π_S ratio. This study provides new insights into understanding the genetic differentiation and adaptation within ancient lineages characterized by multiple reproductive models.

Genomic clines across the species boundary between a hybrid pine and its progenitor in the eastern Tibetan Plateau

Most species have clearly defined distribution ranges and ecological niches. The genetic and ecological causes of species differentiation and the mechanisms that maintain species boundaries between newly evolved taxa and their progenitors are, however, less clearly defined. This study investigated the genetic structure and clines in Pinus densata, a pine of hybrid origin on the southeastern Tibetan Plateau, to gain an understanding of the contemporary dynamics of species barriers. We analyzed genetic diversity in a range-wide collection of P. densata and representative populations of its progenitors, Pinus tabuliformis and Pinus yunnanensis, using exome capture sequencing. We detected four distinct genetic groups within P. densata that reflect its migration history and major gene-flow barriers across the landscape. The demographies of these genetic groups in the Pleistocene were associated with regional glaciation histories. Interestingly, population sizes rebounded rapidly during interglacial periods, suggesting persistence and resilience of the species during the Quaternary ice age. In the contact zone between P. densata and P. yunnanensis, 3.36% of the analyzed loci (57 849) showed exceptional patterns of introgression, suggesting their potential roles in either adaptive introgression or reproductive isolation. These outliers showed strong clines along critical climate gradients and enrichment in a number of biological processes relevant to high-altitude adaptation. This indicates that ecological selection played an important role in generating genomic heterogeneity and a genetic barrier across a zone of species transition. Our study highlights the forces that operate to maintain species boundaries and promote speciation in the Qinghai-Tibetan Plateau and other mountain systems.

The phylogeny and divergence time of Ophiocordyceps sinensis and its host insects based on elongation factor 1 alpha

Ophiocordyceps sinensis Berk. is a fungal parasite that parasitizes the larvae of Hepialidae and is endemic to the Qinghai-Tibet Plateau (QTP). The phylogeny and divergence time of O. sinensis and its host insects were analyzed for 137 individuals from 48 O. sinensis populations based on the elongation factor 1 alpha (EF-1α) gene. Lower nucleotide variation, with only 7 and 16 EF-1α haplotypes, was detected in O. sinensis and its host insects, respectively. The isolated and broad distribution patterns coexisted in both O. sinensis and its host insects on the QTP. The divergence time estimates show that O. sinensis and its host insects originated later than 14.33 million years (Myr) and earlier than 23.60 Myr in the Miocene period, and the major differentiation occurred later than 4 Myr. Their origin and differentiation match well with the second and third uplifts of the QTP, respectively. The host insects from the O. sinensis populations distributed around Qinghai Lake are inferred as an ancient and relict species that has survived various geological events of the QTP. It is suitable to estimate the divergence times of both O. sinensis and its host insects from the same individuals using one gene: EF-1α. Our findings of the origin, phylogeny, and evolution of the endemic species also support the epoch of geological events on the QTP.

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

Introduction

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

Methods

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

Results

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

Conclusion

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

How to survive in the world's third poplar: Insights from the genome of the highest altitude woody plant, Hippophae tibetana (Elaeagnaceae)

Hippophae tibetana (Tibetan sea-buckthorn) is one of the highest distributed woody plants in the world (3,000-5,200 meters a.s.l.). It is characterized by adaptation to extreme environment and important economic values. Here, we combined PacBio Hifi platform and Hi-C technology to assemble a 1,452.75 Mb genome encoding 33,367 genes with a Contig N50 of 74.31 Mb, and inferred its sexual chromosome. Two Hippophae-specific whole-genome duplication events (18.7-21.2 million years ago, Ma; 28.6-32.4 Ma) and long terminal repeats retroelements (LTR-RTs) amplifications were detected. Comparing with related species at lower altitude, Ziziphus jujuba (<1, 700 meters a.s.l.), H. tibetana had some significantly rapid evolving genes involved in adaptation to high altitude habitats. However, comparing with Hippophae rhamnoides (<3, 700 meters a.s.l.), no rapid evolving genes were found except microtubule and microtubule-based process genes, H. tibetana has a larger genome, with extra 2, 503 genes (7.5%) and extra 680.46 Mb transposable elements (TEs) (46.84%). These results suggest that the changes in the copy number and regulatory pattern of genes play a more important role for H. tibetana adapting to more extreme and variable environments at higher altitude by more TEs and more genes increasing genome variability and expression plasticity. This suggestion was supported by two findings: nitrogen-fixing genes of H. tibetana having more copies, and intact TEs being significantly closer genes than fragmentary TEs. This study provided new insights into the evolution of alpine plants.

The significance of recent diversification in the Northern Hemisphere in shaping the modern global flora revealed from the herbaceous tribe of Rubieae (Rubiaceae)

The global herbaceous flora is probably shaped by both ancient and/or recent diversification, companied with the impacts from geographic differences between the Northern and Southern Hemispheres. Therefore, its biogeographic pattern with respect to temporal and spatial divergence is far from full understanding. Tribe Rubieae, the largest herbaceous tribe in the woody-dominant Rubiaceae, provides an excellent opportunity for studying the macroevolution of worldwide colonization. Here, we aim to reconstruct the evolutionary history of Rubieae with regard to climate fluctuation and geological history in the Cenozoic. A total of 204 samples of Rubieae representing all the distribution areas of the tribe were used to infer its phylogenetic and biogeographic histories based on two nrDNA and six cpDNA regions. The ancestral area of Rubieae was reconstructed using a time-calibrated phylogeny in RASP and diversification rates were inferred using Bayesian analysis of macroevolutionary mixtures (BAMM). Our results show Rubieae probably originated in European region during the middle Oligocene, with the two subtribes separating at 26.8 million years ago (Ma). All the genera in Rubieae formed separate clades between 24.79 and 6.23 Ma. The ancestral area of the subtribe Rubiinae was the Madrean-Tethyan plant belt and the North Atlantic land bridge (NALB) provided passage between North America and Europe for Rubiinae. The subtribe Galiinae clade originated in Europe/central Asia during the late Oligocene. Two diversification shifts were detected within Rubieae in the late Neogene. Most extant Rubieae species diverged recently during the Neogene within clades that generally were established during the late Paleogene. The tribe shows complex migration/dispersal patterns within the North Hemisphere combined with multiple recent dispersals into Southern Hemisphere. Our results highlighted the important role of recent biogeographic diversification in the Northern Hemisphere in shaping the modern global herbaceous flora during the latest and rapid worldwide expansion in the Neogene.

Comparative plastomes and phylogenetic analysis of seven Korean endemic Saussurea (Asteraceae)

BACKGROUND

Saussurea is one of the most species-rich genera in the Cardueae, Asteraceae. There are approximately 40 Saussurea species distributed in Korea, with nearly 40% of them endemics. Infrageneric relationships remain uncertain due to insufficient resolutions and low statistical support. In this study, we sequenced the plastid genomes of five Korean endemic Saussurea (S. albifolia, S. calcicola, S. diamantica, S. grandicapitula, and S. seoulensis), and comparative analyses including two other endemics (S. chabyoungsanica and S. polylepis) were conducted.

RESULTS

The plastomes of Korean endemics were highly conserved in gene content, order, and numbers. Exceptionally, S. diamantica had mitochondrial DNA sequences including two tRNAs in SSC region. There were no significant differences of the type and numbers of SSRs among the seven Korean endemics except in S. seoulensis. Nine mutation hotspots with high nucleotide diversity value (Pi > 0.0033) were identified, and phylogenetic analysis suggested that those Korean endemic species most likely evolved several times from diverse lineages within the genus. Moreover, molecular dating estimated that the Korean endemic species diverged since the late Miocene.

CONCLUSIONS

This study provides insight into understanding the plastome evolution and evolutionary relationships of highly complex species of Saussurea in Korean peninsula.

Review of the relationship and underlying mechanisms between the Qinghai-Tibet plateau and host intestinal flora

The intestinal microbial community is the largest ecosystem in the human body, in which the intestinal flora plays a dominant role and has a wide range of biological functions. However, it is vulnerable to a variety of factors, and exposure to extreme environments at high altitudes, as seen on the Qinghai-Tibet plateau, may cause changes in the structure and function of the host intestinal flora. Conversely, the intestinal flora can help the host adapt to the plateau environment through a variety of ways. Herein, we review the relationship and underlying mechanism between the host intestinal flora and the plateau environment by discussing the characteristics of the plateau environment, its influence on the intestinal flora, and the important role of the intestinal flora in host adaptation to the plateau environment. This review aimed to provide a reference for maintaining the health of the plateau population.

Clear phylogeographical structures shed light on the origin and dispersal of the aquatic boreal plant Hippuris vulgaris

Aquatic plants are an important ecological group in the arctic flora; however, their evolutionary histories remain largely unknown. In order to deepen our understanding of the evolution of these plants, we explored the phylogeographical structure of an aquatic boreal plant Hippuris vulgaris in a broad geographical sampling from Eurasia and North America using the chloroplast intergenic spacer psbA-trnH and seven nuclear microsatellite loci. Two closely-related species H. lanceolata and H. tetraphylla were also included because of their taxonomic controversy. Both chloroplast DNA sequences and nuclear microsatellite data revealed three genetic lineages with distinct distribution ranges. Incongruence between nuclear and chloroplast DNA lineages occurred in 14 samples from Russian Far East and Europe caused by inter-lineage hybridization. No private haplotypes or independent genetic clusters were evident in H. lanceolata or H. tetraphylla, suggesting that these two species should be considered conspecific ecotypes of H. vulgaris. Analysis using Approximate Bayesian Computation-Random Forest approach suggests that Hippuris vulgaris originated in China, followed by dispersal into Russia plus Northeast China, then successively westwards into Europe and North America, and finally into the Russian Far East from both North America and Russia plus Northeast China. This study is the first to elucidate the historical dispersal processes of a circumarctic aquatic plant across the entirety of its range.

Lactobacillus johnsonii YH1136 plays a protective role against endogenous pathogenic bacteria induced intestinal dysfunction by reconstructing gut microbiota in mice exposed at high altitude

Background

Intestinal microbiota plays an important role in maintaining the microecological balance of the gastrointestinal tract in various animals. Disturbances in the intestinal microbiota may lead to the proliferation of potentially pathogenic bacteria that become the dominant species, leading to intestinal immune disorders, intestinal inflammation, and other intestinal diseases. Numerous studies have been confirmed that high-altitude exposure affects the normal function of the intestine and the composition of the intestinal microbiota. However, it is still necessary to reveal the changes in intestinal microbiota in high-altitude exposure environments, and clarify the relationship between the proliferation of potentially pathogenic bacteria and intestinal injury in this environment. In addition, explored probiotics that may have preventive effects against intestinal diseases.

Methods and results

C57BL/6 mice were randomly divided into three groups, a high-altitude group (HA), control group (C), and high-altitude probiotic group (HAP). The HA and HAP groups were subjected to hypoxia modeling for 14 days in a low-pressure oxygen chamber with daily gavage of 0.2 mL of normal saline (HA) and Lactobacillus johnsonii YH1136 bacterial fluid (HAP), while the control group was fed normally. L. johnsonii YH1136 was isolated from feces of a healthy Tibetan girl in Baingoin county, the Nagqu region of the Tibet Autonomous Region, at an altitude of 5000 meters. Our observations revealed that gavage of YH1136 was effective in improving the damage to the intestinal barrier caused by high-altitude exposure to hypoxic environments and helped to reduce the likelihood of pathogenic bacteria infection through the intestinal barrier. It also positively regulates the intestinal microbiota to the extent of Lactobacillus being the dominant microbiome and reducing the number of pathogenic bacteria. By analyzing the expression profile of ileal microRNAs and correlation analysis with intestinal microbiota, we found that Staphylococcus and Corynebacterium1 cooperated with miR-196a-1-3p and miR-3060-3p, respectively, to play a regulatory role in the process of high-altitude hypoxia-induced intestinal injury.

Conclusion

These findings revealed the beneficial effect of L. johnsonii YH1136 in preventing potential endogenous pathogenic bacteria-induced intestinal dysfunction in high-altitude environments. The mechanism may be related to the regulation of intestinal injury from the perspective of the gut microbiota as well as miRNAs.

Plastome phylogenomics and biogeography of the subfam. Polygonoideae (Polygonaceae)

Polygonaceae has a complex taxonomic history, although a few studies using plastid or nuclear DNA fragments have explored relationships within this family, intrafamilial relationships remain controversial. Here, we newly sequenced and annotated 17 plastomes representing 12 genera within Polygonaceae. Combined with previously published data, a total of 49 plastomes representing 22/46 Polygonaceae genera and 16/20 Polygonoideae genera were collected to infer the phylogeny of Polygonaceae, with an emphasis on Polygonoideae. Plastome comparisons revealed high conservation within Polygonoideae in structure and gene order. Phylogenetic analyses using both Maximum Likelihood and Bayesian methods revealed two major clades and seven tribes within Polygonoideae. BEAST and S-DIVA analyses suggested a Paleocene origin of Polygonoideae in Asia. While most genera of Polygonoideae originated and further diversified in Asia, a few genera experienced multiple long-distance dispersal events from Eurasia to North America after the Miocene, with a few dispersal events to the Southern Hemisphere also being detected. Both ancient vicariance and long-distance events have played important roles in shaping the current distribution pattern of Polygonoideae.

Comparative and Phylogenetic Analyses of Complete Chloroplast Genomes of Scrophularia incisa Complex (Scrophulariaceae)

The Scrophularia incisa complex is a group of closely related desert and steppe subshrubs that includes S. incisa, S. kiriloviana and S. dentata, which are the only S. sect. Caninae components found in Northwest China. Based on earlier molecular evidence, the species boundaries and phylogenetic relationships within this complex remain poorly resolved. Here, we characterized seven complete chloroplast genomes encompassing the representatives of the three taxa in the complex and one closely related species, S. integrifolia, as well as three other species of Scrophularia. Comparative genomic analyses indicated that the genomic structure, gene order and content were highly conserved among these eleven plastomes. Highly variable plastid regions and simple sequence repeats (SSRs) were identified. The robust and consistent phylogenetic relationships of the S. incisa complex were firstly constructed based on a total of 26 plastid genomes from Scrophulariaceae. Within the monophyletic complex, a S. kiriloviana individual from Pamirs Plateau was identified as the earliest diverging clade, followed by S. dentata from Tibet, while the remaining individuals of S. kiriloviana from the Tianshan Mountains and S. incisa from Qinghai-Gansu were clustered into sister clades. Our results evidently demonstrate the capability of plastid genomes to improve phylogenetic resolution and species delimitation, particularly among closely related species, and will promote the understanding of plastome evolution in Scrophularia.

Adaptive responses drive the success of polyploid yellowcresses (Rorippa, Brassicaceae) in the Hengduan Mountains, a temperate biodiversity hotspot

Polyploids contribute substantially to plant evolution and biodiversity; however, the mechanisms by which they succeed are still unclear. According to the polyploid adaptation hypothesis, successful polyploids spread by repeated adaptive responses to new environments. Here, we tested this hypothesis using two tetraploid yellowcresses (Rorippa), the endemic Rorippa elata and the widespread Rorippa palustris, in the temperate biodiversity hotspot of the Hengduan Mountains. Speciation modes were resolved by phylogenetic modeling using 12 low-copy nuclear loci. Phylogeographical patterns were then examined using haplotypes phased from four plastid and ITS markers, coupled with historical niche reconstruction by ecological niche modeling. We inferred the time of hybrid origins for both species as the mid-Pleistocene, with shared glacial refugia within the southern Hengduan Mountains. Phylogeographic and ecological niche reconstruction indicated recurrent northward colonization by both species after speciation, possibly tracking denuded habitats created by glacial retreat during interglacial periods. Common garden experiment involving perennial R. elata conducted over two years revealed significant changes in fitness-related traits across source latitudes or altitudes, including latitudinal increases in survival rate and compactness of plant architecture, suggesting gradual adaptation during range expansion. These findings support the polyploid adaptation hypothesis and suggest that the spread of polyploids was aided by adaptive responses to environmental changes during the Pleistocene. Our results thus provide insight into the evolutionary success of polyploids in high-altitude environments.

Population genomics reveal deep divergence and strong geographical structure in gentians in the Hengduan Mountains

Understanding the evolutionary and ecological processes driving population differentiation and speciation can provide critical insights into the formation of biodiversity. Here, we examine the link between population genetic processes and biogeographic history underlying the generation of diversity in the Hengduan Mountains (HM), a region harboring a rich and dynamic flora. We used restriction site-associated DNA sequencing to generate 1,907 single-nucleotide polymorphisms (SNPs) and four-kb of plastid sequence in species of the Gentiana hexaphylla complex (Gentianaceae). We performed genetic clustering with spatial and non-spatial models, phylogenetic reconstructions, and ancestral range estimation, with the aim of addressing the processes influencing diversification of G. hexaphylla in the HM. We find the G. hexaphylla complex is characterized by geographic genetic structure with clusters corresponding to the South, North and the central HM. Phylogenetic reconstruction and pairwise F _ST analyses showed deep differentiation between Southern and Northern populations in the HM. The population in Mount Taibai exhibited the highest genetic similarity to the North HM. Ancestral range estimation indicated that the G. hexaphylla complex originated in the central HM and then diverged in the Pliocene and the Early Pleistocene, before dispersing widely, resulting in the current distinct lineages. Overall, we found deep genomic differentiation in the G. hexaphylla complex corresponds to geographic barriers to dispersal in the HM and highlights a critical role of the uplift of the Daxue Mountains and subsequent climatic fluctuations underlying diversification. The colonization of G. hexaphylla in the Mount Taibai region suggests directional dispersal between the alpine flora of the Qinling Mountains and the HM.

Plant Evolution History Overwhelms Current Environment Gradients in Affecting Leaf Chlorophyll Across the Tibetan Plateau

AIMS

Leaf chlorophyll (Chl) is a fundamental component and good proxy for plant photosynthesis. However, we know little about the large-scale patterns of leaf Chl and the relative roles of current environment changes vs. plant evolution in driving leaf Chl variations.

LOCATIONS

The east to west grassland transect of the Tibetan Plateau.

METHODS

We performed a grassland transect over 1,600 km across the Tibetan Plateau, measuring leaf Chl among 677 site-species.

RESULTS

Leaf Chl showed a significantly spatial pattern across the grasslands in the Tibetan Plateau, decreasing with latitude but increasing with longitude. Along with environmental gradient, leaf Chl decreased with photosynthetically active radiation (PAR), but increased with water availability and soil nitrogen availability. Furthermore, leaf Chl also showed significant differences among functional groups (C₄ > C₃ species; legumes < non-legume species), but no difference between annual and perennial species. However, we surprisingly found that plant evolution played a dominant role in shaping leaf Chl variations when comparing the sum and individual effects of all the environmental factors above. Moreover, we revealed that leaf Chl non-linearly decreased with plant evolutionary divergence time. This well-matches the non-linearly increasing trend in PAR or decreasing trend in temperature during the geological time-scale uplift of the Tibetan Plateau.

MAIN CONCLUSION

This study highlights the dominant role of plant evolution in determining leaf Chl variations across the Tibetan Plateau. Given the fundamental role of Chl for photosynthesis, these results provide new insights into reconsidering photosynthesis capacity in alpine plants and the carbon cycle in an evolutionary view.

Phylogeographical Pattern and Population Evolution History of Indigenous Elymus sibiricus L. on Qinghai-Tibetan Plateau

Elymus sibiricus L. is a perennial allotetraploid belonging to Triticeae of Poaceae, Elymus L., as the type species of genus Elymus L. The existing geographical distribution pattern and genetic spatial structure of E. sibiricus on Qinghai-Tibetan Plateau (QTP) are not yet clear. In this study, population genetic structure and demography history of 216 individuals from 44 E. sibiricus populations on QTP were studied used specific-locus amplified fragment sequencing (SLAF-seq). The result of genetic diversity showed that there was no single genetic diversity center was observed across all E. sibiricus populations. The results of genetic variation showed that 44 populations were clearly divided into the following three groups: Qinghai Plateau (Group I), South Tibet (Group II), and Hengduan Mountains (Group III). From the three analyses of AMOVA, Mantel test and Treemix, strong genetic differentiation across all populations and low genetic differentiation among populations within three groups. Molecular dating indicated that E. sibiricus diverged at 16.08 Ma (during the early Miocene) can be linked to the Himalayan Motion stage of QTP uplift. It is speculated that the reasons affecting the current phylogeographical pattern are as follows: (1) The environmental changes due to the uplift of the QTP; (2) The geographic distance between the populations (Groups I and III are close in geographic distance, and gene flow are frequent); (3) Geographical barriers (the Tanggula and Bayangela Mountains between Groups I and II). This study provides new evidence and historical perspective to the future exploration of the evolution and geographic distribution pattern of Elymus L.

Population and Landscape Genetics Provide Insights Into Species Conservation of Two Evergreen Oaks in Qinghai-Tibet Plateau and Adjacent Regions

The combination of population and landscape genetics can facilitate the understanding of conservation strategy under the changing climate. Here, we focused on the two most diverse and ecologically important evergreen oaks: Quercus aquifolioides and Quercus spinosa in Qinghai-Tibetan Plateau (QTP), which is considered as world's biodiversity hotspot. We genotyped 1,657 individuals of 106 populations at 15 nuclear microsatellite loci throughout the species distribution range. Spatial patterns of genetic diversity were identified by mapping the allelic richness (AR) and locally common alleles (LCA) according to the circular neighborhood methodology. Migration routes from QTP were detected by historical gene flow estimation. The response pattern of genetic variation to environmental gradient was assessed by the genotype-environment association (GEA) analysis. The overall genetic structure showed a high level of intra-species genetic divergence of a strong west-east pattern. The West-to-East migration route indicated the complex demographic history of two oak species. We found evidence of isolation by the environment in Q. aqu-East and Q. spi-West lineage but not in Q. aqu-West and Q. spi-East lineage. Furthermore, priority for conservation should be given to populations that retain higher spatial genetic diversity or isolated at the edge of the distribution range. Our findings indicate that knowledge of spatial diversity and migration route can provide valuable information for the conservation of existing populations. This study provides an important guide for species conservation for two oak species by the integration of population and landscape genetic methods.

Extensive Sympatry and Frequent Hybridization of Ecologically Divergent Aquatic Plants on the Qinghai-Tibetan Plateau

Hybridization has fascinated biologists in recent centuries for its evolutionary importance, especially in plants. Hybrid zones are commonly located in regions across environmental gradients due to more opportunities to contact and ecological heterogeneity. For aquatic taxa, intrazonal character makes broad overlapping regions in intermediate environments between related species. However, we have limited information on the hybridization pattern of aquatic taxa in alpines, especially submerged macrophytes. In this study, we aimed to test the hypotheses that niche overlap and hybridization might be extensive in related aquatic plants across an altitudinal gradient. We evaluated the niche overlap in three related species pairs on the Qinghai-Tibetan Plateau and assessed the spatial pattern of hybrid populations. Obvious niche overlap and common hybridization were revealed in all three pairs of related aquatic plants. The plateau edge and river basins were broad areas for the sympatry of divergent taxa, where a large proportion of hybrid populations occurred. Hybrids are also discretely distributed in diverse habitats on the plateau. Differences in the extent of niche overlap, genetic incompatibility and phylogeographic history might lead to variation differences in hybridization patterns among the three species pairs. Our results suggested that plateau areas are a hotspot for ecologically divergent aquatic species to contact and mate and implied that hybridization may be important for the freshwater biodiversity of highlands.