BackgroundPatients with depressive episode and schizophrenia have a high risk of suicide. The sleep electroencephalogram power spectral density characteristics of patients with depressive episode accompanied by suicidal behavior and those with schizophrenia may be different， but there is currently a lack of direct comparative studies on these two groups of patients. ObjectiveTo compare the sleep electroencephalogram power spectral density between depressive episode and schizophrenic patients with suicidal behavior， in order to provide references for exploring predictive indicators of suicidal behavior. MethodsFrom June 2018 to December 2020， 20 patients with depressive episode and 20 patients with schizophrenia who had committed suicide within the past month and were treated at the outpatient department of Shenzhen Kangning Hospital were selected. All of them met the diagnostic criteria for depressive episode or schizophrenia as defined in the International Classification of Diseases， tenth edition （ICD-10）. Using a random sampling method， 20 volunteers with matching gender and age to the patient groups were selected from the Cuiping community in Shenzhen as the control group. The subjective sleep of the patients was evaluated using the Insomnia Severity Index （ISI）， the Dysfunctional Belief and Attitude about Sleep （DBAS）， the Disturbing Dreams and Nightmare Severity Index （DDNSI）， and the Epworth Somnolence Scale （ESS）. The objective sleep of the patients was assessed using polysomnography. The sleep electroencephalogram was filtered and the power spectral density of the brain wave was analyzed and processed for all the subjects. The subjective and objective sleep conditions of the two patient groups were compared， and the sleep electroencephalogram power spectral density of the patient groups and the control group were also compared. ResultsA comparison of subjective and objective sleep conditions between patients with depressive episode accompanied by suicidal behavior and patients with schizophrenia accompanied by suicidal behavior showed no statistically significant differences （P>0.05）. Comparisons of sleep electroencephalogram power spectral density in the W stage （average power of α wave， total power of δ wave， average power of δ wave， average power of θ wave）， N1 stage （average power of β wave， total power of α wave， total power of δ wave）， N2 stage （total power of α wave， average power of α wave， total power of δ wave， average power of δ wave）， N3 stage （average power of α wave， average power of δ wave）， and R stage （total power of α wave， average power of α wave， total power of δ wave， average power of δ wave） between patients with depressive episode accompanied by suicidal behavior， patients with schizophrenia accompanied by suicidal behavior， and the control group showed statistically significant differences （P<0.05 or 0.01）. The total power of δ wave in the W stage and the average power of β wave and δ wave in the N1 stage were higher in two patient groups were higher than those of the control group. The total power of α wave and the average power of α wave in the N2 stage were lower than those of the control group， while the average power of δ wave was higher than that of the control group. The average power of α wave in the N3 stage of both patient groups were lower than that of the control group， while the average power of δ wave was higher than that of the control group. The total power and average power of α wave in the R stage were lower than those of the control group， while the total power and average power of δ wave were higher than those of the control group. All the differences were statistically significant. Patients with depressive episode accompanied by suicidal behavior had higher average powers of α wave， δ wave， and θ wave in the W stage compared with the control group， while the total power of α wave in the N1 stage was lower in the former group. All these differences were statistically significant （P<0.05）. ConclusionThe depressive episode patients accompanied by suicidal behavior have highly overlapping sleep electroencephalogram abnormal patterns with those of schizophrenia patients， mainly manifested as a general decrease in α wave power （N2， N3， R stage） and a general increase in δ wave power （W， N1， N2， N3， R stage） as well as β wave power in N1 stage. At the same time， patients with depressive episode accompanied by suicidal behavior also show specific changes， including an increase in the average power of α and θ waves during the wakefulness period （W stage）， and a decrease in the total power of α wave in N1 stage. ［Funded by Guangdong Province High-level Clinical Key Specialty （with supporting funds from Shenzhen City） （number， SZGSP013）； Shenzhen Key Medical Discipline （number， SZXK041）； Shenzhen Clinical Medicine Research Center Project （number， 20210617155253001）］

ObjectiveTo investigate the effect of mitochondrial calcium uniporter （MCU） on the cytoskeleton of pancreatic ductal epithelial cells in a mouse model of acute pancreatitis （AP） induced by caerulein （CAE）， to analyze the role of MCU in the development of AP， and to provide a theoretical basis for clinical treatment. MethodsIn the in vivo experiment， wild-type male C57BL6/J mice， aged 4 weeks， were randomly divided into control group and AP group， with 6 mice in each group. The mice in the AP group were given intraperitoneal injection of CAE to establish a model of AP， and those in the control group were given intraperitoneal injection of an equal volume of normal saline. Serum and pancreatic tissue samples were collected after 24 hours of modeling. HE staining was used to observe pancreatic histopathological changes； Western Blot was used to measure the expression levels of MCU， glutathione peroxidase 4 （GPX4）， and acyl-CoA synthetase long chain family member 4 （ASCL4）； kits were used to measure the serum level of amylase. In the in vitro experiment， the human pancreatic ductal epithelial cell line HPDE6-C7 was co-cultured with CAE for 24 hours to establish an in vitro AP model， and the cells were divided into control group， CAE group， RR （an MCU activity inhibitor） group， CAE+RR group， Fer-1 （an ferroptosis inhibitor） group， CAE+Fer-1 group， Erastin （an ferroptosis inducer） group， and CAE+Erastin group. CCK-8 assay was used to observe the influence of different agents on cell viability； Western Blot was used to measure the expression levels of MCU， GPX4， and ASCL4； immunofluorescence assay was used to measure reactive oxygen species （ROS）， actin cytoskeleton， and monolayer permeability； kits were used to measure the concentrations of malondialdehyde （MDA）， glutathione （GSH）， Fe²⁺， and total iron. A one-way analysis of variance was used for comparison of continuous data between multiple groups， and the least significant difference t-test was used for comparison between two groups. ResultsIn the in vivo experiment， compared with the control group， the AP group had significant increases in pancreatic histopathological score， the serum level of amylase， and the expression levels of MCU and ASCL4， as well as a significant reduction in the expression of GPX4 （all P<0.05）. In the in vitro experiment， compared with the control group， the CAE group had significant increases in the expression levels of MCU and ASCL4， a significant reduction in the expression of GPX4， and significant increases in the concentrations of Fe²⁺， total iron， and MDA， the green fluorescence intensity of ROS， and monolayer permeability， as well as a significant reduction in the concentration of GSH （all P<0.05）， with the presence of actin cytoskeleton disruption. Compared with the CAE group， the CAE+RR group had a significant increase in the expression level of GPX4， a significant reduction in the expression level of ASCL4， and significant reductions in the concentrations of Fe²⁺， total iron， and MDA， the green fluorescence intensity of ROS， and monolayer permeability and a significant increase in the concentration of GSH （all P<0.05）， with alleviation of actin cytoskeleton disruption. Compared with the CAE group， the CAE+Fer-1 group had significant reductions in the concentrations of Fe²⁺， total iron， and MDA， the green fluorescence intensity of ROS， and monolayer permeability and a significant increase in the concentration of GSH （all P<0.05）， with alleviation of actin cytoskeleton disruption. Compared with the CAE group， the CAE+Erastin group had significant increases in the concentrations of Fe²⁺， total iron， and MDA， the green fluorescence intensity of ROS， and monolayer permeability and a significant reduction in the concentration of GSH （all P<0.05）， with aggravation of actin cytoskeleton disruption. ConclusionDuring the onset of AP， MCU mediates oxidative stress-induced ferroptosis and leads to the disruption of the pancreatic ductal epithelial barrier， which may be one of the possible pathogeneses of AP.

OBJECTIVE To investigate the influence of CYP2C19 gene polymorphism on platelet function and inflammatory cytokines in elderly patients with ischemic stroke， and to analyze potential factors associated with poor prognosis. METHODS A retrospective study was conducted on elderly patients with ischemic stroke admitted to our hospital from June 2024 to June 2025， wh o underwent CYP2C19 genotype testing and received antiplatelet therapy with clopidogrel. The levels of platelet function indicators and inflammatory cytokines before and after treatment were compared among patients with different metabolic phenotypes. Based on the prognosis at 6 months post-treatment， patients were divided into poor prognosis group and good prognosis group. Univariate analysis was performed on general data， metabolic phenotype， the levels of platelet function indicators and inflammatory cytokines. Variables with P ＜0.05 and the levels of inflammatory cytokines before treatment were included in a multivariate Logistic regression analysis to identify independent risk factors for poor prognosis. Multiple linear regression was used to further analyze the relationship between metabolic phenotypes and inflammatory cytokines. RESULTS A total of 448 elderly patients with ischemic stroke were included； among them， 162 cases were normal metabolic phenotype， 218 were intermediate metabolic phenotype， and 68 were poor metabolic phenotype. No rapid or ultrarapid metabolic phenotypes were observed. After treatment， platelet aggregation rate， the levels of P-selectin and platelet activated complex-1 （PAC-1）， high-sensitivity C-reactive Protein （hs-CRP）， interleukin-1β （IL-1β）， IL-6 and tumor necrosis factor-α （TNF-α） in the normal metabolic phenotype group， intermediate metabolic phenotype group， and poor metabolic phenotype group （except for platelet aggregation rate， and the levels of P-selectin and PAC-1 in the poor metabolic phenotype group） were significantly lower than those before treatment in the same group. Moreover， the above indicators in the normal metabolic phenotype group were significantly lower than those in the intermediate and poor metabolic phenotype groups at the corresponding time， and the levels of platelet function indicators in the intermediate metabolic phenotype group were significantly lower than those in the poor metabol ic phenotype group at the corresponding time （ P ＜0.05）. Univariate and multivariate Logistic regression analyses showed that combined with hypertension， combined with diabetes mellitus， and intermediate or poor metabolic genotypes were independent risk factors for poor prognosis in elderly patients with ischemic stroke （ P ＜0.05）. Multiple linear regression analysis showed that serum levels of hs-CRP， IL-1β， IL-6 and TNF-α before treatment were significantly higher in patients with intermediate and poor metabolic genotypes compared to those with normal metabolic genotype （ P ＜0.05）， with a greater magnitude of increase in inflammatory cytokines observed in the patients with poor metabolic genotype. CONCLUSIONS The elderly ischemic stroke patients with CYP2C19 intermediate and poor metabolic genotypes have poor inhibition effect on platelet and higher levels of inflammatory cytokines than normal metabolic genotype； CYP2C19 gene polymorphism， and in combination with hypertension and diabetes， can be used as independent predictors of poor prognosis.

ObjectiveTo describe the disposal process of the first imported case of Zika virus disease in Shanghai, and to provide a reference for the prevention and control of imported infectious diseases in the future. MethodsA retrospective review was conducted of the process by which epidemiological investigation, etiological testing, and case management were performed by the Shanghai municipal-and district-level Centers for Disease Control and Prevention (CDC) after one imported case of Zika virus disease was identified by Shanghai Customs. ResultsOn April 8th, 2025, customs authority at a certain airport in Shanghai identified a febrile inbound case (case A, female, 40 years old). An antecubital venous blood specimen was obtained and analyzed by nucleic acid testing for various vector-borne pathogens, by which Zika virus nucleic acid was found to be positive. On 9 April, the district CDC collected antecubital venous blood specimens again from Case A and her three accompanying travelers (B, C, and D), and nucleic acid testing was conducted for multiple mosquito-borne infectious diseases. Case A again tested positive for Zika virus nucleic acid, while Traveler D (male, aged 6 years) tested positive for dengue virus nucleic acid. The other two travelers tested negative. Case A and Traveler D were subsequently transferred to a designated district hospital for isolation and treatment. After discharge, both cases left China and returned to their overseas residence. ConclusionCase A was the first imported case of Zika virus disease in Shanghai. For cases with a history of living in endemic areas of multiple infectious diseases such as dengue virus, Zika virus and chikungunya virus, it is necessary to carry out rapid testing of multiple pathogens simultaneously to prevent missed diagnosis or misdiagnosis.

Multidrug-resistant (MDR) bacterial infections have emerged as a serious challenge of global public health crisis. The overuse and misuse of conventional antibiotics have dramatically accelerated the emergence, evolution and worldwide spread of drug-resistant bacterial strains, necessitating urgent exploration of novel antibacterial strategies. Bacteriophages serve as natural bacterial predators offering distinct advantages including high host specificity, autonomous self-replication capabilities and cost-effective large-scale production. However, wild-type phages present significant clinical limitations due to their narrow host ranges, susceptibility to rapid immune clearance and poor penetration of bacterial biofilms, which severely restrict their therapeutic applications. The convergence of synthetic biology, nanotechnology and advanced gene editing technologies has accelerated the development of engineered bacteriophage platforms, providing programmable, scalable and clinically translatable pathways to overcome these inherent biological constraints. Here, we systematically delineate four fundamental strategies for engineered bacteriophage development. Chemical modification utilizes reactive functional groups such as amino, carboxyl and thiol moieties on capsid proteins through esterification, amidation or click chemistry reactions to achieve precise drug conjugation and surface functionalization. In vivo editing encompasses ultraviolet or chemical mutagenesis for random mutation induction, homologous recombination for targeted genetic alterations, recombineering methodologies including electroporation-mediated bacteriophage recombination engineering, and CRISPR-Cas systems for precise genome editing to enable exact genetic reconstruction and host range reprogramming. In vitro synthesis leverages genome engineering platforms where intact phage genomes are transferred into yeast or host bacteria to facilitate highly efficient homologous recombination, enabling large DNA fragment assembly and cross-gene host range expansion without bacterial toxicity constraints. Directed evolution combines artificial selection through mutation library screening with rational design approaches involving chimeric receptor binding protein construction or site-specific mutagenesis, effectively balancing the discovery of unknown adaptive pathways with targeted host specificity modification. Moreover, we comprehensively discuss therapeutic applications across diverse clinical scenarios. Engineered bacteriophage effectively disrupt bacterial biofilms through sophisticated functionalized delivery platforms including nanozyme-conjugated phages, phage-liposome nanoconjugates and bio-responsive hydrogels, demonstrating significantly enhanced bactericidal efficiency compared to unmodified free phages. These bioengineered vectors attenuate bacterial virulence and resensitize pathogens to antibiotics by delivering CRISPR-Cas systems or base editors to disrupt critical virulence factors such as pili, capsule synthesis machineries and quorum sensing systems, or by inactivating antibiotic resistance determinants including beta-lactamase genes. As an intelligent nanomedicine delivery platform, engineered bacteriophage enable precise pathogen elimination an through photocatalytic reactive oxygen species generation, immunomodulatory interventions, or controlled release of antibacterial drugs. Furthermore, oral administration of engineered bacteriophage facilitates microbiota modulation, which selectively eliminate intestinal pathogens while preserve beneficial commensal microbiota, thereby restoring microbial community balance and preventing complications associated with dysbiosis. Finally, we critically analyze persistent challenges including host strain matching complexity, evolution of bacterial resistance mechanisms, pharmacokinetic optimization requirements, optimal administration route selection, large-scale production quality control standards and clinical dosing determination protocols. Through multidisciplinary integration of synthetic biology, infectious disease medicine and immunology, future translational medicine studies of bacteriophage should establish comprehensive technical platforms encompassing rapid phage screening, intelligent rational design, rigorous in vivo evaluation and standardized clinical validation processes, ultimately advancing engineered bacteriophage from laboratory innovations to clinically approved therapeutics for effectively combating MDR bacterial infections.

Objective To gain insights into the biological characteristics of different populations of Culex tritaeniorhynchus within and around the Yellow Sea wetland from the perspective of the occurrence of resistance, we investigated the levels of resistance to deltamethrin and kdr gene mutation in the wetland and its peripheral areas. Methods Specimens were collected from Cx. tritaeniorhynchus populations at two monitoring sites in the Rare Bird National Nature Reserve and Tiaozi Ni Wetland Scenic Area, and also from two populations in Yancheng City and the Liuhe District of Nanjing, and the resistance of these mosquitoes to deltamethrin was determined using the CDC biotest bottle method. For each concentration of deltamethrin assessed, a random subset of exposed specimens was selected for amplification of the kdr gene fragment, followed by Sanger sequencing to identify and analyze resistance-associated mutations. Results The LC50 levels of deltamethrin among mosquitoes from the four populations in Luhe, Yancheng, the Rare Bird National Nature Reserve and the Tiaozi Ni Wetland Scenic Area were 2.048 5, 7.798 2, 3.473 3, and 17.695 5 mg/mL, respectively, with corresponding concentrations of deltamethrin ranging from 0.005 to 5.000,0.050 to 50.000,0.050 to 25.000 and 0.050 to 50.000 mg/mL, respectively. Furthermore, the ranges of the KT50 values were 11.76-107.43, 67.05-216.30,29.77-107.43 and 28.40-329.51 min; the 1-h knockdown rates were 34.58%-99.15%, 9.52%-43.80%, 55.09%-73.01%, and 10.09%-68.07%; and the 24-h mortality rates were 12.15%-67.52%,9.52%-79.56%,13.17%-82.21%, and 11.01%-78.99%, respectively. With respect to kdr gene mutation, we assayed a total of 63,70,59, and 57 mosquitoes for the four populations, for which we detected L1014F mutation frequencies of 14.29%, 35.00%, 20.34%, and 31.58%, respectively, with a majority of these mutations being heterozygous for resistance. In addition, five adult mosquitoes were identified has having synonymous mutations at site 1011[i. e. , AAT(asparagine)mutation to AAC(asparagine)]. Conclusions Our findings revealed the clear resistance of Cx. tritaeniorhynchus to deltamethrin in the Yancheng region of the Yellow Sea wetland, and the resistance phenotype and kdr frequency of Cx. tritaeniorhynchus in the wetland environment were comparable to those of Cx. tritaeniorhynchus in the wetland environment, thereby indicating that the resistance of different populations of Cx. tritaeniorhynchus was homogeneous under the pressure of different insecticide selection within and around the wetland. However, the underlying mechanisms need to be further studied.

Background and aims The mechanism of cholestatic liver injury (CLI) is unclear, and effective therapies are lacking. While peroxisome proliferator-activated receptor alpha (PPARα) agonists show potential hepatoprotective effect and pyroptosis is implicated in hepatocellular damage, how PPARα activation mitigates lithocholic acid (LCA)-induced pyroptosis remains unknown. Methods The hepatoprotective effect of PPARα agonists was evaluated in a mouse model of intrahepatic cholestasis induced by LCA. Liver injury was assessed via serum biochemistry, hematoxylin and eosin and TUNEL staining, and electron microscopy. Pyroptosis pathways were analyzed using real-time quantitative polymerase chain reaction, Western blot, and co-immunoprecipitation. Results Combined morphological, histopathological, and biochemical analyses confirmed that PPARα activation protects against CLI. Compared with LCA treatment alone, PPARα activation significantly attenuated the elevation of serum lactate dehydrogenase (LDH), the increased TUNEL-positive cells, and the formation of hepatocyte membrane pores. Mechanistically, PPARα activation suppressed both NOD-like receptor protein 3 (NLRP3) inflammasome-mediated pyroptosis and apoptosis protease-activating factor-1 (APAF-1)/CASPASE-3/GSDME-mediated pyroptosis. Furthermore, PPARα agonist pretreatment inhibited activation of the nuclear factor-kappa B (NF-κB) and forkhead box O1 (FOXO1) signaling pathways. Conclusions PPARα protects against LCA-induced CLI by inhibiting both NLRP3 inflammasome-mediated pyroptosis associated with NF-κB and APAF-1/CASPASE-3/GSDME-mediated pyroptosis associated with the FOXO1 signaling pathway.

ObjectiveThis study aimed to investigate the effects of 4-week high-intensity interval training (HIIT) on synaptic plasticity in the prefrontal cortex (PFC) of rats exposed to chronic unpredictable mild stress (CUMS), and to explore its potential mechanisms. MethodsA total of 48 male Sprague-Dawley rats were randomly divided into 4 groups: control (C), model (M), control plus HIIT (HC), and model plus HIIT (HM). Rats in groups M and HM underwent 8 weeks of CUMS to establish depression-like behaviors, while groups HC and HM received HIIT intervention beginning from the 5th week for 4 consecutive weeks. The HIIT protocol consisted of repeated intervals of 3 min at high speed (85%-90% maximal training speed, S_max) alternated with one minute at low speed (50%-55% S_max), with 3 to 5 sets per session, conducted 5 d per week. Behavioral assessments and tail-vein blood lactate levels were measured at the end of the 4th and 8th weeks. After the intervention, rat PFC tissues were collected for Golgi staining to analyze synaptic morphology. Enzyme-linked immunosorbent assays (ELISA) were employed to detect brain-derived neurotrophic factor (BDNF), monocarboxylate transporter 1 (MCT1), lactate, and glutamate levels in the PFC, as well as serotonin (5-HT) levels in serum. Additionally, Western blot analysis was conducted to quantify the expression of synaptic plasticity-related proteins, including c-Fos, activity-regulated cytoskeleton-associated protein (Arc), and N-methyl-D-aspartate receptor 1 (NMDAR1). ResultsCompared to the control group (C), the CUMS-exposed rats (group M) exhibited significant reductions in sucrose preference rates, number of grid crossings, frequency of upright postures, and entries into and duration spent in open arms of the elevated plus maze, indicating marked depressive-like behaviors. Additionally, the group M showed significantly reduced dendritic spine density in the PFC, along with elevated levels of c-Fos, Arc, NMDAR1 protein expression, and increased concentrations of lactate and glutamate. Conversely, BDNF and MCT1 contents in the PFC and 5-HT levels in serum were significantly decreased. Following HIIT intervention, rats in the group HM displayed considerable improvement in behavioral indicators compared with the group M, accompanied by significant elevations in PFC MCT1 and lactate concentrations. Furthermore, HIIT notably normalized the expression levels of c-Fos, Arc, NMDAR1, as well as glutamate and BDNF contents in the PFC. Synaptic spine density also exhibited significant recovery. ConclusionFour weeks of HIIT intervention may alleviate depressive-like behaviors in CUMS rats by increasing lactate levels and reducing glutamate concentration in the PFC, thereby downregulating the overexpression of NMDAR, attenuating excitotoxicity, and enhancing synaptic plasticity.

ObjectiveTo investigate the awareness and clinical practice of guidelines for the prevention and treatment of chronic hepatitis B （2022 edition） among clinicians. MethodsFrom July 19 to December 31， 2024， a self-designed electronic questionnaire was distributed via the WeChat mini program to collect related data from 1 588 clinicians nationwide， including their awareness and practice based on 18 questions regarding testing and referral， diagnosis and treatment， and follow-up. ResultsAmong all respondents， only 350 clinicians correctly understood all the updated key points of antiviral indications and treatment for special populations in the 2022 edition of guidelines for the prevention and treatment of chronic hepatitis B， with an overall awareness rate of 22.0%. Only 20% ‍—‍ ‍40% of the patients with positive HBV DNA and an age of >30 years receive antiviral therapy， while 80% ‍—‍ ‍100% of the patients with positive HBV DNA and a family history of hepatitis B cirrhosis or hepatocellular carcinoma receive antiviral therapy. The median follow-up rates at 1 year， 3 years， and 5 years were 67.5% 57.5% and 47.5%，respectively， showing a trend of gradual reduction， which might be associated with the influencing factors such as insufficient time for follow-up management by clinicians， insufficient awareness of the disease among patients， and poor adherence to follow-up. ConclusionThere is a gap between the awareness and practice of guidelines for the prevention and treatment of chronic hepatitis B （2022 edition） among clinicians. It is recommended to further strengthen training and focus on the whole process of “detection， diagnosis， treatment， and management” for patients with chronic hepatitis B in healthcare institutions， in order to promote the implementation of the guidelines.

Objective To investigate the suitable habitats of Phlebotomus chinensis in Gansu Province, so as provide insights into effective management of mountain-type zoonotic visceral leishmaniasis (MT-ZVL). Methods The geographical coordinates of locations where MT-ZVL cases were reported were retrieved in Gansu Province from 2015 to 2023, and data pertaining to 26 environmental variables were captured, including 19 climatic variables (annual mean temperature, mean diurnal range, isothermality, temperature seasonality, maximum temperature of the warmest month, minimum temperature of the coldest month, temperature annual range, mean temperature of the wettest quarter, mean temperature of the driest quarter, mean temperature of the warmest quarter, mean temperature of the coldest quarter, annual precipitation, precipitation of the wettest month, precipitation of the driest month, precipitation seasonality, precipitation of the wettest quarter, precipitation of the driest quarter, precipitation of the warmest quarter, and precipitation of the coldest quarter), five geographical variables (elevation, annual normalized difference vegetation index, vegetation type, landform type and land use type), and two population and economic variables (population distribution and gross domestic product). Twelve species distribution models were built using the biomod2 package in R project, including surface range envelope (SRE) model, generalized linear model (GLM), generalized additive model (GAM), multivariate adaptive regression splines (MARS) model, generalized boosted model (GBM), classification tree analysis (CTA) model, flexible discriminant analysis (FDA) model, maximum entropy (MaxEnt) model, optimized maximum entropy (MAXNET) model, artificial neural network (ANN) model, random forest (RF) model, and extreme gradient boosting (XGBOOST) model. The performance of 12 models was evaluated using the area under the receiver operating characteristic curve (AUC), true skill statistics (TSS), and Kappa coefficient, and single models with high performance was selected to build the optimal ensemble models. Factors affecting the survival of Ph. chinensis were identified based on climatic, geographical, population and economic variables. In addition, the suitable distribution areas of Ph. chinensis were predicted in Gansu Province under shared socioeconomic pathway 126 (SSP126), SSP370 and SSP585 scenarios based on climatic data during the period from 1991 to 2020, from 2041 to 2060 (2050s), and from 2081 to 2100 (2090s) . Results A total of 11 species distribution models were successfully built for prediction of potential distribution areas of Ph. chinensis in Gansu Province, and the RF model had the highest predictive accuracy (AUC = 0.998). The ensemble model built based on the RF model, XGBOOST model, GLM, and MARS model had an increased predictive accuracy (AUC = 0.999) relative to single models. Among the 26 environmental factors, precipitation of the wettest quarter (12.00%), maximum temperature of the warmest month (11.58%), and annual normalized difference vegetation index (11.29%) had the greatest contributions to suitable habitats distribution of Ph. sinensis. Under the climatic conditions from 1991 to 2020, the potential suitable habitat area for Ph. chinensis in Gansu Province was approximately 5.80 × 10⁴ km², of which the highly suitable area was 1.42 × 10⁴ km², and primarily concentrated in the southernmost region of Gansu Province. By the 2050s, the unsuitable and lowly suitable areas for Ph. chinensis in Gansu Province had decreased by varying degrees compared to that of 1991 to 2020 period, while the moderately and highly suitable areas exhibited expansion and migration. By the 2090s, under the SSP126 scenario, the suitable habitat area for Ph. chinensis increased significantly, and under the SSP585 scenario, the highly suitable areas transformed into extremely suitable areas, also showing substantial growth. Future global warming is conducive to the survival and reproduction of Ph. chinensis. From the 2050s to the 2090s, the highly suitable areas for Ph. chinensis in Gansu Province will be projected to expand northward. Under the SSP126 scenario, the suitable habitat area for Ph. chinensis in Gansu Province is expected to increase by 194.75% and 204.79% in the 2050s and 2090s, respectively, compared to that of the 1991 to 2020 period. Under the SSP370 scenario, the moderately and highly suitable areas will be projected to increase by 164.40% and 209.03% in the 2050s and 2090s, respectively, while under the SSP585 scenario, they are expected to increase by 195.98% and 211.66%, respectively. Conclusions The distribution of potential suitable habitats of Ph. sinensis gradually shifts with climatic changes. Intensified surveillance and management of Ph. sinensis is recommended in central and eastern parts of Gansu Province to support early warning of MT-ZVL.