Informed species delimitation is crucial in diverse biological fields; however, it can be problematic for species complexes. Showing a peripatric distribution pattern, Stewartia gemmata and S. acutisepala (the S. gemmata complex) provide us with an opportunity to study species boundaries among taxa undergoing nascent speciation. Here, we generated genomic data from representative individuals across the natural distribution ranges of the S. gemmata complex using restriction site-associated DNA sequencing (RAD-seq). Based on the DNA sequence of assembled loci containing 41,436 single-nucleotide polymorphisms (SNPs) and invariant sites, the phylogenetic analysis suggested strong monophyly of both the S. gemmata complex and S. acutisepala, and the latter was nested within the former. Among S. gemmata individuals, the one sampled from Mt. Tianmu (Zhejiang) showed the closest evolutionary affinity with S. acutisepala (which is endemic to southern Zhejiang). Estimated from 2996 high-quality SNPs, the genetic divergence between S. gemmata and S. acutisepala was relatively low (an F_st of 0.073 on a per-site basis). Nevertheless, we observed a proportion of genomic regions showing relatively high genetic differentiation on a windowed basis. Up to 1037 genomic bins showed an F_st value greater than 0.25, accounting for 8.31% of the total. After SNPs subject to linkage disequilibrium were pruned, the principal component analysis (PCA) showed that S. acutisepala diverged from S. gemmata along the first and the second PCs to some extent. By applying phylogenomic analysis, the present study determines that S. acutisepala is a variety of S. gemmata and is diverging from S. gemmata, providing empirical insights into the nascent speciation within a species complex.

As climate change accelerates, plant species largely rely on genetic variation to adapt and survive when they fail to track their ecological niches through range shifts. Predicted genomic vulnerability is able to identify populations lacking the necessary genetic variation for climate change adaptation. However, the role of introgression in genomic vulnerability remains poorly explored. Here, we used the dove tree (Davidia involucrata), a relict species native to southwestern China, to test whether introgression may reduce genomic vulnerability. By integrating population genomics and environmental data collected from 196 individuals of 18 populations, we identified 747 strictly climate-associated loci across the distribution range of D. involucrata, 138 of which were recovered from the genetically admixed populations. We estimated the genomic vulnerability for three genetic lineages and two admixed groups using the gradient forest approach, and found that eastern populations are likely to be at higher risk. The eastern admixed populations exhibited a significant reduction, with introgression from the southern lineage. Cumulative importance analysis showed moderate importance for introgressive loci along environmental gradients. This indicates that the introduction of novel alleles through introgression provides only a partial and insufficient counterbalance to the maladaptation observed in D. involucrata under climate change. Our study highlights the role of intraspecific introgression in response to climate change and emphasizes the importance of genomic vulnerability studies in informing conservation practices for relict and endangered species.