Both temperature fluctuations and East Asian monsoons have driven plant diversification in the karst ecosystems from southern China

Genomic divergence and mutation load in the Begonia masoniana complex from limestone karsts

Understanding genome-wide diversity, inbreeding, and the burden of accumulated deleterious mutations in small and isolated populations is essential for predicting and enhancing population persistence and resilience. However, these effects are rarely studied in limestone karst plants. Here, we re-sequenced the nuclear genomes of 62 individuals of the Begonia masoniana complex (B. liuyanii, B. longgangensis, B. masoniana and B. variegata) and investigated genomic divergence and genetic load for these four species. Our analyses revealed four distinct clusters corresponding to each species within the complex. Notably, there was only limited admixture between B. liuyanii and B. longgangensis occurring in overlapping geographic regions. All species experienced historical bottlenecks during the Pleistocene, which were likely caused by glacial climate fluctuations. We detected an asymmetric historical gene flow between group pairs within this timeframe, highlighting a distinctive pattern of interspecific divergence attributable to karst geographic isolation. We found that isolated populations of B. masoniana have limited gene flow, the smallest recent population size, the highest inbreeding coefficients, and the greatest accumulation of recessive deleterious mutations. These findings underscore the urgency to prioritize conservation efforts for these isolated population. This study is among the first to disentangle the genetic differentiation and specific demographic history of karst Begonia plants at the whole-genome level, shedding light on the potential risks associated with the accumulation of deleterious mutations over generations of inbreeding. Moreover, our findings may facilitate conservation planning by providing critical baseline genetic data and a better understanding of the historical events that have shaped current population structure of rare and endangered karst plants.

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

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

Plastid phylogenomics provides new insights into the systematics, diversification, and biogeography of Cymbidium (Orchidaceae)

Cymbidium (Orchidaceae: Epidendroideae), with around 60 species, is widely-distributed across Southeast Asia, providing a nice system for studying the processes that underlie patterns of biodiversity in the region. However, phylogenetic relationships of Cymbidium have not been well resolved, hampering investigations of species diversification and the biogeographical history of this genus. In this study, we construct a plastome phylogeny of 56 Cymbidium species, with four well-resolved major clades, which provides a framework for biogeographical and diversification rate analyses. Molecular dating and biogeographical analyses show that Cymbidium likely originated in the region spanning northern Indo-Burma to the eastern Himalayas during the early Miocene (∼21.10 Ma). It then rapidly diversified into four major clades in East Asia within approximately a million years during the middle Miocene. Cymbidium spp. migration to the adjacent regions (Borneo, Philippines, and Sulawesi) primarily occurred during the Pliocene-Pleistocene period. Our analyses indicate that the net diversification rate of Cymbidium has decreased since its origin, and is positively associated with changes in temperature and monsoon intensity. Favorable hydrothermal conditions brought by monsoon intensification in the early Miocene possibly contributed to the initial rapid diversification, after which the net diversification rate was reduced with the cooling climate after the middle Miocene. The transition from epiphytic to terrestrial habits may have enabled adaptation to cooler environments and colonization of northern niches, yet without a significant effect on diversification rates. This study provides new insights into how monsoon activity and temperature changes affected the diversification dynamics of plants in Southeast Asia.

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

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

Population differentiation and dynamics of five pioneer species of Gaultheria from the secondary forests in subtropical China

BACKGROUND

The influence of native secondary succession associated with anthropogenic disturbance on the biodiversity of the forests in subtropical China remains uncertain. In particular, the evolutionary response of small understory shrubs, particularly pioneer species inhabiting continuously disturbed habitats, to topographic heterogeneity and climate change is poorly understood. This study aimed to address this knowledge gap by focusing on the Gaultheria crenulata group, a clade of small pioneer shrubs in subtropical China.

RESULTS

We examined the genetic structure and demographic history of all five species of the G. crenulata group with two maternally inherited chloroplast DNA (cpDNA) fragments and two biparentally inherited low-copy nuclear genes (LCG) over 89 natural populations. We found that the genetic differentiation of this group was influenced by the geomorphological boundary between different regions of China in association with Quaternary climatic events. Despite low overall genetic diversity, we observed an isolation-by-distance (IBD) pattern at a regional scale, rather than isolation-by-environment (IBE), which was attributed to ongoing human disturbance in the region.

CONCLUSION

Our findings suggest that the genetic structure of the G. crenulata group reflects the interplay of geological topography, historical climates, and anthropogenic disturbance during the Pliocene-Pleistocene-Holocene periods in subtropical China. The observed IBD pattern, particularly prominent in western China, highlights the role of limited dispersal and gene flow, possibly influenced by physical barriers or decreased connectivity over geographic distance. Furthermore, the east-to-west trend of gene flow, potentially facilitated by the East Asian monsoon system, underscores the complex interplay of biotic and abiotic factors shaping the genetic dynamics of pioneer species in subtropical China's secondary forests. These findings can be used to assess the impact of environmental changes on the adaptation and persistence of biodiversity in subtropical forest ecosystems.

Eco-evolutionary evidence for the global diversity pattern of Cycas (Cycadaceae)

The evolution of the latitudinal diversity gradient (LDG), characterized by a peak in diversity toward the tropics, has captured significant attention in evolutionary biology and ecology. However, the inverse LDG (i-LDG) mechanism, wherein species richness increases toward the poles, remains inadequately explored. Cycads are among one of the oldest lineages of extant seed plants and have undergone extensive diversification in the tropics. Intriguingly, the extant cycad abundance exhibits an i-LDG pattern, and the underlying causes for this phenomenon remain largely elusive. Here, using 1,843 nuclear genes from a nearly complete sampling, we conducted comprehensive phylogenomic analyses to establish a robust species-level phylogeny for Cycas, the largest genus within cycads. We then reconstructed the spatial-temporal dynamics and integrated global environmental data to evaluate the roles of species ages, diversification rates, contemporary environment, and conservatism to ancestral niches in shaping the i-LDG pattern. We found Cycas experienced decreased diversification rates, coupled with the cooling temperature since its origin in the Eocene from continental Asia. Different regions have distinctively contributed to the formation of i-LDG for Cycas, with the northern hemisphere acting as evolutionary museums and the southern hemisphere serving as cradles. Moreover, water-related climate variables, specifically precipitation seasonality and potential evapotranspiration, were identified as paramount factors constraining Cycas species richness in the rainforest biome near the equator. Notably, the adherence to ancestral monsoonal climates emerges as a critical factor in sustaining the diversity pattern. This study underscores the imperative of integrating both evolutionary and ecological approaches to comprehensively unravel the mechanisms underpinning global biodiversity patterns.

Climatic and biogeographic processes underlying the diversification of the pantropical flowering plant family Annonaceae

Tropical forests harbor the richest biodiversity among terrestrial ecosystems, but few studies have addressed the underlying processes of species diversification in these ecosystems. We use the pantropical flowering plant family Annonaceae as a study system to investigate how climate and biogeographic events contribute to diversification. A super-matrix phylogeny comprising 835 taxa (34% of Annonaceae species) based on eight chloroplast regions was used in this study. We show that global temperature may better explain the recent rapid diversification in Annonaceae than time and constant models. Accelerated accumulation of niche divergence (around 15 Ma) lags behind the increase of diversification rate (around 25 Ma), reflecting a heterogeneous transition to recent diversity increases. Biogeographic events are related to only two of the five diversification rate shifts detected. Shifts in niche evolution nevertheless appear to be associated with increasingly seasonal environments. Our results do not support the direct correlation of any particular climatic niche shifts or historical biogeographical event with shifts in diversification rate. Instead, we suggest that Annonaceae diversification can lead to later niche divergence as a result of increasing interspecific competition arising from species accumulation. Shifts in niche evolution appear to be associated with increasingly seasonal environments. Our results highlight the complexity of diversification in taxa with long evolutionary histories.

Phylogeny and diversification of genus Sanicula L. (Apiaceae): novel insights from plastid phylogenomic analyses

BACKGROUND

The genus Sanicula L. is a unique perennial herb that holds important medicinal values. Although the previous studies on Sanicula provided us with a good research basis, its taxonomic system and interspecific relationships have not been satisfactorily resolved, especially for those endemic to China. Moreover, the evolutionary history of this genus also remains inadequately understood. The plastid genomes possessing highly conserved structure and limited evolutionary rate have proved to be an effective tool for studying plant phylogeny and evolution.

RESULTS

In the current study, we newly sequenced and assembled fifteen Sanicula complete plastomes. Combined with two previously reported plastomes, we performed comprehensively plastid phylogenomics analyses to gain novel insights into the evolutionary history of this genus. The comparative results indicated that the seventeen plastomes exhibited a high degree of conservation and similarity in terms of their structure, size, GC content, gene order, IR borders, codon bias patterns and SSRs profiles. Such as all of them displayed a typical quadripartite structure, including a large single copy region (LSC: 85,074-86,197 bp), a small single copy region (SSC: 17,047-17,132 bp) separated by a pair of inverted repeat regions (IRs: 26,176-26,334 bp). And the seventeen plastomes had similar IR boundaries and the adjacent genes were identical. The rps19 gene was located at the junction of the LSC/IRa, the IRa/SSC junction region was located between the trnN gene and ndhF gene, the ycf1 gene appeared in the SSC/IRb junction and the IRb/LSC boundary was located between rpl12 gene and trnH gene. Twelve specific mutation hotspots (atpF, cemA, accD, rpl22, rbcL, matK, ycf1, trnH-psbA, ycf4-cemA, rbcL-accD, trnE-trnT and trnG-trnR) were identified that can serve as potential DNA barcodes for species identification within the genus Sanicula. Furthermore, the plastomes data and Internal Transcribed Spacer (ITS) sequences were performed to reconstruct the phylogeny of Sanicula. Although the tree topologies of them were incongruent, both provided strong evidence supporting the monophyly of Saniculoideae and Apioideae. In addition, the sister groups between Saniculoideae and Apioideae were strongly suggested. The Sanicula species involved in this study were clustered into a clade, and the Eryngium species were also clustered together. However, it was clearly observed that the sections of Sanicula involved in the current study were not respectively recovered as monophyletic group. Molecular dating analysis explored that the origin of this genus was occurred during the late Eocene period, approximately 37.84 Ma (95% HPD: 20.33-52.21 Ma) years ago and the diversification of the genus was occurred in early Miocene 18.38 Ma (95% HPD: 10.68-25.28 Ma).

CONCLUSION

The plastome-based tree and ITS-based tree generated incongruences, which may be attributed to the event of hybridization/introgression, incomplete lineage sorting (ILS) and chloroplast capture. Our study highlighted the power of plastome data to significantly improve the phylogenetic supports and resolutions, and to efficiently explore the evolutionary history of this genus. Molecular dating analysis explored that the diversification of the genus occurred in the early Miocene, which was largely influenced by the prevalence of the East Asian monsoon and the uplift of the Hengduan Mountains (HDM). In summary, our study provides novel insights into the plastome evolution, phylogenetic relationships, taxonomic framework and evolution of genus Sanicula.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

Historical development of karst evergreen broadleaved forests in East Asia has shaped the evolution of a hemiparasitic genus Brandisia (Orobanchaceae)

Brandisia is a shrubby genus of about eight species distributed basically in East Asian evergreen broadleaved forests (EBLFs), with distribution centers in the karst regions of Yunnan, Guizhou, and Guangxi in southwestern China. Based on the hemiparasitic and more or less liana habits of this genus, we hypothesized that its evolution and distribution were shaped by the development of EBLFs there. To test our hypothesis, the most comprehensive phylogenies of Brandisia hitherto were constructed based on plastome and nuclear loci (nrDNA, PHYA and PHYB); then divergence time and ancestral areas were inferred using the combined nuclear loci dataset. Phylogenetic analyses reconfirmed that Brandisia is a member of Orobanchaceae, with unstable placements caused by nuclear-plastid incongruences. Within Brandisia, three major clades were well supported, corresponding to the three subgenera based on morphology. Brandisia was inferred to have originated in the early Oligocene (32.69 Mya) in the Eastern Himalayas-SW China, followed by diversification in the early Miocene (19.45 Mya) in karst EBLFs. The differentiation dates of Brandisia were consistent with the origin of keystone species of EBLFs in this region (e.g., Fagaceae, Lauraceae, Theaceae, and Magnoliaceae) and the colonization of other characteristic groups (e.g., Gesneriaceae and Mahonia). These findings indicate that the distribution and evolution of Brandisia were facilitated by the rise of the karst EBLFs in East Asia. In addition, the woody and parasitic habits, and pollination characteristics of Brandisia may also be the important factors affecting its speciation and dispersal.

Genomic data provide a robust phylogeny backbone for Rhodiola L. (Crassulaceae) and reveal extensive reticulate evolution during its rapid radiation

The Tibetan Plateau and adjacent mountain regions (TP; including the Tibetan Plateau, Himalaya, Hengduan Mountains and Mountains of Central Asia) harbor great biodiversity, some lineages on which may have undergone rapid radiations. However, only a few studies have investigated the evolutionary pattern of such diversification in depth using genomic data. In this study, we reconstructed a robust phylogeny backbone of Rhodiola, a lineage that may have undergone rapid radiation in the TP, using Genotyping-by-sequencing data, and conducted a series of gene flow and diversification analyses. The concatenation and coalescent-based methods yield similar tree topologies, and five well-supported clades were revealed. Potential gene flow and introgression events were detected, both between species from different major clades and closely related species, suggesting pervasive hybridization and introgression. An initial rapid and later slowdown of the diversification rate was revealed, indicating niche filling. Molecular dating and correlation analyses showed that the uplift of TP and global cooling in the mid-Miocene might have played an important role in promoting the rapid radiation of Rhodiola. Our work demonstrates that gene flow and introgression might be an important contributor to rapid radiation possibly by quickly reassembling old genetic variation into new combinations.

Macroevolutionary diversity of traits and genomes in the model yeast genus Saccharomyces

Species is the fundamental unit to quantify biodiversity. In recent years, the model yeast Saccharomyces cerevisiae has seen an increased number of studies related to its geographical distribution, population structure, and phenotypic diversity. However, seven additional species from the same genus have been less thoroughly studied, which has limited our understanding of the macroevolutionary events leading to the diversification of this genus over the last 20 million years. Here, we show the geographies, hosts, substrates, and phylogenetic relationships for approximately 1,800 Saccharomyces strains, covering the complete genus with unprecedented breadth and depth. We generated and analyzed complete genome sequences of 163 strains and phenotyped 128 phylogenetically diverse strains. This dataset provides insights about genetic and phenotypic diversity within and between species and populations, quantifies reticulation and incomplete lineage sorting, and demonstrates how gene flow and selection have affected traits, such as galactose metabolism. These findings elevate the genus Saccharomyces as a model to understand biodiversity and evolution in microbial eukaryotes.

Langbiangia, a new genus of Gesneriaceae endemic to Langbiang Plateau, southern Vietnam and a taxonomic endeavor to achieve key targets of the post-2020 global biodiversity framework

Situated in the southern end of the Annamite Mountain Range, Langbiang Plateau is a major biodiversity hotspot of southern Vietnam known for high species diversity and endemicity. To achieve effective conservation, parts of the plateau were designated as the Langbiang Biosphere Reserve, an UNESCO World Network aiming to improve relationships between inhabitants and their environments. Amongst the rich endemic flora of the plateau are three gesneriads ascribed to Primulina, a calciphilous genus with high species diversity in the vast limestone karsts stretching from southern China to northern Vietnam. However, a recent phylogenetic study questioned the generic placement of the Langbiang Primulina, corroborating with observations on the geographical distribution, habitat preference, and phyllotaxy of the three species. Based on phylogenetic analyses of nuclear ITS and plastid trnL-F DNA sequences of a comprehensive sampling covering nearly all genera of the Old World Gesneriaceae, we demonstrate that the three Langbiang Primulina species form a fully supported clade distantly related to other Primulina. As this clade is biogeographically, ecologically, morphologically, and phylogenetically distinct worthy of generic recognition, we propose to name it Langbiangia gen. nov. to highlight the rich and unique biodiversity of the Langbiang Plateau. By means of this taxonomic endeavor, we are hoping to raise the conservation awareness of this biodiversity heritage of southern Vietnam and promote the importance of Langbiang Biosphere Reserve that is crucial for achieving action-oriented global targets of the post-2020 global biodiversity framework (GBF) of the UN Convention on Biological Diversity (CBD)-effective conservation and management of at least 30% of biodiverse terrestrial, inland water, and costal and marine areas by 2030-that has been agreed at the COP15 in Montréal in December 2022.

Plastid phylogenomics sheds light on divergence time and ecological adaptations of the tribe Persicarieae (Polygonaceae)

Southwestern China, adjacent to the Qinghai-Tibetan Plateau (QTP), is known as a hotspot for plant diversity and endemism, and it is the origin and diversification center of Persicarieae. As one of the major lineages in Polygonaceae, Persicarieae represents a diverse adaptation to various habitats. As a result of morphological plasticity and poorly resolving molecular markers, phylogenetic relationships and infrageneric classification within Persicarieae have long been controversial. In addition, neither plastome phylogenomic studies nor divergence time estimates on a larger sample of Persicarieae species have been made thus far. We sequenced and assembled 74 complete plastomes, including all of the recognized genera within Persicarieae and their relatives. We conducted a comprehensive phylogenetic study of the major clades within Persicarieae and, based on the thus obtained robust phylogeny, also estimated divergence time and the evolution of diagnostic morphological traits. Major relationships found in previous phylogenetic studies were confirmed, including those of the backbone of the tree, which had been a major problem in previous phylogenies of the tribe. Phylogenetic analysis revealed strong support for Koenigia as sister to Bistorta, and together they were sister to the robustly supported Persicaria. Based on the phylogenetic and morphological evidence, we recognize five sections in Persicaria: Persicaria, Amphibia, Tovara, Echinocaulon, and Cephalophilon. It is estimated that the divergence of the Persicarieae began around the late Paleocene, with diversification concentrated in the Eocene and Miocene. In addition, it is suggested that the increasing westerly and monsoon winds in conjunction with the uplift of the QTP may be the driving force for origin and diversification of Persicarieae species. These results provide a valuable evolutionary framework for the study of adaptation in Polygonaceae and insights into plant diversification on the QTP and adjacent areas.

Monsoon intensification in East Asia triggered the evolution of its flora

INTRODUCTION

East Asia (EA), which falls within the region of the Asian monsoon that is composed of the East Asia monsoon (EAM) and the Indian monsoon (IM), is known for its high species diversity and endemism. This has been attributed to extreme physiographical heterogeneity in conjunction with climate and sea-level changes during the Pleistocene, this hypothesis has been widely proven by phylogeographic studies. Recently, dated phylogenies have indicated that the origins (stem age) of the flora occurred after the Oligocene-Miocene boundary and are related to the establishment of the EAM.

METHODS

Hence, this study further examined whether the strengthening of the monsoons triggered floral evolution via a meta-analysis of the tempo-spatial pattern of evolutionary radiation dates (crown ages) of 101 endemic seed plant genera.

RESULTS

Taxonomic diversification began during the late Eocene, whereas the accumulated number of diversifications did not significantly accelerate until the late Miocene. The distribution of the weighted mean and the average divergence times in the EAM, IM, or transitional regions all fall within the mid-late Miocene. Fossils of the Tertiary relict genera are mostly and widely distributed outside EA and only half of the earliest fossils in the EA region are not older than Miocene, while their divergence times are mostly after the late Miocene. The pattern of divergence time of monotypic and polytypic taxa suggest the climatic changes after the late Pliocene exert more influence on monotypic taxa.

DISCUSSION

The two key stages of floral evolution coincide with the intensifications of the EAM and IM, especially the summer monsoon which brings a humid climate. An integrated review of previous studies concerning flora, genus, and species levels further supports our suggestion that monsoon intensification in EA triggered the evolution of its flora.

Genomic insights into inter- and intraspecific mating system shifts in Primulina

The mating system shift from outcrossing to selfing is one of the most frequent evolutionary trends in flowering plants. However, the genomic consequences of this shift remain poorly understood. Specifically, the relative importance of the demographic and genetic processes causing changes in genetic variation and selection efficacy associated with the evolution of selfing is unclear. Here we sequenced the genomes of two Primulina species with contrasting mating systems: P. eburnea (outcrossing) versus P. tabacum (outcrossing, mixed-mating and selfing populations). Whole-genome resequencing data were used to investigate the genomic consequences of mating system shifts within and between species. We found that highly selfing populations of P. tabacum display loss of genetic diversity, increased deleterious mutations, higher genomic burden and fewer adaptive substitutions. However, compared with outcrossing populations, mixed-mating populations did not display loss of genetic diversity and accumulation of genetic load. We find no evidence of population bottlenecks associated with the shift to selfing, which suggests that the genetic effects of selfing on N_e and possibly linked selection, rather than demographic history, are the primary drivers of diversity reduction in highly selfing populations. Our results highlight the importance of distinguishing the relative contribution of mating system and demography on the genomic consequences associated with mating system evolution in plants.

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.

Fingerprints of climatic changes through the late Cenozoic in southern Asian flora: Magnolia section Michelia (Magnoliaceae)

BACKGROUND AND AIMS

Ongoing global warming is a challenge for humankind. A series of drastic climatic changes have been proven to have occurred throughout the Cenozoic based on a variety of geological evidence, which helps to better understand our planet's future climate. Notably, extant biomes have recorded drastic environmental shifts. The climate in southern Asia, which hosts high biodiversity, is deeply impacted by the Asian monsoon. The origins and evolutionary dynamics of biomes occurring between the tropics and sub-tropics in southern Asia have probably been deeply impacted by climatic changes; however, these aspects remain poorly studied. We tested whether the evolutionary dynamics of the above biomes have recorded the drastic, late Cenozoic environmental shifts, by focusing on Magnolia section Michelia of the family Magnoliaceae.

METHODS

We established a fine time-calibrated phylogeny of M. section Michelia based on complete plastid genomes and inferred its ancestral ranges. Finally, we estimated the evolutionary dynamics of this section through time, determining its diversification rate and the dispersal events that occurred between tropical and sub-tropical areas.

KEY RESULTS

The tropical origin of M. section Michelia was dated to the late Oligocene; however, the diversification of its core group (i.e. M. section Michelia subsection Michelia) has occurred mainly from the late Miocene onward. Two key evolutionary shifts (dated approx. 8 and approx. 3 million years ago, respectively) were identified, each of them probably in response to drastic climatic changes.

CONCLUSION

Here, we inferred the underlying evolutionary dynamics of biomes in southern Asia, which probably reflect late Cenozoic climatic changes. The occurrence of modern Asian monsoons was probably fundamental for the origin of M. section Michelia; moreover, the occurrence of asymmetric dispersal events between the tropics and sub-tropics hint at an adaptation strategy of M. section Michelia to global cooling, in agreement with the tropical conservatism hypothesis.

Three new species of Bredia (Sonerileae, Melastomataceae) from the Sino-Vietnamese border area

Brediabullata, B.enchengensis, and B.nitida (Sonerileae, Melastomataceae), three species occurring in Sino-Vietnamese limestone karst regions, are described as new. Molecular phylogenetic analyses and morphological divergence indicate that these species are well separated from their close relatives in Bredia, justifying their recognition as distinct species. Brediabullata is unique in its interveinal areas prominently bullate each with an apical seta, a character otherwise never recorded in the genus. Bredianitida resembles B.malipoensis in habit, leaf shape, and inflorescence morphology, but differs in the glabrescent and nitid adaxial leaf surface (vs. densely pubescent and subvelvety), ovate-elliptic or elliptic calyx lobes (vs. triangular to semiorbicular), and white petals (vs. purplish-red). Brediaenchengensis is closest to B.longiradiosa, but easily recognized by its prostrate habit (vs. erect), the yellowish-green, membranous and fragile leaves (vs. green or dark green, papery), and white anthers (vs. pink to purplish). These new discoveries show that further botanical exploration is warranted in the remote Sino-Vietnamese bordering region.

Comparative plastome analysis of Musaceae and new insights into phylogenetic relationships

Background

Musaceae is an economically important family consisting of 70-80 species. Elucidation of the interspecific relationships of this family is essential for a more efficient conservation and utilization of genetic resources for banana improvement. However, the scarcity of herbarium specimens and quality molecular markers have limited our understanding of the phylogenetic relationships in wild species of Musaceae. Aiming at improving the phylogenetic resolution of Musaceae, we analyzed a comprehensive set of 49 plastomes for 48 species/subspecies representing all three genera of this family.

Results

Musaceae plastomes have a relatively well-conserved genomic size and gene content, with a full length ranging from 166,782 bp to 172,514 bp. Variations in the IR borders were found to show phylogenetic signals to a certain extent in Musa. Codon usage bias analysis showed different preferences for the same codon between species and three genera and a common preference for A/T-ending codons. Among the two genes detected under positive selection (dN/dS > 1), ycf2 was indicated under an intensive positive selection. The divergent hotspot analysis allowed the identification of four regions (ndhF-trnL, ndhF, matK-rps16, and accD) as specific DNA barcodes for Musaceae species.

Bayesian and maximum likelihood phylogenetic analyses using full plastome resulted in nearly identical tree topologies with highly supported relationships between species. The monospecies genus Musella is sister to Ensete, and the genus Musa was divided into two large clades, which corresponded well to the basic number of n = x = 11 and n = x =10/9/7, respectively. Four subclades were divided within the genus Musa. A dating analysis covering the whole Zingiberales indicated that the divergence of Musaceae family originated in the Palaeocene (59.19 Ma), and the genus Musa diverged into two clades in the Eocene (50.70 Ma) and then started to diversify from the late Oligocene (29.92 Ma) to the late Miocene. Two lineages (Rhodochlamys and Australimusa) radiated recently in the Pliocene /Pleistocene periods.

Conclusions

The plastome sequences performed well in resolving the phylogenetic relationships of Musaceae and generated new insights into its evolution. Plastome sequences provided valuable resources for population genetics and phylogenetics at lower taxon.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08454-3.

Plastome phylogenomics and historical biogeography of aquatic plant genus Hydrocharis (Hydrocharitaceae)

BACKGROUND

Hydrocharis L. and Limnobium Rich. are small aquatic genera, including three and two species, respectively. The taxonomic status, phylogenetic relationships and biogeographical history of these genera have remained unclear, owing to the lack of Central African endemic H. chevalieri from all previous studies. We sequenced and assembled plastomes of all three Hydrocharis species and Limnobium laevigatum to explore the phylogenetic and biogeographical history of these aquatic plants.

RESULTS

All four newly generated plastomes were conserved in genome structure, gene content, and gene order. However, they differed in size, the number of repeat sequences, and inverted repeat borders. Our phylogenomic analyses recovered non-monophyletic Hydrocharis. The African species H. chevalieri was fully supported as sister to the rest of the species, and L. laevigatum was nested in Hydrocharis as a sister to H. dubia. Hydrocharis-Limnobium initially diverged from the remaining genera at ca. 53.3 Ma, then began to diversify at ca. 30.9 Ma. The biogeographic analysis suggested that Hydrocharis probably originated in Europe and Central Africa.

CONCLUSION

Based on the phylogenetic results, morphological similarity and small size of the genera, the most reasonable taxonomic solution to the non-monophyly of Hydrocharis is to treat Limnobium as its synonym. The African endemic H. chevalieri is fully supported as a sister to the remaining species. Hydrocharis mainly diversified in the Miocene, during which rapid climate change may have contributed to the speciation and extinctions. The American species of former Limnobium probably dispersed to America through the Bering Land Bridge during the Miocene.

Diversification Slowdown in the Cirrhopetalum Alliance (Bulbophyllum, Orchidaceae): Insights From the Evolutionary Dynamics of Crassulacean Acid Metabolism

Evolutionary slowdowns in diversification have been inferred in various plant and animal lineages. Investigation based on diversification models integrated with environmental factors and key characters could provide critical insights into this diversification trend. We evaluate diversification rates in the Cirrhopetalum alliance (Bulbophyllum, Orchidaceae subfam. Epidendroideae) using a time-calibrated phylogeny and assess the role of Crassulacean acid metabolism (CAM) as a hypothesised key innovation promoting the spectacular diversity of orchids, especially those with an epiphytic habit. An explosive early speciation in the Cirrhopetalum alliance is evident, with the origin of CAM providing a short-term advantage under the low atmospheric CO2 concentrations (pCO2) associated with cooling and aridification in the late Miocene. A subsequent slowdown of diversification in the Cirrhopetalum alliance is possibly explained by a failure to keep pace with pCO2 dynamics. We further demonstrate that extinction rates in strong CAM lineages are ten times higher than those of C3 lineages, with CAM not as evolutionarily labile as previously assumed. These results challenge the role of CAM as a "key innovation" in the diversification of epiphytic orchids.

Immigration dynamics of tropical and subtropical Southeast Asian limestone karst floras

Ex situ origins and dispersal of taxa have played important roles in the assembly of island-like biodiversity hotspots. Insular limestone karsts in Southeast Asia are hotspots of biodiversity and endemism, but the immigration processes of their unique floras are still poorly known. Here, we used Gesneriaceae as a proxy to investigate the immigration dynamics of tropical and subtropical Southeast Asian karst floras. We present the most comprehensive phylogenetic analysis of the Old World gesneriads to date based on twelve loci. By estimating divergence times and reconstructing ancestral states (habitat, soil type and range), we found that immigration into subtropical Southeast Asian karst floras first occurred in the Early Miocene, with two peaks in the Early-Middle Miocene and the Pliocene-Early Pleistocene, whereas immigration into tropical Southeast Asian karsts initiated in the Late Eocene, with two peaks in the Late Oligocene and the Late Miocene. We also discover that Southeast Asian karst biodiversity comprises immigrant pre-adapted lineages and descendants from local acid soil ancestors, although niche shift from acid soil to karst in tropical Southeast Asian islands was lacking. This study advances our understanding of the historical assembly of Southeast Asian karst floras.

When tropical and subtropical congeners met: Multiple ancient hybridization events within Eriobotrya in the Yunnan-Guizhou Plateau, a tropical-subtropical transition area in China

Global climate changes during the Miocene may have created ample opportunities for hybridization between members of tropical and subtropical biomes at the boundary between these zones. Yet, very few studies have explored this possibility. The Yunnan-Guizhou Plateau (YGP) in Southwest China is a biodiversity hotspot for vascular plants, located in a transitional area between the floristic regions of tropical Southeast Asia and subtropical East Asia. The genus Eriobotrya (Rosaceae) comprises both tropical and subtropical taxa, with 12 species recorded in the YGP, making it a suitable basis for testing the hypothesis of between-biome hybridization. Therefore, we surveyed the evolutionary history of Eriobotrya by examining three chloroplast regions and five nuclear genes for 817 individuals (47 populations) of 23 Eriobotrya species (including 19 populations of 12 species in the YGP), plus genome re-sequencing of 33 representative samples. We concluded that: (1) phylogenetic positions for 16 species exhibited strong cytonuclear conflicts, most probably due to ancient hybridization; (2) the YGP is a hotspot for hybridization, with 11 species showing clear evidence of chloroplast capture; and (3) Eriobotrya probably originated in tropical Asia during the Eocene. From the Miocene onwards, the intensification of the Eastern Asia monsoon and global cooling may have shifted the tropical-subtropical boundary and caused secondary contact between species, thus providing ample opportunity for hybridization and diversification of Eriobotrya, especially in the YGP. Our study highlights the significant role that paleoclimate changes probably played in driving hybridization and generating rich species diversity in climate transition zones.

Population dynamics linked to glacial cycles in Cercis chuniana F. P. Metcalf (Fabaceae) endemic to the montane regions of subtropical China

The mountains of subtropical China are an excellent system for investigating the processes driving the geographical distribution of biodiversity and radiation of plant populations in response to Pleistocene climate fluctuations. How the major mountain ranges in subtropical China have affected the evolution of plant species in the subtropical evergreen broadleaved forest is an issue with long-term concern. Here, we focused on Cercis chuniana, a woody species endemic to the southern mountain ranges in subtropical China, to elucidate its population dynamics. We used genotyping by sequencing (GBS) to investigate the spatial pattern of genetic variation among 11 populations. Geographical isolation was detected between the populations located in adjacent mountain ranges, thought to function as geographical barriers due to their complex physiography. Bayesian time estimation revealed that population divergence occurred in the middle Pleistocene, when populations in the Nanling Mts. separated from those to the east. The orientation and physiography of the mountain ranges of subtropical China appear to have contributed to the geographical pattern of genetic variation between the eastern and western populations of C. chuniana. Complex physiography plus long-term stable ecological conditions across glacial cycles facilitated the demographic expansion in the Nanling Mts., from which contemporary migration began. The Nanling Mts. are thus considered as a suitable area for preserving population diversity and large population sizes of C. chuniana compared with other regions. As inferred by ecological niche modeling and coalescent simulations, secondary contact occurred during the warm Lushan-Tali Interglacial period, with intensified East Asia summer monsoon and continuous habitat available for occupation. Our data support the strong influence of both climatic history and topographic characteristics on the high regional phytodiversity of the subtropical evergreen broadleaved forest in subtropical China.

Strong Genetic Structure Observed in Primulina danxiaensis, a Small Herb Endemic to Mount Danxia With Extremely Small Populations

Danxia landform occurring sporadically in southern China is a unique type of petrographic geomorphology. It has nurtured about 400 rare or threatened plant and animal species, whose diversity, endemism, and conservation have called increasing scientific and public attentions. Among them, Primulina danxiaensis (W. B. Liao, S. S. Lin, and R. J. Shen) W. B. Liao and K. F. Chung is a tiny perennial grass species recorded only in Mount Danxia, a natural World Heritage Site as part of China’s Danxia. In this study, restriction site-associated DNA sequencing (RAD-seq) was performed to investigate genetic diversity among these 12 populations of P. danxiaensis. A total of 432,041 variant sites were detected in 84,779 loci across 94 samples. The expected heterozygosity (H_E) ranged from 0.017 to 0.139. Bottleneck signals were detected in most populations, Tajima’s D tests showed that most loci could be under recent positive selection, and one of the six positively selected loci identified by BayeScan was annotated as tRNA^Glu, which may contribute to the species’ adaptation to shady environment. STRUCTURE analysis and phylogenetic tree showed that the 12 populations of P. danxiaensis could be divided into four gene pools (clades) corresponding to their geographic locations, and significant correlation was observed between genetic and geographic distances. Our study demonstrated that P. danxiaensis maintained a middle level of genetic diversity and strong population structure; geographic distance could be an important factor limiting gene flow among populations of P. danxiaensis, which were only sporadically recorded in Mount Danxia.

Phylogenomic and Macroevolutionary Evidence for an Explosive Radiation of a Plant Genus in the Miocene

Mountain systems harbor a substantial fraction of global biodiversity and, thus, provide excellent opportunities to study rapid diversification and to understand the historical processes underlying the assembly of biodiversity hotspots. The rich biodiversity in mountains is widely regarded as having arisen under the influence of geological and climatic processes as well as the complex interactions among them. However, the relative contribution of geology and climate in driving species radiation is seldom explored. Here, we studied the evolutionary radiation of Oreocharis (Gesneriaceae), which has diversified extensively throughout East Asia, especially within the Hengduan Mountains (HDM), using transcriptomic data and a time calibrated phylogeny for 88% (111/126) of all species of the genus. In particular, we applied phylogenetic reconstructions to evaluate the extent of incomplete lineage sorting accompanying the early and rapid radiation in the genus. We then fit macroevolutionary models to explore its spatial and diversification dynamics in Oreocharis and applied explicit birth–death models to investigate the effects of past environmental changes on its diversification. Evidence from 574 orthologous loci suggest that Oreocharis underwent an impressive early burst of speciation starting ca. 12 Ma in the Miocene, followed by a drastic decline in speciation toward the present. Although we found no evidence for a shift in diversification rate across the phylogeny of Oreocharis, we showed a difference in diversification dynamics between the HDM and non-HDM lineages, with higher diversification rates in the HDM. The diversification dynamic of Oreocharis is most likely positively associated with temperature-dependent speciation and dependency on the Asian monsoons. We suggest that the warm and humid climate of the mid-Miocene was probably the primary driver of the rapid diversification in Oreocharis, while mountain building of the HDM might have indirectly affected species diversification of the HDM lineage. This study highlights the importance of past climatic changes, combined with mountain building, in creating strong environmental heterogeneity and driving diversification of mountain plants, and suggests that the biodiversity in the HDM cannot directly be attributed to mountain uplift, contrary to many recent speculations.[East Asian monsoons; environmental heterogeneity; Hengduan Mountains; incomplete lineage sorting; Oreocharis; past climate change; rapid diversification; transcriptome.]

Allopolyploid Speciation Accompanied by Gene Flow in a Tree Fern

Hybridization in plants may result in hybrid speciation or introgression and, thus, is now widely understood to be an important mechanism of species diversity on an evolutionary timescale. Hybridization is particularly common in ferns, as is polyploidy, which often results from hybrid crosses. Nevertheless, hybrid speciation as an evolutionary process in fern lineages remains poorly understood. Here, we employ flow cytometry, phylogeny, genomewide single nucleotide polymorphism data sets, and admixture and coalescent modeling to show that the scaly tree fern, Gymnosphaera metteniana is a naturally occurring allotetraploid species derived from hybridization between the diploids, G. denticulata and G. gigantea. Moreover, we detected ongoing gene flow between the hybrid species and its progenitors, and we found that G. gigantea and G. metteniana inhabit distinct niches, whereas climatic niches of G. denticulata and G. metteniana largely overlap. Taken together, these results suggest that either some degree of intrinsic genetic isolation between the hybrid species and its parental progenitors or ecological isolation over short distances may be playing an important role in the evolution of reproductive barriers. Historical climate change may have facilitated the origin of G. metteniana, with the timing of hybridization coinciding with a period of intensification of the East Asian monsoon during the Pliocene and Pleistocene periods in southern China. Our study of allotetraploid G. metteniana represents the first genomic-level documentation of hybrid speciation in scaly tree ferns and, thus, provides a new perspective on evolution in the lineage.

Phylogenetic signatures of ecological divergence and leapfrog adaptive radiation in Espeletia

PREMISE

Events of accelerated species diversification represent one of Earth's most celebrated evolutionary outcomes. Northern Andean high-elevation ecosystems, or páramos, host some plant lineages that have experienced the fastest diversification rates, likely triggered by ecological opportunities created by mountain uplifts, local climate shifts, and key trait innovations. However, the mechanisms behind rapid speciation into the new adaptive zone provided by these opportunities have long remained unclear.

METHODS

We address this issue by studying the Venezuelan clade of Espeletia, a species-rich group of páramo-endemics showing a dazzling ecological and morphological diversity. We performed several comparative analyses to study both lineage and trait diversification, using an updated molecular phylogeny of this plant group.

RESULTS

We showed that sets of either vegetative or reproductive traits have conjointly diversified in Espeletia along different vegetation belts, leading to adaptive syndromes. Diversification in vegetative traits occurred earlier than in reproductive ones. The rate of species and morphological diversification showed a tendency to slow down over time, probably due to diversity dependence. We also found that closely related species exhibit significantly more overlap in their geographic distributions than distantly related taxa, suggesting that most events of ecological divergence occurred at close geographic proximity within páramos.

CONCLUSIONS

These results provide compelling support for a scenario of small-scale ecological divergence along multiple ecological niche dimensions, possibly driven by competitive interactions between species, and acting sequentially over time in a leapfrog pattern.

Parallel ddRAD and Genome Skimming Analyses Reveal a Radiative and Reticulate Evolutionary History of the Temperate Bamboos

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.]

The East Asian Winter Monsoon Acts as a Major Selective Factor in the Intraspecific Differentiation of Drought-Tolerant Nitraria tangutorum in Northwest China

The influence of Quaternary climate fluctuation on the geographical structure and genetic diversity of species distributed in the regions of the Qinghai–Tibet Plateau (QTP) has been well established. However, the underlying role of the East Asian monsoon system (EAMS) in shaping the genetic structure of the population and the demography of plants located in the arid northwest of China has not been explored. In the present study, Nitraria tangutorum, a drought-tolerant desert shrub that is distributed in the EAMS zone and has substantial ecological and economic value, was profiled to better understand the influence of EAMS evolution on its biogeographical patterns and demographic history. Thus, the phylogeographical structure and historical dynamics of this plant species were elucidated using its five chloroplast DNA (cpDNA) fragments. Hierarchical structure analysis revealed three distinct, divergent lineages: West, East-A, and East-B. The molecular dating was carried out using a Bayesian approach to estimate the time of intraspecies divergence. Notably, the eastern region, which included East-A and East-B lineages, was revealed to be the original center of distribution and was characterized by a high level of genetic diversity, with the intraspecific divergence time dated to be around 2.53 million years ago (Ma). These findings, combined with the data obtained by ecological niche modeling analysis, indicated that the East lineages have undergone population expansion and differentiation, which were closely correlated with the development of the EAMS, especially the East Asian winter monsoon (EAWM). The West lineage appears to have originated from the migration of N. tangutorum across the Hexi corridor at around 1.85 Ma, and subsequent colonization of the western region. These results suggest that the EAWM accelerated the population expansion of N. tangutorum and subsequent intraspecific differentiation. These findings collectively provide new information on the impact of the evolution of the EAMS on intraspecific diversification and population demography of drought-tolerant plant species in northwest China.

Two new combinations in Oreocharis (Gesneriaceae) based on morphological, molecular and cytological evidence

The newly-circumscribed genus Oreocharis is recently enlarged by incorporating ten other genera with high floral diversity. In this study, our morphological, molecular and cytological evidence supports our adding two species from other two different genera (Boeica and Beccarinda) to Oreocharis. The special corolla shape (campanulate or flat-faced) and related short filament of these two new combinations, Oreocharis guileana and O. baolianis, further enrich the diversity of floral characters of the enlarged Oreocharis. Meanwhile, some supplementary and amended descriptions of these two species are made here. Our morphological, molecular and geographical data indicate that O. guileana is related to O. pilosopetiolata to a certain extent. For O. baolianis, however, our current dataset does not allow conclusions on the species relationship within Oreocharis.

Biogeography and evolution of Asian Gesneriaceae based on updated taxonomy

Based on an updated taxonomy of Gesneriaceae, the biogeography and evolution of the Asian Gesneriaceae are outlined and discussed. Most of the Asian Gesneriaceae belongs to Didymocarpoideae, except Titanotrichum was recently moved into Gesnerioideae. Most basal taxa of the Asian Gesneriaceae are found in the Indian subcontinent and Indo-China Peninsula, suggesting Didymocarpoideae might originate in these regions. Four species diversification centers were recognized, i.e. Sino-Vietnam regions, Malay Peninsula, North Borneo and Northwest Yunnan (Hengduan Mountains). The first three regions are dominated by limestone landscapes, while the Northwest Yunnan is well-known for its numerous deep gorges and high mountains. The places with at least 25% species are neoendemics (newly evolved and narrowly endemic) which were determined as evolutionary hotspots, including Hengduan Mountains, boundary areas of Yunnan-Guizhou-Guangxi in Southwest China, North Borneo, Pahang and Terengganu in Malay Peninsula, and mountainous areas in North Thailand, North Sulawesi Island. Finally, the underlying mechanisms for biogeographical patterns and species diversification of the Asian Gesneriaceae are discussed.

An ancient tropical origin, dispersals via land bridges and Miocene diversification explain the subcosmopolitan disjunctions of the liverwort genus Lejeunea

Understanding the biogeographical and diversification processes explaining current diversity patterns of subcosmopolitan-distributed groups is challenging. We aimed at disentangling the historical biogeography of the subcosmopolitan liverwort genus Lejeunea with estimation of ancestral areas of origin and testing if sexual system and palaeotemperature variations can be factors of diversification. We assembled a dense taxon sampling for 120 species sampled throughout the geographical distribution of the genus. Lejeunea diverged from its sister group after the Paleocene-Eocene boundary (52.2 Ma, 95% credibility intervals 50.1-54.2 Ma), and the initial diversification of the crown group occurred in the early to middle Eocene (44.5 Ma, 95% credibility intervals 38.5-50.8 Ma). The DEC model indicated that (1) Lejeunea likely originated in an area composed of the Neotropics and the Nearctic, (2) dispersals through terrestrial land bridges in the late Oligocene and Miocene allowed Lejeunea to colonize the Old World, (3) the Boreotropical forest covering the northern regions until the late Eocene did not facilitate Lejeunea dispersals, and (4) a single long-distance dispersal event was inferred between the Neotropics and Africa. Biogeographical and diversification analyses show the Miocene was an important period when Lejeunea diversified globally. We found slight support for higher diversification rates of species with both male and female reproductive organs on the same individual (monoicy), and a moderate positive influence of palaeotemperatures on diversification. Our study shows that an ancient origin associated with a dispersal history facilitated by terrestrial land bridges and not long-distance dispersals are likely to explain the subcosmopolitan distribution of Lejeunea. By enhancing the diversification rates, monoicy likely favoured the colonisations of new areas, especially in the Miocene that was a key epoch shaping the worldwide distribution.

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.

Phylogenomics of the aquatic plant genus Ottelia (Hydrocharitaceae): Implications for historical biogeography

Ottelia Pers. is the second largest genus of the family Hydrocharitaceae, including approximately 23 extant species. The genus exhibits a diversity of both bisexual and unisexual flowers, and complex reproductive system comprising cross-pollinated to cleistogamous flowers. Ottelia has been regarded as a pivotal group to study the evolution of Hydrocharitaceae, but the phylogenic relationships and evolutionary history of the genus remain unresolved. Here, we reconstructed a robust phylogenetic framework for Ottelia using 40 newly assembled complete plastomes. Our results resolved Ottelia as a monophyletic genus consisting of two major clades, which correspond to the main two centers of diversity in Asia and Africa. According to the divergence time estimation analysis, the crown group Ottelia began to diversify around 13.09 Ma during the middle Miocene. The biogeographical analysis indicated the existence of the most recent common ancestor somewhere in Africa/Australasia/Asia. Basing on further insights from the morphological evolution of Ottelia, we hypothesized that the ancestral center of origin was in Africa, from where the range expanded by transoceanic dispersal to South America and Australasia, and further from Australasia to Asia. We suggested that the climatic change and global cooling since the mid-Miocene, such as the development of East Asian monsoon climate and tectonic movement of the Yunnan-Guizhou Plateau (YGP), might have played a crucial role in the evolution of Ottelia in China.

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.

Primulina hochiensis var. ochroleuca (Gesneriaceae), a new variety from a limestone area of Guangxi, China, and errata on five new species of Primulina

Primulina hochiensis var. ochroleuca, a new variety from a limestone hill of karst areas, Guangxi, China is described with color photographs. It resembles P. hochiensis var. hochiensis, P. hochiensis var. ovata and P. hochiensis var. rosulata, but can be easily distinguished by a combination of characteristics, especially by its corolla color. We found only one population with approx. 3000 mature individuals at the type locality. This variety is provisionally assessed as vulnerable [VU C1] using IUCN criteria.

The genetic basis of hybrid male sterility in sympatric Primulina species

BACKGROUND

Sympatric sister species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain species boundaries. The persistence of morphologically and genetically distinct populations in sympatry can only occur if some degree of reproductive isolation exists. A pair of sympatric sister species of Primulina (P. depressa and P. danxiaensis) was used to explore the genetic architecture of hybrid male sterility.

RESULTS

We mapped one major- and seven minor-effect quantitative trait loci (QTLs) that underlie pollen fertility rate (PFR). These loci jointly explained 55.4% of the phenotypic variation in the F2 population. A Bateson-Dobzhansky-Muller (BDM) model involving three loci was observed in this system. We found genotypic correlations between hybrid male sterility and flower morphology, consistent with the weak but significant phenotypic correlations between PFR and floral traits.

CONCLUSIONS

Hybrid male sterility in Primulina is controlled by a polygenic genetic basis with a complex pattern. The genetic incompatibility involves a three-locus BDM model. Hybrid male sterility is genetically correlated with floral morphology and divergence hitchhiking may occur between them.

Analyses of Plastome Sequences Improve Phylogenetic Resolution and Provide New Insight Into the Evolutionary History of Asian Sonerileae/Dissochaeteae

Sonerileae/Dissochaeteae (Melastomataceae) comprises ca. 50 genera, two thirds of which occur in Southeast Asia. Phylogenetic relationships within this clade remain largely unclear, which hampers our understanding of its origin, evolution, and biogeography. Here, we explored the use of chloroplast genomes in phylogenetic reconstruction of Sonerileae/Dissochaeteae, by sampling 138 species and 23 genera in this clade. A total of 151 complete plastid genomes were assembled for this study. Plastid genomic data provided better support for the backbone of the Sonerileae/Dissochaeteae phylogeny, and also for relationships among most closely related species, but failed to resolve the short internodes likely resulted from rapid radiation. Trees inferred from plastid genome and nrITS sequences were largely congruent regarding the major lineages of Sonerileae/Dissochaeteae. The present analyses recovered 15 major lineages well recognized in both nrITS and plastid phylogeny. Molecular dating and biogeographical analyses indicated a South American origin for Sonerileae/Dissochaeteae during late Eocene (stem age: 34.78 Mya). Two dispersal events from South America to the Old World were detected in late Eocene (33.96 Mya) and Mid Oligocene (28.33 Mya) respectively. The core Asian clade began to diversify around early Miocene in Indo-Burma and dispersed subsequently to Malesia and Sino-Japanese regions, possibly promoted by global temperature changes and East Asian monsoon activity. Our analyses supported previous hypothesis that Medinilla reached Madagascar by transoceanic dispersal in Miocene. In addition, generic limits of some genera concerned were discussed.

A Preliminary Evaluation of The Karst Flora of Brazil Using Collections Data

Karst is defined as landscapes that are underlain by soluble rock in which there is appreciable water movement arising from a combination of high rock solubility and well-developed secondary (fracture) porosity. Karsts occupy approximately 20% of the planet’s dry ice-free land and are of great socioeconomic importance, as they supply water to up to 25% of the world’s population and represent landscapes of cultural and touristic importance. In Southeast Asia karst is associated with high species-richness and endemism in plants and seen as priority areas for the conservation of biodiversity. There has been little research into the floras associated with karst in South America, most of which occurs in Brazil. We therefore sought to evaluate the importance of Brazilian karst with respect to its species-richness and endemism. We sought to do so using curated plant specimen data in the Botanical Information and Ecology Network (BIEN) dataset. We show that, except for Amazonia, the BIEN dataset is representative of the Brazilian flora with respect to the total number of species and overall patterns of species richness. We found that karst is under-sampled, as is the case for much of Brazil. We also found that whilst karst represent an important source of plant diversity for Brazil, including populations of approximately 1/3 of the Brazilian flora, it is not significantly more species-rich or richer in small-range and endemic species than surrounding landscapes. Similarly, whilst important for conservation, comprising populations of 26.5–37.4% of all Brazilian species evaluated as of conservation concern by International Union for Nature Conservation (IUCN), karst is no more so than the surrounding areas. Whilst experimental error, including map resolution and the precision and accuracy of point data may have under-estimated the species-richness of Brazilian karst, it likely represents an important biodiversity resource for Brazil and one that can play a valuable role in conservation. Our findings are in sharp contrast to those for Southeast Asia where karst represents a more important source of species-richness and endemism. We also show that although BIEN represents a comprehensive and curated source of point data, discrepancies in the application of names compared to current more comprehensive taxonomic backbones, can have profound impacts on estimates of species-richness, distribution ranges and estimates of endemism.

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.

Landscape Features and Climatic Forces Shape the Genetic Structure and Evolutionary History of an Oak Species (Quercus chenii) in East China

Major topographic features facilitate intraspecific divergence through geographic isolation. This process may be enhanced by environmental isolation along climatic gradients, but also may be reduced by range shifts under rapid climatic changes. In this study, we examined how topography and climate have interacted over time and space to influence the genetic structure and evolutionary history of Quercus chenii, a deciduous oak species representative of the East China flora. Based on the nuclear microsatellite variation at 14 loci, we identified multiple genetic boundaries that were well associated with persistent landscape barriers of East China. Redundancy analysis indicated that both geography and climate explained similar amounts of intraspecific variation. Ecological differences along altitudinal gradients may have driven the divergence between highlands and lowlands. However, range expansions during the Last Interglacial as inferred from approximate Bayesian computation (ABC) may have increased the genetic diversity and eliminated the differentiation of lowland populations via admixture. Chloroplast (cp) DNA analysis of four intergenic spacers (2,866 bp in length) identified a total of 18 haplotypes, 15 of which were private to a single population, probably a result of long-term isolation among multiple montane habitats. A time-calibrated phylogeny suggested that palaeoclimatic changes of the Miocene underlay the lineage divergence of three major clades. In combination with ecological niche modeling (ENM), we concluded that mountainous areas with higher climatic stability are more likely to be glacial refugia that preserved higher phylogenetic diversity, while plains and basins may have acted as dispersal corridors for the post-glacial south-to-north migration. Our findings provide compelling evidence that both topography and climate have shaped the pattern of genetic variation of Q. chenii. Mountains as barriers facilitated differentiation through both geographic and environmental isolation, whereas lowlands as corridors increased the population connectivity especially when the species experienced range expansions.

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

The rediscovery and delimitation of Elatostemasetulosum W.T.Wang (Urticaceae)

Of the 280 species of Elatostema documented in China, 189 are known only from a single collection. Elatostemasetulosum is one such species, having been known only from the type collection for nearly half a century, until recent field investigations in Guangxi. Due to its morphological similarity to E.huanjiangense and E.tetracephalum, we undertook a critical review of all three species using morphological and molecular evidence. Our results suggest that all three names refer to the same species, which based on priority should be known as Elatostemasetulosum. We recognize E.huanjiangense and E.tetracephalum as synonyms. A distribution map of E.setulosum and the extinction risk according to the IUCN criteria is provided. After recircumscription, the taxon must be considered as Least Concern (LC).

Phylogeographic Analyses of the East Asian Endemic Genus Prinsepia and the Role of the East Asian Monsoon System in Shaping a North-South Divergence Pattern in China

Prinsepia Royle (Rosaceae) is a genus native to China and the Himalayan region. In order to explain its current fragmented distribution pattern and to compare the impact of relatively recent climate changes on the genetic structure of Prinsepia species in different regions of China, a total of 66 populations and 617 individuals of four species of Prinsepia were genotyped, using three cpDNA markers. Meanwhile, phylogenetic reconstructions and divergence dating were conducted using the cpDNA haplotypes dataset and the nuclear ribosomal internal transcribed spacer (ITS) dataset, respectively. Ecological niche modeling (ENM) was performed to predict the potential distribution of each species of Prinsepia at present and during the Last Glacial Maximum. Both ITS and cpDNA gene trees support a north-south divergence of Prinsepia species in China. The divergence time of the northern and southern Clades occurred around the late Oligocene epoch. Combining the present distribution of Prinsepia species and their habitats, we inferred that the transition to a monsoon climate system in East Asia during the late Oligocene epoch, created a humid forest vegetation zone from central to East China, which potentially gave rise to the north-south divergence of Prinsepia species. Both regional climates and allopatric divergence may have played an important role in the speciation of P. sinensis and P. uniflora. P. sinensis had the lowest genetic diversity and a putative northward post-glacial colonization. The distribution range of P. uniflora was also extremely sensitive to interglacial-glacial cycles. P. utilis from southwestern China preserved more haplotypes than P. sinensis and P. uniflora due to its multiple and isolated refugia.

Checklist to the Elatostema (Urticaceae) of Vietnam including 19 new records, ten new combinations, two new names and four new synonyms

Elatostema (Urticaceae) comprises several hundred herbaceous species distributed in tropical and subtropical Africa, Asia, Australia and Oceania. The greatest species richness occurs on limestone karst in Southeast Asia. Taxonomic revisions of Elatostema are largely out of date and contradict each other with respect to the delimitation of Elatostema and Pellionia. Most herbaria in SE Asia and worldwide contain significant amounts of unidentified material. As part of a broader revision of Elatostema in SE Asia, we present an updated checklist for Vietnam based on field visits, a review of specimens in herbaria worldwide, a review of type material and nomenclature. We recognize 77 taxa (75 species and two infraspecific taxa) of Elatostema in Vietnam, 23 of which were previously ascribed to Pellionia. Nineteen of these are new records for the country, i.e., E. attenuatoides, E. austrosinense, E. backeri, E. brunneinerve, E. crassiusculum, E. crenatum, E. fengshanense, E. glochidioides, E. malacotrichum, E. nanchuanense, E. oblongifolium, E. obtusum, E. oppositum, E. pergameneum, E. prunifolium, E. pseudolongipes, E. pycnodontum, E. salvinioides and E. xichouense. We place E. baviensis in synonymy of E. platyphyllum, E. colaniae in synonymy of E. myrtillus, P. macroceras in synonymy of E. hookerianum, and P. tetramera in synonymy of E. dissectum for the first time. Fourteen taxa (18% of all the recognized taxa) are endemic to Vietnam, which makes Elatostema one of the richest genera for endemic species in this country; this level of endemism is comparable to levels observed in Orchidaceae. Our checklist suggests that the highest diversity and endemism of Elatostema occurs in northern Vietnam, and that there is the greatest floristic similarity of northern Vietnam to SW China. The relationship among floristic regions is also investigated. We could find no records of Elatostema for 33 out of 63 provincial units of Vietnam, including all the southernmost provinces. We propose that further studies on the diversity of Elatostema in central and southern Vietnam are severely needed.