ObjectiveTo evaluate the clinical efficacy of Shenqi Yangxin decoction in the treatment of ischemic cardiomyopathy （ICM） with Qi and Yin deficiency and blood stasis syndrome and its effect on serum hydrogen sulfide （H₂S） and calcium ion （Ca²⁺）. MethodsA total of 64 ICM patients with Qi and Yin deficiency and blood stasis syndrome who met the inclusion criteria were randomly divided into a control group （n=32） and a treatment group （n=32）. All patients received conventional Western medicine treatment. The treatment group was additionally given Shenqi Yangxin decoction. The TCM syndrome score， Minnesota Living with Heart Failure Questionnaire （MLHFQ） score， left ventricular ejection fraction （LVEF）， N-terminal pro-B-type natriuretic peptide （NT-proBNP）， 6-minute walk test （6MWT）， New York Heart Association （NYHA） cardiac function classification， and serum H₂S and Ca²⁺ levels were compared between the two groups pre- and post-treatment. ResultsTwo cases dropped out from each group during the study. Finally， 30 patients in each group were included in the analysis. There were no significant differences in age， gender， course of coronary heart disease， underlying diseases， and laboratory tests between the two groups. Compared with baseline， the TCM syndrome score， MLHFQ score， and NT-proBNP in both treatment group and control group decreased significantly （P<0.01）， LVEF， 6MWT， and H₂S increased significantly （P<0.01）， and serum Ca²⁺ increased （P<0.05）. Compared with the control group after treatment， the MLHFQ score and NT-proBNP in the treatment group decreased （P<0.05）， the TCM syndrome score decreased significantly （P<0.01）， LVEF， 6MWT， and serum Ca²⁺ increased （P<0.05）， and H₂S increased significantly （P<0.01）. The improvement degree of the NYHA cardiac function classification in the treatment group was higher than that in the control group， but there was no significant difference. ConclusionShenqi Yangxin decoction is effective in treating ICM patients with Qi and Yin deficiency and blood stasis， which could significantly improve cardiac function and quality of life， and its therapeutic effect may be related to the regulation of serum H₂S and Ca²⁺ levels.

Objective To identify factors affecting and key environmental factors of the Oncomelania hupensis snail density in the Yangtze River Delta region using machine learning methods. Methods Administrative village-level O. hupensis snail survey data in the Yangtze River Delta (including Shanghai Municipality, Jiangsu Province, Zhejiang Province and Anhui Province) from 2011 to 2021 were retrieved from the Information Management System for Parasitic Disease Control of Chinese Center for Disease Control and Prevention. Environmental factor data were captured from the Google Earth Engine platform, including elevation, slope, terrain, normalized difference vegetation index (NDVI), vegetation type, soil type, total petroleum hydrocarbon (TPH), ammonium nitrogen, inorganic nitrogen, dissolved oxygen, pH of water, chemical oxygen demand (COD) and inorganic phosphorus, and climatic factor data in the study region were retrieved from the Copernicus Climate Data Store, including annual precipitation, aridity index and annual mean temperature (AMT). O. hupensis snail survey data in the Yangtze River Delta region from 2011 to 2021 were randomly divided into a training set (70%) and a test set (30%), and five machine learning models were selected for machine learning model construction and comparative analysis of the O. hupensis snail density using the software R 4.3.0, including random forest (RF), eXtreme gradient boosting (XGBoost), support vector machine (SVM), gradient boosting machine (GBM) and neural network (NN). The XGBoost model was employed to construct a predictive model for the O. hupensis snail density, and the impact of each environmental factor on O. hupensis snail distribution was quantified. The SHapley Additive exPlanations (SHAPs) values were calculated to estimate the average contribution of each variable to the model prediction, and the core environmental factors affecting the O. hupensis snail population density were screened. Results Among the five machine learning models, the XGBoost model exhibited the optimal comprehensive performance, with the coefficient of determination (R²) of 0.855, mean squared error (MSE) of 0.188, root mean squared error (RMSE) of 0.434 and mean absolute error (MAE) of 0.155, respectively. Analysis of factors affecting the O. hupensis snail density with the XGBoost model showed that among the 16 environmental factors, the top four high-impact factors ranked by SHAPs values included annual precipitation, elevation, aridity index and NDVI, with cumulative SHAPs contributions of 75%, which was higher than that of other environmental factors. If NDVI was higher than 0.6, the O. hupensis snail density increased with NDVI and peaked if NDVI was 0.8 (1.60 snails/0.1 m²). The O. hupensis snail density increased with elevation if the elevation ranged from 14 to 40 m, and slowly rose if the annual precipitation ranged from 900 to 1 300 mm, and then increased rapidly to the peak (1.52 snails/0.1 m²) if the annual precipitation ranged from 1 300 to 1 500 mm. In addition, the O. hupensis snail density increased rapidly to the maximum (1.60 snails/0.1 m²) if the aridity index ranged from 0.8 to 1.1, and decreased gradually if the aridity index exceeded 1.1. Conclusions The XGBoost model shows excellent performance in prediction of the O. hupensis snail density and identification of key environmental factors in the Yangtze River Delta region. Annual precipitation, elevation, aridity index and NDVI are key environmental factors affecting the distribution and density of O. hupensis snails in the Yangtze River Delta region.

The frequent extreme climatic events post multifaceted impacts on the distribution of Oncomelania hupensis, the intermediate host of Schistosoma japonicum in the context of global climate change. This article systematically reviews the effects of four types of extreme climatic events, including floods, droughts, heat waves, and cold waves, on the survival, reproduction, and distribution of Oncomelania hupensis. Floods may expand suitable snail habitats, and increase both emerging and re-emerging snail habitats; however, the impact of floods on O. hupensis density is characterized by a lag effect of a decline followed by a rise. Droughts may cause fragmentation of suitable O. hupensis snail habitats, reduced O. hupensis snail egg production rates, and increased O. hupensis snail mortality, and heat waves may cause an increase in O. hupensis snail mortality, a reduction in numbers of O. hupensis snail populations and shrinking of O. hupensis snail distribution, while cold waves may cause a reduction in O. hupensis snail density and a rise in O. hupensis snail mortality. Extreme climate events pose both shortand long-term effects on the distribution of O. hupensis. Intensified surveillance of O. hupensis snails is required in high-risk environments.

ObjectiveTo observe the effects of Yangjing Zhongyutang （YJZYT） on mitochondrial biogenesis and oxidative stress damage mediated by the silent information regulator 1 （SIRT1）/peroxisome proliferator-activated receptor gamma coactivator-1alpha （PGC-1α） signaling pathway in cyclophosphamide （CTX）-induced rats with diminished ovarian reserve （DOR）， and to explore its mechanism in improving ovarian reserve function and follicular development. MethodsForty-two 8-week-old female SD rats with normal estrous cycles were randomly divided into a blank control group （n=7） and a model group （n=35）. Rats in the model group received a single intraperitoneal injection of CTX （90 mg·kg^-1） to establish the DOR model. After modeling， estrous cycles were monitored for 7 consecutive days， and model success was confirmed based on criteria for estrous cycle disruption. After successful modeling， rats were divided into groups for intervention： estradiol valerate group （0.09 mg·kg^-1）， and YJZYT high-， medium-， and low-dose groups （19.98， 9.99， 5.00 g·kg^-1）. The blank control group and model group were given an equal volume of distilled water by gavage. All groups received daily gavage once for 4 consecutive weeks. The general state， body weight， and ovarian wet weight of rats were observed and recorded， and the ovarian organ index was calculated. Enzyme-linked immunosorbent assay （ELISA） was used to measure serum levels of follicle-stimulating hormone （FSH）， luteinizing hormone （LH）， estradiol （E₂）， anti-Müllerian hormone （AMH）， superoxide dismutase （SOD）， and glutathione peroxidase （GSH-Px）. Hematoxylin-eosin （HE） staining was performed to observe ovarian histomorphological changes and follicular development status. Immunofluorescence was used to detect reactive oxygen species （ROS） expression levels. Colorimetric assays were employed to measure adenosine triphosphate （ATP） and malondialdehyde （MDA） content in ovarian tissues. Quantitative Real-time polymerase chain reaction （Real-time PCR） was used to detect mitochondrial DNA （mtDNA） copy number and the mRNA expression levels of key genes including SIRT1， PGC-1α， nuclear respiratory factor 1 （NRF1）， and mitochondrial transcription factor A （TFAM）. Western blot was performed to detect the protein expression levels of SIRT1， PGC-1α， NRF1， and TFAM. ResultsCompared with the blank group， rats in the model group exhibited disrupted estrous cycles， obviously reduced body weight， and decreased ovarian index （P<0.05）. Ovarian histopathology revealed cortical thinning， loose structure， and a significant reduction in both primordial and growing follicles （P<0.01）. Serum FSH and LH levels were significantly elevated （P<0.01）， while E₂ and AMH levels were obviously reduced （P<0.05， P<0.01）. ATP content and mtDNA copy number decreased in ovarian tissue （P<0.01）， ROS expression increased， MDA levels rose， while SOD and GSH-Px activities obviously decreased （P<0.05， P<0.01）， mRNA and protein expression levels of SIRT1， PGC-1α， NRF1， and TFAM were obviously downregulated （P<0.05， P<0.01）. After treatment， compared with the model group， body weight and ovarian index obviously recovered in rats administered various doses of YJZYT （P<0.05）， serum E₂ and AMH levels increased， while FSH and LH levels obviously decreased （P<0.05， P<0.01）， ovarian tissue ATP content and mtDNA copy number were up-regulated， ROS and MDA levels decreased， and antioxidant enzymes SOD and GSH-Px activity obviously increased （P<0.05， P<0.01）， Gene and protein expression levels related to the SIRT1/PGC-1α /NRF1/TFAM signaling pathway were obviously up-regulated compared to the model group （P<0.05， P<0.01）， HE staining revealed that ovarian structure gradually recovered to integrity in all treatment groups， with a obviously increase in the number of primordial and growing follicles （P<0.05， P<0.01）. Granulosa cells were neatly arranged， indicating marked improvement in ovarian function. ConclusionYJZYT may improve ovarian function and follicular development in rats with diminished ovarian reserve by activating the SIRT1/PGC-1α signaling pathway， promoting mitochondrial biogenesis， enhancing mitochondrial function， and alleviating oxidative stress damage.

Background: Influenza A virus (IAV) is a negative-sense RNA virus of the Orthomyxoviridae family and is the etiological agent of a highly contagious acute respiratory disease that can lead to acute lung injury. Objective: To elucidate the molecular mechanisms of IAV infection, an integrative research approach combining gene expression profiling, multinetwork analysis, and in vivo experimental validations was employed. Methods: First, a series of network-based analyses were performed, including protein-protein interaction network construction, weighted gene co-expression network analysis, and subsequent gene set enrichment analysis, to identify the major underlying mechanisms of IAV infection. Following gene expression analysis, core targets, both direct and indirect regulators, were screened. An IAV (H1N1) strain A/PR/8/34-induced acute lung injury mouse model was constructed for in vivo validations. Batch one included two groups to evaluate findings from the multi-network analysis: Mock (n = 10; 5 males and 5 females) and IAV (n = 10; 5 males and 5 females). Batch two included three groups to assess the role of toll-like receptor 3 (TLR3) in IAV infection: Mock (n = 6; 3 males and 3 females), IAV (n = 6; 3 males and 3 females), and TLR3 inhibitor (n = 6; 3 males and 3 females). Body weight was measured on days 0, 3, and 5 after infection. On day 5, lung tissues were collected to assess viral load and histopathological changes. Key targets were examined using enzyme-linked immunosorbent assay, Western blotting, and immunofluorescence staining, both in sera and lung tissues. Results: IAV infection was significantly associated with dysregulation of the immune-inflammation system, such as the LTR, nucle-otide-binding oligomerization domain-(NOD) like receptor, retinoic acid-inducible gene I-like receptor, and nuclear factor kappa-B signaling pathways. Gene set enrichment analysis further indicated that the TLR and necroptosis signaling pathways played crucial roles in the progression of IAV infection (TLR signaling pathway normalized enrichment score = 2.3941, P = 1.00 × 10 ⁻¹⁰; necroptosis normalized enrichment score = 1.9421, P = 6.21 × 10 ⁻⁷). Among the core targets, TLR3 and mixed lineage kinase domain-like protein (MLKL) may regulate gene expression at the transcriptional level (all P < 0.05). In vivo validation using an IAV (PR8) infected acute lung injury mouse model demonstrated increased viral load and lung index, alveolar structural damage, and inflammatory cell infiltration. Immunofluorescence staining exhibited large gaps in Lamin B1 staining and breaches in Emerin signals following IAV-PR8 infection. Expression levels of TLR3, p-receptor-interacting serine/threonine-protein kinase 3 (RIPK3)/RIPK3, and p-mixed lineage kinase domain-like protein (MLKL)/MLKL proteins in lung tissues, as well as proinflammatory factors and mediators in sera, were significantly elevated after IAV infection. Moreover, enhanced neutrophil infiltration (myeloperoxidase) and citrullinated histone H3 (a neutrophil extracellular trap-specific marker), both established indicators of neutrophil extracellular trap formation, were observed. Notably, treatment with a TLR3 inhibitor significantly ameliorated IAV-induced acute lung injury by regulating necroptosis-related targets. Conclusion: Our study provides network-based in vivo evidence that TLR3-receptor-interacting serine/threonine-protein kinase 3-MLKL-mediated necroptosis may underlie IAV-induced acute lung injury and could serve as a potential therapeutic target in severe influenza cases.

Objective To analyze the epidemiological characteristics of hand, foot, and mouth disease (HFMD) in Shanxi Province from 2012 to 2024, and predict the incidence trend for 2025, and to provide a scientific basis for the formulation of prevention and control strategies. Methods Based on the surveillance data of HFMD in Shanxi Province from 2012 to 2024, the spatial and temporal distribution characteristics and time trends of the disease were analyzed. The ARIMA model was constructed and used to predict the incidence trend in 2025. Results From 2012 to 2024, a total of 254 028 HFMD cases were reported in Shanxi Province, with an average annual incidence rate of 54.17 per 100 000 population, a severe case rate of 0.56%, and a case fatality rate (CFR) of 12.60 per 100 000 population. Joinpoint regression analysis showed that the incidence rate, severe case rate, mortality rate, and case fatality rate all presented a downward trend. The epidemic exhibited obvious seasonal distribution characteristics, with the peak period from April to November, and two incidence peaks in June-July and October-November. The male-to-female incidence ratio was 1.41:1. Children aged 1-5 years accounted for 89.24% of the total cases, among which scattered children (58.48%) and nursery children (33.54%) were the high-risk groups. Linfen City (96.06 /100 000) and Taiyuan City (88.54 /100 000) had relatively high incidence rates. After 2017, the proportion of enterovirus A71 (EV-A71) decreased, while coxsackievirus A16 (Cox-A16) and other enteroviruses became the main epidemic strains. The ARIMA(1,0,1)(0,1,1)₁₂ model predicted that the incidence of HFMD in 2025 would remain at the level of 2023-2024, and the dual-peak characteristic would continue. Conclusion From 2012 to 2024, the overall HFMD epidemic in Shanxi Province generally shows a significant downward trend. The high-risk population includes scattered children and nursery children under 5 years old, with high-incidence areas concentrated in the central and southern regions, requiring focused attention. The seasonal ARIMA model can effectively fit the evolutionary trend of HFMD incidence in Shanxi Province and possesses short-term predictive capability.

ObjectiveTo investigate the optimal termination time for acupuncture in treating patients with acute peripheral facial paralysis and its cost-effectiveness. MethodsA total of 120 eligible patients with acute-stage peri-pheral facial paralysis were randomly assigned to either the mild dysfunction termination group and the complete recovery termination group， with 60 patients in each group. Both groups received the standard acupuncture treatment protocol. Treatment in the mild dysfunction termination group was terminated when the Sunnybrook facial grade scale （SFGS） score first reached or exceeded 83 points， while that in the complete recovery termination group was terminated when the SFGS score first reached or exceeded 95 points. Assessments were conducted before treatment， 6 and 12 months after onset. SFGS， facial disability index （FDI） including physical function （FDIp） and social function （FDIs）， self-rating anxiety scale （SAS）， and self-rating depression scale （SDS） scores were assessed before treatment， and 6 and 12 months after onset. Any acupuncture-related adverse events during treatment were recorded for safety evaluation. Treatment sessions and medical costs including direct costs， indirect costs， insurance coverage， total societal costs， and patient out-of-pocket expenses were also recorded， and an economic evaluation was conducted including cost-effectiveness ratio （CER） and incremental cost-effectiveness ratio （ICER）. ResultsUltimately， 56 patients in the mild dysfunction termination group and 55 in the complete recovery termination group completed the follow-up. At 6 and 12 months after onset， SFGS and FDIp scores in both groups improved significantly while FDIs， SAS and SDS scores decreased （P＜0.05）. Comparison of scores between groups 6 months and 12 months after onset showed no statistically significant differences （P＞0.05）. During the trial， the incidence of adverse events was 13.3% （8/60） in the mild dysfunction termination group and 18.3% （11/60） in the complete recovery termination group， with no statistically significant difference （P＞0.05）. The number of treatment sessions， total social costs， and out-of-pocket expenses in the mild dysfunction termination group were significantly lower than those in the complete recovery termination group （P＜0.05）. The CER of the mild dysfunction termination group in SFGS， FDIp， FDIs， SAS， and SDS scores was lower than that of the complete recovery termination group. The ICER analysis showed that continuing treatment until full recovery incurred an additional cost of 573.30 CNY/point in SFGS improvement， whereas 1-point improvement in FDIp， FDIs， SAS， and SDS required 21,355.25 CNY， 1779.60 CNY， 3713.96 CNY， and 2755.52 CNY， respectively. ConclusionFor acupuncture in treating acute peripheral facial palsy， terminating treatment when mild dysfunction is achieved yields long-term efficacy comparable to that of continuing treatment until complete recovery， while significantly reducing medical costs and socioeconomic burden.

ObjectiveTo explore the efficacy and mechanism of Euphorbia humifusa on acute kidney injury （AKI） based on network pharmacology， molecular docking and experimental verification. MethodsThe active components and targets of E. humifusa were retrieved from TCMSP and SwissTargetPrediction database， and the AKI targets were screened by GeneCards and Online Mendelian Inheritance in Man（OMIM） databases. The drug targets and disease targets were intersected to construct a protein-protein interaction network， and the intersection targets were subjected to gene ontology （GO） and Kyoto encyclopedia of genes and genomes （KEGG） enrichment analysis. Discover Studio software was used to verify the molecular docking of key components and core targets. Gentamicin （GM） was used to induce AKI rat model. Control group， model group， verapamil （16 mg·kg^-1） group， E. humifusa extract （18， 54， 162 mg·kg^-1·d^-1） group and E. humifusa 70% ethanol extract （423 mg·kg^-1） group were continuously administered for 14 days. Urine volume was detected 24 h after modeling and administration. Serum creatinine （SCr）， Blood urea nitrogen （BUN）， 24-hour urine protein （24 hUTP） and uric acid （UA） content； the contents of malondialdehyde （MDA）， glutathione （GSH）， superoxide dismutase （SOD）， carbon monoxide synthase （NOS） and lactate dehydrogenase （LDH） in kidney were measured. The levels of interleukin （IL）-6 and tumor necrosis factor （TNF）-α in serum were detected by enzyme linked immunosorbent assay（ELISA） kit. The pathological changes of renal tissue were detected by hematoxylin-eosin （HE） and Masson staining. Western blot was used to detect the expression of PI3K/protein kinase B（Akt）/NF-κB signaling pathway-related proteins. ResultsIn this study， 13 active components such as kaempferol， luteolin， apigenin， gallic acid and quercetin were screened and identified from E. humifusa. Through bioinformatics analysis， these components and AKI have a total of 289 targets， of which 62 are core targets， including Akt1， TNF， tumor protein p53（TP53） and IL-1β. These targets are mainly involved in the regulation of biological processes such as NF-κB signaling pathway， HIF-1 signaling pathway， TNF signaling pathway， PI3K/Akt signaling pathway and mitogen-activated protein kinase（MAPK） signaling pathway. In animal experiments， we successfully constructed a GM-induced AKI model in rats. Compared with the model group， E. humifusa extract could significantly reduce the levels of 24 hUTP， BUN and SCr in rats （P<0.01）， indicating its improvement effect on renal function. In addition， the extract of E. humifusa also significantly reduced LDH activity and MDA content in rat kidney tissue （P<0.05， P<0.01）， and significantly increased SOD， NOS activity and GSH content （P<0.05）， indicating that the extract of E. humifusa has the potential of anti-oxidation and protection of renal function. Further analysis of inflammatory factors showed that the levels of IL-6 and TNF-α in serum of rats treated with E. humifusa extract were significantly decreased （P<0.01）， indicating that E. humifusa extract had anti-inflammatory effects. In addition， the extract of E. humifusa can also regulate the protein expression of PI3K/Akt/NF-κB signaling pathway， which further confirmed its mechanism of reducing GM-induced AKI. ConclusionThe extract of E. humifusa has a significant therapeutic effect on acute kidney injury through its multi-component and multi-target mechanism. Its effect is reflected in improving renal function， anti-oxidation， anti-inflammation and regulating immune response. These findings provide a scientific basis for the application of E. humifusa in the treatment of acute kidney injury， and point out the direction for future drug development and clinical research.

ObjectiveTo analyze the incidence and spatiotemporal distribution of typhus fever in Baoshan City, Yunnan Province from 2005 to 2023, to identify high-risk populations and regions, so as to provide a scientific basis for optimizing the allocation of local prevention and control resources and developing targeted intervention measures. MethodsData of typhus fever cases in Baoshan City from 2005 to 2023 were obtained from the Infectious Disease Information Management System of the Chinese Center for Disease Control and Prevention. Descriptive epidemiological methods were used to analyze the temporal, spatial and demographic distribution of typhus fever cases. Spatial clustering was assessed using spatial dynamic window scan statistics (circular and elliptical windows), flexible spatial scan statistics, and local spatial autocorrelation methods (including local Moran’s I, local Geary’s C, and Getis-Ord G_i*). Retrospective spatiotemporal scan statistics were employed to detect spatiotemporal clusters. ResultsA total of 1 099 typhus fever cases were reported in Baoshan City from 2005 to 2023. The incidence rate peaked at 6.31/ 100 000 in 2007, followed by a decline until reaching its lowest level at 0.21/100 000 in 2015 , and subsequently rebounded during 2016‒2023. The highest proportion of cases was among children under 10 years of age (31.12%), and the top three occupations of cases were farmers, students, and children, accounting for 88.62% of all cases. Cases occurred predominantly between June and September each year. The incidence was relatively high in Jiucheng Town (62.58/100 000), Yaoguan Town (57.15/100 000), and Dianyang Town (46.81/100 000) of Shidian County. Spatial clustering analyses indicated that high-risk areas were mainly located in the southern part of Baoshan City, showing a south-to-north trend. Spatiotemporal scan analyses identified five clusters, with the most likely cluster centered around Yaoguan Town, covering ten towns (subdistricts) during the period 2007‒2010. ConclusionThe incidence of typhus fever in Baoshan City exhibits a clear seasonal and spatial clustering pattern, with peak incidence occurring in summer and autumn. Spatially, cases are primarily distributed in the southern part of Baoshan City, and high-risk clusters exhibit a south-to-north trend. Farmers, students, and children are the high-risk groups.

Objective To assess the survival vulnerability of Oncomelania hupensis in Jiangxi Province under future climate scenarios, and to identify low-vulnerability areas for its survival in this province. Methods Village-level O. hupensis snail survey and O. hupensis snail control with chemical treatments in Jiangxi Province from 2016 to 2024 were captured from the Parasitic Disease Prevention and Control Information Management System of China Disease Prevention and Control Information System. Climatic data were primarily sourced from the Resource and Environmental Science Data Platform, Chinese Academy of Sciences (http://www.resdc.cn/), including annual average temperature, annual average precipitation, annual accumulated temperature above 10 °C, annual accumulated temperature above 0 °C, annual maximum temperature, annual minimum temperature, and annual average relative humidity, and nineteen bioclimatic variables were downloaded from the WorldClim website (https://www.worldclim.org/), including mean diurnal range, isothermality, temperature seasonality, and so on. Elevation and normalized difference vegetation index were catprued from the Resource and Environmental Science Data Platform, Chinese Academy of Sciences (http://www.resdc.cn/), and distance to rivers was downloaded from the WorldPop website (http://www.worldpop.org), and land use and land cover (LULC) data were downloaded from the Big Earth Data Center, Chinese Academy of Sciences (https://data.casearth.cn/), and nature reserve data were obtained from the China Nature Reserve Specimen Resource Sharing Platform (http://www.papc.cn/). Three Shared Socioeconomic Pathways (SSPs) from the Beijing Climate Center-Climate System Model version 2-Medium Resolution (BCC-CSM2-MR) global climate model were employed as future climate scenarios, including SSP126, SSP245, SSP585, and the biomod2 ensemble model in R package was used to simulate suitable habitats for O. hupensis snails in Jiangxi Province in 2050 and 2070 under these scenarios. A snail survival vulnerability index was constructed based on the area of suitable snail habitats, area covered by snail control through chemical treatment, area covered by nature reserves, and changes in snail habitat fragmentation, and a map of snail survival vulnerability distribution was plotted. Results The real area of snail habitats ranged from 78 486.76 to 85 309.47 hm², and the area of snail control with chemical treatment ranged from 10 138.98 to 13 240.16 hm² in Jiangxi Province from 2016 to 2024. There were 429 to 531 villages detected with snails during the nine-year period, and the number of actually snail-infested villages ranged from 645 to 686. A total of 818 snail-present points and 1 996 snail-absent points were obtained from snail survey records. The best performance of the biomod2 ensemble model was achieved if a weighted mean approach was used as the ensemble strategy, with a true skill statistic value of 0.799 and an area under the receiver operating characteristic curve of 0.957, and modeling identified annual average relative humidity and annual average precipitation as two most influencing climatic variables for snail distribution. Relative to the current areas of suitable snail habitats under present climate conditions, the area of suitable snail habitats was projected to expand by 24.49% to 46.28% in Jiangxi Province under future climate scenarios, and the proportion of nature reserves areas in the areas of suitable snail habitats was projected to decrease slightly from the current 2.77% to approximately 2.52%, while the proportion of areas of snail control through chemical treatment in areas of suitable snail habitats varied from 0.64% to 19.57%, and the percentage of changes in snail habitat fragmentation ranged from 3.86% to 12.23%. Based on these four indicators, the snail survival vulnerability index was estimated to range from –1.96 to 0.62 in Jiangxi Province. The arithmetic mean of the snail survival vulnerability index differed under three SSP scenarios (SSP126, SSP245 and SSP585), with the highest mean value (–0.69) in 2070 under SSP126, and the lowest mean value (–0.78) in 2070 under SSP585. Conclusions The snail survival vulnerability index ranges from –1.96 to 0.62 in Jiangxi Province under future climate scenarios, and the suitable habitats for O. hupensis snails appear an overall tendency towards expansion. Low-vulnerability snail habitats are mainly distributed along the shores of Poyang Lake and the Yangtze River in Jiangxi Province, partially overlapping with nature reserves. Intensified surveillance of O. hupensis snails is recommended in these areas in the future.