Skip to main content

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

Abstract

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.

Peer Review reports

Background

Understanding the drivers of the global pattern in biodiversity is a main task for ecologists and biogeographers [1]. The latitudinal diversity gradient (LDG), which refers to the decrease in species richness from the equator to the poles, stands as one of the most pervasive and widely debated patterns in biogeography [2]. However, there are also some taxonomic lineages that exhibit inverse LDGs, in which species richness reaches its zenith in extratropical regions and decreases toward the tropical regions. Such inverse LDGs have been observed in insects [3, 4], amphibians [5], birds [6, 7], mammals [6], angiosperms [8,9,10], and gymnosperms [11, 12]. For example, the stenammine ants species are most diverse in the 35°N to 40°N range and decrease towards the equator [3], and the Rhamnaceae species diversity is markedly higher in temperate than in tropical regions, resulting in a bimodal latitudinal diversity gradient [10]. Determining the driving forces behind these less common richness patterns is helpful in deepening our understanding of the regular LDGs.

To elucidate the evolutionary mechanisms of species richness patterns, several hypotheses have been proposed in previous studies, including the time-for-speciation hypothesis [13], the diversification rate hypothesis [14], and the niche conservatism hypothesis [15]. The time-for-speciation hypothesis suggests that regions with high species richness were colonized earlier and hence had more time to increase their diversity via ongoing speciation [13]. In contrast, the diversification rate hypothesis assumes that elevated regional species richness might arise from high diversification rates (speciation rate minus extinction rate), driven by one or more ecological variables [14, 16]. Moreover, the contemporary environment could also affect species distribution [2, 17]. The niche conservatism hypothesis states that most species tend to maintain their ancestral niches [18]. Consequently, the species richness of a group is typically high in regions with environmental conditions similar to its ancestral niche, whereas in regions with environmental conditions different from its ancestral niche, diversity tends to be low due to niche conservatism limitations on dispersal [15, 19]. Importantly, these three hypotheses are not necessarily exclusive and might act together to explain richness patterns [2, 14].

In this study, we evaluate the roles of evolutionary time, diversification rates, and niche conservatism in shaping the distribution pattern of the species-rich genus Potentilla L. (Rosaceae). With an approximate count of 500 species, this genus is widely distributed in Europe, temperate Asia and North America. Only a few species can also penetrate into tropical regions, a pattern that strongly suggests an inverse LDG [20,21,22,23,24,25]. Most species of this genus occur in the open habitats situated within mountainous and arctic terrains, as well as in the drylands of western North America [26]. The divergence times estimation and ancestral area reconstruction suggested that the genus originated in Asia during the Eocene, and subsequently spread to Europe and North America [26, 27]. However, the inadequate sampling and the absence of distribution data make it difficult to answer why Potentilla exhibits high diversity in temperate regions.

Additionally, rapid diversification may hold significance in the evolutionary history of Potentilla, and the information on the drivers of accelerated diversification will provide valuable insights into the underlying mechanism towards inverse LDG within the genus. Most of the extant Potentilla species appeared during the late Miocene to Quaternary, a period considered to have experienced rapid radiation events in the genus [26]. In particular, the Argentea clade, which is the most diverse clade of Potentilla, predominantly diverged after ca. 5 Ma, yet it contains more than 300 species [28, 29]. Evolutionary radiation can be triggered by abiotic factors, such as abrupt climate or tectonic changes [30,31,32], as well as biotic factors, such as species intrinsic properties [30, 33]. Based on the considerable proportion of polyploids in Potentilla, Dobeš and Paule [26] hypothesized that polyploidization might be an important element to explain its rapid speciation, yet this hypothesis has not been confirmed to date. In addition, many adaptive morphological traits within Potentilla have been formed in order to adapt to extreme environments. For example, Potentilla species in Arctic region are typically dwarf herbs covered with dense hairs to resist severe cold [20]. Unfortunately, the contributions of these abiotic and biotic factors to the high-level diversity in Potentilla remain unclear due to the lack of modelling diversification rates.

In order to conduct comprehensive macroevolutionary analyses, a robust phylogenetic framework is a key necessity. Although previous studies have provided a basic phylogenetic disentanglement of Potentilla based on a few short DNA markers [26, 28, 34,35,36], in-depth infrageneric relationships remain unclear, thereby hampering a proper study on spatiotemporal evolutionary trends within the genus. Recent studies have shown that plastomes can fully resolve the phylogenetic relationships between different clades and within most clades [37, 38]. Therefore, this study also utilizes plastome sequences to reconstruct a solid phylogenetic backbone with more complete sampling. We aim to: (1) estimate the divergence times and reconstruct the ancestral areas of Potentilla; (2) carry out diversification rate analyses, and explore the roles of abiotic and biotic factors in the accelerated diversification of Potentilla; (3) reconstruct its ancestral thermal niche; and (4) reveal the relationships between species richness and evolutionary time, diversification rates, deviation from ancestral niche in latitudinal gradients, and test the three hypotheses on the inverse LDG of Potentilla.

Results

Phylogenomic analysis of Potentilla

Based on whole plastome sequences, Potentilla is strongly supported as monophyletic (Fig. 1; Fig. S1). A total of eight clades were identified with high bootstrap values (BS = 100% in Fig. S1). The Anserina clade is the earliest-diverging lineage and forms a sister group to Potentilla s.s. Within Potentilla s.s., the Centigrana clade is the basalmost, followed by Alba, Reptans, Ancistrifolia, Fragarioides, Ivesioid, and Argentea. Except for the Argentea clade, the internal relationships were fully or basally resolved within other clades. However, for the Argentea, which is the most diverse clade in Potentilla, the internal relationships were less clear. Along with three early derived species, six subclades could be retrieved in the Argentea clade, but all subclades only received low support values.

Fig. 1
figure 1

Historical biogeography of Potentilla. (a) Five biogeographic regions and potential dispersal routes of Potentilla. (b) Ancestral range estimation using BioGeoBEARS with the BAYAREALIKE + J model. Nodes 1–7 refer to the stem ages of the eight clades, nodes 8–15 refer to the crown ages of the eight clades, and nodes 16–17 refer to the dispersal events. (c) Temporal patterns in the net diversification rate of Potentilla across the globe and five biogeographic regions estimated by the empirical time-calibrated tree; red curve shows global temperature differences over the last 37.44 Ma as compared to current temperature and is modified from Westerhold et al. [88]. The map used in this study was downloaded from DIVA-GIS (http://www.diva-gis.org/Data). Qua. = Quaternary

Divergence time and ancestral range estimation

The divergence times estimated by the penalized likelihood method in treePL indicated that the stem age of Potentilla was 47.64 Ma (95% highest posterior density [HPD]: 45.85–50.65 Ma, Fig. S2). Potentilla split into the Anserina clade and Potentilla s.s. at 37.44 Ma (95% HPD: 35.63–40.13 Ma; node 1) (Fig. 1B). The Anserina clade began to diversify at 20.07 Ma (95% HPD: 19.13–21.05; node 8). Within Potentilla s.s., the stem node ages of the extant clades ranged from 18.39 to 7.95 Ma (nodes 2–7). However, except for the Centigrana and Alba clades, the other clades began to diversify after 7 Ma (nodes 10–15). All species of the Fragarioides clade (node 12), and the Argentea clade (node 15), as well as most species of the Ivesioid clade (node 14) originated after 3.08 Ma. Among the 154 taxa sampled in this study, 86% of them originated after ca. 3 Ma, and 52% originated after 1 Ma (Table. S10). For the 302 species without plastome, their simulated ages were presented in Fig. S3 and Table S11.

Ancestral area reconstruction showed that the genus Potentilla originated on the Qinghai-Tibet Plateau during the late Eocene (Fig. 1B, node 1). The Anserina clade was mainly distributed in the Qinghai-Tibet Plateau. Within Potentilla s.s., several dispersal events could be inferred (Fig. 1A, B). The early two dispersal events occurred between the Qinghai-Tibet Plateau and East Asia (nodes 2 and 9). Subsequently, three dispersal events were inferred from Asia to Europe and North America (nodes 7, 16, and 17). Among these, the first dispersal event occurred in the Alba clade, of which the most recent common ancestor originated on the Qinghai-Tibet Plateau and from where members of the clade subsequently dispersed into Europe during the middle Miocene (node 16). A second dispersal event occurred in the Reptans clade, of which the most recent common ancestor originated in East Asia and from where members of the clade subsequently dispersed towards Europe, North Asia, and North America during the Pliocene (node 17). A third dispersal event took place in the late Miocene (node 7) from East Asia to North America, thereby forming the Ivesioid clade being confined to eastern North America. Finally, Potentilla recolonized North Asia in the late Miocene from North America (node 15) after which some of its member dispersed again to the other four biogeographic regions. In accordance with the ancestral area reconstruction, multiple dispersal events between North America and Eurasia occurred in the Argentea clade. In addition, during a very short period (ca. 3 Ma), each biogeographic region is characterized by a rapid radiation, thereby forming the highly diverse Argentea clade (Fig. 1B).

Diversification rate analyses

The phylorate plot of the speciation rate revealed two shifts in diversification within Potentilla having high posterior probability (Fig. 2A; Fig. S4): one in the Ivesioid clade and the other in the Argentea clade. Diversification rates were low in early diverging clades, such as Anserina, Alba, and Reptans, but high in more recently diverged clades such as Fragarioides, Ivesioid, and Argentea. The Argentea clade showed rapid diversification after roughly 3 Ma with a diversification rate of 1.6 species per million years (Fig. 2A). In addition, the rate-through-time plots showed that the speciation and net diversification rates of Potentilla increased slowly prior to the Pliocene and then significantly accelerated, while extinction rates remained almost always constant (Fig. 2B). For the net diversification rates of Potentilla in the five biogeographical regions, whether based on the empirical or the simulated time-calibrated tree, their temporal trends are almost synchronized (Fig. 1C; Fig. S5). The early diversification within the Potentilla is shown to be rather slow, mainly occurring on the Qinghai-Tibet Plateau. Then, diversification rates throughout the whole Potentilla drastically increased in all biogeographical regions nearly simultaneously after the late Pliocene, a period when global temperatures were strongly decreasing (Fig. 1C). Interestingly, diversification rates in regions at lower latitudes (Qinghai-Tibet Plateau and East Asia) became gradually exceeded by those in regions at higher latitudes (North America, Europe, and North Asia) (Fig. 1C; Fig. S5).

Fig. 2
figure 2

Diversification rate analyses of Potentilla. (a) Phylorate plot showing the speciation rates along each branch. Arrows, shifts in speciation rates. The four traits used for trait-dependent analyses (A–D) are plotted at right—A: ploidy; B: leaves hair; C: root; D: basal leaves. The colors of the circles in A–D correspond to the boxplots c–f. (b) Speciation, extinction, and net diversification rates over time according to BAMM analysis. (c–e) Binary trait-dependent diversification inferred by HiSSE analyses. (f) Multistate trait-dependent diversification estimated by MuSSE analysis. Asterisks in (c-e) indicate significant differences according to t tests. ****P < 0.0001; N.S., not significant

The paleotemperature-dependent analysis employed in RPANDA v.1.9 show that the linear speciation without extinction with temperature dependence is the best model (AICc = 637.606) (Table S5), and the diversification rate increased with the decrease of temperature (α=-0.055), which was consistent with the results of BAMM analysis (Fig. 1C). For the trait-dependent analyses, the best model for the three binary traits was the full HiSSE model, which had unique speciation, extinction, and transition rates between the two observed character states and hidden states (AICc = 622.863, 746.362, and 622.863, respectively; Table S8). Polyploid species, as well as species with densely hairy leaves and robust root with few branches have higher net diversification rates (Fig. 2C–E). For the MuSSE analysis, the full model was selected as best best-fitting model, i.e., it had the lowest AICc (AICc = 930.716) (Table S9). Species with palmate shaped (≥ 5) basal leaves had higher net diversification rates than other basal leaf types (Fig. 2F). Interestingly, this palmate leaf type is generally present among the species of the Argentea clade.

Distribution pattern of species richness, evolutionary time, diversification rate, and deviation from ancestral niche

Distribution pattern analyses showed that Potentilla mainly occurred in the northern temperate zone, with only a few species extending into tropical regions (Fig. 3A, B). Species richness of Potentilla peaked around 45 °N, and decreased toward both high and low latitudes (Fig. 4). Most Potentilla species are confined to montane areas at middle latitudes, such as the Alpine Mountains in Europe, the Tianshan and Altai Mountains in North Asia, the Hengduan and Qilian Mountains within Qinghai-Tibet Plateau, the Qinling and Taihang Mountains in East Asia, and the Nevada Mountains in western North America (Fig. 3A, B). The distribution pattern of species richness based on the 149-dataset resembled that of the 451-dataset, except for the latter exhibiting higher Potentilla richness in western North America. Sampling bias analysis showed that the waterbodies have more obvious effect on sampling intensity than the airports, cities and roads (Fig. S6). Northeast and Central Asia showed a low number of collection records (Fig. S7).

Fig. 3
figure 3

Distribution of species richness (a, b), evolutionary time (c, d), mean diversification rates (e, f), and level of deviation from ancestral climate (g, h) of Potentilla calculated based on the 149-dataset and 451-dataset, respectively. The resolution of grid cells was 100 km × 100 km

Fig. 4
figure 4

Latitudinal trends in species richness, evolutionary time, diversification rates, and level of deviation from ancestral niche calculated based on the 149-dataset (a), and the 451-dataset (b)

Regarding the distribution pattern of evolutionary time, Potentilla species in Qinghai-Tibet Plateau and East Asia are much older than those in Europe, North America and North Asia (Fig. 3C, D). In the latitudinal gradient, whether based on empirical or simulated datasets, the evolutionary time at ca. 20 °N is significantly higher than that at other latitudes (Fig. 4). These results are consistent with the dispersal histories of Potentilla as revealed by the ancestral area reconstruction (Fig. 1).

The mean diversification rates in Europe, western North Asia and western North America were significantly higher than in all other regions, corroborating the higher net diversification rates estimated by BAMM analyses in these regions (Fig. 1C; Fig. S5). In the latitudinal gradient, diversification rates generally increase with increasing latitude, punctuated by a slight trough around 65 °N (Fig. 4). However, at high latitudes (from 60 °N to 80 °N), the diversification rates calculated based on the 451-dataset were noticeably higher than those based on the 149-dataset. This is possibly because the 451-dataset includes more species from North Asia, North America, and Europe, most of which appeared explosively after 3 Ma (Fig. S3; Table S3).

The ancestral thermal niches of Potentilla reconstructed for MAT, MTWM, and MTDQ were 6.66 ℃, 23.68 ℃, and 1.22 ℃, respectively (Fig. S8), implying that Potentilla likely originated in a relatively cool environment. Notably, ancestral temperature inferred from the simulated time-calibrated tree was slightly higher than that from the empirical time-calibrated tree (Figs. S8 and S9). This could be attributed to the presence of some inserted Anserina species, which are native to low-latitude regions, particularly New Guinea (Table S2), and consequently elevated the estimated temperature of the root node. The level of deviation from ancestral niche, as well as the ΔMAT, ΔMTWM, and ΔMTDQ, showed lower values at the diversity hotspots of Potentilla, such as Europe, western North America, East Asia, and the surrounding area of Qinghai-Tibet Plateau (Fig. 3G, H; Fig. S10). Both datasets showed that the contemporary environment at ca. 45 °N is most similar to the ancestral niche of Potentilla, with increased deviation towards to lower and higher latitudes (Fig. 4).

Discussion

Phylogeny of Potentilla

We present a robust phylogenetic framework for Potentilla based on whole plastome data and a global sampling (Fig. S1). The infrageneric relationships among Potentilla were generally congruent with prior studies, except for some minor differences [26, 28, 35, 37, 38]. Based on three plastid DNA makers, Dobeš and Paule [26] placed Argentea as sister to Fragarioides, and their relationships with Ivesioid clade are unclear. Our plastome-based phylogeny supported Ivesioid as a sister clade of Argentea with high support. For the Alba clade, Feng et al. [28] indicated that the Himalayan clade (including P. purpurea, P. tetrandra, and P. suavis) formed a sister group to the remaining species of Potentilla s.s. in plastid gene tree. However, we resolved the Himalayan clade in the Alba clade as an early derived lineage, which was supported by other plastid-based studies [26, 35,36,37,38]. In general, compared with previous studies, most clades had good sampling representativeness with the internal relationships being well resolved, thereby facilitating further biogeographic, evolutionary, and ecological studies.

Biogeographic history of Potentilla in the Northern Hemisphere

The current divergence time estimation suggested a crown group age of 37.44 Ma (95% HPD: 35.63–40.13) for Potentilla. This is in agreement with the results of Dobeš and Paule [26] (95% HPD: 34.3–45.2 Ma). Additionally, the divergence times of major nodes within the genus Potentilla are in congruence with the above-mentioned study. Our ancestral area reconstruction analysis indicated that the most recent common ancestor of Potentilla originated in the Qinghai-Tibet Plateau around the Eocene-Oligocene boundary (Fig. 1B, node 1), which was more accurate than the previously suggested Asian origin [26]. Then, the Potentilla diverged into two groups: the Anserina clade and Potentilla s.s. The former was mainly distributed in high-altitude areas of the Qinghai-Tibet Plateau, especially in the Hengduan Mts (Fig. S11H). The Anserina clade started to diverge from the early Miocene to the Quaternary, which corresponds with the uplift time of the southeastern Qinghai-Tibet Plateau [39]. As such, the uplift of the Plateau might potentially have triggered the in situ diversification of the Anserina clade, similar to what happened in Rhodiola [40] and Saxifraga [41].

From the middle Miocene to the early Pliocene, there were three out-of-Asia dispersal events. The first dispersal event occurred from the Qinghai-Tibet Plateau towards Europe, and involved the Alba clade (Fig. 1B, node 16). Except for P. biflora, which was also found in North America [24], the other species of the Alba clade are confined to high mountains of Eurasia, particularly the Himalayas and the Alps (Fig. S11F). Therefore, we hypothesize that the Alba clade spread over the Himalayas towards the Caucasus (P. divina and P. brachypetala) and the Alps (P. nitida, P. caulescens, and P. alba). Although the Himalayas Mts have uplifted to a considerable elevation during this period [42], this dispersal path was still feasible given that most species of this clade were distributed at high altitudes (Fig. S12). Similar dispersal routes were reported for Gentiana [43]. The other two out-of-Asia dispersal events took place in the Reptans clade (node 17) during the early Pliocene and in the Ivesioid-Argentea clade (node 7) during the late Miocene. We postulate that the Reptans clade spread from East Asia towards North Asia and Europe (P. reptans and P. erecta), and then further migrated towards eastern North America (P. canadensis) via the North Atlantic land bridge (NALB). Although it was previously assumed that the NALB became inexistent after the middle Miocene [44], recent studies suggested that a discontinuous land bridge could still have acted as a corridor for plant migration until the onset of the Pliocene [45, 46]. For the latter, given that the Bering Land Bridge (BLB) served as a corridor for plant migrations from the late Cretaceous to late Neogene [46], we inferred that the Ivesioid-Argentea clade spread from East Asia to North America via the BLB. Subsequently, the Ivesioid-Argentea clade split into two clades: the Ivesioid clade, which diversified in west North America, and the Argentea clade, which returned to North Asia. The Ivesioid clade, except P. biennis, is characterized by many drought-tolerant traits enabling it to adapt to extremely xeric conditions, indicating that aridification in western North America triggered diversification of the Ivesioid clade [47]. The species-rich clade, Argentea, probably spread from North America to North Asia via the BLB and underwent multiple transcontinental migrations and evolutionary radiation in the Northern Hemisphere over a very short period beginning in the late Pliocene.

Rapid diversification of Potentilla from the late Pliocene to Quaternary

The diversification rates in the genus Potentilla accelerated during the late Pliocene to Quaternary, when the global temperature markedly dropped (Fig. 1C). Furthermore, paleotemperature and trait-dependent analyses showed that a decrease in temperature, polyploidization, and some adaptive traits may have promoted an increased diversification in Potentilla (Fig. 2). Abrupt environmental changes, such as cooling, may increase the frequency of ameiotic gamete formation, thereby promoting polyploid formation [48]. Compared with diploid taxa, polyploids tend to be more adaptable to extreme environments, which has been suggested to be important for the adaptation to lower temperatures [49]. In China, the polyploid frequency of angiosperms increased from low to high latitudes [50]. In this study, polyploid species in Potentilla are characterized by higher diversification rates than for diploids. In addition, polyploid Potentilla’s diverged recently, and are distributed in colder environments (Fig. S13). The Argentea clade contains a large number of polyploids and had the highest species richness diversification rates (Fig. 2). All species of this clade rapidly diverged after the late Pliocene, with many of them adapted to extremely low temperatures (Figs. S8 and S9). Moreover, some morphological traits in Potentilla, such as densely hairy leaves and robust root were considered to help plants protect against cold conditions [51,52,53]. These adaptive traits are typically observed in the cold-adapted Argentea clade, as well as in the high-altitude adapted Alba and Anserina clades. By contrast, species of the Reptans and Fragarioides clades prefer warmth and occur in lower altitudes, and typically have sparsely hairy leaves and slender root (Fig. 2; Figs. S8 and S9). Therefore, we postulate that the rapid cooling after the late Pliocene can be regarded as an important driver of polyploidization in Potentilla. The polyploidization, accompanied by innovations of adaptive traits, including densely hairy leaves and robust root with few branches enhanced the adaptation of Potentilla species to the cold environment. Although paleobotanists believe that Quaternary speciation was rare due to climatic oscillations [54], a recent review of molecular phylogenetic studies found that plant speciation and radiation were ubiquitous during the Quaternary [55]. Our findings provide an example of how plants adapted to low temperatures during the Quaternary and underwent rapid radiation in high latitude regions.

Evaluation of alternative LDG hypotheses

The LDG can be explained by latitudinal variations in environmental factors, diversification rates, and the time for species accumulation [2]. We found that Potentilla has a strong inverse diversity gradient, with most species inhabiting the northern temperate zone between 35 and 55°N latitude. In contrast, evolutionary time generally decreased from low (ca. 20 °N) to high latitudes, possibly because the oldest clade of Potentilla, Anserina, was mainly distributed in the Qinghai-Tibet Plateau. Ancestral area reconstruction showed that Potentilla originated in the Qinghai-Tibet Plateau during the late Eocene (ca. 37 Ma) and gradually spread to other regions of the Northern Hemisphere after the late Miocene (ca. 10 Ma). Although being colonized later, regions at higher latitudes (North America, Europe, and North Asia) have accumulated more Potentilla species compared to regions at lower latitudes (Qinghai-Tibet Plateau and East Asia) (Table S3). Consequently, the time-for-speciation hypothesis is not applicable to the inverse LDG of Potentilla.

The rapid radiation after the late Pliocene contributed 90% of extant Potentilla diversity, in sharp contrast to the relatively ancient origin of Potentilla, indicating that rapid radiation contributed to the overall diversity in Potentilla. Despite the differences between the 149- and 451-datasets, it is undeniable that diversification rates in the North America, Europe, and North Asia are significantly higher than those in the Qinghai-Tibet Plateau and East Asia (Figs. 1C and 3E-F, and 4). This trend aligns with the higher Potentilla richness observed in the former three regions (Table S3). Thus, the high diversity of Potentilla in middle latitudinal regions of the Northern Hemisphere likely results from the elevated diversification rates in these regions.

Ancestral niche reconstructions indicate that most of the Potentilla species have existed in temperate environmental conditions since their origin (Figs. S8 and S9), suggesting a strong tendency to remain in their ancestral climatic niche. With increased deviation from the ancestral climate, the diversity of overall Potentilla species decreased (Fig. 4). Previous studies have also shown that plants originating in temperate climates are less adaptable to warmer conditions but more cold-tolerant than those originating in tropical climates [56]. Therefore, we speculate that the conservatism of Potentilla species in cool environments restricted their distribution to the northern temperate zone, and imposed climatic constraints on dispersal into tropical regions. Meanwhile, for species that are conservative in their climatic niche, altering their existing distribution range may prove challenging [8]. Given that many Potentilla plants are common herbal medicines used in folk medicine [57], human activities such as overharvesting may cause some species to become threatened or even extinct. So, it is necessary to carry out in situ or ex situ conservation for these species.

In conclusion, our results suggest that the radiating diversification and niche conservatism jointly shaped the inverse LDG of Potentilla. It is pervasive for species to retain the climatical preferences of their ancestors [18, 19]. An increasing number of studies have found that niche conservatism play a major role in shaping various distribution patterns, including LDGs [15, 58], inverse LDGs [3, 5, 8, 10], and others [59, 60]. However, the influence of evolutionary time and diversification rates on driving species distribution patterns varies greatly across lineages, regardless of whether these distribution patterns conform to LDGs or inverse LDGs [3, 6, 9, 10, 12]. Generally, the lineages with ancient origins tend to favor the time-for-speciation hypothesis, such as Cycadaceae (stem age: ca. 260 Ma) [12], Pinus (ca.155 Ma) [11], stenammine ants (ca. 90 Ma) [3], Caudata (ca. 198 Ma) [6], while lineages with relatively young origins usually support the diversification rate hypothesis, such as Zygophyllaceace (ca. 69 Ma) [60], Rhododendron (ca. 68 Ma) [17], Anseriformes (ca. 71 Ma) [6], Procellariiformes (ca. 61 Ma) [6]. This may be because the effect of the diversification rate on the distribution pattern of a taxon can manifest over a short period, while the effect of evolutionary time is relatively delayed [2]. Although a consensus on the mechanisms behind LDGs or inverse LDGs has not yet been reached, it is evident that the species distribution patterns are driven by multiple factors, and different hypotheses are not necessarily exclusive [2, 6, 60], Therefore, it is crucial to integrate ecological and evolutionary approaches to investigate the mechanisms underlying the global patterns of biological groups.

Conclusion

Through interdisciplinary evidence from evolutionary ecology, we have illuminated the underlying mechanism driving the inverse LDG of Potentilla. Our results indicate that Potentilla originated in the Qinghai-Tibet Plateau during the late Eocene and then spread to other regions of the Northern Hemisphere after the late Miocene. It is postulated that the decreased temperature after the late Pliocene significantly accelerated the diversification rates of Potentilla. The reconstruction of thermal niche suggests that the ancestor of Potentilla preferred relatively cool temperatures, and the regions where environmental conditions similar to its ancestral niche usually have high species diversity. The middle latitude regions, been colonized more recently but harbored higher diversity, suggest that the time-for-speciation effect may not fully explain the inverse LDG of Potentilla. While the diversification rates hypotheses and niche conservatism hypotheses jointly contributed to the inverse LDG of Potentilla. The integrative framework incorporating both ecological and evolutionary approaches demonstrated here can lead to a better understanding of the factors shaping global patterns of plant diversity.

Although this study reconstructed a time-calibrated phylogeny for Potentilla with most complete sampling to date, our sampling for species in North America, Europe, and North Asia remains limited. Most of these unsampled species are regional endemic, posing a huge challenge for our sampling. Future studies should increase the sampling of Potentilla species from middle-high latitude regions. Furthermore, field surveys in Northeast and Central Asia should be intensified to compile a more complete geographical distribution database for Potentilla.

Materials and methods

Taxon sampling, genome-skimming sequencing, and plastome assembly

In this study, we adopted a broad circumscription for Potentilla, as in prior reports [26, 35, 61], including the genera Horkelia, Horkeliella, and Ivesia of the Ivesioid clade, and Argentina of the Anserina clade. A total of 184 plastomes were analyzed. Among these, 11 were newly sequenced and 8 were assembled from genome skimming data downloaded from the GenBank SRA database; the others were downloaded from GenBank (Table S1). Samples investigated comprised 179 species of Rosaceae, of which 149 are ingroup species and 30 are outgroup species. Our sampling covered all eight proposed clades of Potentilla [35, 37]. For the most species-rich Argentea clade, we sampled 94 taxa, comprising 92 species, 1 variety, and 1 subspecies. The newly generated plastomes have been uploaded to GenBank. The formal identification of these samples was undertaken by Tiantian Xue and Shengxiang Yu (Institute of Botany, Chinese Academy of Sciences). All voucher specimens were deposited in the herbarium of the Institute of Botany, Chinese Academy of Sciences (PE). Taxon names, voucher information, and GenBank accession numbers are listed in Table S1. To avoid under sampling, we also included species not covered in phylogeny that met both of the following criteria: (1) species name is accepted in the Plant of the World Online (POWO, https://powo.science.kew.org/), (2) species can be determined which phylogenetic lineage it belongs to by previous molecular analyses, Flora, type specimens or other literatures, and (3) occurrence points can be obtained from Global Biodiversity Information Facility (GBIF) or Chinese Virtual Herbarium (CVH). Finally, a total of 302 species were selected for supplemental analyses (Table S2). The sample fraction was calculated at global and region scales (Table S3).

We extracted total genomic DNA from silica-dried leaf materials or herbarium specimens using a modified CTAB method [62]. All of the samples were sequenced at Novogene Corporation (Beijing, China), which generated sequencing libraries using the NEB Next® Ultra™ DNA Library Prep Kit for Illumina following the manufacturer’s recommendations. The DNA libraries were sequenced on the Illumina HiSeq X-Ten platform and 150-bp paired-end reads were generated. We obtained the genome skimming data for each sample (ca. 4 Gb). Before assembling the plastomes, the raw sequence reads were subjected to quality control by Novogene Corporation using fastq software [63] to obtain high-quality clean reads. De novo plastome assembly was performed in the GetOrganelle [64] pipeline with default parameters.

Phylogenetic analysis and divergence time estimation

Whole plastomes with one inverted repeat region excluded were used for phylogenetic reconstruction and molecular clock analysis. We first aligned the plastomes using MAFFT v. 7.490 [65] with the auto strategy, and poorly aligned regions were removed by trimAL v. 1.4 [66] using the command ‘-automated1.’ Subsequently, we performed a maximum-likelihood (ML) analysis in RAxML-NG (RAxML Next Generation) [67], with 1,000 rapid bootstrap (BS) replicates. The best-fit substitution model (GTR + F + I + G4) was selected by ModelFinder [68] according to the corrected Akaike information criterion (AICc).

Divergence times were estimated using the penalized likelihood method in treePL, which is suitable for dealing with large datasets with hundreds of species [69]. The ML tree generated from RAxML-NG as the input tree. First, we identified the best optimization parameters using the prime command. Next, we carried out a cross-validation analysis with test values of 10− 20 to 1020 to determine the best smoothing value, which was 10− 11. Finally, 1,000 ML bootstrap trees with branch lengths generated from RAxML-NG were used to estimate the age of each internal node, and the 95% HPD intervals were calculated using TreeAnnotator v. 2.6.7 in the BEAST package [70]. Although a number of fossils have been reported for Potentilla, most of which lacked key identified traits, and their reliability needs to be further evaluated (http://ifpni.org/genus.htm?id=CFE35DDB-0C1E-4248-B033-1D88273E7B95, accessed 15 July 2022). Therefore, two calibration points were utilized within outgroups. The crown age of Rosaceae was constrained with maximum and minimum ages between 96.36 and 94.46 Ma based on the estimate of Zhang et al. [71]. The oldest fossil record of the genus Rosa (Rosa germerensis, 55.8–48.6 Ma) was used to calibrate the crown age of the Rosa-Potentilleae clade (cf. Palaeobiology database https://paleobiodb.org) [72]. To assess the impact of missing species on subsequent analyses, the 302 Potentilla species without plastomes were added randomly to their respective clades in the empirical time-calibrated tree, using the add.species function in the R package phytools v. 0.7-00 [73]. This approach has been commonly employed in previous macroevolutionary studies of highly diverse groups [74,75,76]. Ultimately, two time-calibrated trees were generated: an empirical tree comprising 149 species with plastomes, and a simulated tree comprising 451 species.

Ancestral area reconstruction

To trace the biogeographic history of Potentilla, ancestral area reconstruction was conducted using BioGeoBEARS [77] implemented in RASP v. 4.3 [78]. The time-calibrated phylogenetic tree was obtained from treePL, and outgroups were pruned. According to the endemism centers and extant distribution patterns of Potentilla, five geographic regions were defined: North America, Europe, North Asia, Qinghai-Tibet Plateau, and East Asia. Six models were tested in RASP v. 4.3: the dispersal-extinction-cladogenesis (DEC), the likelihood version of dispersal-vicariance (DIVALIKE), the BayArea likelihood version of the range evolution model (BAYAREA), and the + J versions of these models, which include founder-event speciation. The BAYAREALIKE + J model with the highest log-likelihood (LnL) value and low AICc value (LnL = -406.3, AICc = 818.7, Table S4) was selected as the best-fit model.

Diversification rate analyses

The diversification rate of Potentilla was estimated using Bayesian Analysis of Macroevolutionary Mixtures (BAMM) (v.2.5) [79]. Initially, the prior values were selected using the setBAMMpriors function of BAMMtools in R v.4.2.2 [80]. The parameter globalSamplingFraction was set to be 0.3 following the sampling coverage in the current phylogeny to account for incomplete sampling, and the Markov chain Monte Carlo (MCMC) was run for 10 million generations and sampled every 5,000 generations. The initial 10% of samples of the MCMC run were discarded as burn-in, and the remaining data were assessed for convergence to ensure that the ESS values were > 200. We used the addBAMMshifts function to infer possible shifts of speciation across the phylogeny. Subsequently, we used the PlotRateThroughTime function in BAMMtools to plot the net diversification rates through time for the genera at the global scale and in the five biogeographic regions. Given that the uneven sampling fraction may introduce bias into the diversification rate estimation in these five biogeographical regions. We also plotted the net diversification rates through time for the five biogeographic regions based on the simulated time-calibrated tree, with globalSamplingFraction set to 0.9.

To explore the effects of abiotic and biotic factors on the diversification rates of Potentilla, we carried out paleoenvironment- and trait-dependent analyses, respectively. We used RPANDA v. 1.9 [81] to evaluate seven paleotemperature-dependent birth–death models (Table S5). The relationships of speciation and extinction rates with paleotemperature were assumed to be constant, linear, or exponential. The global paleotemperature data set was obtained from Sun et al. [82]. We compared the likelihood values, and the model with the smallest AICc was considered the best diversification model.

For the trait-dependent analyses, four traits were selected, including three binary traits (diploidy vs. polyploidy; leaves densely hairy vs. sparsely hairy; and roots robust with less branching vs. slender with much branching) and one multi-state trait (basal leaves ternate vs. palmate ≥ 5 vs. pinnate ≤ 10 pairs vs. pinnate > 10 pairs) (Table S6). Chromosome ploidy data of Potentilla were extracted from the Chromosome Counts Database (http://ccdb.tau.ac.il/, accessed 7 July 2022) and Index to Plant Chromosome Numbers (http://legacy.tropicos.org/Project/IPCN, accessed 7 July 2022). The ploidy of a species was divided into diploidy and polyploidy. For species with multiple cytotypes, we scored ploidy as the most frequently reported ploidy level. The three morphological traits were selected and coded based on the Flora of China [23], Flora of North America [24, 25], Flora Europaea [22], and Flora URSS [21] and the herbarium specimens in JSTOR (https://plants.jstor.org/) and the Chinese Virtual Herbarium (http://www.cvh.ac.cn/). The effects of the three binary states on diversification rates were estimated using the hidden state speciation and extinction (HiSSE) model, which enables testing of hypotheses related to the effects of observed traits and incorporates unmeasured factors [83]. Twenty-five models were tested in the R package HISSE v. 1.9.10 and the model with the lowest AICc was selected as the best-fit model (Tables S7 and S8). The extinction, speciation, and net diversification rates of each state were calculated. For the basal leaf type, multi-state speciation and extinction (MuSSE) analysis was performed in DIVERSITREE v. 0.9.10 [84]. Four models were tested: a null model with fully constrained variables; a full model allowing all variables to change independently; a model constraining the extinction rate (µ) to be equal and allowing the speciation rate (λ) to vary (free λ); and a model constraining the λ to be equal and allowing µ to vary (free µ) (Table S9). Next, the net diversification rate for each state was obtained by a Bayesian approach (MCMC analysis) with an exponential prior with 5000 generations.

Ancestral niche reconstruction

In order to estimate the ancestral niche of Potentilla, we first compiled a database including the occurrence data and climate data, and conducted principal component analysis (PCA) to determine the main driving climatic factors for the distribution pattern of Potentilla. The occurrence data were obtained from CVH (http://www.cvh.ac.cn/) and GBIF (downloaded on 01 January 2021, from https://doi.org/10.15468/dl.fwp545). Then, we removed records meeting one of the following criteria: (1) without longitude or latitude; (2) equal latitude and longitude; (3) zero and/or integer latitude and longitude; (4) coordinates falling within sea; (5) coordinates were outside the species’ native range recorded in POWO. A total of 214,231 distribution records for 451 Potentilla species passed the above filtering criteria. For each record, we extracted the values of 19 bioclimatic variables from WorldClim (http://worldclim.org) using ArcGIS 10.6. Highly correlated variables (absolute Pearson’s correlation coefficient r ≥ 0.75) were removed to reduce collinearity. The remaining eight bioclimatic variables were: mean annual temperature (MAT), isothermality, max temperature of warmest month (MTWM), mean temperature of wettest quarter (MTWQ), mean temperature of driest quarter (MTDQ), annual precipitation (AP), precipitation seasonality (PS), precipitation of coldest quarter (PCQ). Then, a PCA was performed using the R packages FactoMineR v.2.6 [85] and factoextra v.1.0.7 [86]. The first two principal components (PCs) of bioclimatic variables explained 70.2% of the total climate variation in Potentilla, and the MAT, MTWM, and MTDQ had the largest contribution to the first two PCs (Fig. S14). Therefore, we reconstructed the ancestral states of these three bioclimatic variables using the fastAnc function in phytools. Mean values of the MAT, MTWM, and MTDQ for each species were used in the reconstruction.

Distribution pattern of species richness, evolutionary time, diversification rate, and deviation from ancestral climate

Considering that the Potentilla comprises at least 500 species and it is difficult to obtain a complete sampling at the plastome level, here we evaluated the effect of missing species on macroevolutionary analyses. Two datasets were employed to investigate the distribution pattern of Potentilla: one dataset comprised solely of the 149 species with genomic data (referred to as the “149-dataset”), while the other dataset encompassed species with and without genomic data (referred to as the “451-dataset”). To mitigate the effect of area on diversity estimation and potential sampling incompleteness, we divided the map of world into 31,203 grid cells with a 100 km × 100 km (1 × 104 km2) resolution utilizing ArcGIS 10.6. Both datasets were spatially aligned with these grid cells using ArcGIS 10.6. Subsequently, species richness was computed within each grid cell. In addition, we also used the R package ‘sampbias’ [87] to quantify the specific and combined effects of roads, rivers, airports and cities on recorded occurrence data.

Evolutionary time was represented by the age of the oldest species present within a given grid cell. Diversification rate was represented by the mean diversification rate of all species in the grid cell. Species ages were extracted from the time-calibrated tree, while diversification rate of each species was calculated using the DR_statistic function in phytools (Table S10). For the level of deviation from ancestral niche, we first calculate the mean value of MAT, MTWM, and MTDQ in each grid cell. The difference between the present and ancestral states in each grid cell was calculated and expressed by ΔMAT, ΔMTWM, and ΔMTDQ, respectively. Just as the equations below:

ΔMAT = |MAT present – MAT ancestor|

ΔMTWM = |MTWM present – MTWM ancestor|

ΔMTDQ = |MTDQ present – MTDQ ancestor|

Then, the three metrics (ΔMAT, ΔMTWM, and ΔMTDQ) were normalized across a range of 0 to 1, with their summation serving as an indicator to measure the deviation between ancestral and contemporary niche. Grid cells with small value of the metrics and their sum indicate their environmental conditions are similar to the ancestral niche of Potentilla. Ultimately, the distributions of species richness, evolutionary time, mean diversification rate, and deviation from ancestral climate were visualized using ArcGIS 10.6. Moreover, their latitudinal trends were modeled via local polynomial regression.

Data availability

New sequenced and other published plastomes in this study can be found in GenBank (https://www.ncbi.nlm.nih.gov/genbank/), and the accession numbers showed in Table S1.

Abbreviations

AICc:

Corrected Akaike Information Criterion

BAMM:

Bayesian analysis of macro-evolutionary mixtures

BAYAREA:

BayArea likelihood version of the range evolution model

BS:

Bootstrap

ca.:

Circa

HiSSE:

Hidden state speciation and extinction

HPD:

Highest posterior density

LDG:

Latitudinal diversity gradient

Ma:

Million years ago

MAT:

Mean annual temperature

MCMC:

Markov chain Monte Carlo

ML:

Maximum likelihood

MTDQ:

Mean temperature of driest quarter

MTWM:

Max temperature of warmest month

MuSSE:

Multi-state speciation and extinction

PCA:

Principal component analysis

References

  1. Schluter D, Pennell MW. Speciation gradients and the distribution of biodiversity. Nature. 2017;546(7656):48–55.

    Article  CAS  PubMed  Google Scholar 

  2. Pontarp M, Bunnefeld L, Cabral JS, Etienne RS, Fritz SA, Gillespie R, et al. The latitudinal diversity gradient: novel understanding through mechanistic eco-evolutionary. Trends Ecol Evol. 2019;34(3):211–23.

    Article  PubMed  Google Scholar 

  3. Branstetter MG, Longino JT, Reyes-Lopez JL, Brady SG, Schultz TR. Out of the temperate zone: a phylogenomic test of the biogeographical conservatism hypothesis in a contrarian clade of ants. J Biogeogr. 2022;49(9):1640–53.

    Article  Google Scholar 

  4. Orr MC, Hughes AC, Chesters D, Pickering J, Zhu C-D, Ascher JS. Global patterns and drivers of Bee distribution. Curr Biol. 2021;31(3):451–.

    Article  CAS  PubMed  Google Scholar 

  5. Pyron RA, Burbrink FT. Can the tropical conservatism hypothesis explain temperate species richness patterns? An inverse latitudinal biodiversity gradient in the New World snake tribe Lampropeltini. Global Ecol Biogeogr. 2009;18(4):406–15.

    Article  Google Scholar 

  6. Cerezer FO, Machac A, Rangel TF, Dambros CS. Exceptions to the rule: relative roles of time, diversification rates and regional energy in shaping the inverse latitudinal diversity gradient. Global Ecol Biogeogr. 2022;31(9):1794–809.

    Article  Google Scholar 

  7. Chown SL, Gaston KJ, Williams PH. Global patterns in species richness of pelagic seabirds: the Procellariiformes. Ecography. 1998;21(4):342–50.

    Article  Google Scholar 

  8. Yang Y, Bian Z, Ren G, Liu J, Shrestha N. Niche conservatism limits the distribution of Medicago in the tropics. Ecography. 2022;2022(9):e06085.

    Article  Google Scholar 

  9. Yan H-F, Zhang C-Y, Anderberg AA, Hao G, Ge X-J, Wiens JJ. What explains high plant richness in East Asia? Time and diversification in the tribe Lysimachieae (Primulaceae). New Phytol. 2018;219(1):436–48.

    Article  PubMed  Google Scholar 

  10. Tian Q, Stull GW, Kellermann J, Medan D, Nge FJ, Liu S-Y, et al. Rapid in situ diversification rates in Rhamnaceae explain the parallel evolution of high diversity in temperate biomes from global to local scales. New Phytol. 2024;241(4):1851–65.

    Article  PubMed  Google Scholar 

  11. Jin WT, Gernandt DS, Wehenkel C, Xia XM, Wei XX, Wang XQ. Phylogenomic and ecological analyses reveal the spatiotemporal evolution of global pines. Proc Natl Acad Sci U S A. 2021;118(20):e2022302118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Liu J, Lindstrom AJ, Nagalingum NS, Wiens JJ, Gong X. Testing the causes of richness patterns in the paleotropics: time and diversification in cycads (Cycadaceae). Ecography. 2021;44(11):1606–18.

    Article  Google Scholar 

  13. Stephens PR, Wiens JJ. Explaining species richness from continents to communities: the time-for-speciation effect in emydid turtles. Am Nat. 2003;161(1):112–28.

    Article  PubMed  Google Scholar 

  14. Mittelbach GG, Schemske DW, Cornell HV, Allen AP, Brown JM, Bush MB, et al. Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography. Ecol Lett. 2007;10(4):315–31.

    Article  PubMed  Google Scholar 

  15. Wiens JJ, Donoghue MJ. Historical biogeography, ecology and species richness. Trends Ecol Evol. 2004;19(12):639–44.

    Article  PubMed  Google Scholar 

  16. Li H, Wiens JJ. Time explains regional richness patterns within glades more often than diversification rates or area. Am Nat. 2019;193(4):514–29.

    Article  PubMed  Google Scholar 

  17. Shrestha N, Wang Z, Su X, Xu X, Lyu L, Liu Y, et al. Global patterns of Rhododendron diversity: the role of evolutionary time and diversification rates. Global Ecol Biogeogr. 2018;27(8):913–24.

    Article  Google Scholar 

  18. Wiens JJ, Graham CH. Niche conservatism: integrating evolution, ecology, and conservation biology. Annu Rev Ecol Evol Syst. 2005;36:519–39.

    Article  Google Scholar 

  19. Wiens JJ, Ackerly DD, Allen AP, Anacker BL, Buckley LB, Cornell HV, et al. Niche conservatism as an emerging principle in ecology and conservation biology. Ecol Lett. 2010;13(10):1310–24.

    Article  PubMed  Google Scholar 

  20. Wolf T. Monographie Der Gattung Potentilla. Bibl Bot. 1908;16(71):1–715.

    Google Scholar 

  21. Juzepczuk SV. Potentilla L. In: Komarov VL, editor. Flora URSS. Volume 10. & Leningrad: Izdatel’stvo Akademii Nauk SSSR; 1941. pp. 78–223. Moscow.

    Google Scholar 

  22. Ball PW, Pawlowski B, Walters SM. Potentilla L. In: Tutin TG, Heywood VH, Burges NA, Moore DM, Valentine DH, Walters SM, et al. editors. Flora Europaea. Volume 2. Cambridge: Cambridge University Press; 1968. pp. 36–47.

    Google Scholar 

  23. Li CL, Ikeda H, Ohba H. Potentilla Linnaeus. In: Wu ZY, Raven PH, Hong DY, editors. Flora of China. Volume 9. Beijing: Science; 2003. pp. 291–327.

    Google Scholar 

  24. Ertter B, Elven R, Reveal JL, Murray DF. Potentilla Linnaeus. Flora of North America. Volume 9. Oxford: Oxford University Press; 2014. pp. 121–218. Flora of North America Editorial Committee.

    Google Scholar 

  25. Ertter B, Reveal JL. Ivesia Torrey and A. Gray, Horkelia Chamisso and Schlechtendal, Horkeliella Rydberg, Duchesnea Smith. In: Flora of North America Editorial Committee, editor. Flora of North America (vol. 9). Oxford: Oxford University Press; 2014. pp.219–274.

  26. Dobeš C, Paule J. A comprehensive chloroplast DNA-based phylogeny of the genus Potentilla (Rosaceae): implications for its geographic origin, phylogeography and generic circumscription. Mol Phylogen Evol. 2010;56(1):156–75.

    Article  Google Scholar 

  27. Feng T. Systematics and Biogeography of Potentilleae (Rosaceae). Ph.D. Thesis. China: University of Chinese Academy of Sciences; 2017.

  28. Feng T, Moore MJ, Yan MH, Sun YX, Zhang HJ, Meng AP, et al. Phylogenetic study of the tribe Potentilleae (Rosaceae), with further insight into the disintegration of Sibbaldia. J Syst Evol. 2017;55(3):177–91.

    Article  Google Scholar 

  29. Eriksson T, Persson NL, Smedmark JEE. What is Potentilla? A phylogeny-based taxonomy for Potentillinae (Rosaceae). Taxon. 2022;71(3):493–505.

    Article  Google Scholar 

  30. Zhang X, Landis JB, Sun Y, Zhang H, Lin N, Kuang T et al. Macroevolutionary pattern of Saussurea (Asteraceae) provides insights into the drivers of radiating diversification. Proc R Soc B-Biol Sci. 2021;288(1962):20211575.

  31. Ding WN, Ree RH, Spicer RA, Xing YW. Ancient orogenic and monsoon-driven assembly of the world’s richest temperate alpine flora. Science. 2020;369(6503):578–81.

    Article  CAS  PubMed  Google Scholar 

  32. Wang W, Xiang XG, Xiang KL, Ortiz RC, Jabbour F, Chen ZD. A dated phylogeny of Lardizabalaceae reveals an unusual long-distance dispersal across the Pacific Ocean and the rapid rise of east Asian subtropical evergreen broadleaved forests in the late Miocene. Cladistics. 2020;36(5):447–57.

    Article  PubMed  Google Scholar 

  33. Han TS, Zheng QJ, Onstein RE, Rojas-Andres BM, Hauenschild F, Muellner-Riehl AN, et al. Polyploidy promotes species diversification of Allium through ecological shifts. New Phytol. 2020;225(1):571–83.

    Article  CAS  PubMed  Google Scholar 

  34. Eriksson T, Donoghue MJ, Hibbs MS. Phylogenetic analysis of Potentilla using DNA sequences of nuclear ribosomal internal transcribed spacers (ITS), and implications for the classification of Rosoideae (Rosaceae). Plant Syst Evol. 1998;211(3–4):155–79.

    Article  CAS  Google Scholar 

  35. Töpel M, Lundberg M, Eriksson T, Eriksen B. Molecular data and ploidal levels indicate several putative allopolyploidization events in the genus Potentilla (Rosaceae). PLoS Curr. 2011;3:RRN1237.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Eriksson T, Lundberg M, Töpel M, Östensson P, Smedmark JEE. Sibbaldia: a molecular phylogenetic study of a remarkably polyphyletic genus in Rosaceae. Plant Syst Evol. 2015;301(1):171–84.

    Article  Google Scholar 

  37. Xue T, Janssens SB, Liu B, Yu S. Phylogenomic conflict analyses of the plastid and mitochondrial genomes via deep genome skimming highlight their independent evolutionary histories: a case study in the cinquefoil genus Potentilla Sensu lato (Potentilleae, Rosaceae). Mol Phylogen Evol. 2024;190.

  38. Li Q-Q, Khasbagan, Zhang Z-P, Wen J, Yu Y. Plastid phylogenomics of the tribe potentilleae (Rosaceae). Mol Phylogen Evol. 2024;190.

  39. Spicer RA, Farnsworth A, Su T. Cenozoic topography, monsoons and biodiversity conservation within the Tibetan Region: an evolving story. Plant Divers. 2020;42(4):229–54.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Ren CQ, Zhang DQ, Liu XY, Zhang JQ. Genomic data provide a robust phylogeny backbone for Rhodiola L. (Crassulaceae) and reveal extensive reticulate evolution during its rapid radiation. Mol Phylogen Evol. 2023;186:107863.

    Article  CAS  Google Scholar 

  41. Ebersbach J, Muellner-Riehl AN, Michalak I, Tkach N, Hoffmann MH, Roeser M, et al. In and out of the Qinghai-Tibet Plateau: divergence time estimation and historical biogeography of the large arctic-alpine genus Saxifraga L. J Biogeogr. 2017;44(4):900–10.

    Article  Google Scholar 

  42. Ding L, Kapp P, Cai F, Garzione CN, Xiong Z, Wang H, et al. Timing and mechanisms of Tibetan Plateau uplift. Nat Rev Earth Environ. 2022;3(10):652–67.

    Article  Google Scholar 

  43. Favre A, Michalak I, Chen CH, Wang JC, Pringle JS, Matuszak S, et al. Out-of-Tibet: the spatio-temporal evolution of Gentiana (Gentianaceae). J Biogeogr. 2016;43(10):1967–78.

    Article  Google Scholar 

  44. Tiffney BH. The Eocene North Atlantic land bridge: its importance in Tertiary and modern phytogeography of the Northern Hemisphere. J Arnold Arbor. 1985;66(2):243–73.

    Article  Google Scholar 

  45. Denk T, Grimsson F, Zetter R. Episodic migration of oaks to Iceland: evidence for a North Atlantic land bridge in the latest Miocene. Am J Bot. 2010;97(2):276–87.

    Article  PubMed  Google Scholar 

  46. Graham A. The role of land bridges, ancient environments, and migrations in the assembly of the north American flora. J Syst Evol. 2018;56(5):405–29.

    Article  Google Scholar 

  47. Töpel M, Antonelli A, Yesson C, Eriksen B. Past climate change and plant evolution in western North America: a case study in Rosaceae. PLoS ONE. 2012;7(12):e50358.

    Article  PubMed  PubMed Central  Google Scholar 

  48. De Storme N, Geelen D. The impact of environmental stress on male reproductive development in plants: biological processes and molecular mechanisms. Plant Cell Environ. 2014;37(1):1–18.

    Article  PubMed  Google Scholar 

  49. Cai L, Xi Z, Amorim AM, Sugumaran M, Rest JS, Liu L, et al. Widespread ancient whole-genome duplications in Malpighiales coincide with Eocene global climatic upheaval. New Phytol. 2019;221(1):565–76.

    Article  CAS  PubMed  Google Scholar 

  50. Wang KL, Deng PR, Yao Z, Dong JY, He Z, Yang P, et al. Biogeographic patterns of polyploid species for the angiosperm flora in China. J Syst Evol. 2022;61:776–89.

    Article  Google Scholar 

  51. Jonasson S, Callaghan TV. Root mechanical-properties related to disturbed and stressed habitats in the Arctic. New Phytol. 1992;122(1):179–86.

    Article  PubMed  Google Scholar 

  52. Körner C. Alpine plant life. Cham, Switzerland: Springer Nature Switzerland; 2021.

    Google Scholar 

  53. Thomas DN. Arctic Ecology. West Sussex, UK: John Wiley & Sons Ltd.; 2021.

    Book  Google Scholar 

  54. Bennett KD. Is the number of species on earth increasing or decreasing? Time, chaos and the origin of species. Palaeontology. 2013;56(6):1305–25.

    Article  Google Scholar 

  55. Kadereit JW, Abbott RJ. Plant speciation in the Quaternary. Plant Ecol Divers. 2021;14(3–4):105–42.

    Article  Google Scholar 

  56. Bennett JM, Sunday J, Calosi P, Villalobos F, Martinez B, Molina-Venegas R et al. The evolution of critical thermal limits of life on Earth. Nat Commun. 2021;12(1).

  57. Augustynowicz D, Latte KP, Tomczyk M. Recent phytochemical and pharmacological advances in the genus Potentilla L. Sensu lato - an update covering the period from 2009 to 2020. J Ethnopharmacol. 2021;266.

  58. Kerkhoff AJ, Moriarty PE, Weiser MD. The latitudinal species richness gradient in New World woody angiosperms is consistent with the tropical conservatism hypothesis. Proc Natl Acad Sci U S A. 2014;111(22):8125–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Su XY, Shrestha N, Xu XT, Sandanov D, Wang QG, Wang SY, et al. Phylogenetic conservatism and biogeographic affinity influence woody plant species richness-climate relationships in eastern Eurasia. Ecography. 2020;43(7):1027–40.

    Article  Google Scholar 

  60. Wang Q, Wu S, Su X, Zhang L, Xu X, Lyu L, et al. Niche conservatism and elevated diversification shape species diversity in drylands: evidence from Zygophyllaceae. Proc R Soc B-Biol Sci. 2018;285(1890):20181742.

    Article  Google Scholar 

  61. Persson N, Toresen I, Andersen HL, Smedmark JEE, Eriksson T. Detecting destabilizing species in the phylogenetic backbone of Potentilla (Rosaceae) using low-copy nuclear markers. Aob Plants. 2020;12(3):plaa017.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987;19:11–5.

    Google Scholar 

  63. Chen SF, Zhou YQ, Chen YR, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):884–90.

    Article  Google Scholar 

  64. Jin JJ, Yu WB, Yang JB, Song Y, dePamphilis CW, Yi TS, et al. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020;21:241.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Capella-Gutiérrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25(15):1972–3.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics. 2019;35(21):4453–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Kalyaanamoorthy S, Bui Quang M, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Smith SA, O’Meara BC. treePL: divergence time estimation using penalized likelihood for large phylogenies. Bioinformatics. 2012;28(20):2689–90.

    Article  CAS  PubMed  Google Scholar 

  70. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Zhang SD, Jin JJ, Chen SY, Chase MW, Soltis DE, Li HT, et al. Diversification of Rosaceae since the late cretaceous based on plastid phylogenomics. New Phytol. 2017;214(3):1355–67.

    Article  CAS  PubMed  Google Scholar 

  72. Hollick CA. The tertiary floras of Alaska. Washington DC: US Goverment Printing Office; 1936.

    Google Scholar 

  73. Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3(2):217–23.

    Article  Google Scholar 

  74. Rabosky DL, Chang J, Title PO, Cowman PF, Sallan L, Friedman M, et al. An inverse latitudinal gradient in speciation rate for marine fishes. Nature. 2018;559(7714):392–5.

    Article  CAS  PubMed  Google Scholar 

  75. Li XQ, Xiang XG, Zhang Q, Jabbour F, Ortiz RdC, Erst AS, et al. Immigration dynamics of tropical and subtropical southeast Asian limestone karst floras. Proc R Soc B-Biol Sci. 2022;289(1966):20211308.

    Article  CAS  Google Scholar 

  76. Zhou P, Li JH, Liu YZ, Zhu ZW, Luo Y, Xiang XG. Species richness disparity in tropical terrestrial herbaceous floras: evolutionary insight from Collabieae (Orchidaceae). Mol Phylogen Evol. 2023;186:107860.

    Article  Google Scholar 

  77. Matzke NJ. Probabilistic historical biogeography: new models for founder-event speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Front Biogeogr. 2013;5(4):242–8.

    Article  Google Scholar 

  78. Yu Y, Blair C, He XJ. RASP 4: ancestral state reconstruction tool for multiple genes and characters. Mol Biol Evol. 2020;37(2):604–6.

    Article  CAS  PubMed  Google Scholar 

  79. Rabosky DL, Grundler M, Anderson C, Title P, Shi JJ, Brown JW, et al. BAMMtools: an R package for the analysis of evolutionary dynamics on phylogenetic trees. Methods Ecol Evol. 2014;5(7):701–7.

    Article  Google Scholar 

  80. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2022. https://www.R-project.org.

    Google Scholar 

  81. Morlon H, Lewitus E, Condamine FL, Manceau M, Clavel J, Drury J. RPANDA: an R package for macroevolutionary analyses on phylogenetic trees. Methods Ecol Evol. 2016;7(5):589–97.

    Article  Google Scholar 

  82. Sun M, Folk RA, Gitzendanner MA, Soltis PS, Chen Z, Soltis DE, et al. Recent accelerated diversification in rosids occurred outside the tropics. Nat Commun. 2020;11(1):3333.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Beaulieu JM, O’Meara BC. Detecting hidden diversification shifts in models of trait-dependent speciation and extinction. Syst Biol. 2016;65(4):583–601.

    Article  PubMed  Google Scholar 

  84. FitzJohn RG. Diversitree: comparative phylogenetic analyses of diversification in R. Methods Ecol Evol. 2012;3(6):1084–92.

    Article  Google Scholar 

  85. Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25(1):1–18.

    Article  Google Scholar 

  86. Kassambara A, Mundt F. Factoextra: extract and visualize the results of multivariate data analyses. https://CRAN.R-project.org/package=factoextra. 2020.

  87. Zizka A, Antonelli A, Silvestro D. Sampbias, a method for quantifying geographic sampling biases in species distribution data. Ecography. 2021;44(1):25–32.

    Article  Google Scholar 

  88. Westerhold T, Marwan N, Drury AJ, Liebrand D, Agnini C, Anagnostou E, et al. An astronomically dated record of Earth’s climate and its predictability over the last 66 million years. Science. 2020;369(6509):1383–7.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Weibin Xu (Guangxi Institute of Botany) for field work assistance and sample collection, Ran Wei, Xiaoqian Li, Guanlong Cao (Institute of Botany, Chinese Academy of Sciences), and Peng Zhou (Nanchang University) for assistance with data analysis.

Funding

This work was supported by the National Natural Science Foundation of China [32071654 to S.X.Y].

Author information

Authors and Affiliations

Authors

Contributions

S.Y. and S.J. conceived the project and designed the research. S.Y, T.X., T.F. and X.Y. collected the samples. T.X., Y.L., X.Y., F.Q. and J.Y. conducted the analyses. T.X. drafted the manuscript. S.Y., S.J. and T.F. revised the manuscript. All authors read and improved the manuscript.

Corresponding authors

Correspondence to Steven B. Janssens or Shengxiang Yu.

Ethics declarations

Ethics approval and consent to participate

The collection of all samples completely complies with institutional, national, and local guidelines and legislation. All plant species are not rare species for protection. We have the permission to collect the plant samples by the local government. The collecting of all samples in this study followed the Regulations on the Protection of Wild Plants of China, the IUCN Policy Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xue, T., Feng, T., Liang, Y. et al. Radiating diversification and niche conservatism jointly shape the inverse latitudinal diversity gradient of Potentilla L. (Rosaceae). BMC Plant Biol 24, 443 (2024). https://doi.org/10.1186/s12870-024-05083-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05083-8

Keywords