Skip to main content

Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense

Abstract

Background

Relict species are important for enhancing the understanding of modern biogeographic distribution patterns. Although both geological and climatic changes since the Cenozoic have affected the relict flora in East Asia, the contributions of geographical processes remain unclear. In this study, we employed restriction-site associated DNA sequencing (RAD-seq) and shallow genome sequencing data, in conjunction with ecological niche modeling (ENM), to investigate the spatial genetic patterns and population differentiation history of the relict species Rehderodendron kwangtungense Chun.

Results

A total of 138 individuals from 16 populations were collected, largely covering the natural distribution of R. kwangtungense. The genetic diversity within the R. kwangtungense populations was extremely low (HO = 0.048 ± 0.019; HE = 0.033 ± 0.011). Mantel tests revealed isolation-by-distance pattern (R2 = 0.38, P < 0.001), and AMOVA analysis showed that the genetic variation of R. kwangtungense occurs mainly between populations (86.88%, K = 7). Between 23 and 21 Ma, R. kwangtungense underwent a period of rapid differentiation that coincided with the rise of the Himalayas and the establishment of the East Asian monsoon. According to ENM and population demographic history, the suitable area and effective population size of R. kwangtungense decreased sharply during the glacial period and expanded after the last glacial maximum (LGM).

Conclusion

Our study shows that the distribution pattern of southern China mountain relict flora may have developed during the panplain stage between the middle Oligocene and the early Miocene. Then, the flora later fragmented under the force of orogenesis, including intermittent uplift during the Cenozoic Himalayan orogeny and the formation of abundant rainfall associated with the East Asian monsoon. The findings emphasized the predominant role of geographical processes in shaping relict plant distribution patterns.

Peer Review reports

Introduction

Relict species are considered fascinating ‘living fossils’ or remnants of ancient clades and biotas and can be divided into geographical relicts and phylogenetic relicts [1, 2]. These organisms are often viewed as providing compelling evidence for the conservation of ancestral character states in terms of morphology, ecology [3, 4], and spatial distribution [5, 6]. Relict species are often confined to ecological refugia and have fragmented distributions, and they are commonly high priority targets for biodiversity conservation [7]. They must adjust their physiological adaptability [8] and purge the accumulation of deleterious mutations within the population to quickly respond to the pressure of environmental changes [9]. Thus, relict species serve as a “window to the past”, allowing us to understand the conditions that enabled them to survive for so long [1, 10, 11].

East Asia contains the most abundant Tertiary relict flora in the Northern Hemisphere [12,13,14]. Southeast China is an important host of such relict flora and is considered a floristic museum characterised by phylogenetic and distribution centres of ancient relict species [15, 16]. In general, climatic fluctuations during the Cenozoic prompted the southward migration of Northern Hemisphere flora [17]. High topographic heterogeneity and climate stability has allowed East Asia to serve as a long-term stable habitat, and many species remained in situ during the Pleistocene glaciation [18,19,20,21,22]. Many phylogeographic studies have indicated that plants can survive the climatic fluctuations associated with an ice age and take refuge in mountains in situ [23,24,25,26]. It is clear that the establishment of relict flora in East Asia occurred much earlier than the Pleistocene. Previous studies have shown that past geographic events and climate changes have jointly influenced the lineage differentiation of relict plants in East Asia [13, 27,28,29]. However, the distribution ranges of relict species prior to fragmentation have not been explored, and the independent effects of climate and geography on relict patterns still need further study.

In fact, a topographic reversal has occurred in South China during the Cenozoic era, and the implications on biota of this phenomenon may have been previously overlooked. Before the late Eocene, the terrain of South China was high along the eastern coast of Cathaysia and low in the western Yangtze interior [30,31,32]. Several lacustrine sedimentary basin, such as the Danxia Basin, Nanxiong Basin, Heyuan Basin, and Guiping Basin, formed in the interior under hot and arid climates [26, 33]. During the early Oligocene, the uplift of the SE Tibetan Plateau triggered a topographic reversal in the South China Block [34,35,36,37,38]. This period was also the key period for the transformation of the palaeo-Pearl River, which flowed from east to west, to the present Pearl River, which flows from west to east [39, 40]. Although, the exact time frame during which the modern Pearl River formed, such as during the early Miocene [41], late Oligocene [39], or  30 Ma [42], is debated. There is no doubt that a peneplain existed between the Cathaysian coast (current coastal mountains of South China) and the Yangtze interior between the Oligocene and Miocene [43]. This peneplain can still be observed as the present mountain planation surface in the Nanling Mountains, specifically at elevations of approximately 1000–1350 m in northern Guangdong [44, 45] and 1416–1780 m in northern Guangxi and southern Hunan [46]. A large number of relict plants appeared in the high altitudes of the Nanling Mountains, e.g. Fokienia hodginsii [22], Bretschneidera sinensis [25]. The occurrence of this continuous relict flora is highly improbable to be attributed to random chance. Thus, we suspect that the peneplain stage between the Oligocene and Miocene in South China was important for establishing continuous warm East Asian relict flora.

Rehderodendron Hu (Styracaceae) is a typical relict genus, whose fossil fruits and pollen are widely known in Europe from the early Eocene to the Pliocene [47,48,49,50,51]. However, the modern species of Rehderodendron are restricted mostly to warm tropical and subtropical climates in the mountains of East Asia [52]. Our previous phylogenetic evidence indicated that R. kwangtungense Chun is the most ancient species within Rehderodendron (unpublished data). This fact is further demonstrated by the similarities between modern pollen and pollen discovered in the lower Eocene strata of England [48]. R. kwangtungense is sporadically distributed in fog-dependent forests in the subtropical mountains of southern China, exhibiting an island-like (near-mountaintop, mountain planation surface) distribution pattern (Fig. 1a). This island-like distribution pattern is widely observed in East Asian relict species but remains poorly studied. R. kwangtungense is an ancient species with low dispersal capacity, making it an ideal species for tracing the spatiotemporal dynamic history of relict flora in East Asia.

Here, we examined the genetic diversity and population divergence of R. kwangtungense by employing RAD-seq data [53]. Then, we used ENM to identify the climatic niche of this species under current environmental conditions and to predict its potential distribution during the LGM and mid-Holocene. We combined this information to assess the formation of the modern distribution pattern of R. kwangtungense. Our specific objectives were to: (1) whether the distribution patterns of R. kwangtungense formed in the panplain stage in the Oligocene, (2) reveal the population structure and population demographic history, and (3) conservation genetics research on R. kwangtungense. This study is expected to provide new understanding of East Asian relict flora.

Methods

The studied species, sample collection and DNA extraction

Rehderodendron kwangtungense is a 5–15 m tall arbor with a diameter up to 30 cm at breast height. Its corolla is white, bell-shaped and fragrant (Fig. 1b), and its fruits are cylindrical ellipsoid in shape (Fig. 1c). The woody exocarp is approximately 1 mm thick, the fibrous corky mesocarp is 8–12 mm thick, and the woody endocarp is hard and radiates into the mesocarp (Fig. 1d). The tillering ability of R. kwangtungense is relatively strong. On the one hand, when the trunk of a tree is broken, the base can branch and grow into a larger tree; On the other hand, after the main trunk falls down and touches the soil, it is easy to take root and clone many trees. During field investigation, we found two major distribution areas of R. kwangtungense with obvious phenological differences, one on the eastern margin of the Yunnan-Guizhou Plateau (northwestern populations, P1-7; Fig. 1a) and the other in the Nanling and southern mountains (southeastern populations, P8-16; Fig. 1a). In March and April, when the southeastern populations are blooming, the northwestern populations are still in the early bud stage.

Fresh leaves were collected from a total of 138 individuals (16 populations) of R. kwangtungense, covering most of its geographic range (Table 1; Fig. 1a). The geographical information of the populations was recorded using a Garmin GPS unit (GPSMAP 62sc, Shanghai). Voucher specimens were deposited at the Herbarium of Sun Yat-sen University (SYS). The fresh leaves were dried and deposited with silica gel in sealed bags. Genomic DNA was isolated using the modified cetyl trimethylammonium bromide (CTAB) method [54].

Fig. 1
figure 1

Geographic distribution and morphology of R. kwangtungense. (a) The white dots indicate the 16 sampling sites in this study. (b) Flowering branch. (c) Fruiting branch. (d) Transverse section of a fruit

RAD-seq data processing

RAD-seq library preparation and Illumina sequencing were performed by Novogene Bioinformatics Technology Co. Ltd. (Tianjin, China). The standard protocol was followed. Stacks was developed to work with RAD-seq for the purpose of building genetic maps and conducting population genomics and phylogeographic analyses. Therefore, we processed the RAD-Seq data using Stacks v2.55 [55]. The program “process_radtags” was used to decomplex reads for each sample, and “denovo_map” was used to determine the optimum values for M (number of mismatches allowed between stacks within individuals) and n (number of mismatches allowed between stacks among individuals) with subset data containing 12 samples. Then, “denovo_map” was reused to process all 138 samples with the optimum M and n values (M = n = 4). Finally, the results were filtered using the program “populations”, in which polymorphic RAD loci that were present in at least 40% of the individuals across populations, had a minor allele frequency > 0.05 and had an observed heterozygosity < 0.7 were retained (-R 0.4 --min-maf 0.05 --max-obs-het 0.7). And for each locus, unlinked SNP were retained using the command “--write-single-snp”. To detect non-neutral loci under selection, Tajima’s D was calculated using VCFtools [56], in which the window size was set as 3,000 bp. Loci with Tajima’s D value < − 1.771 or > 2.086 were deemed as non-neutral [57]. The produced vcf files (Supplementary Material S1) were used for downstream analyses, and converted to other formats using VCFtools.

Genetic diversity and population structure analysis

The program “populations” of Stacks was used to calculate genetic diversity parameters (HO, HE, π and FIS) and pairwise fixation index (FST) [55]. Population structure was examined using Admixture v1.3.0 [58]. The value of the hypothetical ancestral population (K) was set from 1 to 10, the optimal K was determined by the smallest cross-validation (CV) error value, and Q estimates were used as a proxy for ancestry fractions. Principal component analysis (PCA) was performed using Plink v1.90b to evaluate population genetic variation [59]. Analysis of molecular variance (AMOVA) was performed to estimate standardized genetic differentiation using Arlequin v3.5 [60], in which the 16 populations of R. kwangtungense were divided into two regions (northwestern and southeastern regions according to variation in flowering time), four regions according to PCA, and seven regions according to Admixture analysis. The python script “vcf2phylip.py” (https://github.com/edgardomortiz/vcf2phylip) was used to convert the vcf file to the phylip format, and the maximum likelihood phylogenetic tree was produced with IQ-TREE v2.2 [61] and plotted with iTOL v6 [62].

Population historical dynamics inference

The site frequency spectra (SFS) were generated using easySFS (https://github.com/isaacovercast/easySFS), and the value for projecting down samples was set to 112 to maximize the number of segregating sites. Stairway Plot v2.1.1 [63] was subsequently used to infer demographic history based on the folded SFS. Stairway Plot is based on SFS, which does not require whole-genome sequence data or a reference genome to infer recent population size changes [64]. The mutation rate of R. kwangtungense was set to 2.29 e-9 per site per year according to the estimate for walnut by Luo et al. [65]. The generation time was set to ten years, which was observed in the field investigation.

Shallow genome sequencing and plastid genome, ribosomal cistron assembly

For each population, we randomly selected one individual for shallow genome sequencing. In addition, 5 other species of Rehderodendron have undergone shallow sequencing and were utilized in the construction of the phylogenetic tree. The extracted DNA was sent to JieRui BioScience Co. Ltd. (Guangzhou, China) for library construction and Illumina sequencing on the Illumina 2000 platform. The program Novoplasty v2.7.2 [66] was used to assemble the chloroplast genome, using the R. macrocarpum chloroplast genome (GenBank accession MG719844) as a reference and its chloroplast rbcl gene sequence as a seed. Meanwhile, GetOrganellle v1.7.7 [67] was used to assemble the ribosomal cistron, using the Melliodendron xylocarpum ribosomal cistron (GenBank accession MF171073) as a seed.

Phylogenetic tree construction and divergence time estimation

According to Yan et al. [68], the plastid genomes of 6 related genera were downloaded from the NCBI nucleotide database (Supplementary Table 1). Together with the 16 newly assembled plastid genomes of R. kwangtungense and 5 other species of Rehderodendron in this study, 31 plastid genome sequences were aligned using MAFFT v7.508 [69]. Gaps and ambiguous bases (N) were removed using MEGA 7 [70]. A similar method was used to standardize ribosomal DNA data (Supplementary Table 2). The processed sequences were subsequently used to construct a phylogenetic tree, and the divergence times were inferred via Beast v1.8.4 [71]. The tree was calibrated with three fossils, Rehderodendron stonei (52 − 49 Ma) [49], Pterostyrax coronatus (33.9–28.1 Ma) [72] and Halesia reticulata (37.2–33.9 Ma) [73]. The lognormal prior distributions were enforced with a mean = 1.5, a Stdev = 2.0 and an offset set so that 95% of the distribution would fit the age of the upper (youngest) stratum from which each fossil was described. A GTR + I + F + G4 substitution model and lognormal relaxed-clock model were used to allow rate variation among branches. Additionally, a constant size was chosen as the tree prior. The analysis was run for 400 million Markov Chain Monte Carlo (MCMC) steps. MCMC samples were imported into TRACER v1.5 (available from http://beast.bio.ed.ac.uk/Tracer) to inspect the sampling adequacy and convergence of the chains to a stationary distribution. The final plot was visualized using Figtree v1.4.4.

Ecological niche modelling

The ecological niche model (ENM) was used to predict the potential distribution of R. kwangtungense. In addition to our sample locations, geographic distributions of the species were also obtained from the Chinese virtual herbarium portal (CVH, https://www.cvh.ac.cn/). After filtering duplicate records and invalid records without a specific location, a total of 51 records were used in the following step. Climate layer information for 19 bioclimatic variables was downloaded from WorldClim [74, 75] at a resolution of 2.5 arc-minutes. Based on geographical distribution data from R. kwangtungense and 19 climatic factors, a MaxEnt model was initially established to evaluate the contributions of different variables to the distribution of climate data [76]. To avoid overfitting of the simulation results due to the mutual influence between highly correlated climatic variables, Pearson correlation coefficient analysis was performed among the 19 climatic variables using SPSS v26.0. A pair of climatic variables with correlation coefficients greater than |0.8| was considered non-independent. Finally, for each pair of significantly correlated variables, only one variable with a large contribution was used for projection. After the selection procedure, 7 variables were retained for modelling the distribution of R. kwangtungense: the mean diurnal range (Bio2), isothermality (Bio3), minimum temperature of coldest month (Bio6), temperature annual range (Bio7), annual precipitation (Bio12), precipitation of driest quarter (Bio17), and precipitation of coldest quarter (Bio19). We employed MaxEnt 3.4.4 to simulate and predict the evolution of potentially suitable areas of R. kwangtungense in different periods, using default settings and specific parameters: random test (25%), training (75%), regularization multiplier (1), maximum iterations (5000), convergence threshold (0.00001), maximum background points (10,000), and 10 bootstrap replications. Model accuracy was gauged using the AUC of receiver operating characteristic (ROC). The output file was imported into ArcGIS 10.8, and the suitability was manually divided at 0.05 intervals, and areas with suitability values less than 0.1 were considered unsuitable areas.

Correlation of genetic structure with geographic isolation and environment factor

The isolation-by‐distance (IBD) reveals that genetic differentiation among populations increases with geographic distance, while the isolation‐by‐environment (IBE) assumes a linear relationship between genetic variation and environmental differences among populations [77]. To explore the existence of IBD and IBE, we applied Mantel test to assess the relationships between genetic distance and geographic or environmental distance with R v4.3.2 (https://www.r-project.org/). The pairwise genetic distance among populations was calculated by FST/(1 − FST), the geographic distance matrix was transformed by the geographic coordinates of populations, and 7 climatic variables used in ENM were matrix transformed using PCA to obtain the environmental distance. A partial mantel test was conducted to control the influence of geographical and environmental distance, respectively.

Results

RAD data processing and genetic diversity

A total of 11,190,163 loci were genotyped by the “denovo_map” program, and the mean, minimum and maximum values of effective per-sample coverage were 17.1×, 10.1× and 23.7×, respectively. After data filtration, a total of 1,943 SNPs remained for all subsequent analyses. The Tajima’s D values of all 1943 variant sites ranged from − 0.394 to 2.029, all sites are neutral. The observed heterozygosity (HO) for the 16 populations ranged from 0.022 to 0.078, with an average value of 0.048 ± 0.019; the expected heterozygosity (HE) ranged from 0.019 to 0.052, with an average value of 0.033 ± 0.011; and the inbreeding coefficient (FIS) was close to 0 (-0.0009 to 0.0005) (Table 1).

Table 1 Geographic information and genetic diversity parameters of R. kwangtungense

Population structure

Based on the results of the Admixture program, the CV error decreases and becomes stable at K = 7 (Supplementary Fig. 1). When K = 2, the populations of R. kwangtungense were divided into northwestern group (P1-7) and southeastern group (P8-16) (Fig. 2a), which was consistent with phenological differences. When K = 7, the populations were subdivided gradually, with the northwestern cluster divided into four groups (G1-4) and the southeastern cluster divided into three groups (G5-7) (Fig. 2b).

Fig. 2
figure 2

Delimitation of 16 populations of R. kwangtungense based on Admixture analysis. (a) Genetic structure of R. kwangtungense based on K = 2. (b) Genetic structure of R. kwangtungense based on K = 7. Each individual (indicated as columns along the X-axis) is probabilistically assigned (probability of assignment q on the Y-axis) to one of the inferred genetic clusters

PCA showed that the contribution rates of the first and second principal components were 21.57% and 16.64%, respectively. The two-dimensional plot of PC1 and PC2 divided all the 138 individuals into four large clusters (Fig. 3a): the seven populations located in the northwestern area were ascribed to one cluster, including the four groups G1-G4 according to Admixture analysis (Fig. 2b), while the nine populations in the southeastern area were divided into three clusters, corresponding to the three groups G5-G7 according to Admixture analysis (Fig. 2b). The unrooted phylogenetic tree based on population data also confirmed that individuals of most group (G1-G6) clustered into a monophyletic group with a high support value (Fig. 3b). All terminal branches are very short in length (Fig. 3b), indicating little genetic variation within the population.

Fig. 3
figure 3

Genetic and phylogenetic analysis of R. kwangtungense. (a) Two-dimensional clustering of individuals obtained via principal component analysis. (b) Phylogenetic tree constructed from RAD-seq data

AMOVA showed that the genetic variation among regions increased from 27.41 to 59.72% and that the variation among populations within groups decreased from 60.41 to 27.16% as the K value increased from 2 to 7, and the variation within populations remained stable (12.18-13.12%, Table 2).

Table 2 Analysis of molecular variance (AMOVA) of R. kwangtungense

IBD and IBE

The mantel test showed that genetic distance [FST/(1–FST)] was significantly correlated with geographical distance and conformed to an isolation-by-distance pattern (R2 = 0.38, P < 0.001, Fig. 4a). Meanwhile, environmental factors also affected the genetic variation (R2 = 0.21, P < 0.001, Fig. 4b). The partial mantel test showed that after controlling the influence of environment matrix, geographical distance and genetic distance were still related (R2 = 0.24, P < 0.001). However, after controlling for the influence of geographical matrix, there was no correlation between environmental factors and genetic distance (R2 = 0.018, P = 0.881). Therefore, geographical isolation rather than environmental factors was the main factor that caused the genetic differences of R. kwangtungense populations.

Fig. 4
figure 4

Isolation-by-distance (IBD) and Isolation-by-environment (IBE) patterns for populations of Rehderodendron kwangtungense. (a) Mantel test result of genetic distance and geographical. (b) Mantel test result of genetic distance and environment distance. Each dot represents a pair of populations

Population demographic history

A demographic history investigation showed that the effective population size of R. kwangtungense increased quickly from approximately 16,000 years ago (YA) to 10,000 YA and then remained stable above 2000 until 8,00 YA. It subsequently gradually declined to the very low level observed at present (approximately 200) (Fig. 5).

Fig. 5
figure 5

Estimates of the effective population size of R. kwangtungense

Divergence history of populations inferred by phylogenetic tree

The plastid and ribosomal cistron phylogenetic tree yielded similar topologies and intra-genus divergence time. However, ribosomal phylogenetic tree showed lower posterior probability and larger time estimation range (Supplementary Fig. 2), and ribosomal DNA has limited credibility as one of the nuclear genes. Thus, we mainly referred to the plastid phylogenetic tree. Furthermore, the inclusion of other species within Rehderodendron genus exerts minimal impact on the estimation of populations divergence time, given that R. kwangtungense represents the most ancestral clade (Supplementary Fig. 3). Therefore, only the phylogenetic tree constructed from plastid genomes excluding other Rehderodendron species was presented (Fig. 6a and b), which showed that existing R. kwangtungense populations differentiated at ca. 26.25 Ma. In addition, populations P14, P15, P16 and P8 were placed on the root of the tree; these populations diverged from the remaining populations at ca. 22–23 Ma. The populations P1-P7 formed a monophyletic clade, and the crown age was estimated to be ca. 12.92 Ma. Additionally, populations P9-P13 formed another monophyletic clade that diverged only recently (ca. 1.73 Ma).

Fig. 6
figure 6

(a) The presumed migration path of R. kwangtungense. (b) Phylogenetic tree and divergent time estimation (below the bars) of R. kwangtungense based on the plastid genome. The fossil points are marked with triangles. Maximum likelihood standard bootstrap support values (BS) and Bayesian posterior probabilities (PP) are shown above the bars (*: BS ≥ 95% or PP ≥ 0.99). The median ages of the nodes are shown below the branches (those younger than 1 Ma are not displayed), with blue bars indicating the 95% highest posterior density intervals

Suitable distribution based on ENM

Among the 7 bioclimatic variables selected for ENM simulation, the bioclimatic variables with greater contribution percent are Bio17, Bio7, Bio6; the bioclimatic variables with greater permutation importance are Bio7, Bio2, Bio6; and the bioclimatic variables with greater regularized training gain of the jackknife test are Bio2, Bio17, Bio19 (Supplementary Tables 3 & Supplementary Fig. 4). In conclusion, the bioclimatic variables that dominate the potential geographical distribution of R. kwangtungense are mean diurnal range (bio2), min temperature of coldest month (Bio6), temperature annual range (Bio7) and precipitation of driest quarter (Bio17). The suitable areas for R. kwangtungense were few and scattered in southern China during the LGM (5.01 km2, Fig. 7a), increased obviously at ca. 6000 YA (26.10 km2, Fig. 7b), and reached a maximum at present (43.51 km2, Fig. 7c). The predicted areas at ca. 6000 YA and at present are larger than the actual distribution area of the species, mostly in southeastern China.

Fig. 7
figure 7

Presence records (dark circles) and predicted distribution probability according to climatic conditions during the (a) last glacial maximum (22 kya), (b) mid-Holocene (6 kya), and (c) current (0 kya)

Discussion

Origin and early dispersal of R. kwangtungense

Although there have been some phylogenetic studies of family Styracaceae, our study is the first to infer the internal divergence time of family Styracaceae, so there is no intra-family divergence time reference. Based on large-scale phylogenetic tree, Rose et al. [78] showed that the divergence time between the family Styracaceae and the family Diapensiaceae was about 93.2 Ma, and Li et al. [79] showed that the formation time of Ericales was between 103 and 107 Ma, which is close to our estimated time. Rehderodendron is suggested to be a member of “boreotropical flora” [48] because of the existence of a continuous zone of maritime influenced vegetation along the Tethys during the Cenozoic [80]. Our phylogenetic results showed that populations P14 and P15, which are located in the coastal mountains of Southeast China, harbor the most ancient plastid genomes (Fig. 6b). These findings suggest that the ancestor species of R. kwangtungense may have initially spread to the Cathaysian coast through the lowland region along the Tethys Ocean during the early Eocene, when the Tibetan Plateau had not yet fully uplifted [34]. This early diffusion event is supported by the similarity of R. kwangtungense pollen and pollen from the early Eocene [48].

Most of the eastern populations (P8, P14, P15, P16) and early clades (P1-P7, P9-P11) of R. kwangtungense quickly diverged between 26.25 Ma and 21.48 Ma (Fig. 6b). This result implies that the modern distribution area of R. kwangtungense was already established before the early Miocene.

Field observation revealed that the fruits of R. kwangtungense are large and cylindrical with thick spongy mesocarps (Fig. 1c-d), and long-distance transport by animals is difficult (only rodents eat them). Thus, water and gravity are the main forces driving lateral fruit distribution [48]. The Mantel test also confirmed the weak dispersal ability of R. kwangtungense (Fig. 4). Therefore, the fruits of R. kwangtungense do not have the ability to spread between mountains in different river basins. The current island-like distribution pattern most likely formed before mountain uplift and erosion. Palynological evidence also suggests that the genus Rehderodendron seems to be more adapted to braided river system transmission [80]. The existence of a panplain stage between the Oligocene and Miocene with many nondirectional water systems [39], seems to be the most plausible explanation for the early dispersal of R. kwangtungense.

Differentiation of R. kwangtungense since the establishment of the east Asian monsoon

The results of population structure analysis suggest that R. kwangtungense populations should be divided into seven groups (Fig. 2). Group 1–4 are located on the eastern side of the Yunnan-Guizhou Plateau on the second level of the terrain ladder in China, and the flowering period of these populations is late. Group 5–7 are located in the Nanling Mountains and the southern region on the third level of the terrain ladder, and the flowering period of these populations is early. The phylogenetic tree shows that the crown age of extant R. kwangtungense was estimated to be 26.25 Ma, and this species underwent a period of rapid populations divergence at 23 − 21 Ma (Fig. 6b), coinciding with the rapid rise of the Himalayas and subsequent establishment of the East Asian monsoon. Ding et al. proposed that the Himalayas experienced a rapid increase in elevation to 4000 m from 23 to 19 Ma [81], while the Yunnan-Guizhou Plateau was uplifted at a much greater rate than the Nanling Mountains, thus forming the three-step ladder of China’s topography and restructuring the monsoon circulation system. Sun & Wang also suggested that the topographic inversion that occurred at ca. 25 Ma, resulted in the formation of the East Asian monsoon, which had a large impact on the dry and wet patterns in China [31]. The geographical isolation caused by geological movements hindered gene flow, and monsoonal rains intensified mountain erosion, which together promoted the differentiation of R. kwangtungense into two groups from 26 − 21 Ma. Moreover, a differentiated population structure between the Yunnan-Guizhou Plateau and South China is common among relict species, such as Bretschneidera sinensis [25], Eomecon chionantha [82] and Eurycorymbus cavaleriei [24].

Subsequently, continuous geographical processes have led to the formation of more complex mountain systems, and several phases of intensification in the Asian monsoon have led to climate change [83]. The R. kwangtungense in the Yunnan-Guizhou Plateau has been continuously differentiated since 12.69 Ma. Similar divergence processes were also observed in Cyclocarya paliurus and Eriobotrya during this period [84, 85].

Effects of Quaternary glaciation on the R. kwangtungense and implications for conservation

Southern China, which has numerous mountains and valleys, provides long-term stable habitats for many species, allowing them to maintain in situ, and harbours a great abundance of Tertiary relict flora [18, 19]. Under the influence of repeated ice ages in the Quaternary, the population of R. kwangtungense contracted dramatically and took refuge in mountainous areas. The geographic isolation caused by population contraction restricted gene flow, which combined with a changing climate, facilitated the second phase of rapid divergence in P1-2, P4-6 and P9-13. Similarly, the lineage differentiation has also occurred in Cercdiphyllum japonicum, Kalopanax septemlobus and other species [86, 87].

Based on ENM and historical dynamic analysis, when the LGM occurred, the area suitable for R. kwangtungense decreased sharply, and the effective population size reached an extremely low level, indicating that R. kwangtungense was not adapted to the cold climate and that a large number of populations died out, with only a small minority taking refuge in the mountains. As the climate became warm and humid after the LGM, the suitable habitat area for R. kwangtungense continued to expand, and the historical dynamic simulation confirmed that the population size increased simultaneously, reaching a peak at approximately 10,000 YA (Figs. 5 and 7). However, since the modern landscape had been basically formed at this point, the dispersal of R. kwangtungense was limited, and the species only expanded in situ, forming an island-like distribution pattern. AMOVA indicated that genetic variation within populations was minimal (Table 2), as was that of HE and HO in various groups (Table 1), and the phylogenetic tree constructed from RAD data showed a very short branch length for all individuals (Fig. 3b). This evidence points to the limited genetic diversity of R. kwangtungense populations, and a population may be the result of the cloning of a few individuals. Thus, we surmised that the dispersal type of fruits and pollination mode determined whether population expansion occurred after the ice age. Similarly, Cathaya argyrophylla did not experience long-distance dispersal or population expansion after the glacial period [23]. In contrast, Cyclocarya paliurus with samara experienced substantial spatial expansion [85].

As temperatures and precipitation become more favorable, the effective population size has shrunk over the past 1,000 years. The main reason for the inconsistency may be the weak diffusion ability and competition ability of R. kwangtungense. According to the field observation, the genus Rehderodendron can quickly grow many small trees after the forest gap period, while in the original forest with good vegetation, the seedlings are basically invisible (because the fruit will soon become mildew if it is not dried in time). In other words, we assume that the suitable habitat for R. kwangtungense is the plain, with rivers and flat land, and the current mountain habitat is not suitable for the diffusion and competition of R. kwangtungense. The ENM result is mainly based on temperature and precipitation, but in reality, temperature and precipitation are not the only limiting factors.

Another reason for the reduction in effective population size could be human activities. To our knowledge, many populations are currently located in natural reserves where the species may be well protected (P8, P10 and P13), although in some regions, such as P5, the species are susceptible to the impacts of nearby villages. The fragile genetic diversity of the populations suggests that all populations should be protected, and priority should be given to representative populations in each group. Field observations have shown that seedlings are rare, possibly due to failure to compete with other species, while manual cultivation has shown that successful germination of seeds can be as high as 80% [88]. Therefore, artificial breeding and other measures can be adopted for protection.

Conclusion

The understanding of the origin of the contemporary distribution pattern of relict plants is crucial for the scientific conservation of rare and endangered plant species. In this study, we revealed the population differentiation history of the relict species Rehderodendron kwangtungense Chun during the Cenozoic. The modern distribution area of R. kwangtungense was already established during the Oligocene-Miocene transition period, and subsequently underwent fragmentation in Miocene as a result of allometric plate uplift and mountain erosion. During the Quaternary glacial period, R. kwangtungense experienced a significant population decline and sought refuge in situ. The findings of our study emphasize the predominant role of geographical processes in shaping relict plant distribution patterns, while the impact of climate fluctuations primarily manifests through their influence on local population size. This means that scenario of population expansion and contraction under climate oscillations may be overestimated, especially for those species with low dispersal capacity. Furthermore, we propose a perspective on the existence of a panplain in south China during the Oligocene to Miocene epochs, which has significant implications for understanding the formation process of relict flora in East Asia.

Data availability

The RAD-seq data has been uploaded to NCBI (https://www.ncbi.nlm.nih.gov/) Sequence Read Archive under the BioProject PRJNA1056300. The assembled chloroplast genomes data has been uploaded to NCBI GenBank under the accessions PP048903-PP048918, PP761266-PP761270. And the ribosomal DNA data has been uploaded to NCBI GenBank under the accessions PP765768-PP765788.

Abbreviations

RAD-seq:

Restriction-site associated DNA sequencing

ENM:

Ecological niche modelling

LGM:

Last glacial maximum

H O :

Observed heterozygosity

H E :

Expected heterozygosity

π:

Nucleotide diversity

F IS :

Inbreeding coefficient

F ST :

Fixation index

K:

Hypothetical ancestral population

CV:

Cross-validation

PCA:

Principal component analysis

AMOVA:

Analysis of molecular variance

SFS:

Site frequency spectra

IBD:

Isolation-by-distance

IBE:

Isolation-by-environment

YA:

Years ago

References

  1. Habel JC, Assmann T, Schmitt T, Avise JC. Relict species: From Past to Future. In: Habel JC, Assmann T, editors. Relict species. Berlin, Heidelberg: Springer; 2010. pp. 1–5. https://doi.org/10.1007/978-3-540-92160-8_1

    Chapter  Google Scholar 

  2. Grandcolas P, Nattier R, Trewick S. Relict species: a relict concept? Trends Ecol Evol. 2014;29(12):655–63. https://doi.org/10.1016/j.tree.2014.10.002

    Article  PubMed  Google Scholar 

  3. Friedman WE, Ryerson KC. Reconstructing the ancestral female gametophyte of angiosperms: insights from Amborella and other ancient lineages of flowering plants. Am J Bot. 2009;96(1):129–43. https://doi.org/10.3732/ajb.0800311

    Article  PubMed  Google Scholar 

  4. Soltis DE, Soltis PS, Albert VA, Oppenheimer DG, DePamphilis CW, Ma H, Frohlich MW, Theißen G. Missing links: the genetic architecture of flower and floral diversification. Trends Plant Sci. 2002;7(1):22–31. https://doi.org/10.1016/S1360-1385(01)02098-2

    Article  CAS  PubMed  Google Scholar 

  5. Jones ME, Tennyson AJ, Worthy JP, Evans SE, Worthy TH. A sphenodontine (Rhynchocephalia) from the Miocene of New Zealand and palaeobiogeography of the tuatara (Sphenodon). Proc Biol Sci. 2009;276(1660):1385–90. https://doi.org/10.1098/rspb.2008.1785

    Article  PubMed  PubMed Central  Google Scholar 

  6. Ladiges PY, Cantrill D. New Caledonia–Australian connections: biogeographic patterns and geology. Aust Syst Bot. 2007;20(5):383–9. https://doi.org/10.1071/SB07018

    Article  Google Scholar 

  7. Bartolini F, Giomi F. Microclimate drives intraspecific thermal specialization: conservation perspectives in freshwater habitats. Conserv Physiol. 2021;9(1):coab006. https://doi.org/10.1093/conphys/coab006

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Cotado A, Munné-Bosch S, Pintó‐Marijuan M. Strategies for severe drought survival and recovery in a Pyrenean relict species. Physiol Plant. 2020;169(2):276–90. https://doi.org/10.1111/ppl.13072

    Article  CAS  PubMed  Google Scholar 

  9. Feng Y, Comes HP, Chen J, Zhu S, Lu R, Zhang X, Li P, Qiu J, Olsen KM. Qiu Y.Genome sequences and population genomics provide insights into the demographic history, inbreeding, and mutation load of two ‘living fossil’ tree species of Dipteronia. Plant J. 2023. https://doi.org/10.1111/tpj.16486

    Article  PubMed  Google Scholar 

  10. Tan J, Kelly CK, Jiang L. Temporal niche promotes biodiversity during adaptive radiation. Nat Commun. 2013;4(1):2102. https://doi.org/10.1038/ncomms3102

    Article  CAS  PubMed  Google Scholar 

  11. Yoder J, Clancey E, Des Roches S, Eastman J, Gentry L, Godsoe W, Hagey T, Jochimsen D, Oswald B, Robertson J. Ecological opportunity and the origin of adaptive radiations. J Evol Biol. 2010;23(8):1581–96. https://doi.org/10.1111/j.1420-9101.2010.02029.x

    Article  CAS  PubMed  Google Scholar 

  12. Milne RI, Abbott RJ. The origin and evolution of Tertiary relict floras. Adv Bot Res. 2002;38(4):281–314. https://doi.org/10.1016/S0065-2296(02)38033-9

    Article  Google Scholar 

  13. Qiu Y, Lu Q, Zhang Y, Cao YJBS. Phylogeography of East Asia’s Tertiary relict plants: current progress and future prospects. Biodivers Sci. 2017;25(2):136. https://doi.org/10.17520/biods.2016292

    Article  Google Scholar 

  14. Tiffney BH. The Eocene North Atlantic land bridge: its importance in Tertiary and modern phytogeography of the Northern Hemisphere. J Arnold Arboretum. 1985;66(2):243–73. https://www.jstor.org/stable/43782164

    Article  Google Scholar 

  15. Lu L-M, Mao L-F, Yang T, Ye J-F, Liu B, Li H-L, Sun M, Miller JT, Mathews S, Hu H-H. Evolutionary history of the angiosperm flora of China. Nat Plants. 2018;554(7691):234–8. https://doi.org/10.1038/nature25485

    Article  CAS  Google Scholar 

  16. Qian H, Deng T, Jin Y, Mao L, Zhao D, Ricklefs RE. Phylogenetic dispersion and diversity in regional assemblages of seed plants in China. Proceedings of the National Academy of Sciences. 2019;116(46):23192–23201. https://doi.org/10.1073/pnas.1822153116

  17. Qian H, Ricklefs RE. A comparison of the taxonomic richness of vascular plants in China and the United States. Am Nat. 1999;154(2):160–81. https://doi.org/10.1086/303230

    Article  PubMed  Google Scholar 

  18. López-Pujol J, Zhang FM, Sun HQ, Ying TS, Ge S. Centres of plant endemism in China: places for survival or for speciation? J Biogeogr. 2011;38(7):1267–80. https://doi.org/10.1111/j.1365-2699.2011.02504.x

    Article  Google Scholar 

  19. Qian H, Ricklefs RE. Large-scale processes and the Asian bias in species diversity of temperate plants. Nat Plants. 2000;407(6801):180–2. https://doi.org/10.1038/35025052

    Article  CAS  Google Scholar 

  20. Tang CQ, Matsui T, Ohashi H, Dong Y-F, Momohara A, Herrando-Moraira S, Qian S, Yang Y, Ohsawa M, Luu HT. Identifying long-term stable refugia for relict plant species in East Asia. Nat Commun. 2018;9(1):4488. https://doi.org/10.1038/s41467-018-06837-3

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Woolbright SA, Whitham TG, Gehring CA, Allan GJ, Bailey JK. Climate relicts and their associated communities as natural ecology and evolution laboratories. Trends Ecol Evol. 2014;29(7):406–16. https://doi.org/10.1016/j.tree.2014.05.003

    Article  PubMed  Google Scholar 

  22. Yin QY, Fan Q, Li P, Truong D, Zhao WY, Zhou RC, Chen SF, Liao WB. Neogene and Quaternary climate changes shaped the lineage differentiation and demographic history of Fokienia Hodginsii (Cupressaceae Sl), a Tertiary relict in East Asia. J Syst Evol. 2021;59(5):1081–99. https://doi.org/10.1111/jse.12582

    Article  Google Scholar 

  23. WANG HW, Ge S. Phylogeography of the endangered cathaya argyrophylla (Pinaceae) inferred from sequence variation of mitochondrial and nuclear DNA. Mol Ecol. 2006;15(13):4109–22. https://doi.org/10.1111/j.1365-294X.2006.03086.x

    Article  CAS  PubMed  Google Scholar 

  24. Wang J, Gao P, Kang M, Lowe AJ, Huang H. Refugia within refugia: the case study of a canopy tree (Eurycorymbus Cavaleriei) in subtropical China. J Biogeogr. 2009;36(11):2156–64. https://doi.org/10.1111/j.1365-2699.2009.02165.x

    Article  Google Scholar 

  25. Wang M-N, Duan L, Qiao Q, Wang Z-F, Zimmer EA, Li Z-C, Chen H-F. Phylogeography and conservation genetics of the rare and relict Bretschneidera sinensis (Akaniaceae). PLoS ONE. 2018;13(1):e0189034. https://doi.org/10.1371/journal.pone.0189034

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Zhang J, Li Z, Fritsch PW, Tian H, Yang A, Yao X. Phylogeography and genetic structure of a Tertiary relict tree species, Tapiscia Sinensis (Tapisciaceae): implications for conservation. Ann Botany. 2015;116(5):727–37. https://doi.org/10.1093/aob/mcv112

    Article  CAS  Google Scholar 

  27. Cao YN, Zhu SS, Chen J, Comes HP, Wang IJ, Chen LY, Sakaguchi S, Qiu YX. Genomic insights into historical population dynamics, local adaptation, and climate change vulnerability of the east Asian tertiary relict Euptelea (Eupteleaceae). Evol Appl. 2020;13(8):2038–55. https://doi.org/10.1111/eva.12960

    Article  PubMed  PubMed Central  Google Scholar 

  28. Wan Q, Huang K, Chen S-F, Boyer F, Taberlet P, Li H, Chen C, Tang Y, Zheng Z, Zhang X. Fagus diversification in China in relation to east Asian monsoon evolution. Q Sci Rev. 2023;320:108350. https://doi.org/10.1016/j.quascirev.2023.108350

    Article  Google Scholar 

  29. Xu W, Dong W-J, Fu T-T, Gao W, Lu C-Q, Yan F, Wu Y-H, Jiang K, Jin J-Q, Chen H-M, et al. Herpetological phylogeographic analyses support a Miocene focal point of himalayan uplift and biological diversification. Natl Sci Rev. 2020;8(9). https://doi.org/10.1093/nsr/nwaa263

  30. Chen P. Coastal mountains of SE China, desertization and saliniferous lakes of Central China during the Upper cretaceous. J Stratigraphy. 1997;21(3):203–13.

    Google Scholar 

  31. Sun X, Wang P. How old is the Asian monsoon system?—Palaeobotanical records from China. Palaeogeography. Palaeoclimatology Palaeoecology. 2005;222(3–4):181–222. https://doi.org/10.1016/j.palaeo.2005.03.005

    Article  Google Scholar 

  32. Wang P. Cenozoic deformation and the history of sea-land interactions in Asia. Geophys Monogr Ser. 2004;149:1–22. https://doi.org/10.1029/149GM01

    Article  CAS  Google Scholar 

  33. Min ZHANG, Hai-long LI, Ling TANG, Shi-hua QI, Jian KUANG, Wen WANG. Jin-long Z: deep structure, sedimentary sequence and tectonic evolution of Heyuan Basin, Eastern Guangdong Province. Acta Geoscientica Sinica. 2023;44(4):1–19.

    Google Scholar 

  34. Ding L, Kapp P, Cai F, Garzione CN, Xiong Z, Wang H, Wang C. Timing and mechanisms of Tibetan Plateau uplift. Nat Reviews Earth Environ. 2022;3(10):652–67. https://doi.org/10.1038/s43017-022-00318-4

    Article  Google Scholar 

  35. Ding W-N, Ree RH, Spicer RA, Xing Y-W. Ancient orogenic and monsoon-driven assembly of the world’s richest temperate alpine flora. Science. 2020;369(6503):578–81. https://doi.org/10.1126/science.abb4484

    Article  CAS  PubMed  Google Scholar 

  36. Hoke GD, Liu-Zeng J, Hren MT, Wissink GK, Garzione CN. Stable isotopes reveal high southeast Tibetan Plateau margin since the Paleogene. Earth Planet Sci Lett. 2014;394:270–8. https://doi.org/10.1016/j.epsl.2014.03.007

    Article  CAS  Google Scholar 

  37. Li S-F, Valdes PJ, Farnsworth A, Davies-Barnard T, Su T, Lunt DJ, Spicer RA, Liu J, Deng W-Y-D, Huang J. Orographic evolution of northern Tibet shaped vegetation and plant diversity in eastern Asia. Sci Adv. 2021;7(5):eabc7741. https://doi.org/10.1126/sciadv.abc7741

    Article  PubMed  PubMed Central  Google Scholar 

  38. Su T, Spicer RA, Li S-H, Xu H, Huang J, Sherlock S, Huang Y-J, Li S-F, Wang L, Jia L-B. Uplift, climate and biotic changes at the eocene–oligocene transition in south-eastern Tibet. Natl Sci Rev. 2019;6(3):495–504. https://doi.org/10.1093/nsr/nwy062

    Article  CAS  PubMed  Google Scholar 

  39. He J, Garzanti E, Cao L, Wang H. The zircon story of the Pearl River (China) from cretaceous to present. Earth Sci Rev. 2020;201:103078. https://doi.org/10.1016/j.earscirev.2019.103078

    Article  CAS  Google Scholar 

  40. Yu Y, Zhang C, Li S, Zhu R, Liu J, Qin C, Zhang Z. Influences of Tibetan Plateau uplift on provenance evolution of the paleo-Pearl River. Chin J Geochem. 2015;34:208–18. https://doi.org/10.1007/s11631-015-0032-z

    Article  CAS  Google Scholar 

  41. Cao L, Shao L, Qiao P, Zhao Z, van Hinsbergen DJ. Early miocene birth of modern Pearl River recorded low-relief, high-elevation surface formation of SE Tibetan Plateau. Earth Planet Sci Lett. 2018;496:120–31. https://doi.org/10.1016/j.epsl.2018.05.039

    Article  CAS  Google Scholar 

  42. Jin H, Wan S, Clift PD, Liu C, Huang J, Jiang S, Li M, Qin L, Shi X, Li A. Birth of the Pearl River at 30 ma: evidence from sedimentary records in the northern South China Sea. Earth Planet Sci Lett. 2022;600:117872. https://doi.org/10.1016/j.epsl.2022.117872

    Article  CAS  Google Scholar 

  43. Zhang L, Fang X. Palaeogeography of China: the formation of China’s natural environment. Beijing: In.: Science; 2012.

    Google Scholar 

  44. Zhang K, Huang Y. Researches on the planation surfaces in north Guangdong. Trop Geogr. 1995;15(4):295–305.

    Google Scholar 

  45. Zhou S, Liu J, Xi Z, Zhou L. The residual top planation surfaces in the Nanling mountains and their adjacent regions. J Glaciology Geocryology. 2008;30(6):938–45.

    Google Scholar 

  46. Liu F, Xue Y, CJ H, Xi Z, Xiong L, Yang Y. Research on the geomorphological characteristics of the planation surfaces at the junction of Hunan and Guangxi based on DEM. J Hunan City Unive Rsity (Natural Science). 2023;32(4):19–26. https://doi.org/10.3969/j.issn.1672-7304.2023.04.0004

    Article  Google Scholar 

  47. Geissert F, Gregor H. Einige interessante und neue sommergrüne Pflanzenelemente (Fruktifikationen) aus dem Elsäßer Pliozän (Genera Sabia Colebr., Wikstroemia Endl., Alangium Lam., Nyssa L., Halesia Ellis, Rehderodendron Hu). Mitt Badischen Landesver Naturk Naturs eV Freiburg Breisgau 1981, 12:233–239.

  48. Hofmann C-C, Zhao W-Y. Unravelling the palaeo-biogeographical history of the living fossil genus Rehderodendron (Styracaceae)–an Asian relic element of the eurasian cenozoic flora. 2022. https://doi.org/10.1186/s12862-022-02097-4

  49. Mai D. Subtropische elemente im europäischen Tertiär I. Paläont Abh, B 1970, 3(3–4):441–503.

  50. Mai DH, Petrescu I. Eine Neue Rehderodendron-Art (Styracaceae) aus dem oberen Pliozän Des Baraolt-Beckens (SR Rumänien). Z für Geologische Wissenschaften. 1983;11(7):915–25.

    Google Scholar 

  51. Martinetto E. East Asian elements in the Plio-Pleistocene floras of Italy. In: Proceedings of the International Symposium on Floristic Character Diversity of East Asian Plants: 1998. Springer-Verlag Berlin, Germany: 71–87.

  52. Zhao W-Y, Fritsch PW, Liu Z-C, Fan Q, Jin J-H, Liao W-B. New combinations and synonyms in Rehderodendron. (Styracaceae) PhytoKeys. 2020;161:79. https://doi.org/10.3897/phytokeys.161.54970

    Article  PubMed  PubMed Central  Google Scholar 

  53. Rochette NC, Catchen JM. Deriving genotypes from RAD-seq short-read data using Stacks. Nat Protoc. 2017;12(12):2640–59. https://doi.org/10.1038/nprot.2017.123

    Article  CAS  PubMed  Google Scholar 

  54. Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bull. 1987;19(1):11–5.

    Google Scholar 

  55. Catchen J, Hohenlohe PA, Bassham S, Amores A, Cresko WA. Stacks: an analysis tool set for population genomics. Mol Ecol. 2013;22(11):3124–40. https://doi.org/10.1111/mec.12354

    Article  PubMed  PubMed Central  Google Scholar 

  56. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8. https://doi.org/10.1093/bioinformatics/btr330

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123(3):585–95. https://doi.org/10.1093/genetics/123.3.585

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19(9):1655–64. https://doi.org/10.1101/gr.094052.109

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4(1):s13742. https://doi.org/10.1186/s13742-015-0047-8. -13015-10047-13748.

    Article  CAS  Google Scholar 

  60. Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10(3):564–7. https://doi.org/10.1111/j.1755-0998.2010.02847.x

    Article  PubMed  Google Scholar 

  61. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, Von Haeseler A, Lanfear R. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biology Evol. 2020;37(5):1530–4. https://doi.org/10.1093/molbev/msaa015

    Article  CAS  Google Scholar 

  62. Letunic I, Bork P. Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6. https://doi.org/10.1093/nar/gkab301

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Liu X, Fu Y-X. Stairway plot 2: demographic history inference with folded SNP frequency spectra. Genome Biol. 2020;21(1):1–9. https://doi.org/10.1186/s13059-020-02196-9

    Article  Google Scholar 

  64. Liu X, Fu Y-X. Exploring population size changes using SNP frequency spectra. Nat Genet. 2015;47(5):555–9. https://doi.org/10.1038/ng.3254

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Luo M-C, You FM, Li P, Wang J-R, Zhu T, Dandekar AM, Leslie CA, Aradhya M, McGuire PE, Dvorak J. Synteny analysis in Rosids with a walnut physical map reveals slow genome evolution in long-lived woody perennials. BMC Genomics. 2015;16:1–17. https://doi.org/10.1186/s12864-015-1906-5

    Article  CAS  Google Scholar 

  66. Dierckxsens N, Mardulyn P, Smits G. NOVOPlasty: de novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2017;45(4):e18–18. https://doi.org/10.1093/nar/gkw955

    Article  CAS  PubMed  Google Scholar 

  67. Jin J-J, Yu W-B, Yang J-B, Song Y, DePamphilis CW, Yi T-S, Li D-Z. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020;21:1–31. https://doi.org/10.1186/s13059-020-02154-5

    Article  Google Scholar 

  68. Yan M-H, Li C-Y, Fritsch PW, Cai J, Wang H-C. Phylogeny of the Styracaceae Revisited based on whole plastome sequences, including Novel Plastome Data from Parastyrax. Syst Bot. 2021;46(1):162–74. https://doi.org/10.1600/036364421X16128061189576

    Article  Google Scholar 

  69. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biology Evol. 2013;30(4):772–80. https://doi.org/10.1093/molbev/mst010

    Article  CAS  Google Scholar 

  70. Kumar S, Stecher G, Tamura K. MEGA7. Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biology Evol. 2016;33(7):1870–4. https://doi.org/10.1093/molbev/msw054

    Article  CAS  Google Scholar 

  71. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biology Evolutionary Appl. 2012;29(8):1969–73. https://doi.org/10.1093/molbev/mss075

    Article  CAS  Google Scholar 

  72. Mai DH. Contribution to the flora of the middle Oligocene Calau beds in Brandenburg, Germany. Rev Palaeobotany Palynology. 1998;101(1–4):43–70. https://doi.org/10.1016/S0034-6667(97)00069-9

    Article  Google Scholar 

  73. MacGinitie HD. Fossil plants of the Florissant beds, Colorado. (No Title). 1953.

  74. Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int J Climatol. 2017;37(12):4302–15. https://doi.org/10.1002/joc.5086

    Article  Google Scholar 

  75. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatology: J Royal Meteorological Soc. 2005;25(15):1965–78. https://doi.org/10.1002/joc.1276

    Article  Google Scholar 

  76. Phillips SJ, Dudík M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008;31(2):161–75. https://doi.org/10.1111/j.0906-7590.2008.5203.x

    Article  Google Scholar 

  77. Liu H, Wang Z, Zhang Y, Li M, Wang T, Su Y. Geographic isolation and environmental heterogeneity contribute to genetic differentiation in Cephalotaxus oliveri. Ecol Evol. 2023;13(3):e9869. https://doi.org/10.1002/ece3.9869

    Article  PubMed  PubMed Central  Google Scholar 

  78. Rose JP, Kleist TJ, Löfstrand SD, Drew BT, Schönenberger J, Sytsma KJ, Evolution. Phylogeny, historical biogeography, and diversification of angiosperm order Ericales suggest ancient neotropical and east Asian connections. Mol Phylogenetics. 2018;122:59–79. https://doi.org/10.1016/j.ympev.2018.01.014

    Article  Google Scholar 

  79. Li H-T, Yi T-S, Gao L-M, Ma P-F, Zhang T, Yang J-B, Gitzendanner MA, Fritsch PW, Cai J, Luo Y. Origin of angiosperms and the puzzle of the jurassic gap. Nat Plants. 2019;5(5):461–70. https://doi.org/10.5061/dryad.bq091cg

    Article  PubMed  Google Scholar 

  80. Hofmann C-C, Sachse M. SEM pollen analysis of Miocene deposits of Entrischenbrunn (Bavaria, Germany) reveal considerable amounts of pollen of subhumid and sclerophyllous together with azonal water plants reflecting the vegetation mosaic of a braided river system. Rev Palaeobotany Palynology. 2023;308:104787. https://doi.org/10.1016/j.revpalbo.2022.104787

    Article  Google Scholar 

  81. Ding L, Spicer R, Yang J, Xu Q, Cai F, Li S, Lai Q, Wang H, Spicer T, Yue Y. Quantifying the rise of the Himalaya orogen and implications for the south Asian monsoon. Geology. 2017;45(3):215–8. https://doi.org/10.1130/G38583.1

    Article  Google Scholar 

  82. Tian S, Kou Y, Zhang Z, Yuan L, Li D, López-Pujol J, Fan D, Zhang Z. Phylogeography of Eomecon Chionantha in subtropical China: the dual roles of the Nanling Mountains as a glacial refugium and a dispersal corridor. BMC Evol Biol. 2018;18(1):1–12. https://doi.org/10.1186/s12862-017-1093-x

    Article  Google Scholar 

  83. Wan S, Li A, Clift PD, Stuut J-BW. Development of the east Asian monsoon: mineralogical and sedimentologic records in the northern South China Sea since 20 ma. Palaeogeography. Palaeoclimatology Palaeoecology. 2007;254(3–4):561–82. https://doi.org/10.1016/j.palaeo.2007.07.009

    Article  Google Scholar 

  84. Chen S, Milne R, Zhou R, Meng K, Yin Q, Guo W, Ma Y, Mao K, Xu K, Kim YD. When tropical and subtropical congeners met: multiple ancient hybridization events within Eriobotrya in the Yunnan-Guizhou Plateau, a tropical‐subtropical transition area in China. Mol Ecol. 2022;31(5):1543–61. https://doi.org/10.1111/mec.16325

    Article  CAS  PubMed  Google Scholar 

  85. Kou Y, Cheng S, Tian S, Li B, Fan D, Chen Y, Soltis DE, Soltis PS, Zhang Z. The antiquity of Cyclocarya paliurus (Juglandaceae) provides new insights into the evolution of relict plants in subtropical China since the late early miocene. J Biogeogr. 2016;43(2):351–60. https://doi.org/10.1111/jbi.12635

    Article  Google Scholar 

  86. Qi XS, Chen C, Comes HP, Sakaguchi S, Liu YH, Tanaka N, Sakio H, Qiu YX. Molecular data and ecological niche modelling reveal a highly dynamic evolutionary history of the East Asian Tertiary relict Cercidiphyllum (Cercidiphyllaceae). New Phytologist 2012, 196(2):617–630. https://doi.org/10.1111/j.1469-8137.2012.04242.x

  87. Sakaguchi S, QIU YX, LIU YH, QI XS, KIM SH, Han J, Takeuchi Y, Worth JR, Yamasaki M, Sakurai S. Climate oscillation during the quaternary associated with landscape heterogeneity promoted allopatric lineage divergence of a temperate tree Kalopanax septemlobus (Araliaceae) in East Asia. Mol Ecol. 2012;21(15):3823–38. https://doi.org/10.1111/j.1365-294X.2012.05652.x

    Article  PubMed  Google Scholar 

  88. Wan C. Germination methods of three kinds of rare and endangered plant seeds. Seed. 1991;5:65–7.

    Google Scholar 

Download references

Acknowledgements

We thank Fan Ye, Su-Juan Wei for their help in the field work.

Funding

This work was supported by National Natural Science Foundation of China (32300172), Natural Science Foundation of Guangdong Province (2021A1515110425), Guangdong Special Fund for Natural Resources Management and Ecological Forestry Construction (2021GJGY001), Guangdong Province Ecological Forestry Development Project, and the project of the Research on Coevolution of Vegetation and Geological Environment of Shenzhen Dapeng Peninsula Geopark (2020F36).

Author information

Authors and Affiliations

Authors

Contributions

Jiehao Jin: study design, fieldwork, data analysis, writing; Wanyi Zhao: study design, fieldwork, writing; Sufang Chen: data analysis, writing; Chao Gu, Zhihui Chen, Zhongcheng Liu: fieldwork; Wenbo Liao, Qiang Fan: study design.All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wanyi Zhao or Qiang Fan.

Ethics declarations

Ethics approval and consent to participate

This study did not involve any human tissue materials or animal tissue materials. It did not require ethical approval. We declare that all plant materials used in this study were collected and sequenced in accordance with relevant institutional, national, and international guidelines and legislation. Relevant permits, permissions, and licenses were obtained prior to the collection and use of plant materials. All steps in the study were performed in accordance with the relevant guidelines and regulations. W.Z. undertook the formal identification of the plant material used in this study. Voucher specimens were deposited at the Herbarium of Sun Yat-sen University (SYS) under the voucher numbers ZWY-1728 (population NKS), ZWY-1762 (population DYS), ZWY-1831 (population SSP), ZWY-1850 (population SD), ZWY-1868 (population SHS) and ZWY-1900 (population TJS).

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

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

Jin, J., Zhao, W., Chen, S. et al. Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense. BMC Plant Biol 24, 459 (2024). https://doi.org/10.1186/s12870-024-05181-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05181-7

Keywords