1.Comparison of sleep EEG power spectral density between depressive episode patients and schizophrenia patients with suicidal behavior
Jingwen LIU ; Yunfei ZHOU ; Jingchu HU ; Jiaoyan ZHOU ; Junwei YANG ; Jie LIANG ; Hong XU ; Yu CANG ; Shimeng MA
Sichuan Mental Health 2026;39(1):50-57
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)]
2.Engineered Bacteriophages for The Treatment of Multidrug-resistant Bacterial Infections
Yu-Ying CHEN ; Chun-Mei HUANG ; Jin-Zhi PAN ; De-Liang LIU ; Yang ZHOU ; Gui-Qin DAI ; Peng-Fei ZHAO ; Hong-Zhou LU ; Ming-Bin ZHENG
Progress in Biochemistry and Biophysics 2026;53(6):1581-1596
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.
3.Differences in deltamethrin resistance and kdr gene mutation in Culex tritaeniorhynchus population in and outside the Yellow Sea wetland
Xiao-er ZHANG ; Zhi-ming WU ; Ye TIAN ; Qian CUI ; Yu-qian JI ; Huan WANG ; Shu-juan YANG ; Yi-chao ZHAO ; Yu WANG ; Hua-yu YIN ; Yu DING ; Guo-jin YAN ; Min-sen ZHAO ; Shou-gang ZHANG ; Bing-dong SONG ; Hong-na CHEN ; Jian GAO ; Wei-fang YANG ; Yu-fu ZHANG ; Hui LIU ; Hong-liang CHU
Acta Parasitologica et Medica Entomologica Sinica 2026;33(2):101-107
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.
4.Construction of a Dual-dimensional Outcome System for Musculoskeletal Pain:A Common Core Outcome Set and A Specific Outcome Set for Traditional Chinese Medicine External Treatment
Yuanwen LIANG ; Hong CHEN ; Yanan SUN ; Changhe YU
Journal of Traditional Chinese Medicine 2026;67(14):1503-1508
ObjectiveTo develop a core outcome set for musculoskeletal pain (MKP-COS) and a specific outcome set for traditional Chinese medicine external treatment (MKP-SOS-ETCM), providing a standardized basis for clinical research in musculoskeletal pain and serving as a reference for core outcome set (COS) development in related fields. MethodsExisting COS studies on musculoskeletal pain (MKP) were systematically searched, and core outcome indicators were extracted to form an initial item pool. The Delphi method was used, in which experts were invited to reach consensus on MKP-COS and MKP-SOS-ETCM. The inclusion criteria for consensus are set as an importance score of 7-9 points with ≥75% of experts reaching agreement. ResultsA total of 47 COS studies were included, and after discussion and consolidation, a preliminary item pool containing 34 outcome indicators was established. After three rounds of Delphi expert consensus, MKP-COS was identified, including five outcomes including pain, motor function, joint range of motion, activities of daily living, and adverse events/ adverse reactions; MKP-SOS-ETCM included six outcomes involving stiffness (including morning stiffness), joint tenderness count, work status/ productivity, overall physician assessment, structural damage, and muscle strength. ConclusionThis study has established MKP-COS and MKP-SOS-ETCM, forming a dual-dimensional outcome assessment system integrating both disease commonality and intervention specificity.
5.Construction of a Dual-dimensional Outcome System for Musculoskeletal Pain:A Common Core Outcome Set and A Specific Outcome Set for Traditional Chinese Medicine External Treatment
Yuanwen LIANG ; Hong CHEN ; Yanan SUN ; Changhe YU
Journal of Traditional Chinese Medicine 2026;67(14):1503-1508
ObjectiveTo develop a core outcome set for musculoskeletal pain (MKP-COS) and a specific outcome set for traditional Chinese medicine external treatment (MKP-SOS-ETCM), providing a standardized basis for clinical research in musculoskeletal pain and serving as a reference for core outcome set (COS) development in related fields. MethodsExisting COS studies on musculoskeletal pain (MKP) were systematically searched, and core outcome indicators were extracted to form an initial item pool. The Delphi method was used, in which experts were invited to reach consensus on MKP-COS and MKP-SOS-ETCM. The inclusion criteria for consensus are set as an importance score of 7-9 points with ≥75% of experts reaching agreement. ResultsA total of 47 COS studies were included, and after discussion and consolidation, a preliminary item pool containing 34 outcome indicators was established. After three rounds of Delphi expert consensus, MKP-COS was identified, including five outcomes including pain, motor function, joint range of motion, activities of daily living, and adverse events/ adverse reactions; MKP-SOS-ETCM included six outcomes involving stiffness (including morning stiffness), joint tenderness count, work status/ productivity, overall physician assessment, structural damage, and muscle strength. ConclusionThis study has established MKP-COS and MKP-SOS-ETCM, forming a dual-dimensional outcome assessment system integrating both disease commonality and intervention specificity.
6.Current Status and Prospects of Artificial Intelligence Technologyin Minimally Invasive Gastric Cancer Surgery
Tao ZHANG ; Boer SU ; Guanxing LIANG ; Shiman DAI ; Jiawei CHEN ; Zhengjie LIU ; Cheng PENG ; Rong LIU ; Qinglan LIN ; Yidan WU ; Yuhui WU ; Jiaming WEN ; Hong WANG ; Hao CHEN ; Jiang YU
Medical Journal of Peking Union Medical College Hospital 2026;17(4):933-942
Gastric cancer remains a highly prevalent malignancy worldwide, with surgical resection currently constituting the cornerstone of treatment aimed at improving long-term patient survival. Owing to their notable advantages, including reduced surgical trauma and accelerated postoperative recovery, minimally invasive procedures are progressively supplanting conventional open surgery and have become the mainstream approach in gastric cancer management. Concurrently, the rapid advancement of artificial intelligence (AI) technologies has enabled real-time intraoperative monitoring of surgical scenes, thereby furnishing novel technical support for adjunctive decision-making, surgical navigation, and skill assessment during gastrectomy. This article provides a systematic review of the current status of AI applications in minimally invasive gastric cancer surgery, with a particular focus on research progress pertaining to instrument recognition, surgical phase identification, delineation of normal anatomical structures, detection of metastatic foci, and early warning of intraoperative adverse events. Furthermore, we discuss the potential value of AI in enhancing surgical efficiency, ensuring patient safety, and optimizing surgical education. On this basis, we further analyze the principal challenges and inherent risks confronting current AI systems, with the aim of informing future technological innovation and facilitating clinical translation.
7.Analysis of The Characteristics of Brain Functional Activity in Gross Motor Tasks in Children With Autism Based on Functional Near-infrared Spectroscopy Technology
Wen-Hao ZONG ; Qi LIANG ; Shi-Yu YANG ; Feng-Jiao WANG ; Meng-Zhao WEI ; Hong LEI ; Gui-Jun DONG ; Ke-Feng LI
Progress in Biochemistry and Biophysics 2025;52(8):2146-2162
ObjectiveBased on functional near-infrared spectroscopy (fNIRS), we investigated the brain activity characteristics of gross motor tasks in children with autism spectrum disorder (ASD) and motor dysfunctions (MDs) to provide a theoretical basis for further understanding the mechanism of MDs in children with ASD and designing targeted intervention programs from a central perspective. MethodsAccording to the inclusion and exclusion criteria, 48 children with ASD accompanied by MDs were recruited into the ASD group and 40 children with typically developing (TD) into the TD group. The fNIRS device was used to collect the information of blood oxygen changes in the cortical motor-related brain regions during single-handed bag throwing and tiptoe walking, and the differences in brain activation and functional connectivity between the two groups of children were analyzed from the perspective of brain activation and functional connectivity. ResultsCompared to the TD group, in the object manipulative motor task (one-handed bag throwing), the ASD group showed significantly reduced activation in both left sensorimotor cortex (SMC) and right secondary visual cortex (V2) (P<0.05), whereas the right pre-motor and supplementary motor cortex (PMC&SMA) had significantly higher activation (P<0.01) and showed bilateral brain region activity; in terms of brain functional integration, there was a significant decrease in the strength of brain functional connectivity (P<0.05) and was mainly associated with dorsolateral prefrontal cortex (DLPFC) and V2. In the body stability motor task (tiptoe walking), the ASD group had significantly higher activation in motor-related brain regions such as the DLPFC, SMC, and PMC&SMA (P<0.05) and showed bilateral brain region activity; in terms of brain functional integration, the ASD group had lower strength of brain functional connectivity (P<0.05) and was mainly associated with PMC&SMA and V2. ConclusionChildren with ASD exhibit abnormal brain functional activity characteristics specific to different gross motor tasks in object manipulative and body stability, reflecting insufficient or excessive compensatory activation of local brain regions and impaired cross-regions integration, which may be a potential reason for the poorer gross motor performance of children with ASD, and meanwhile provides data support for further unraveling the mechanisms underlying the occurrence of MDs in the context of ASD and designing targeted intervention programs from a central perspective.
8.Processing technology of calcined Magnetitum based on concept of QbD and its XRD characteristic spectra.
De-Wen ZENG ; Jing-Wei ZHOU ; Tian-Xing HE ; Yu-Mei CHEN ; Huan-Huan XU ; Jian FENG ; Yue YANG ; Xin CHEN ; Jia-Liang ZOU ; Lin CHEN ; Hong-Ping CHEN ; Shi-Lin CHEN ; Yuan HU ; You-Ping LIU
China Journal of Chinese Materia Medica 2025;50(9):2391-2403
Guided by the concept of quality by design(QbD), this study optimizes the calcination and quenching process of calcined Magnetitum and establishes the XRD characteristic spectra of calcined Magnetitum, providing a scientific basis for the formulation of quality standards. Based on the processing methods and quality requirements of Magnetitum in the Chinese Pharmacopoeia, the critical process parameters(CPPs) identified were calcination temperature, calcination time, particle size, laying thickness, and the number of vinegar quenching cycles. The critical quality attributes(CQAs) included Fe mass fraction, Fe~(2+) dissolution, and surface color. The weight coefficients were determined by combining Analytic Hierarchy Process(AHP) and the criteria importance though intercrieria correlation(CRITIC) method, and the calcination process was optimized using orthogonal experimentation. Surface color was selected as a CQA, and based on the principle of color value, the surface color of calcined Magnetitum was objectively quantified. The vinegar quenching process was then optimized to determine the best processing conditions. X-ray diffraction(XRD) was used to establish the characteristic spectra of calcined Magnetitum, and methods such as similarity evaluation, cluster analysis, and orthogonal partial least squares-discriminant analysis(OPLS-DA) were used to evaluate the quality of the spectra. The optimized calcined Magnetitum preparation process was found to be calcination at 750 ℃ for 1 h, with a laying thickness of 4 cm, a particle size of 0.4-0.8 cm, and one vinegar quenching cycle(Magnetitum-vinegar ratio 10∶3), which was stable and feasible. The XRD characteristic spectra analysis method, featuring 9 common peaks as fingerprint information, was established. The average correlation coefficient ranged from 0.839 5-0.988 1, and the average angle cosine ranged from 0.914 4 to 0.995 6, indicating good similarity. Cluster analysis results showed that Magnetitum and calcined Magnetitum could be grouped together, with similar compositions. OPLS-DA discriminant analysis identified three key characteristic peaks, with Fe_2O_3 being the distinguishing component between the two. The final optimized processing method is stable and feasible, and the XRD characteristic spectra of calcined Magnetitum was initially established, providing a reference for subsequent quality control and the formulation of quality standards for calcined Magnetitum.
X-Ray Diffraction/methods*
;
Drugs, Chinese Herbal/chemistry*
;
Quality Control
;
Particle Size
9.Four Weeks of HIIT Modulates Lactate-mediated Synaptic Plasticity to Improve Depressive-like Behavior in CUMS Rats
Yu-Mei HAN ; Zi-Wei ZHANG ; Jia-Ren LIANG ; Chun-Hui BAO ; Jun-Sheng TIAN ; Shi ZHOU ; Huan XIANG ; Yong-Hong YANG
Progress in Biochemistry and Biophysics 2025;52(6):1499-1510
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, Smax) alternated with one minute at low speed (50%-55% Smax), 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.
10.Prediction of suitable habitats of Phlebotomus chinensis in Gansu Province based on the Biomod2 ensemble model
Dawei YU ; Yandong HOU ; Aiwei HE ; Yu FENG ; Guobing YANG ; Chengming YANG ; Hong LIANG ; Hailiang ZHANG ; Fan LI
Chinese Journal of Schistosomiasis Control 2025;37(3):276-283
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 × 104 km2, of which the highly suitable area was 1.42 × 104 km2, 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.


Result Analysis
Print
Save
E-mail