What explains high plant richness in East Asia? Time and diversification in the tribe Lysimachieae (Primulaceae)

Radiating diversification and niche conservatism jointly shape the inverse latitudinal diversity gradient of Potentilla L. (Rosaceae)

BACKGROUND

The latitudinal diversity gradient (LDG), characterized by an increase in species richness from the poles to the equator, is one of the most pervasive biological patterns. However, inverse LDGs, in which species richness peaks in extratropical regions, are also found in some lineages and their causes remain unclear. Here, we test the roles of evolutionary time, diversification rates, and niche conservatism in explaining the inverse LDG of Potentilla (ca. 500 species). We compiled the global distributions of ~ 90% of Potentilla species, and reconstructed a robust phylogenetic framework based on whole-plastome sequences. Next, we analyzed the divergence time, ancestral area, diversification rate, and ancestral niche to investigate the macroevolutionary history of Potentilla.

RESULTS

The genus originated in the Qinghai-Tibet Plateau during the late Eocene and gradually spread to other regions of the Northern Hemisphere posterior to the late Miocene. Rapid cooling after the late Pliocene promoted the radiating diversification of Potentilla. The polyploidization, as well as some cold-adaptive morphological innovations, enhanced the adaptation of Potentilla species to the cold environment. Ancestral niche reconstruction suggests that Potentilla likely originated in a relatively cool environment. The species richness peaks at approximately 45 °N, a region characterized by high diversification rates, and the environmental conditions are similar to the ancestral climate niche. Evolutionary time was not significantly correlated with species richness in the latitudinal gradient.

CONCLUSIONS

Our results suggest that the elevated diversification rates in middle latitude regions and the conservatism in thermal niches jointly determined the inverse LDG in Potentilla. This study highlights the importance of integrating evolutionary and ecological approaches to explain the diversity pattern of biological groups on a global scale.

Rapid in situ diversification rates in Rhamnaceae explain the parallel evolution of high diversity in temperate biomes from global to local scales

The macroevolutionary processes that have shaped biodiversity across the temperate realm remain poorly understood and may have resulted from evolutionary dynamics related to diversification rates, dispersal rates, and colonization times, closely coupled with Cenozoic climate change. We integrated phylogenomic, environmental ordination, and macroevolutionary analyses for the cosmopolitan angiosperm family Rhamnaceae to disentangle the evolutionary processes that have contributed to high species diversity within and across temperate biomes. Our results show independent colonization of environmentally similar but geographically separated temperate regions mainly during the Oligocene, consistent with the global expansion of temperate biomes. High global, regional, and local temperate diversity was the result of high in situ diversification rates, rather than high immigration rates or accumulation time, except for Southern China, which was colonized much earlier than the other regions. The relatively common lineage dispersals out of temperate hotspots highlight strong source-sink dynamics across the cosmopolitan distribution of Rhamnaceae. The proliferation of temperate environments since the Oligocene may have provided the ecological opportunity for rapid in situ diversification of Rhamnaceae across the temperate realm. Our study illustrates the importance of high in situ diversification rates for the establishment of modern temperate biomes and biodiversity hotspots across spatial scales.

Lysimachiadanxiashanensis, a new species of Primulaceae from Guangdong, China

Lysimachiadanxiashanensis, a new Primulaceae species, endemic to the Danxia landscape in Guangdong Province, China, is described and illustrated. This new species is morphologically similar to L.pseudohenryi, L.phyllocephala, L.congestiflora and L.kwangtungensis, but it differs from the similar species by its purplish-red plants, petiole without wings, calyx with orange glandular and the corolla margin serrated on upper half with orange-red glandular punctates. This new species belongs to Lysimachiasubgen.Lysimachiasect.Nummularia. Phylogenetic analysis confirmed that L.danxiashanensis is a distinct clade, based on the combined data of ITS and rbcL sequences. The conservation status of the new species was evaluated as Endangered (EN) according to IUCN Red List Categories and Criteria.

Plastid genome and its phylogenetic implications of Asiatic Spiraea (Rosaceae)

BACKGROUND

Spiraea L. is a genus comprising approximately 90 species that are distributed throughout the northern temperate regions. China is recognized as the center of species diversity for this genus, hosting more than 70 species, including 47 endemic species. While Spiraea is well-known for its ornamental value, its taxonomic and phylogenetic studies have been insufficient.

RESULTS

In this study, we conducted sequencing and assembly of the plastid genomes (plastomes) of 34 Asiatic Spiraea accessions (representing 27 Asiatic Spiraea species) from China and neighboring regions. The Spiraea plastid genome exhibits typical quadripartite structures and encodes 113-114 genes, including 78-79 protein-coding genes (PCGs), 30 tRNA genes, and 4 rRNA genes. Linear regression analysis revealed a significant correlation between genome size and the length of the SC region. By the sliding windows method, we identified several hypervariable hotspots within the Spiraea plastome, all of which were localized in the SC regions. Our phylogenomic analysis successfully established a robust phylogenetic framework for Spiraea, but it did not support the current defined section boundaries. Additionally, we discovered that the genus underwent diversification after the Early Oligocene (~ 30 Ma), followed by a rapid speciation process during the Pliocene and Pleistocene periods.

CONCLUSIONS

The plastomes of Spiraea provided us invaluable insights into its phylogenetic relationships and evolutionary history. In conjunction with plastome data, further investigations utilizing other genomes, such as the nuclear genome, are urgently needed to enhance our understanding of the evolutionary history of this genus.

The evolutionary responses of life-history strategies to climatic variability in flowering plants

The evolution of annual or perennial strategies in flowering plants likely depends on a broad array of temperature and precipitation variables. Previous documented climate life-history correlations in explicit phylogenetic frameworks have been limited to certain clades and geographic regions. To gain insights which generalize to multiple lineages we employ a multi-clade approach analyzing 32 groups of angiosperms across eight climatic variables. We utilize a recently developed method that accounts for the joint evolution of continuous and discrete traits to evaluate two hypotheses: annuals tend to evolve in highly seasonal regions prone to extreme heat and drought; and annuals tend to have faster rates of climatic niche evolution than perennials. We find that temperature, particularly highest temperature of the warmest month, is the most consistent climatic factor influencing the evolution of annual strategy in flowering plants. Unexpectedly, we do not find significant differences in rates of climatic niche evolution between perennial and annual lineages. We propose that annuals are consistently favored in areas prone to extreme heat due to their ability to escape heat stress as seeds, but they tend to be outcompeted by perennials in regions where extreme heat is uncommon or nonexistent.

Plastome evolution and phylogenomic insights into the evolution of Lysimachia (Primulaceae: Myrsinoideae)

BACKGROUND

Lysimachia L., the second largest genus within the subfamily Myrsinoideae of Primulaceae, comprises approximately 250 species worldwide. China is the species diversity center of Lysimachia, containing approximately 150 species. Despite advances in the backbone phylogeny of Lysimachia, species-level relationships remain poorly understood due to limited genomic information. This study analyzed 50 complete plastomes for 46 Lysimachia species. We aimed to identify the plastome structure features and hypervariable loci of Lysimachia. Additionally, the phylogenetic relationships and phylogenetic conflict signals in Lysimachia were examined.

RESULTS

These fifty plastomes within Lysimachia had the typical quadripartite structure, with lengths varying from 152,691 to 155,784 bp. Plastome size was positively correlated with IR and intron length. Thirteen highly variable regions in Lysimachia plastomes were identified. Additionally, ndhB, petB and ycf2 were found to be under positive selection. Plastid ML trees and species tree strongly supported that L. maritima as sister to subg. Palladia + subg. Lysimachia (Christinae clade), while the nrDNA ML tree clearly placed L. maritima and subg. Palladia as a sister group.

CONCLUSIONS

The structures of these plastomes of Lysimachia were generally conserved, but potential plastid markers and signatures of positive selection were detected. These genomic data provided new insights into the interspecific relationships of Lysimachia, including the cytonuclear discordance of the position of L. maritima, which may be the result of ghost introgression in the past. Our findings have established a basis for further exploration of the taxonomy, phylogeny and evolutionary history within Lysimachia.

Species richness disparity in tropical terrestrial herbaceous floras: evolutionary insight from Collabieae (Orchidaceae)

Species richness is spatially heterogeneous even in the hyperdiverse tropical floras. The main cause of uneven species richness among the four tropical regions are hot debated. To date, higher net diversification rates and/or longer colonization time have been usually proposed to contribute to this pattern. However, there are few studies to clarify the species richness patterns in tropical terrestrial floras. The terrestrial tribe Collabieae (Orchidaceae) unevenly distributes in the tropical regions with a diverse and endemic center in Asia. Twenty-one genera 127 species of Collabieae and 26 DNA regions were used to reconstruct the phylogeny and infer the biogeographical processes. We compared the topologies, diversification rates and niche rates of Collabieae and regional lineages on empirical samplings and different simulated samplings fractions respectively. Our results suggested that the Collabieae originated in Asia at the earliest Oligocene, and then independently spread to Africa, Central America, and Oceania since the Miocene via long-distance dispersal. These results based on empirical data and simulated data were similar. BAMM, GeoSSE and niche analyses inferred that the Asian lineages had higher net diversification and niche rates than those of Oceanian and African lineages on the empirical and simulated analyses. Precipitation is the most important factor for Collabieae, and the Asian lineage has experienced more stable and humid climate, which may promote the higher net diversification rate. Besides, the longer colonization time may also be associated with the Asian lineages' diversity. These findings provided a better understanding of the regional diversity heterogeneity in tropical terrestrial herbaceous floras.

Lysimachiafenghwaiana (Primulaceae), a new species from Hunan Province, China

A new species, Lysimachiafenghwaiana G.Hao & H.F.Yan (Primulaceae), from Hunan Province, China, is described and illustrated. This new species belongs to Lysimachiasubgen.Lysimachiasect.Nummularia and is morphologically similar to L.crista-galli and L.carinata, but is distinctive in its leaf shape and arrangement of flowers. It can be further distinguished from L.crista-galli by the absence of calyx lobule spur, and from L.carinata by the black glandular striates in the corolla lobes, rather than punctate.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

Ecological and biogeographic drivers of biodiversity cannot be resolved using clade age-richness data

Estimates of evolutionary diversification rates - speciation and extinction - have been used extensively to explain global biodiversity patterns. Many studies have analyzed diversification rates derived from just two pieces of information: a clade's age and its extant species richness. This "age-richness rate" (ARR) estimator provides a convenient shortcut for comparative studies, but makes strong assumptions about the dynamics of species richness through time. Here we demonstrate that use of the ARR estimator in comparative studies is problematic on both theoretical and empirical grounds. We prove mathematically that ARR estimates are non-identifiable: there is no information in the data for a single clade that can distinguish a process with positive net diversification from one where net diversification is zero. Using paleontological time series, we demonstrate that the ARR estimator has no predictive ability for real datasets. These pathologies arise because the ARR inference procedure yields "point estimates" that have been computed under a saturated statistical model with zero degrees of freedom. Although ARR estimates remain useful in some contexts, they should be avoided for comparative studies of diversification and species richness.

Characterization of the complete plastid genome of Lysimachia christinae Hance (Primulaceae)

Lysimachia christinae Hance is widely distributed in subtropical China at the elevational range from 500–2300 m. The species is an important medicinal herb for treating jaundice, urinary disorders, and the liver. Here, we sequenced and characterized the whole plastid genome of L. christinae. It is 154,810 bp in length, containing two copies of inverted repeat (IR) regions (26,034 bp, each), a large single-copy (LSC) region (84,809 bp), and a small single-copy (SSC) region (17,933 bp). It has 114 genes, of which 80 are protein-coding, 30 are tRNA, and 4 are rRNA genes. The ML tree indicates L. christinae is closely related to Lysimachia congestiflora Hemsl. This genome information can help us better construct a backbone phylogeny of Lysimachia in the future.

RASP 4: Ancestral State Reconstruction Tool for Multiple Genes and Characters

With the continual progress of sequencing techniques, genome-scale data are increasingly used in phylogenetic studies. With more data from throughout the genome, the relationship between genes and different kinds of characters is receiving more attention. Here, we present version 4 of RASP, a software to reconstruct ancestral states through phylogenetic trees. RASP can apply generalized statistical ancestral reconstruction methods to phylogenies, explore the phylogenetic signal of characters to particular trees, calculate distances between trees, and cluster trees into groups. RASP 4 has an improved graphic user interface and is freely available from http://mnh.scu.edu.cn/soft/blog/RASP (program) and https://github.com/sculab/RASP (source code).

Temporal patterns of diversification in Brassicaceae demonstrate decoupling of rate shifts and mesopolyploidization events

BACKGROUND AND AIMS

Whole-genome duplication (WGD) events are considered important driving forces of diversification. At least 11 out of 52 Brassicaceae tribes had independent mesopolyploid WGDs followed by diploidization processes. However, the association between mesopolyploidy and subsequent diversification is equivocal. Herein we show the results from a family-wide diversification analysis on Brassicaceae, and elaborate on the hypothesis that polyploidization per se is a fundamental driver in Brassicaceae evolution.

METHODS

We established a time-calibrated chronogram based on whole plastid genomes comprising representative Brassicaceae taxa and published data spanning the entire Rosidae clade. This allowed us to set multiple calibration points and anchored various Brassicaceae taxa for subsequent downstream analyses. All major splits among Brassicaceae lineages were used in BEAST analyses of 48 individually analysed tribes comprising 2101 taxa in total using the internal transcribed spacers of nuclear ribosomal DNA. Diversification patterns were investigated on these tribe-wide chronograms using BAMM and were compared with family-wide data on genome size variation and species richness.

KEY RESULTS

Brassicaceae diverged 29.9 million years ago (Mya) during the Oligocene, and the majority of tribes started diversification in the Miocene with an average crown group age of about 12.5 Mya. This matches the cooling phase right after the Mid Miocene climatic optimum. Significant rate shifts were detected in 12 out of 52 tribes during the Mio- and Pliocene, decoupled from preceding mesopolyploid WGDs. Among the various factors analysed, the combined effect of tribal crown group age and net diversification rate (speciation minus extinction) is likely to explain sufficiently species richness across Brassicaceae tribes.

CONCLUSIONS

The onset of the evolutionary splits among tribes took place under cooler and drier conditions. Pleistocene glacial cycles may have contributed to the maintenance of high diversification rates. Rate shifts are not consistently associated with mesopolyploid WGD. We propose, therefore, that WGDs in general serve as a constant 'pump' for continuous and high species diversification.

Accelerated diversification correlated with functional traits shapes extant diversity of the early divergent angiosperm family Annonaceae

A major goal of phylogenetic systematics is to understand both the patterns of diversification and the processes by which these patterns are formed. Few studies have focused on the ancient, species-rich Magnoliales clade and its diversification pattern. Within Magnoliales, the pantropically distributed Annonaceae are by far the most genus-rich and species-rich family-level clade, with c. 110 genera and c. 2,400 species. We investigated the diversification patterns across Annonaceae and identified traits that show varied associations with diversification rates using a time-calibrated phylogeny of 835 species (34.6% sampling) and 11,211 aligned bases from eight regions of the plastid genome (rbcL, matK, ndhF, psbA-trnH, trnL-F, atpB-rbcL, trnS-G, and ycf1). Twelve rate shifts were identified using BAMM: in Annona, Artabotrys, Asimina, Drepananthus, Duguetia, Goniothalamus, Guatteria, Uvaria, Xylopia, the tribes Miliuseae and Malmeeae, and the Desmos-Dasymaschalon-Friesodielsia-Monanthotaxis clade. TurboMEDUSA and method-of-moments estimator analyses showed largely congruent results. A positive relationship between species richness and diversification rate is revealed using PGLS. Our results show that the high species richness in Annonaceae is likely the result of recent increased diversification rather than the steady accumulation of species via the 'museum model'. We further explore the possible role of selected traits (habit, pollinator trapping, floral sex expression, pollen dispersal unit, anther septation, and seed dispersal unit) in shaping diversification patterns, based on inferences of BiSSE, MuSSE, HiSSE, and FiSSE analyses. Our results suggest that the liana habit, the presence of circadian pollinator trapping, androdioecy, and the dispersal of seeds as single-seeded monocarp fragments are closely correlated with higher diversification rates; pollen aggregation and anther septation, in contrast, are associated with lower diversification rates.

Chloroplast Genomes and Comparative Analyses among Thirteen Taxa within Myrsinaceae s.str. Clade (Myrsinoideae, Primulaceae)

The Myrsinaceae s.str. clade is a tropical woody representative in Myrsinoideae of Primulaceae and has ca. 1300 species. The generic limits and alignments of this clade are unclear due to the limited number of genetic markers and/or taxon samplings in previous studies. Here, the chloroplast (cp) genomes of 13 taxa within the Myrsinaceae s.str. clade are sequenced and characterized. These cp genomes are typical quadripartite circle molecules and are highly conserved in size and gene content. Three pseudogenes are identified, of which ycf15 is totally absent from five taxa. Noncoding and large single copy region (LSC) exhibit higher levels of nucleotide diversity (Pi) than other regions. A total of ten hotspot fragments and 796 chloroplast simple sequence repeats (SSR) loci are found across all cp genomes. The results of phylogenetic analysis support the notion that the monophyletic Myrsinaceae s.str. clade has two subclades. Non-synonymous substitution rates (d_N) are higher in housekeeping (HK) genes than photosynthetic (PS) genes, but both groups have a nearly identical synonymous substitution rate (d_S). The results indicate that the PS genes are under stronger functional constraints compared with the HK genes. Overall, the study provides hypervariable molecular markers for phylogenetic reconstruction and contributes to a better understanding of plastid gene evolution in Myrsinaceae s.str. clade.

Lysimachia fanii, a new species of Primulaceae from limestone area of Guangxi, China

Lysimachia fanii, a new species of Lysimachia (Subgen. Idiophyton, Primulaceae), is described and illustrated from Guangxi, China based on morphological and molecular data. Lysimachia fanii differs from L. verbascifolia, L. rupestris and L. alpestris mainly by the habit being nearly rosulate, leaves congested at the apex of the rhizome, leaf blades spatulate to narrowly oblanceolate and flowers solitary. Phylogenetic analyses supported L. verbascifolia as sister to L. fanii. This new species is endemic to limestone areas in Liucheng county of Guangxi, China.

Origin and macroevolution of micro-moths on sunken Hawaiian Islands

The origins and evolution of Hawaiian biodiversity are a matter of controversy, and the mechanisms of lineage diversification for many organisms on this remote archipelago remain unclear. Here we focus on the poorly known endemic leaf-mining moth genus Philodoria (Lepidoptera, Gracillariidae), whose species feed on a diversity of Hawaiian plant lineages, many of which are critically endangered. We use anchored hybrid enrichment to assemble the first phylogenomic dataset (507 loci) for any Hawaiian animal taxon. To uncover the timing and pattern of diversification of these moths, we apply two frequently used dating calibration strategies, biogeographic calibrations and secondary calibrations. Island calibrations on their own resulted in much younger and unrealistic dates compared to strategies that relied on secondary calibrations. Philodoria probably originated on the now partially sunken islands of Laysan or Lisianski, approximately 21 Ma, and were associated with host plants in the families Ebenaceae, Malvaceae or Primulaceae. Major feeding groups associated with specific host-plant families originated soon after the plants colonized the islands. Allopatric isolation and host shifts, in concert and independently, probably play major roles in the diversification of Philodoria Our dating results indicate that Philodoria is among the oldest known Hawaiian arthropod lineages, and that island calibrations alone can lead to unrealistically young dates.