Potential Suitable Habitat Prediction and Distribution Patterns of Primula L. in China Under Climate Change

Climate change is increasingly reshaping species habitat suitability worldwide. Primula L., the largest genus in Primulaceae, comprises 404 species in China (including 296 endemic species) and is characterized by high endemism and numerous rare and endangered taxa. However, global warming has intensified habitat fragmentation and loss, while its distribution patterns and key environmental drivers remain insufficiently understood. We compiled 7647 occurrence records of 404 wild Primula species in China and integrated 60 environmental variables (climatic, topographic, and soil factors). Using the MaxEnt model combined with ArcGIS spatial analysis, we assessed current and future habitat suitability, identified dominant environmental drivers, and quantified conservation gaps under multiple climate scenarios. Species richness is highly concentrated in Sichuan (186 species), Yunnan (177 species), and Xizang (165 species), with the Hengduan Mountains and eastern Himalayas representing the core distribution area and showing clear peripheral differentiation. The optimized MaxEnt model performed well (AUC = 0.858), identifying temperature seasonality (bio4, 39.8%) and elevation (27.1%) as the main limiting factors. The total suitable habitat area is 268.52 × 10⁴ km², with high-suitability areas mainly distributed in the Hengduan Mountains, southeastern Qinghai-Xizang Plateau, and the Central Mountain Range of Taiwan. Under three shared socioeconomic pathway (SSP) scenarios (SSP126, SSP245, and SSP585), suitable habitat shows a persistent decline, most pronounced under SSP585 in the 2090s (-20.73%), accompanied by a 25.86% reduction in low-suitability areas. Localized expansion of high-suitability habitats suggests that the Hengduan Mountains and southeastern Qinghai-Xizang Plateau may act as potential climatic refugia. Habitat loss consistently exceeds habitat gain, while the distribution centroid shifts westward and northwestward, with migration distances increasing under higher-emission scenarios. Conservation gap analysis indicates that 90.01% of high-suitability habitats lie outside the current protected area network, revealing a strong mismatch between biodiversity hotspots and conservation coverage. These findings highlight the urgent need to expand protected areas and establish micro-reserves in key gap regions (southwestern Sichuan, northwestern Yunnan, southeastern Xizang, and southern Gansu), and to integrate climate-driven migration corridors into conservation planning to support long-term alpine plant persistence under climate change.