Speciation history of a species complex of Primulina eburnea (Gesneriaceae) from limestone karsts of southern China, a biodiversity hot spot

Hiptage yangshuoensis (Malpighiaceae), a new species on karst hills close to Lijiang River, Guangxi, China, based on molecular and morphological data

Hiptage yangshuoensis K.Tan & K.S.Nguyen, a new species of Hiptage collected from a karst cliff close to the Lijiang River, Northeast of Guangxi Zhuang Autonomous Region, China, is described and illustrated based on molecular and morphological data. Hiptage yangshuoensis shares some morphological similarities with the H. multiflora F.N.Wei, but easily distinguished by its long pedicels with articulate at top, one large calyx gland, oblanceolate middle wing and lanceolate lateral wings of samara, and young branch covered rusty sericeous. The new species status is also supported by molecular phylogenetic analyses based on nuclear ribosome internal transcribed spacer (nrITS), which showed distinct systematic distinctiveness from the most morphologically similar species, H. multiflora.

Comprehensive Comparative Analyses of Aspidistra Chloroplast Genomes: Insights into Interspecific Plastid Diversity and Phylogeny

Limestone karsts are renowned for extremely high species richness and endemism. Aspidistra (Asparagaceae) is among the highly diversified genera distributed in karst areas, making it an ideal group for studying the evolutionary mechanisms of karst plants. The taxonomy and identification of Aspidistra species are mainly based on their specialized and diverse floral structures. Aspidistra plants have inconspicuous flowers, and the similarity in vegetative morphology often leads to difficulties in species discrimination. Chloroplast genomes possess variable genetic information and offer the potential for interspecies identification. However, as yet there is little information about the interspecific diversity and evolution of the plastid genomes of Aspidistra. In this study, we reported chloroplast (cp) genomes of seven Aspidistra species (A. crassifila, A. dolichanthera, A. erecta, A. longgangensis, A. minutiflora, A. nankunshanensis, and A. retusa). These seven highly-conserved plastid genomes all have a typical quartile structure and include a total of 113 unique genes, comprising 79 protein-coding genes, 4 rRNA genes and 30 tRNA genes. Additionally, we conducted a comprehensive comparative analysis of Aspidistra cp genomes. We identified eight divergent hotspot regions (trnC-GCA-petN, trnE-UUC-psbD, accD-psaI, petA-psbJ, rpl20-rps12, rpl36-rps8, ccsA-ndhD and rps15-ycf1) that serve as potential molecular markers. Our newly generated Aspidistra plastomes enrich the resources of plastid genomes of karst plants, and an investigation into the plastome diversity offers novel perspectives on the taxonomy, phylogeny and evolution of Aspidistra species.

Genome-wide analysis of the cellulose toolbox of Primulina eburnea, a calcium-rich vegetable

BACKGROUND

Human-guided crop domestication has lasted for more than 10,000 years. In terms of the domestication and breeding of vegetables, cellulose content in edible tissues is one of the most important traits. Primulina eburnea is a recently developed calcium-rich vegetable with a high soluble and bioavailable calcium content in its leaves. However, the high cellulose content in the leaves hampers the taste, and no research has been reported on the genetic basis of cellulose biosynthesis in this calcium-rich vegetable.

RESULTS

We identified 36 cellulose biosynthesis-involved genes belonging to eight gene families in the P. eburnea genome. The cellulose accumulated decreasingly throughout leaf development. Nineteen genes were considered core genes in cellulose biosynthesis, which were highly expressed in buds but lowly expressed in mature leaves. In the nitrogen fertilization experiment, exogenous nitrogen decreased the cellulose content in the buds. The expressing pattern of 14 genes were consistent with phenotypic variation in the nitrogen fertilization experiment, and thus they were proposed as cellulose toolbox genes.

CONCLUSIONS

The present study provides a strong basis for the subsequent functional research of cellulose biosynthesis-involved genes in P. eburnea, and provides a reference for breeding and/or engineering this calcium-rich vegetable with decreased leaf cellulose content to improve the taste.

Genetic and ecophysiological evidence that hybridization facilitated lineage diversification in yellow Camellia (Theaceae) species: a case study of natural hybridization between C. micrantha and C. flavida

BACKGROUND

Hybridization is generally considered an important creative evolutionary force, yet this evolutionary process is still poorly characterized in karst plants. In this study, we focus on natural hybridization in yellow Camellia species, a group of habitat specialists confined to karst/non-karst habitats in southwestern China.

RESULTS

Based on population genome data obtain from double digest restriction-site associated DNA (ddRAD) sequencing, we found evidence for natural hybridization and introgression between C. micrantha and C. flavida, and specifically confirmed their hybrid population, C. "ptilosperma". Ecophysiological results suggested that extreme hydraulic traits were fixed in C. "ptilosperma", these being consistent with its distinct ecological niche, which lies outside its parental ranges.

CONCLUSION

The identified hybridization event is expected to have played a role in generating novel variation during, in which the hybrid population displays different phenological characteristics and novel ecophysiological traits associated with the colonization of a new niche in limestone karst.

Comparative transcriptome analyses of three Gentiana species provides signals for the molecular footprints of selection effects and the phylogenetic relationships

Species in Gentiana section Cruciata are important alpine plants with a center of diversity and speciation in Qinghai-Tibet Plateau (QTP), and some of these species are sympatrically distributed in northeastern QTP. Studies on genome features and natural selection signatures of sympatric species in section Crucata have been impeded by a lack of genomic resources. Here, we showed transcript characterizations and molecular footprints of selection effects on G. straminea, G. dahurica and G. officinalis based on the comparative transcriptome. A total of 62.97 Gb clean reads were obtained with unigene numbers per species ranging from 141,819 to 236,408 after assembly. We found that these three species had similar distribution of functional categories in different databases, and key enzyme-encoding genes involved in the iridoids biosynthesis were also obtained. The selective pressure analyses indicated that most paired orthologs between these three species were subject to negative selection, and only a low proportion of the orthologs that underwent positive selection were detected. We found that some positive selected genes were involved in "catalytic activity", "metabolic process", "response to stimulus" and "response to stress". Besides, large numbers of SSR primer pairs with transferabilities were successfully designed based on the available transcriptome datasets of three Gentiana species. The phylogenetic relationships reconstructed based on 352 single-copy nuclear genes provided a rough phylogenetic framework for this genus and confirmed the monophyly of section Cruciata. Our study not only provides insights for the natural selection effects on sympatric Gentiana species, but also enhances future genetic breeding or evolutionary studies on Qinjiao species.

Phylogeographic structure of Heteroplexis (Asteraceae), an endangered endemic genus in the limestone karst regions of southern China

Though the karst regions in south and southwest China are plant diversity hotspots, our understanding of the phylogeography and evolutionary history of the plants there remains limited. The genus Heteroplexis (Asteraceae) is one of the typical representative plants isolated by karst habitat islands, and is also an endangered and endemic plant to China. In this study, species-level phylogeographic analysis of the genus Heteroplexis was conducted using restriction site-associated DNA sequencing (RADseq). The genetic structure showed a clear phylogeographic structure consistent with the current species boundaries in the H. microcephala, H. incana, H. vernonioides, H. sericophylla, and H. impressinervia. The significant global (R = 0.37, P < 0.01) and regional (R = 0.650.95, P < 0.05) isolation by distance (IBD) signals among species indicate strong geographic isolation in the karst mountains, which may result in chronically restricted gene flow and increased genetic drift and differentiation. Furthermore, the phylogeographic structure of Heteroplexis suggested a southward migration since the last glacial period. Demographic analysis revealed the karst mountains as a refuge for Heteroplexis species. Finally, both Treemix and ABBA-BABA statistic detected significant historical gene flow between species. Significant historical gene flow and long-term stability of effective population size (Ne) together explain the high genome-wide genetic diversity among species (π = 0.05370.0838). However, the recent collapse of Ne, widespread inbreeding within populations, and restricted contemporary gene flow suggest that Heteroplexis species are probably facing a high risk of genetic diversity loss. Our results help to understand the evolutionary history of karst plants and guide conservation.

Plastomes of limestone karst gesneriad genera Petrocodon and Primulina, and the comparative plastid phylogenomics of Gesneriaceae

Petrocodon and Primulina are two characteristic genera of Gesneriaceae that exhibit remarkable species and floral diversity, and high endemism across the Sino-Vietnamese Limestone Karsts. To better understand the evolution of limestone gesneriad plastomes, we report nine complete plastomes of seven Primulina and two Petrocodon which have never been assembled before. The newly generated plastomes range from 152,323 to 153,786 bp in size and display a typical quadripartite structure. To further explore the plastome evolution across Gesneriaceae, we assembled five additional plastomes from public reads data and incorporated 38 complete Gesneriaceae plastomes available online into comparative and phylogenomic analyses. The comparison of 52 Gesneriaceae plastomes reveals that not only Primulina and Petrocodon but all gesneriad genera analyzed are highly conserved in genome size, genome structure, gene contents, IR boundary configurations, and codon usage bias. Additionally, sliding window analyses were implemented across alignments of Primulina and Petrocodon for identifying highly variable regions, providing informative markers for future studies. Meanwhile, the SSRs and long repeats of Gesneriaceae plastomes were characterized, serving as useful data in studying population and repetitive sequence evolutions. The results of plastome phylogenetics represent a preliminary but highly resolved maternal backbone genealogy of Primulina and the Old World subtribes of Gesneriaceae.

Reticulate evolution in the Pteris fauriei group (Pteridaceae)

The Pteris fauriei group (Pteridaceae) has a wide distribution in Eastern Asia and includes 18 species with similar but varied morphology. We collected more than 300 specimens of the P. fauriei group and determined ploidy by flow cytometry and inferred phylogenies by molecular analyses of chloroplast and nuclear DNA markers. Our results reveal a complicated reticulate evolution, consisting of seven parental taxa and 58 hybrids. The large number of hybrid taxa have added significant morphological complexity to the group leading to difficult taxonomic issues. The hybrids generally had broader ranges and more populations than their parental taxa. Genetic combination of different pairs of parental species created divergent phenotypes of hybrids, exhibited by both morphological characteristics and ecological fidelities. Niche novelty could facilitate hybrid speciation. Apogamy is common in this group and potentially contributes to the sustainability of the whole group. We propose that frequent hybridizations among members of the P. fauriei group generate and maintain genetic diversity, via novel genetic combinations, niche differentiation, and apogamy.

Gene flow, linked selection, and divergent sorting of ancient polymorphism shape genomic divergence landscape in a group of edaphic specialists

Interpreting the formation of genomic variation landscape, especially genomic regions with elevated differentiation (i.e. islands), is fundamental to a better understanding of the genomic consequences of adaptation and speciation. Edaphic islands provide excellent systems for understanding the interplay of gene flow and selection in driving population divergence and speciation. However, discerning the relative contribution of these factors that modify patterns of genomic variation remains difficult. We analysed 132 genomes from five recently divergent species in Primulina genus, with four species distributed in Karst limestone habitats and the fifth one growing in Danxia habitats. We demonstrated that both gene flow and linked selection have contributed to genome-wide variation landscape, where genomic regions with elevated differentiation (i.e., islands) were largely derived by divergent sorting of ancient polymorphism. Specifically, we identified several lineage-specific genomic islands that might have facilitated adaptation of P. suichuanensis to Danxia habitats. Our study is amongst the first cases disentangling evolutionary processes that shape genomic variation of plant specialists, and demonstrates the important role of ancient polymorphism in the formation of genomic islands that potentially mediate adaptation and speciation of endemic plants in special soil habitats.

The identification of an R2R3-MYB transcription factor involved in regulating anthocyanin biosynthesis in Primulina swinglei flowers

Primulina genus is an ideal wild ornamental flower and emerging model for studying biosynthesis, diversity, and evolution of flower pigment. However, the molecular mechanism underlying anthocyanin biosynthesis and regulation in Primulina remains unknown. Here, changes in anthocyanin content and the expression profiles of anthocyanin biosynthetic structural genes were examined in developing Primulina swinglei flowers and three other organs. Seventy-three R2R3-MYB transcription factor genes were identified from transcriptome of P. swinglei flowers, two of which, PsMYB1 and PsMYB2, are candidate regulators of anthocyanin biosynthesis according to clustering analysis. Furthermore, transient over-expression studies using tobacco leaves showed distinct pigment accumulation following co-infection with PsMYB1 and MrbHLH1 (a previously confirmed anthocyanin regulator from Morella rubra). Additionally, dual luciferase assays showed that PsMYB1 trans-activated the PsANS promoter, with the addition of MrbHLH1 resulting in a 5-fold increase in the intensity of this interaction. PsMYB1 did not, however, have any effect on the PsF3H promoter. The expression profile and dual luciferase assays showed that PsMYB2 plays no roles in anthocyanin regulation. Therefore, PsMYB1 is proposed to be the transcription factor gene regulating anthocyanin biosynthesis in P. swinglei.

The genome of a cave plant, Primulina huaijiensis, provides insights into adaptation to limestone karst habitats

Although whole genome duplication (WGD) has been suggested to facilitate adaptive evolution and diversification, the role of specific WGD events in promoting diversification and adaptation in angiosperms remains poorly understood. Primulina, a species-rich genus with > 180 species associated with limestone karst habitat, constitutes an ideal system for studying the impact of WGD events on speciation and evolutionary adaptation. We sequenced and assembled a chromosome-level genome of the cave-dwelling species P. huaijiensis to study gene family expansion and gene retention following WGDs. We provide evidence that P. huaijiensis has undergone two WGDs since the γ triplication event shared by all eudicots. In addition to a WGD shared by almost all Lamiales (L event), we identified a lineage-specific WGD (D event) that occurred in the early Miocene around 20.6-24.2 Myr ago and that is shared by almost the entire subtribe Didymocarpinae. We found that gene retentions following the D event led to gene family proliferation (e.g. WRKYs) that probably facilitated adaptation to the high salinity and drought stress in limestone karst. Our study highlights the role of lineage-specific WGD in species diversification and adaptation of plants from special habitats.

Population Differentiation and Demographic History of the Cycas taiwaniana Complex (Cycadaceae) Endemic to South China as Indicated by DNA Sequences and Microsatellite Markers

Historical geology, climatic oscillations, and seed dispersal capabilities are thought to influence the population dynamics and genetics of plants, especially for distribution-restricted and threatened species. Investigating the genetic resources within and among taxa is a prerequisite for conservation management. The Cycas taiwaniana complex consists of six endangered species that are endemic to South China. In this study, we investigated the relationship between phylogeographic history and the genetic structure of the C. taiwaniana complex. To estimate the phylogeographic history of the complex, we assessed the genetic structure and divergence time, and performed phylogenetic and demographic historical analyses. Two chloroplast DNA intergenic regions (cpDNA), two single-copy nuclear genes (SCNGs), and six microsatellite loci (SSR) were sequenced for 18 populations. The SCNG data indicated a high genetic diversity within populations, a low genetic diversity among populations, and significant genetic differentiation among populations. Significant phylogeographical structure was detected. Structure and phylogenetic analyses both revealed that the 18 populations of the C. taiwaniana complex have two main lineages, which were estimated to diverge in the Middle Pleistocene. We propose that Cycas fairylakea was incorporated into Cycas szechuanensis and that the other populations, which are mainly located on Hainan Island, merged into one lineage. Bayesian skyline plot analyses revealed that the C. taiwaniana complex experienced a recent decline, suggesting that the complex probably experienced a bottleneck event. We infer that the genetic structure of the C. taiwaniana complex has been affected by Pleistocene climate shifts, sea-level oscillations, and human activities. In addition to providing new insights into the evolutionary legacy of the genus, the genetic characterizations will be useful for the conservation of Cycas species.

Primulina serrulata (Gesneriaceae), a new species from southeastern Guizhou, China

Primulina serrulata R.B.Zhang & F. Wen, a new species from a limestone area in southeastern Guizhou, China, is described and illustrated here. The new species is morphologically related to P. fimbrisepala (Hand.-Mazz.) Y.Z.Wang. We examined the morphological differences between these congeners and provide illustrations and photographs of this new species in this paper.

Molecular systematics of Chiritopsis-like Primulina (Gesneriaceae): one new species, one new name, two new combinations, and new synonyms

BACKGROUND

The Gesneriaceae genus Chiritopsis, confined almost exclusively to cave or cave-like microhabitats of limestone karsts of southern China, was described to distinguish it from Chirita by much smaller flowers and generally miniature plant sizes in the former genus. However, molecular phylogenetic analyses showed that Chiritopsis is polyphyletic and its species delimitation has been problematic. To understand how many times Chiritopsis-like species have evolved from within the recircumscribed Primulina and to further clarify their species identification, we sampled all but two recently described species of Chiritopsis-like Primulina and reconstructed their phylogenetic relationship based on DNA sequences of nuclear ITS and chloroplast trnL-F and trnH-psbA.

RESULTS

With 182 accessions of 165 taxa of Primulina sampled, our analyses placed the 40 accessions of 25 taxa of Chiritopsis-like Primulina in 17 unrelated positions, indicating at least 17 independent origins of the traits associated with caves or cave-like microhabitats. Of the 17 clades containing Chiritopsis-like Primulina, Clade 1 is composed of P. bipinnatifida, P. cangwuensis, P. jianghuaensis, P. lingchuanensis, and P. zhoui, as well as additional samples that show variable and overlapping morphology in leaf shapes. Clade 10 includes P. cordifolia, P. huangii, and P. repanda, while Primulina repanda var. guilinensis is not placed within Clade 10. Primulina glandulosa var. yangshuoensis is not placed in the same clade of P. glandulosa.

CONCLUSIONS

Based on our data, P. cangwuensis, P. jianghuaensis, and P. lingchuanensis are proposed to synonymize under P. bipinnatifida, with P. zhoui treated as a variety of P. bipinnatifida. Primulina repanda var. guilinensis is transferred as P. subulata var. guilinensis comb. nov. and Primulina pseudoglandulosa nom. nov. is proposed for P. glandulosa var. yangshuoensis. One new species is named P. chingipengii to honor the late Dr. Ching-I Peng (1950-2018).

Primulina malipoensis (Gesneriaceae), a new species from Sino-Vietnamese border area

Primulina malipoensis, a new species from limestone areas around the Sino-Vietnamese border, is described and illustrated. This new species is morphologically similar to P. maguanensis and P. lungzhouensis, but obviously differs from the latter two species by its pale greenish-yellow flowers (vs. purple, with different colour patterns). The phylogenetic affinity, illustration and photographs of this new species are provided in this paper.

Genomic Signature of Adaptive Divergence despite Strong Nonadaptive Forces on Edaphic Islands: A Case Study of Primulina juliae

Both genetic drift and divergent selection are expected to be strong evolutionary forces driving population differentiation on edaphic habitat islands. However, the relative contribution of genetic drift and divergent selection to population divergence has rarely been tested simultaneously. In this study, restriction-site associated DNA-based population genomic analyses were applied to assess the relative importance of drift and divergent selection on population divergence of Primulina juliae, an edaphic specialist from southern China. All populations were found with low standing genetic variation, small effective population size (NE), and signatures of bottlenecks. Populations with the lowest genetic variation were most genetically differentiated from other populations and the extent of genetic drift increased with geographic distance from other populations. Together with evidence of isolation by distance, these results support neutral drift as a critical evolutionary driver. Nonetheless, redundancy analysis revealed that genomic variation is significantly associated with both edaphic habitats and climatic factors independently of spatial effects. Moreover, more genomic variation was explained by environmental factors than by geographic variables, suggesting that local adaptation might have played an important role in driving population divergence. Finally, outlier tests and environment association analyses identified 31 single-nucleotide polymorphisms as candidates for adaptive divergence. Among these candidates, 26 single-nucleotide polymorphisms occur in/near genes that potentially play a role in adaptation to edaphic specialization. This study has important implications that improve our understanding of the joint roles of genetic drift and adaptation in generating population divergence and diversity of edaphic specialists.

Speciation history of a species complex of Primulina eburnea (Gesneriaceae) from limestone karsts of southern China, a biodiversity hot spot

Limestone karsts in southern China are characterized by high edaphic and topographic heterogeneity and host high levels of species richness and endemism. However, the evolutionary mechanisms for generating such biodiversity remain poorly understood. Here, we performed species delimitation, population genetic analyses, simulations of gene flow and analyses of floral morphological traits to infer the geographic history of speciation in a species complex of Primulina eburnea from limestone karsts of southern China. Using Bayesian species delimitation, we determined that there are seven distinct species that correspond well to the putative morphological species. Species tree reconstruction, Structure and Neighbor‐Net analyses all recovered four lineages in agreement with currently species geographic boundaries. High levels of genetic differentiation were observed both within and among species. Isolation–migration coalescent analysis provides evidence for significant but low gene flow among species. Approximate Bayesian computation (ABC) analysis supports a scenario of historical gene flow rather than recent contemporary gene flow for most species divergences. Finally, we found no evidence of divergent selection contributing to population differentiation of a suite of flower traits. These results support the prevalence of allopatric speciation and highlight the role of geographic isolation in the diversification process. At small geographic scales, limited hybridization occurred in the past between proximate populations but did not eliminate species boundaries. We conclude that limited gene flow might have been the predominant evolutionary force in promoting population differentiation and speciation.