Analysis of the population genetic structure and demographic history of Tilia amurensis and Tilia japonica in China using SSR markers

Tilia amurensis, a deciduous broad-leaved tree, is distributed in regions neighboring those of Tilia japonica populations in China; these plants belong to the Malvaceae family. To date, the genetic status and demographic history of T. japonica remain unknown, and the demographic history of T. amurensis is not well defined. This study uses standard population genetic and approximate Bayesian computation (ABC) analyses of SSR data to determine the genetic status and divergence times of these two lime taxa. The results revealed that the genetic diversity of T. japonica (I = 1.181, H_O = 0.484, H_E = 0.602) was high at the population level. The results of the genetic structure revealed that the genetic variation was dominated primarily by within-population variation. Additionally, there was significant genetic differentiation and bidirectional introgression between T. amurensis and T. japonica. The ABC analysis suggests a (Middle) Pleistocene divergence. These findings have important implications for the formulation of appropriate conservation strategies. Specifically, the clarified divergence time and evidence of genetic exchange indicate that the two species should be protected based on their respective population conditions.

1 Introduction

Deforestation and forest fragmentation are the primary threats to forest plants (Lindenmayer and Fischer, 2013; Noh et al., 2022). Owing to the threat posed by various factors (such as human interference, habitat destruction and reproductive challenges), the population of Tilia has been declining annually, its distribution range has been shrinking and becoming fragmented. Among them, Tilia amurensis has been listed as a Class II Protected Wild Plant in China (Cui and Mu, 2016; Chen et al., 2023; Zhang et al., 2024). The evolutionary potential of a species in nature depends on its level of genetic variation and differentiation. Studies of threatened species can improve the understanding of genetic distribution within and among populations to establish effective strategies for species conservation and management (Hedrick and Miller, 1992).

The earliest Tilia fossils were found at the end of the Late Cretaceous and are based on morphological genera established from incomplete leaf fossils, whereas more reliable linden fossils occur in post-Palaeocene strata (Engler, 1909). Tilia fossils have been recorded for all periods from the Tertiary to the Pleistocene and are particularly well represented in East Asia (Engler, 1909; Berry, 1918). This suggests that the ancestral genealogy of T. amurensis and Tilia japonica in China may date back tens of millions of years. With the expectation of continued climate change and further range shifts, habitat fragmentation caused by human disturbance and overexploitation is threatening the survival of T. amurensis and T. japonica (Logan et al., 2018; Zhang et al., 2024; Mu and Liu, 2007). T. amurensis, a deciduous broad-leaved tree, grows mainly in northern China and North Korea, west of the upper reaches of the Yenisei River and east of the Korean Peninsula, where its exhibits an overlapping distribution with the closely related T. japonica in the Shandong region of China (Ya and Ren, 1996). There are subtle differences in morphology between T. amurensis and T. japonica in terms of leaves, bracts and inflorescences, but individuals in the ecotone are difficult to distinguish morphologically (Chen, 1990; Committee, 2007). Consequently, elucidating the spatiotemporal changes across their distribution ranges and the genetic profiles of the closely related species T. amurensis and T. japonica populations is of significant importance. Previously, we performed a preliminary study of the genetic status of T. amurensis and detected possible introgression between T. amurensis and T. japonica (Wu et al., 2024). However, the genetic background and demographic history of T. amurensis and its closely related species have not been thoroughly investigated. Therefore, in this study, we used simple sequence repeat (SSR) markers based on previous studies to determine the genetic structure of natural populations of T. amurensis and T. japonica in different distribution regions in China and inferred their divergence time from a common ancestor, thereby elucidating the evolutionary history of T. amurensis and T. japonica and providing a theoretical basis for the effective conservation and rational utilization of Tilia genetic resources.

2 Materialw and methods

2.1 Sample collection and DNA isolation

The targeted sites were in Heilongjiang, Jilin, Liaoning, Beijing, Hebei, Shandong, Anhui and Zhejiang Provinces in China (four sites in Shandong and one site in each of the remaining provinces). A total of 242 individuals were collected with eight T. amurensis populations (TS, LU, LS, ZS, YDS, DGL, MPS, and DXG) and three T. japonica populations (HS, MS and TMS) (Table 1; Figure 1). The sampling work was carried out from May to July 2023. The selection of target plants strictly followed the principle of random sampling, with a minimum distance of 30 meters between any two sampled individuals (the target sample size for each population was 30 plants; if the actual number of individuals was fewer than 30, all were sampled). Fresh young leaves of linden trees were collected, identified by Prof. Dekui Zang, and preserved in silica gel for subsequent DNA extraction. Voucher specimens were also preserved in the herbarium of Shandong Agricultural University.

Table 1

Table 1. Information for the Tilia sampling sites.

Figure 1

Figure 1. Geographic locations of the populations of Tilia included in this study. Descriptions of abbreviations for the study sites can be found in Table 1.

We used the SteadyPure Plant Genomic DNA Extraction Kit (CTAB method) for DNA extraction, and the experimental procedures were performed according to the methods of Wu’s study (Wu et al., 2024). To evaluate the quality and concentration of the extracted DNA, we used the NanoDrop™ One/OneC Microvolume UV-Vis Spectrophotometer (Thermo Fisher Scientific, USA) for measurement. In accordance with the experimental requirements, the DNA solution was diluted to 30 ng·μl^-1 for subsequent PCR amplification experiments.

2.2 PCR amplification of SSRs

We obtained SSR primers reported for the study of Tilia species and selected 27 microsatellite primers with relatively simple repeat motifs and high polymorphism rates for screening. The primers were synthesized by Shanghai Personalbio Technology Co., Ltd. Fluorescent dyes were added to the 5’ end of each primer. The addition of a fluorescent label enhances the functionality, stability, and detection efficiency of the primers, facilitating the rapid identification of amplification products in subsequent high-efficiency capillary electrophoresis. A total of 15 SSR primers from Yue’s study were used for subsequent experiments (Supplementary Table S1) (Yue et al., 2022). The reaction system for amplification and the conditions and parameters for PCR were as described in our previous study (Wu et al., 2024). The PCR amplification products were subsequently detected using 1.5% agarose gel electrophoresis. The PCR products that met the experimental requirements were purified and sent to Shanghai Personalbio Technology Co., Ltd., for high-efficiency capillary electrophoresis analysis.

2.3 Genetic diversity and population structure

The genetic diversity was estimated on the basis of the SSR data and analyzed using GenAlExv6.51 software to obtain the observed allele number (Na), efficient allele number (Ne), Shannon’s information index (I), observed heterozygosity (H_O), expected heterozygosity (H_E) and inbreeding coefficient (F_IS) (Peakall and Smouse, 2006). The obtained data were subsequently subjected to AMOVA and PCoA via GenAlExv6.51 (Smouse and Peakall, 2012). The genetic differentiation among populations was calculated via the genetic differentiation coefficient (F_ST), and the gene flow between populations was determined using the F_ST value: N_m = (1/F_ST − 1)/4 (Weir and Cockerham, 1984). UPGMA clustering analysis was performed using PopGen32, and a phylogenetic tree was constructed via MEGA7 to determine the genetic relationships among populations (Kumar et al., 2016). We used the Bayesian clustering method in STRUCTUREv.2.3.4 to investigate population structure (Pritchard et al., 2000). The optimal genetic group number K was determined by setting the “length of burn-in period” to 1,000, the “number of MCMC reps after burn-in” to 20,000, and calculating the ΔK value (Evanno et al., 2005).

2.4 Demographic history and effective population size

The approximate Bayesian computation (ABC) method, implemented in DIYABC v 2.1.0, was used to infer the demographic history of T. amurensis and T. japonica (Cornuet et al., 2014). On the basis of the findings from three approaches (STRUCTURE analysis, phylogenetic tree and PCoA) (Figures 2, 3) and the geographic distribution of the Tilia populations (Figure 1), four population groups were defined (Figure 4). Pop1 was composed of three Northeast China sites (DGL, MPS and DXG), Pop2 was composed of two North China sites (YDS and ZS), Pop3 was composed of three Shandong sites (TS, LU and LS), and Pop4 included the remaining three T. japonica sites (HS, MS and TMS). In scenario 1, Pop2 and Pop3 merged at time t1, Pop1 was present at time t2, and Pop4 was present at t3. The ancestors (NA) reached sufficient population size at t4. In scenario 2, Pop3 was assumed to have originated from an admixture of Pop2 and Pop4 at time t1. The rate of admixture of Pop2 and Pop3 was set as “ra”, and that of Pop3 and Pop4 was set as “1-ra”. Pop2 merged with Pop1 at time t2 and finally with Pop4 at time t3. NA reached sufficient population size at t4. In both scenarios 1 and 2, T. amurensis and T. japonica completed their divergence at time t3. Finally, in scenario 3, the four population groups merged simultaneously at t1. The prior values were set for the effective population size and divergence time estimated with a uniform distribution for all the parameters (Supplementary Table S2). The effective population size for each population group (N1, N2, N3, and N4) ranged from 10 to 20,000, with the constraint that all values were greater than or equal to the ancestral population size (NA). The divergence time (td) was set to a minimum of 10 generations and a maximum of 10,000 generations. Tilia plants can first flower when they are 12–40 years old, and the trees can live for >450 years (Pigott, 2012).With this in mind, we used a conservative generation time of 100 years. For each scenario, a mutation model was assumed with Mean mutation rate ranging from 10⁻⁴ to 10⁻³, Individual locus mutation rate ranging from 10⁻⁵ to 10⁻², Mean coefficient P ranging from 10⁻¹ to 3×10⁻¹, Individual locus coefficient P ranging from 10⁻² to 9, Mean SNI rate ranging from 10⁻⁸ to 10⁻⁴, and Individual locus SNI rate ranging from 10⁻⁹ to 10⁻³. Single-nucleotide insertion/deletion (SNI) rates were derived from 10,000 simulations generated for each putative scenario. The posterior probabilities were determined following 1,000,000 simulations for each scenario. The goodness-of-fit of the three scenarios was also evaluated by principal component analysis (PCA) using the option “model checking” in DIYABC.

Figure 2

Figure 2. Bar plots for eleven Tilia populations at K = 2 and K = 3 according to STRUCTURE analysis. Descriptions of abbreviations for the study sites can be found in Table 1.

Figure 3

Figure 3. (A) Principal coordinate analysis (PCoA) of eleven Tilia populations. (B) Phylogenetic trees of eleven Tilia populations. Descriptions of abbreviations for the study sites can be found in Table 1.

Figure 4

Figure 4. Three scenarios of population group demographic history examined via DIYABC analysis of T. amurensis and T. japonica; t# is the time scale measured for the generation. N1–4 and NA are the effective population sizes of the corresponding population groups, that is, Pop1 (DGL, MPS and DXG), Pop2 (YDS and ZS), Pop3 (TS, LU and LS), Pop4 (HS, MS and TMS) and the ancestral population.

3 Results

3.1 Genetic diversity

The results of the SSR data analysis revealed the level of genetic diversity in the T. amurensis and T. japonica samples (Table 2). The level of genetic diversity in T. japonica was determined from Shannon’s information index (I), the observed heterozygosity (H_O), and the expected heterozygosity (H_E). The Shannon’s information index (I) of the T. amurensis populations varied between 0.889 and 1.231, whereas that of T. japonica populations ranged from 0.915 to 1.381. The observed heterozygosity (H_O) of T. amurensis populations ranged from 0.329 to 0.647, and the expected heterozygosity (H_E) ranged from 0.440 to 0.588. The observed heterozygosity (H_O) of T. japonica populations ranged from 0.426 to 0.606, and the expected heterozygosity (H_E) ranged from 0.520 to 0.657. The fixation indices (F_IS) were positive for all three T. japonica populations (Table 2).

Table 2

Table 2. Genetic diversity parameters among the three populations of T. amurensis and T. japonica.

3.2 Genetic structure

The analysis of molecular variance (AMOVA) was implemented on the basis of SSR data and revealed the molecular variation attributable to differentiation among and within the populations for T. amurensis and T. japonica (Table 3). For T. amurensis, 12% of genetic variation exists among populations, while 88% exists within populations. In contrast, for T. japonica, 14% of the genetic variation was observed among populations, with the remainder (84%) existing within populations. All the results were highly significant (P < 0.001). The genetic differentiation coefficient (F_ST) among the populations of T. amurensis was 0.073, and the gene flow (N_m) was 3.191 (Table 3). The T. japonica populations, on the other hand, presented relatively high genetic differentiation (F_ST = 0.094) coefficients and relatively low gene flow (N_m = 2.403) (Table 3).

Table 3

Table 3. Hierarchical partitioning of genetic variance in Tilia using AMOVA.

To better understand the relationships between T. amurensis and T. japonica, the gene flow (N_m) and genetic differentiation coefficient (F_ST) were estimated for all pairs of populations (Table 4). The highest gene flow (N_m = 14.568) and smallest genetic differentiation coefficient (F_ST = 0.017) were observed between ZS and YDS. The lowest gene flow (N_m = 0.670) and greatest genetic differentiation coefficient (F_ST = 0.272) were detected between DGL and TMS.

Table 4

Table 4. Estimates of pairwise N_m (above diagonal) and F_ST (below diagonal) values between Tilia populations.

The Bayesian analysis implemented in STRUCTURE revealed a peak at K = 2. The K = 2 model was supported by the highest ΔK value (Figure 5), indicating that the most likely number of genetic clusters was 2 (ΔK = 250.88). In contrast, a lower ΔK value was observed at K = 3 (132.47). The bar plot for Tilia at K = 2 and 3 is shown in Figure 2. At K = 2, the 8 T. amurensis populations and the 3 T. japonica populations were clearly divided into two groups, with all T. amurensis populations categorized in the first group (red) and all T. japonica populations categorized in the second group (green). At K = 3, the populations show admixture between different genetic clusters, this being most evident within the MS population.

Figure 5

Figure 5. The ΔK value calculated by Structure 2.3.4 shows the maximum at K = 2.

To validate the results of the population structure analysis, phylogenetic tree were constructed for eleven Tilia populations, and PCoA was performed (Figure 3). Both analyses showed that populations from TS, LU and LS clustered into one taxon, those from YDS and ZS clustered into one taxon, and those from DGL, MPS and DXG clustered into another taxon. The three T. japonica populations (MS, HS and TMS) are clustered together and clearly separated from the T. amurensis population. The results were in close agreement with the structural analysis (Figure 2).

3.3 Demographic history

The highest posterior probability after evaluation using the DIYABC analysis of Tilia was scenario 1, with p=0.8544 (Table 5; Figure 6). In Scenario 1, the ancestral population first diverged to give rise to the T. amurensis and T. japonica groups, with the T. amurensis group subsequently differentiating into three distinct population groups. The DIYABC model was evaluated via principal component analysis (PCA) in scenario 1 on the basis of 1000,000 simulations showing the fit between the observed and simulated datasets (Supplementary Figure S1). The posterior distribution of each effective population size was also analyzed for the four populations, Pop1, Pop2, Pop3 and Pop4, with median values of 3810, 5840, 9340 and 7610, respectively. The median value of NA was 75.8 for the ancestral population. The mean values of the divergence times t1, t2 and t3 were 287, 416, and 2090 generations (approximately 28.7, 41.6, 209 ka BP), respectively. The time scale of the ancestral population size change, t4, was 4880 generations (approximately 488 ka BP). The prior and posterior distribution patterns are presented in Supplementary Figure S2. The mean mutation rates of microsatellites and single-nucleotide insertions/deletions (SNIs) at the examined loci were 8.13×10–⁴ and 8.12×10^-6, respectively.

Table 5

Table 5. Posterior probability of three putative scenarios in the DIYABC analysis of Tilia.

Figure 6

Figure 6. Logistic regression plot for the simulated scenarios in the DIYABC analysis of Tilia.

4 Discussion

4.1 Genetic diversity

Based on the statistical analysis of hundreds of research papers on genetic variation and population genetic differentiation in different types of plants by Hamrick and Godt (1990), the average expected heterozygosity (H_O) for dicotyledonous plants is 0.136. For long-lived perennial woody plants, the average H_O is 0.177. For plants with an allogamous breeding system and animal-dependent pollination, the average H_O is 0.167. Combined with our previous findings, T. japonica (I = 1.181, H_O = 0.484, H_E = 0.602) and T. amurensis (I = 0.982, H_O = 0.516, H_E = 0.515) exhibit comparable levels of genetic diversity, both showing relatively high genetic variation at the population level. It is noteworthy that the conclusions drawn by Hamrick and Godt were based on allozyme data. However, when compared with genetic diversity revealed by molecular markers in species such as the ecologically similar Picea obovata (H_O = 0.408; H_E = 0.423) (Kravchenko et al., 2016) and congeneric species including Tilia sibirica (H_E = 0.318) (Logan et al., 2018)and Tilia tomentosa (I = 0.50) (Gabur et al., 2019), T. japonica and T. amurensis still demonstrate the relatively high level of genetic diversity. According to the statistical analysis by Hamrick and Godt (1990) of hundreds of research papers on genetic variation and population genetic differentiation in different types of plants, the factors influencing the magnitude of genetic variation at the species level are, in order: taxonomic status, distribution range, life form, breeding system, and seed dispersal mechanism. At the population level, the influencing factors are: breeding system, distribution range, life history, taxonomic status, and seed dispersal mechanism. Woody tree species with a wide distribution range, long lifespan, predominantly allogamous breeding system, entomophilous pollination, frequent gene exchange, and abundant seeds that occur in late successional communities exhibit high levels of genetic variation (Hamrick and Godt, 1990). Woody tree species with a wide distribution range, long lifespan, predominantly allogamous breeding system, entomophilous pollination, and frequent gene exchange exhibit high levels of genetic variation (Hamrick and Godt, 1990).

The high genetic diversity of T. japonica may be determined by its life history and reproductive system characteristics. Tóth et al. (2022) also suggested that species of the genus Tilia are able to maintain a high level of genetic diversity due to their abundance during warmer periods, outcrossing mating system, and increased lifespan. Entomophilous pollination of Tilia and the intertwining of intraspecific differentiation and germplasm infiltration phenomena provide opportunities for gene recombination, which leads to high levels of gene flow in T. japonica (Ren and Ya, 1995). The reproductive system of T. japonica, which is predominantly characterized by sexual reproduction via entomophily, ensures a high level of genetic diversity within its populations (Boo and Park, 2016; Barker et al., 2022). Compared with the narrow distribution ranges of its congeners T. sibirica and T. argentea, T. japonica has a relatively wide distribution in eastern China and Japan. This broader distribution range may have led to higher levels of genetic diversity in T. japonica than in other plants of the same genus. Even though the Ne value was lower than the Na value in all T. japonica populations, genetic drift did not lead to a reduction in genetic diversity levels. Our results revealed that the genetic diversity of the TMS population was significantly lower than that of the other T. japonica populations. The inbreeding coefficients of T. japonica populations were greater than 0, and that of the TMS (F_IS = 0.221) population was the largest (Table 2), suggesting the possibility of selfing in T. japonica. Selfing reduces genetic diversity by decreasing genetic variation and weakening gene flow among populations, which may account for the relatively low genetic diversity observed in TMS populations (Jullien et al., 2019; Clo et al., 2019).

4.2 Genetic differentiation and introgression

According to the clear rules of genetic differentiation between populations (F_ST < 0.05, low; 0.05< F_ST < 0.15, medium; F_ST >0.15, high) defined by Wright (1978) and the results of the current study, there is a moderate degree of genetic differentiation between populations of both T. amurensis (F_ST = 0.073) and T. japonica (F_ST = 0.094). Moreover, AMOVA revealed that most of the genetic variation was distributed within the populations of T. amurensis and T. japonica. Genetic differentiation among populations is strongly influenced by gene flow and genetic drift (Andrews, 2010). When N_m > 1, there is sufficient gene flow between populations to exert a homogenizing effect, counteracting genetic drift and preventing differentiation. When N_m > 4, there is high gene flow between populations, and the populations can be regarded as a single unit of random mating (Wright, 1931). Both populations of T. amurensis (N_m = 3.191) and T. japonica (N_m = 2.403) exhibit strong gene flow, which can reduce local variation and prevent adaptive differentiation, thereby leading to lower genetic differentiation among their populations (Xiang et al., 2015; Wang et al., 2025). According to the study by Hamrick and Loveless (2019), the level of gene flow in plants is significantly influenced by their reproductive systems and seed dispersal mechanisms. Specifically, the mode of plant reproduction (such as selfing or outcrossing) and the mode of seed dispersal (such as wind dispersal or animal dispersal) are key factors in determining the intensity of gene flow. For instance, plants that predominantly rely on outcrossing typically exhibit higher levels of gene flow, whereas selfing plants tend to have relatively lower gene flow (Ballesteros-Mejia et al., 2016; Noreen et al., 2016; Gamba and Muchhala, 2023). Tilia species primarily reproduce through outcrossing via insect-mediated pollination, which can be considered the main factor contributing to the gene flow and population structure of populations of T. amurensis and T. japonica.

The results of the STRUCTURE analysis of 242 individuals revealed that the T. amurensis populations and the T. japonica populations could be well distinguished when the K value was 2, consistent with the two clustering results revealed by the phylogenetic tree and PCoA analysis. All individuals of T. amurensis have genetic components from T. japonica, and all individuals of the MS population in T. japonica, as well as several individuals of HS and TMS, also experienced genetic infiltration from T. amurensis. Relative to the HS and TMS populations, due to geographical proximity to the distribution area of T. amurensis, the MS population and T. amurensis populations exhibit greater gene flow and a lower degree of differentiation. This may be a factor leading to the MS population experiencing greater more infiltration of genes from T. amurensis. When the value of K is 3, all individuals of T. amurensis have genetic components from T. japonica, and all individuals of T. japonica have genetic components from T. amurensis. Among them, DGL, MPS, and DXG populations had the lowest degree of gradual infiltration, while easily distinguishable from other T. amurensis populations. The main reason is also that DGL, MPS, and DXG populations are geographically furthest away from T. japonica distribution area. The presence of significant gene infiltration between the two species implies that the two species may have originated from the same ancestor (Sork, 2016; Ribicoff et al., 2025).

4.3 Pleistocene split between T. amurensis and T. japonica

The ABC analysis in DIYABC revealed that T. amurensis and T. japonica share the same ancestral origin in this scenario, after which T. amurensis diverged into three population groups. Assuming a generation time of 100 years for Tilia and considering the summary statistics, the parameter set and the SSR markers used, we can infer a Middle Pleistocene split between T. amurensis and T. japonica, approximately 2090 generations ago (approximately 209 ka BP). In the Late Pleistocene, the T. amurensis population group in Northeast China (DGL, MPS and DXG) diverged from those in North China (ZS and YDS), approximately 416 generations ago (approximately 41.6 ka BP). Additionally, divergence from the Shandong population group (TS, LU and LS) occurred approximately 287 generations ago (approximately 28.7 ka BP). Some of the earliest evidence for putative Tilia species dates back to the Tertiary (Wolfe and Wehr, 1987; Pigott, 2012), and the approximate age of Tilia was inferred to be 17 million years (Richardson et al., 2015). A recently approximated age for the divergence of Tilia was inferred to be approximately 447 ka BP, also during the Middle Pleistocene (Logan et al., 2018). Some studies suggest that climatic fluctuations from the Middle to Late Pleistocene strongly affected tree species, causing lineage divergence. For example, Quercus spp diverged during the Early Pleistocene, Populus alba and P. davidiana diverged in the Middle Pleistocene as they spread southward from the Qinling Mountains, and the North American tree species Populus balsamifera and P. trichocarpa diverged during the Late Pleistocene at approximately ∼75 ka BP (Levsen et al., 2012; Richardson et al., 2015; Hou and Li, 2022). During the Middle and Late Pleistocene, the pronounced alternation of glacial and interglacial periods and the expansion and contraction of the ice sheet greatly influenced the population dynamics of Tilia (Groisman et al., 2013). During the Ice Age, when the climate is cold and dry, the distribution of linden populations shrinks to lower altitudes and latitudes, while during the interglacial period, when the climate warms, Tilia populations expand rapidly (Bolikhovskaya and Shunkov, 2014). T. amurensis taxa that had been migrating southward since the Miocene were strongly affected in their genetic structure during this repeated expansion and contraction, leading to genetic isolation and differences, which may have subsequently led to a Middle/Late Pleistocene split. The formation of today’s Tilia range is the result of redispersal after the end of the most recent ice age (approximately 15 ka BP).

4.4 Conservation implications

Both T. amurensis and T. japonica are vulnerable to habitat loss and overexploitation. Habitat fragmentation caused by anthropogenic disturbances threatens the survival of Tilia (Mu and Liu, 2007; Zhang et al., 2024). Currently, there are no effective conservation measures in place for T. amurensis and T. japonica. Our results indicate that the level of genetic diversity is high in T. japonica and that genetic variation occurs mainly within populations of T. amurensis and T. japonica. Therefore, in situ conservation of these populations is the most effective strategy to preserve their species diversity. The HS and MS populations of T. japonica maintain high genetic diversity and should be prioritized for in situ conservation. The complete biological conservation strategy should combine in situ and translocated conservation (Volis and Blecher, 2010; Dulloo et al., 2010; Duman et al., 2024). Therefore, to avoid inbreeding decline and population degradation, relocation protection and seed collection should also be strengthened. The shared ancestral background and the presence of bidirectional gene introgression between T. amurensis and T. japonica suggest that genetic connectivity between these two species should be considered in their conservation and management. To prevent population degradation and loss of genetic resources due to excessive hybridization between the two species, the protection of the three T. amurensis populations in Northeast China (DGL, MPS and DXG) and the TMS population of T. japonica, which have a purer genetic structure with less genetic drift, should be strengthened. Furthermore, conserving the genetic diversity of these two species requires attention not only to their respective core distribution areas but also to their zones of distributional overlap, which represent critical regions for gene flow.

5 Conclusion

We observed significant genetic differences and bidirectional gene infiltration through genetic studies of T. amurensis and T. japonica. Considering a generation time of 100 years, we estimate by DIYABC that the splitting of T. amurensis and T. japonica likely occurred in the Middle Pleistocene. The differentiation of the three T. amurensis population groups also occurred in the Late Pleistocene. Combined with the results of previous studies, our findings revealed that T. amurensis and T. japonica presented high levels of genetic diversity. However, given that both T. amurensis and T. japonica are facing dual pressures of habitat loss and overexploitation, we still strongly recommend a focus on conservation work for both species. The natural distribution areas of T. amurensis and T. japonica should be restored, and protection and management of these species should continue.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

YZ: Conceptualization, Methodology, Software, Writing – original draft, Writing – review & editing. FZ: Conceptualization, Methodology, Software, Writing – original draft. DL: Methodology, Writing – review & editing. YM: Software, Writing – review & editing, Methodology. YL: Writing – review & editing, Software. QW: Conceptualization, Data curation, Formal Analysis, Supervision, Writing – review & editing. DZ: Conceptualization, Funding acquisition, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This study was supported by the Subject of Key R & D Program of Shandong Province (2024LZGC003), the Subject of Key R & D Plan of Shandong Province (Major Scientific and Technological Innovation Project) “Mining and Accurate Identification of Forest Tree Germplasm Resources” (2021LZGC023), 2021 Shandong Province Higher Education Institutions “Qingchuang Science and Technology Support Program” Project: Species Conservation and Biological Invasion Innovation Team (LuKeJiaoHan [2021] No. 51).

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1651814/full#supplementary-material

References

Andrews, C. A. (2010). Natural selection, genetic drift, and gene flow do not act in isolation in natural populations. Nat. Educ. Knowledge 3, 5.

Google Scholar

Ballesteros-Mejia, L., Lima, N. E., Lima-Ribeiro, M. S., and Collevatti, R. G. (2016). Pollination mode and mating system explain patterns in genetic differentiation in neotropical plants. PloS One 11, e0158660. doi: 10.1371/journal.pone.0158660

PubMed Abstract | Crossref Full Text | Google Scholar

Barker, C., Davis, M. L., and Ashton, P. (2022). Reproductive strategy of a temperate canopy tree Tilia cordata Mill.(Malvaceae) is related to temperature during flowering and density of recent recruits. Tree Genet. Genomes 18, 22.

Google Scholar

Berry, E. W. (1918). The history of the linden and ash. Plant World 21, 163–175.

Google Scholar

Bolikhovskaya, N. and Shunkov, M. (2014). Pleistocene environments of northwestern Altai: Vegetation and climate. Archaeology Ethnology Anthropology Eurasia 42, 2–17. doi: 10.1016/j.aeae.2015.01.001

Crossref Full Text | Google Scholar

Boo, D. and Park, S. J. (2016). Molecular phylogenetic study of Korean Tilia L. Korean J. Plant Resour. 29, 547–554. doi: 10.7732/kjpr.2016.29.5.547

Crossref Full Text | Google Scholar

Chen, H. (1990). Flora of Shandong (Science Press).

Google Scholar

Chen, B., Zou, H., Zhang, B., Zhang, X., Wang, C., Zhang, X., et al. (2023). Distribution change and protected area planning of Tilia amurensis in China: A study of integrating the climate change and present habitat landscape pattern. Global Ecol. Conserv. 43, e02438. doi: 10.1016/j.gecco.2023.e02438

Crossref Full Text | Google Scholar

Clo, J., Gay, L., and Ronfort, J. (2019). How does selfing affect the genetic variance of quantitative traits? An updated meta-analysis on empirical results in angiosperm species. Evolution 73, 1578–1590. doi: 10.1111/evo.13789

PubMed Abstract | Crossref Full Text | Google Scholar

Committee F o C E (2007). Flora of China (Hippocastanaceae through Theaceae) Vol. 12. Eds. Wu, C. Y., Raven, P. H., and Hong, D. Y. (Beijing & St. Louis: Science Press & Missouri Botanical Garden Press), 1–302. Fl. China.

Google Scholar

Cornuet, J. M., Pudlo, P., Veyssier, J., Dehne-Garcia, A., Gautier, M., Leblois, R., et al. (2014). DIYABC v2. 0: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics 30, 1187–1189. doi: 10.1093/bioinformatics/btt763

PubMed Abstract | Crossref Full Text | Google Scholar

Cui, L. and Mu, L. Q. (2016). Ectomycorrhizal communities associated with Tilia amurensis trees in natural versus urban forests of Heilongjiang in Northeast China. J. Forestry Res. 27, 401–406. doi: 10.1007/s11676-015-0158-1

Crossref Full Text | Google Scholar

Dulloo, M. E., Hunter, D., and Borelli, T. (2010). Ex situ and in situ conservation of agricultural biodiversity: major advances and research needs. Notulae Botanicae Horti Agrobotanici Cluj-Napoca 38, 123–135.

Google Scholar

Duman, H., Doğan, M., Atlı, Ö., and Celep, F. (2024). Ex situ and in situ conservation approaches in species-rich Anatolian steppe ecosystem: A case study from Ankara, Türkiye. Ecologies 5, 664–678. doi: 10.3390/ecologies5040039

Crossref Full Text | Google Scholar

Engler, V. (1909). Monographie der Gattung Tilia: aus d. Botan. Ed. Korn, D. v. W. G. (Garten d. Univ. Breslau).

Google Scholar

Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x

PubMed Abstract | Crossref Full Text | Google Scholar

Gabur, I., Lipsa, F. D., Adumitresei, L., Tanase, C., and Simioniuc, D. P. (2019). Assessment of genetic variation of Tilia tomentosa by RAPD markers. J. Plant Dev. 26, 85. doi: 10.33628/jpd.2019.26.1.85

Crossref Full Text | Google Scholar

Gamba, D. and Muchhala, N. (2023). Pollinator type strongly impacts gene flow within and among plant populations for six Neotropical species. Ecology 104, e3845. doi: 10.1002/ecy.3845

PubMed Abstract | Crossref Full Text | Google Scholar

Groisman, P. Y., Blyakharchuk, T. A., Chernokulsky, A. V., Arzhanov, M. M., Marchesini, L. B., Bogdanova, E. G., et al. (2013). “Climate changes in siberia,” in Regional environmental changes in Siberia and their global consequences, 57–109.

Google Scholar

Hamrick, J. L. and Godt, M. W. (1990). Allozyme diversity in plant species.

Google Scholar

Hamrick, J. and Loveless, M. (2019). “The genetic structure of tropical tree populations: associations with reproductive biology,” in The evolutionary ecology of plants (CRC Press), 129–146.

Google Scholar

Hedrick, P. W. and Miller, P. S. (1992). Conservation genetics: techniques and fundamentals. Ecol. Appl. 2, 30–46. doi: 10.2307/1941887

PubMed Abstract | Crossref Full Text | Google Scholar

Hou, Z. and Li, A. (2022). Genomic differentiation and demographic histories of two closely related Salicaceae species. Front. Plant Sci. 13, 911467. doi: 10.3389/fpls.2022.911467

PubMed Abstract | Crossref Full Text | Google Scholar

Jullien, M., Navascués, M., Ronfort, J., Loridon, K., and Gay, L. (2019). Structure of multilocus genetic diversity in predominantly selfing populations. Heredity 123, 176–191. doi: 10.1038/s41437-019-0182-6

PubMed Abstract | Crossref Full Text | Google Scholar

Kravchenko, A. N., Ekart, A. K., and Larionova, A. Y. (2016). Genetic diversity and differentiation of Siberian spruce populations at nuclear microsatellite loci. Russian J. Genet. 52, 1142–1148. doi: 10.1134/S1022795416090088

PubMed Abstract | Crossref Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | Crossref Full Text | Google Scholar

Levsen, N. D., Tiffin, P., and Olson, M. S. (2012). Pleistocene speciation in the genus Populus (Salicaceae). Systematic Biol. 61, 401. doi: 10.1093/sysbio/syr120

PubMed Abstract | Crossref Full Text | Google Scholar

Lindenmayer, D. B. and Fischer, J. (2013). Habitat fragmentation and landscape change: an ecological and conservation synthesis (Island Press).

Google Scholar

Logan, S. A., Chytrý, M., and Wolff, K. (2018). Genetic diversity and demographic history of the Siberian lime (Tilia sibirica). Perspect. Plant Ecology Evol. Systematics 33, 9–17. doi: 10.1016/j.ppees.2018.04.005

Crossref Full Text | Google Scholar

Mu, L. Q. and Liu, Y. N. (2007). Genetic diversity of Tilia amurensis populations in different geographical distribution regions. Chin. J. Plant Ecol. 31, 1190.

Google Scholar

Noh, J. K., Echeverria, C., Gaona, G., Kleemann, J., Koo, H., Fürst, C., et al. (2022). Forest ecosystem fragmentation in Ecuador: Challenges for sustainable land use in the Tropical Andean. Land 11, 287. doi: 10.3390/land11020287

Crossref Full Text | Google Scholar

Noreen, A., Niissalo, M., Lum, S., et al. (2016). Persistence of long-distance, insect-mediated pollen movement for a tropical canopy tree species in remnant forest patches in an urban landscape. Heredity 117, 472–480. doi: 10.1038/hdy.2016.64

PubMed Abstract | Crossref Full Text | Google Scholar

Peakall, R. and Smouse, P. E. (2006). GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 6, 288–295. doi: 10.1111/j.1471-8286.2005.01155.x

PubMed Abstract | Crossref Full Text | Google Scholar

Pigott, D. (2012). Lime-trees and basswoods: a biological monograph of the genus Tilia (Cambridge University Press).

Google Scholar

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959. doi: 10.1093/genetics/155.2.945

PubMed Abstract | Crossref Full Text | Google Scholar

Ren, Z. and Ya, T. (1995). Morphological evolution and biogeography of Tilia. J. Southwest Forestry Univ. 15, 1–14.

Google Scholar

Ribicoff, G., Garner, M., Pham, K., Althaus, K. N., Cavender‐Bares, J., Crowl, A. A., et al. (2025). Introgression, Phylogeography, and Genomic Species Cohesion in the Eastern North American White Oak Syngameon. Mol Ecol.[J], 34, e17822.

PubMed Abstract | Google Scholar

Richardson, J. E., Whitlock, B. A., Meerow, A. W., and Madriñán, S. (2015). The age of chocolate: a diversification history of Theobroma and Malvaceae. Front. Ecol. Evol. 3, 120. doi: 10.3389/fevo.2015.00120

Crossref Full Text | Google Scholar

Smouse, R. P. P. and Peakall, R. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 28, 2537–2539.

PubMed Abstract | Google Scholar

Sork, V. L. (2016). Gene flow and natural selection shape spatial patterns of genes in tree populations: implications for evolutionary processes and applications. Evolutionary Appl. 9, 291–310. doi: 10.1111/eva.12316

PubMed Abstract | Crossref Full Text | Google Scholar

Tóth, E. G., Szilágyi, K., Patyi, A., and György, Z. (2022). Genetic diversity in a historic lime tree allée of Széchenyi Castle in Nagycenk, Hungary. Genet. Resour. Crop Evol. 69, 1407–1418.

Google Scholar

Volis, S. and Blecher, M. (2010). Quasi in situ: a bridge between ex situ and in situ conservation of plants. Biodiversity Conserv. 19, 2441–2454. doi: 10.1007/s10531-010-9849-2

Crossref Full Text | Google Scholar

Wang, J., Tang, L., Yu, X., Yu, C., Tang, X., Wang, D., et al. (2025). Gene flow enhances genetic diversity and local adaptation in Pyropia yezoensis populations. Water Biol. Secur., 100411. doi: 10.1016/j.watbs.2025.100411

Crossref Full Text | Google Scholar

Weir, B. S. and Cockerham, C. C. (1984). Estimating F-statistics for the analysis of population structure. Evolution, 1358–1370. doi: 10.2307/2408641

PubMed Abstract | Crossref Full Text | Google Scholar

Wolfe, J. A. and Wehr, W. (1987). Middle Eocene dicotyledonous plants from Republic, northeastern Washington (USGPO).

Google Scholar

Wright, S. (1931). Evolution in Mendelian populations. Genetics 16, 97. doi: 10.1093/genetics/16.2.97

PubMed Abstract | Crossref Full Text | Google Scholar

Wright, S. (1978). Evolution and the genetics of populations: a treatise in four volumes: Vol. 4: variability within and among natural populations (University of Chicago Press).

Google Scholar

Wu, Q., Zhang, Y., Xie, X., Tong, B., Liu, D., Ma, Y., et al. (2024). Analysis of the genetic diversity and population structure of Tilia amurensis from China using SSR markers: Implications for conservation. Global Ecol. Conserv. 54, e03173. doi: 10.1016/j.gecco.2024.e03173

Crossref Full Text | Google Scholar

Xiang, X. Y., Zhang, Z. X., Duan, R. Y., Zhang, X. P., and Wu, G. L. (2015). Genetic diversity and structure of Pinus dabeshanensis revealed by expressed sequence tag-simple sequence repeat (EST-SSR) markers. Biochem. Systematics Ecol. 61, 70–77. doi: 10.1016/j.bse.2015.06.001

Crossref Full Text | Google Scholar

Ya, T. and Ren, Z. (1996). Geographical distribution of Tilia Linn. Acta Phytotaxonomica Sin. 34, 254–264.

Google Scholar

Yue, Y., Yan, L., Huang, X., Wang, H., and Tang, S. (2022). Development EST-SSR markers based on transcriptome sequences of Tilia miqueliana maxim. Mol. Plant Breed., 1–8.

Google Scholar

Zhang, F. G., Zhang, S., Wu, K., Zhao, R., Zhao, G., Wang, Y., et al. (2024). Potential habitat areas and priority protected areas of Tilia amurensis Rupr in China under the context of climate change. Front. Plant Sci. 15, 1365264. doi: 10.3389/fpls.2024.1365264

PubMed Abstract | Crossref Full Text | Google Scholar