Search

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

Parnassia L., a perennial herbaceous genus in the family Celastraceae, consists of about 60 species and is mainly distributed in the Pan-Himalayan and surrounding mountainous regions. The taxonomic position and phylogenetic relationships of the genus are still controversial. Herein, we reassessed the taxonomic status of Parnassia and its intra- and inter-generic phylogeny within Celastraceae. To that end, we sequenced and assembled the whole plastid genomes and nuclear ribosomal DNA (nrDNA) of 48 species (74 individuals), including 25 species of Parnassia and 23 species from other genera of Celastraceae. We integrated high throughput sequence data with advanced statistical toolkits and performed the analyses. Our results supported the Angiosperm Phylogeny Group IV (APG IV) taxonomy which kept the genus to the family Celastraceae. Although there were topological conflicts between plastid and nrDNA phylogenetic trees, Parnassia was fully supported as a monophyletic group in all cases. We presented a first attempt to estimate the divergence of Parnassia, and molecular clock analysis indicated that the diversification occurred during the Eocene. The molecular phylogenetic results confirmed numerous taxonomic revisions, revealing that the morphological characters used in Parnassia taxonomy and systematics might have evolved multiple times. In addition, we speculated that hybridization/introgression might exist during genus evolution, which needs to be further studied. Similarly, more in-depth studies will clarify the diversification of characters and species evolution models of this genus.

The chloroplast (cp) genome is useful in the study of phylogenomics, molecular dating, and molecular evolution. Gentiana sect. Kudoa is a predominantly alpine flowering plant that is valued for its contributions to medicine, ecology, and horticulture. Previous evolutionary studies showed that the plastid gene loss pattern and intra-sectional phylogenetics in sect. Kudoa are still unclear. In this study, we compared 11 Gentiana plastomes, including 7 newly sequenced plastomes from sect. Kudoa, to represent its three serious: ser. Ornatae, ser. Verticillatae, and ser. Monanthae. The cp genome sizes of the seven species ranged from 137,278 to 147,156 bp. The plastome size variation mainly occurred in the small single-copy and long single-copy regions rather than the inverted repeat regions. Compared with sect. Cruciata, the plastomes in ser. Ornatae and ser. Verticillatae had lost approximately 11 kb of sequences containing 11 ndh genes. Conversely, far fewer losses were observed in ser. Monanthae. The phylogenetic tree revealed that sect. Kudoa was not monophyletic and that ser. Monanthae was more closely related to other sections rather than sect. Kudoa. The molecular dating analysis indicated that ser. Monanthae and sect. Kudoa diverged around 8.23 Ma. In ser. Ornatae and ser. Verticillatae, the divergence occurred at around 0.07–1.78 Ma. The nucleotide diversity analysis indicated that the intergenic regions trnH-psbA, trnK-trnQ, ycf3-trnS and rpl32-trnL constituted divergence hotspots in both sect. Kudoa and Gentiana, and would be useful for future phylogenetic and population genetic studies.

The varying topography and environment that resulted from paleoorogeny and climate fluctuations of the Himalaya–Hengduan Mountains (HHM) areas had a considerable impact on the evolution of biota during the Quaternary. To understand the phylogeographic pattern and historical dynamics of Triosteum himalayanum (Caprifoliaceae), we sequenced three chloroplast DNA fragments (rbcL-accD, rps15-ycf1, and trnH-psbA) from 238 individuals representing 20 populations. Nineteen haplotypes (H1–H19) were identified based on 23 single-site mutations and eight indels. Most haplotypes were restricted to a single population or neighboring populations. Analysis of molecular variance revealed that variations among populations were much higher than that within populations for the overall gene pool, as well as for the East Himalayan group (EH group) and the North Hengduan group (NHM group), but not for the Hengduan Mountains group (HM group). Ecoregions representing relatively high genetic diversity or high frequencies of private haplotypes were discovered, suggesting that this alpine herbaceous plant underwent enhanced allopatric divergence in isolated and fragmented locations during the Quaternary glaciations. The current phylogeographic structure of T. himalayanum might be due to heterogeneous habitats and Quaternary climatic oscillations. Based on the phylogeographic structure of T. himalayanum populations, the phylogenetic relationship of identified haplotypes and palaeodistributional reconstruction, we postulated both westwards and northwards expansion from the HM group for this species. The westwards dispersal corridor could be long, narrow mountain areas and/or the Yarlung Zangbo Valley, while the northwards movement path could be south–north oriented mountains and low-elevation valleys.

The effects of rapid, recent uplift of the Hengduan Mountains on evolution and diversification of young floristic lineages still remain unclear. Here, we investigate diversification of three closely related Saxifraga species with a distribution restricted to the Hengduan Mountains (HM) and southern Tibet, and comment on their taxonomy based on molecular evidence. Three chloroplast DNA fragments (rbcL, trnL-F, trnS-G) and the nuclear ribosomal DNA internal transcribed spacer (ITS) were employed to study genetic structure across 104 individuals from 12 populations of Saxifraga umbellulata, S. pasumensis, and S. banmaensis. Chloroplast DNA (cpDNA) phylogenies revealed two well supported clades, corresponding to S. umbellulata and S. pasumensis plus S. banmaensis. Topology of the ITS phylogeny was largely congruent with that generated from cpDNA haplotypes, but with minor conflicts which might be caused by incomplete lineage sorting. Analyses of molecular variance of both cpDNA and ITS datasets revealed that most variation was held between S. pasumensis s.l. (with S. banmaensis) and S. umbellulata (92.31% for cpDNA; 69.78% for ITS), suggesting a high degree of genetic divergence between them. Molecular clock analysis based on ITS dataset suggested that the divergence between S. pasumensis s.l. and S. umbellulata can be dated to 8.50 Ma, probably a result of vicariant allopatric diversification associated with the uplift events of the HM. Vicariance associated with HM uplifts may also have been responsible for infraspecific differentiation in S. pasumensis. In contrast, infraspecific differentiation in S. umbellulata was most likely triggered by Quaternary glaciations. The much lower levels of gene diversity within populations of S. pasumensis compared with S. umbellulata could have resulted from both range contractions and human collection on account of its putative medicinal properties. Combining evidence from morphology, geographical distributions and molecular phylogenetic data, we recommend that S. banmaensis should be treated as a synonym of S. pasumensis which in turn, and based on the same sources of evidence, should be treated as a separate species rather than as a variety of S. umbellulata.

The effects of rapid, recent uplift of the Hengduan Mountains on evolution and diversification of young floristic lineages still remain unclear. Here, we investigate diversification of three closely related Saxifraga species with a distribution restricted to the Hengduan Mountains (HM) and southern Tibet, and comment on their taxonomy based on molecular evidence. Three chloroplast DNA fragments (rbcL, trnL-F, trnS-G) and the nuclear ribosomal DNA internal transcribed spacer (ITS) were employed to study genetic structure across 104 individuals from 12 populations of Saxifraga umbellulata, S. pasumensis, and S. banmaensis. Chloroplast DNA (cpDNA) phylogenies revealed two well supported clades, corresponding to S. umbellulata and S. pasumensis plus S. banmaensis. Topology of the ITS phylogeny was largely congruent with that generated from cpDNA haplotypes, but with minor conflicts which might be caused by incomplete lineage sorting. Analyses of molecular variance of both cpDNA and ITS datasets revealed that most variation was held between S. pasumensis s.l. (with S. banmaensis) and S. umbellulata (92.31% for cpDNA; 69.78% for ITS), suggesting a high degree of genetic divergence between them. Molecular clock analysis based on ITS dataset suggested that the divergence between S. pasumensis s.l. and S. umbellulata can be dated to 8.50 Ma, probably a result of vicariant allopatric diversification associated with the uplift events of the HM. Vicariance associated with HM uplifts may also have been responsible for infraspecific differentiation in S. pasumensis. In contrast, infraspecific differentiation in S. umbellulata was most likely triggered by Quaternary glaciations. The much lower levels of gene diversity within populations of S. pasumensis compared with S. umbellulata could have resulted from both range contractions and human collection on account of its putative medicinal properties. Combining evidence from morphology, geographical distributions and molecular phylogenetic data, we recommend that S. banmaensis should be treated as a synonym of S. pasumensis which in turn, and based on the same sources of evidence, should be treated as a separate species rather than as a variety of S. umbellulata.

Studying closely related species and divergent populations provides insight into the process of speciation. Previous studies showed that the Sibiraea complex's evolutionary history on the Qinghai-Tibetan Plateau (QTP) was confusing and could not be distinguishable on the molecular level. In this study, the genetic structure and gene flow of Sibiraea laevigata and Sibiraea angustata on the QTP was examined across 45 populations using 8 microsatellite loci. Microsatellites revealed high genetic diversity in Sibiraea populations. Most of the variance was detected within populations (87.45%) rather than between species (4.39%). We found no significant correlations between genetic and geographical distances among populations. Bayesian cluster analysis grouped all individuals in the sympatric area of Sibiraea into one cluster and other individuals of S. angustata into another. Divergence history analysis based on the approximate Bayesian computation method indicated that the populations of S. angustata at the sympatric area derived from the admixture of the 2 species. The assignment test assigned all individuals to populations of their own species rather than its congeneric species. Consistently, intraspecies were detected rather than interspecies first-generation migrants. The bidirectional gene flow in long-term patterns between the 2 species was asymmetric, with more from S. angustata to S. laevigata. In conclusion, the Sibiraea complex was distinguishable on the molecular level using microsatellite loci. We found that the high genetic similarity of the complex resulted from huge bidirectional gene flow, especially on the sympatric area where population admixtures occurred. This study sheds light on speciation with gene flow in the QTP.

Studying closely related species and divergent populations provides insight into the process of speciation. Previous studies showed that the Sibiraea complex's evolutionary history on the Qinghai-Tibetan Plateau (QTP) was confusing and could not be distinguishable on the molecular level. In this study, the genetic structure and gene flow of Sibiraea laevigata and Sibiraea angustata on the QTP was examined across 45 populations using 8 microsatellite loci. Microsatellites revealed high genetic diversity in Sibiraea populations. Most of the variance was detected within populations (87.45%) rather than between species (4.39%). We found no significant correlations between genetic and geographical distances among populations. Bayesian cluster analysis grouped all individuals in the sympatric area of Sibiraea into one cluster and other individuals of S. angustata into another. Divergence history analysis based on the approximate Bayesian computation method indicated that the populations of S. angustata at the sympatric area derived from the admixture of the 2 species. The assignment test assigned all individuals to populations of their own species rather than its congeneric species. Consistently, intraspecies were detected rather than interspecies first-generation migrants. The bidirectional gene flow in long-term patterns between the 2 species was asymmetric, with more from S. angustata to S. laevigata. In conclusion, the Sibiraea complex was distinguishable on the molecular level using microsatellite loci. We found that the high genetic similarity of the complex resulted from huge bidirectional gene flow, especially on the sympatric area where population admixtures occurred. This study sheds light on speciation with gene flow in the QTP.

Swertia subuniflora is described from Fujian Province, China. Morphological and molecular evidences indicate its close relationship with S. franchetiana, S. mussotii and S. punicea, three species from Swertia series Japonica (further grouped in section Swertopsis, subgenus Ophelia). The new species can be recognized by its solitary flowers and ovate-cordate stem leaves.