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.

Objective To examine the impact of extreme temperature and drought stress on the survival of Bulinus globosus, so as to provide the theoretical evidence for the genomic research of Bulinus in absence of reference genes. Methods B. globosus snail samples were collected from Kiwani Shehia in Pemba Island, Zanzibar, Tanzania, and offspring snails were obtained through laboratory breeding and reproduction. A total of 120 10-week-old B. globosus snails from the same generation were selected and randomly assigned into four groups, including the high-temperature drought (HD) group, normal temperature drought (D) group, low-temperature drought (LD) group, and the control (C) group, of 30 snails in each group. Snails in HD, D, and LD groups were placed in beakers containing dry soil at the bottom and subsequently housed in climate chambers at 35, 26 ℃, and 10 ℃, respectively, while snails in Group C were maintained in 500 mL petri dishes containing dechlorinated tap water at 26 ℃. Following 3 days of breeding, living snails in each group were collected, and soft tissues were dissected and isolated. Total RNA was extracted from snail soft tissues for library construction, followed by high-throughput sequencing on the Illumina HiSeq 4000 sequencing system. De novo transcriptome assembly was performed using the Trinity software, and the longest transcripts were selected as unigenes. Gene functional annotations of unigenes were conducted using the Diamond software against Gene Ontology (GO) knowledgebase, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, NCBI non-redundant (NR) protein sequences database, Protein Family (Pfam) database, and UniProtKB/Swiss-Prot (Swiss-Prot) knowledgebase. GO and KEGG enrichment analyses of differentially expressed genes (DEGs) were performed using the topGO and clusterProfiler software, respectively. In addition, four relevant genes were selected for validation using a real-time quantitative PCR (qRT-PCR) assay to verify the reliability of transcriptome sequencing results. Results Following 3 days of breeding, there were 7, 20, 28, and 30 survival B. globosus snails in HD, LD, D, and C groups, with corresponding survival rates of 23.33% (7/30), 66.67% (20/30), 93.33% (28/30), and 100.00% (30/30), respectively (χ² = 52.72, P < 0.001). De novo transcriptome assembly generated 176 942 unigenes, with annotation rates of 0.98%, 13.49%, 26.46%, 12.48%, and 14.39% against GO knowledgebase, KEGG pathway database, NR protein sequences database, Pfam database, and Swiss-Prot knowledgebase, respectively. There were 33 up-regulated and 72 down-regulated genes in Group D, 483 up-regulated and 815 down-regulated genes in Group HD, and 245 up-regulated and 172 down-regulated genes in Group LD relative to in Group C. Following removal of overlapping genes across groups and unmatched genes, 11 candidate genes were identified. GO and KEGG analyses revealed 3 heat shock protein (HSP)-related DEGs in these 11 candidate genes, which were annotated as HSP12.2, HSP70, and HSP20 genes and were all significantly up-regulated in each treatment group. Three immune and nervous system-related DEGs were identified, and were all significantly down-regulated in each treatment group, which were involved in the neural cell adhesion molecule L1-like protein pathway, fibrinogen binding protein pathway, and leukocyte elastase inhibitor-like protein pathway. qRT-PCR assay quantified that the expression trends of four genes related to temperature and drought stress across different treatment groups were highly consistent with transcriptome sequencing data. Conclusion The survival rate of B. globosus significantly reduces under combined stresses of extreme temperature and drought, possibly due to an imbalance in its cellular homeostasis regulatory system.

Objective Objective To analyze the potential spatial factors affecting the genetic differentiation of Oncomelania hupensis robertsoni in Yunnan Province. Methods A total of 13 administrative villages were selected from schistosomiasis-endemic areas of Yunnan Province as O. hupensis snail sampling sites. At least 200 snails were collected in each site, and the spatial variable data of each site were recorded, including longitude, latitude and altitude. Thirty active and Schistosoma japonicum uninfected O. hupensis snails were selected from each sampling site by means of the crawling method and the cercarial shedding method. Genomic DNA was extracted from O. hupensis snails. Following PCR amplification, purification of PCR amplification products and sequencing, the gene sequences of O. hupensis snail samples were spliced and edited using the DNAstar software and the NCBI database to yield the complete mitochondrial sequences of O. hupensis snails at each sampling site, and the mitochondrial genetic distance matrix of O. hupensis robertsoni was calculated at each sampling site. The geographical coordinates of each sampling site were marked using the software ArcGIS 10.2, and the straight-line geographical distance between each sampling site was calculated. The altitude difference, longitude difference and latitude difference between each sampling site were calculated using the Excel software, and the correlation between the mitochondrial genetic distance matrix of O. hupensis robertsoni and each spatial variable matrix was examined by using the Mantel test at 13 sampling sites in Yunnan Province. Results Among the 13 O. hupensis snail sampling sites in Yunnan Province, the largest mitochondrial genetic distance of O. hupensis robertsoni snail populations was seen between Anding Village, Nanjian Yi Autonomous County and Caizhuang Village, Midu County (26.244 2), and the largest geographical distance was seen between Dongyuan Village, Gucheng District and Cangling Village, Chuxiong County (272.64 km). The highest altitude difference was seen between Anding Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1 086.10 m), and the largest longitude difference was found between Qiandian Village, Eryuan County and Cangling Village, Chuxiong County (1.86°), while the largest latitude difference was measured between Leqiu Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1.81°). In addition, the mitochondrial genetic distance of O. hupensis robertsoni snail populations was positively correlated with altitude at 13 snail sampling sites in Yunnan Province (r = 0.542 8, P < 0.001), and showed no significant correlations with geographical distance (r = 0.093 4, P > 0.05), longitude (r = −0.199 5, P > 0.05) or latitude (r = 0.205 7, P > 0.05). Conclusion Altitude may be a potential spatial factor affecting the genetic differentiation of O. hupensis robertsoni in Yunnan Province.

Based on review of global and Chinese schistosomiasis control progress and the evolution of control strategies, this article focuses on Chinese practical experiences of schistosomiasis control and systematically summarizes five key determinants for schistosomiasis transmission risk, including source of infections, intermediate host snails, high-risk populations, natural environments, and social factors. To address these risks and challenges associated with these determinants, the article reviews the advances in techniques for assessment of schistosomiasis transmission risk and their applications, including conventional risk assessment approaches, mathematical model-based tools for prediction of schistosomiasis transmission risk, and indicator-systembased techniques for assessment of schistosomiasis transmission risk. This review underscores the essential role of interdisciplinary integration and dynamic management in precision schistosomiasis control and recommends the intensification of verification of field adaptation and dynamic updates of indicator systems to promote the widespread application of assessment tools across diverse regions and contexts, so as to provide strategic guidance and methodological support to achieve the target for elimination of schistosomiasis across China in 2030.

ObjectiveTo analyze the research hotspots and development trends in the field of spinal cord stimulation (SCS) for spinal cord injury (SCI). MethodsLiterature about SCS for SCI was retrieve from the Web of Science (WOS) Core Collection database, with a time range from January, 1999 to July, 2025. VOSviewer 1.6.20 and CiteSpace 6.4.R2 were used to analyze the annual publication volume, countries, authors, institutions, journals and keywords. ResultsA total of 636 literatures were included. From 1999 to 2025, the overall publication trend in this field showed an upward trajectory, with recent years fluctuating but tending to stabilize. The country with the most publications was the United States (429 papers), followed by Russia (98 papers) and China (70 papers). The institution with the highest number of publications was the University of California, Los Angeles (76 papers), the author with the most publications was V. Reggie Edgerton (70 papers), and the journal with the most publications was Journal of Clinical Medicine (31 papers). The most frequently cited study focused on exploring the combination of epidural spinal cord stimulation with task-specific training to restore motor function in patients with complete SCI. Keyword analysis showed that the research hotspots in this field were mainly focused on neuroregulation mechanisms, recovery of motor and autonomic nervous dysfunction, artificial intelligence, closed-loop stimulation and brain-computer interface technology innovations. In recent years, the research focus gradually shifted from basic mechanisms to personalized and precise multifunctional rehabilitation strategies. ConclusionThe field of SCS for SCI has undergone phases of basic mechanism exploration and clinical application expansion. Current research hotspots and future trends focus primarily on the development of new stimulation paradigms and combined innovative technologies.

ObjectiveTo delineate imaging findings using an imaging platform and investigate the correlation between MRI characteristics of spinal cord atrophy and clinical diagnosis in children with spinal cord injury (SCI). MethodsImaging data of 150 children with SCI admitted to Beijing Bo'ai Hospital, China Rehabilitation Research Center, from January, 2002 to March, 2024 were collected and imported into the imaging platform. The anteroposterior and transverse diameters of the middle part of the spinal cord at the cross-section with the most severe atrophy were measured, and the relevant indicators of the previous normal spinal cord segment were measured as controls; the radiomic features were extracted. Clinical data of the children including gender, age, cause of injury, sensory level, motor level, spinal cord injury level, injury severity and disease course were collected. ResultsSpinal cord atrophy was identified in 81 cases (54%), among which 78 cases (96%) were American Spinal Injury Association Impairment Scale (AIS) grade A and 3 cases (4%) were AIS grade C. The upper boundary of the spinal cord atrophy site strongly correlated with the injury level, motor level and sensory level (r > 0.8, P < 0.001). ConclusionMore than half of children with SCI may develop secondary spinal cord atrophy, the vast majority of whom suffer from complete spinal cord injury; the upper boundary of spinal cord atrophy is correlated with the injury level.

ObjectiveTo explore the effect of piercing electroacupuncture on ankle muscle strength, dynamic balance and subjective stability in patients with functional ankle instability (FAI). MethodsA total of 48 patients with FAI admitted to the Affiliated Hospital of Traditional Chinese Medicine, Shandong Second Medical University from April, 2024 to July, 2025 were enrolled and randomly divided into control group (n = 24) and experimental group (n = 24). The control group received multimodal functional training, while the experimental group received piercing electroacupuncture in addition, for six weeks. Relative peak torque of ankle plantarflexion, dorsiflexion, inversion and eversion were measured with isokinetic muscle strength testing system before and after treatment, while Y-Balance Test (YBT) and Cumberland Ankle Instability Tool (CAIT) were adopted for evaluation. ResultsOne case dropped out in each group. Relative peak torque of ankle plantarflexion, dorsiflexion, inversion and eversion at each angular velocity, the normalized distance in each direction and composite score of YBT, as well as the score of CAIT increased in both groups after treatment (|t| > 2.891, P < 0.01), while all the indicators were better in the experimental group than in the control group (|t| > 2.104, P < 0.05), except dorsiflexion muscle strength. ConclusionCombination of piercing electroacupuncture can further improve perimalleolar muscle strength, dynamic balance and subjective stability in patients with FAI.

ObjectiveTo explore the application of virtual reality biofeedback training combined with oral localization therapy in dysphagia after oral cancer surgery. MethodsFrom May, 2023 to July, 2024, 86 patients with dysphagia after oral cancer surgery in Zhejiang Provincial People's Hospital were randomly divided into control group (n = 43) and experimental group (n = 43). The control group received conventional swallowing function training, while the experimental group added virtual reality biofeedback training combined with oral positioning therapy, for four weeks. The Standardized Swallowing Function Assessment Scale (SSA), Functional Oral Intake Scale (FOIS) and M.D.Anderson Dysphagia Inventory (MDADI) were used for evaluation before intervention, and two weeks, four weeks and eight weeks after intervention. ResultsFor scores of SSA , the main effects of group (F = 150.190, P < 0.001, η²_p = 0.641) and time (F = 230.870, P < 0.001, η²_p = 0.733), as well as the interaction effect (F = 16.910, P < 0.001, η²_p = 0.168) were all significant. For scores of FOIS, the main effects of group (F = 59.601, P < 0.001, η²_p = 0.415) and time (F = 89.464, P < 0.001, η²_p = 0.516), as well as the interaction effect (F = 7.990, P < 0.001, η²_p = 0.087) were all significant. For scores of MDADI, the main effects of group (F = 33.133, P < 0.001, η²_p = 0.283) and time (F = 49.650, P < 0.001, η²_p = 0.371), as well as the interaction effect (F = 3.224, P = 0.023, η²_p = 0.037) were all significant. ConclusionVirtual reality biofeedback training combined with oral localization therapy could improve the swallowing function, oral feeding ability and overall quality of life of patients with dysphagia after oral cancer surgery.

ObjectiveTo explore the effect of presentation-assimilation-discussion (PAD) classroom model-based classroom-community collaborative teaching on course of sports rehabilitation therapy under the background of sports-medicine integration. MethodsFrom Febrary, 2024 to June, 2025, a total of 46 undergraduate students majoring in sports rehabilitation (Grade 2022, Capital University of Physical Education and Sports) were randomly divided into control group (n = 23) and experimental group (n = 23). The control group received routine lectures combined with laboratory practice, while the experimental group was taught using a classroom-community collaborative teaching based on PAD classroom structure. The two groups were compared in terms of excellent rate in practical assessment, achievement in higher-order dimensions of Bloom's Taxonomy, real-case rehabilitation scenario transformation ability and teaching satisfaction at the end of the semester. ResultsAt the end of the semester, the excellent rate in practical assessment was 60.87% in the experimental group, higher than 21.74% in the control group (χ² = 7.725, P < 0.05), while the scores in the analyse, evaluation and creation dimensions of Bloom's Taxonomy were significantly higher (t > 8.490, P < 0.001), real-case rehabilitation scenario transformation ability was better (t > 3.873, P < 0.001) and teaching satisfaction was higher (χ² = 8.576, P = 0.003). ConclusionPAD-based classroom-community collaborative teaching can effectively improve practical skills, higher-order cognitive abilities and practical application abilities among students majoring in sports rehabilitation.

The Golgi body, a core organelle in eukaryotic cells, plays a critical role in protein modification, sorting, vesicular transport, and serves as a key site for lipid synthesis and glycosylation. Glucose and lipid metabolism are central processes for cellular energy maintenance and biosynthesis, and are closely linked to Golgi function. Recent studies have revealed the extensive involvement of the Golgi body in regulating glucose and lipid metabolism, where maintaining its structural and functional homeostasis is crucial for normal physiological activity. Under various stress conditions such as acidosis, hypoxia, and nutrient deficiency, the Golgi body undergoes structural and functional disruption, leading to Golgi stress. This in turn activates specific signaling pathways, such as those mediated by the cAMP-responsive element binding protein 3 (CREB3) and proteoglycans, to alleviate Golgi stress and enhance Golgi function. Golgi stress contributes to glucose and lipid metabolic disorders by affecting the activity of insulin receptors, glucose transporters, and lipid metabolism-related enzymes. For example, Golgi stress triggers the cleavage and release of the active fragment of CREB3, which enters the nucleus and upregulates the transcription of ADP-ribosylation factor 4 (ARF4) and key gluconeogenic enzymes, including phosphoenolpyruvate carboxykinase (PEPCK) and glucose-6-phosphatase (G6Pase). ARF4 promotes vesicle retrograde transport between the Golgi and endoplasmic reticulum, maintains secretory capacity, and enhances hepatic glucose output. This pathway is particularly active under high-fat or lipotoxic stress, leading to fasting hyperglycemia. When damaged Golgi components accumulate beyond a tolerable threshold, the cell initiates an autophagic response, selectively encapsulating the damaged Golgi into autophagosomes, which then fuse with lysosomes to form autolysosomes, leading to Golgiphagy. This process results in the degradation and clearance of damaged Golgi, thereby regulating Golgi quantity, quality, and function. Golgiphagy also plays a significant role in regulating glucose and lipid metabolism. For instance, under high-glucose conditions, autophagic flux may be suppressed, impairing the timely clearance and renewal of damaged Golgi, compromising its normal function, and further exacerbating glucose metabolism disorders. Additionally, Golgiphagy may participate in lipid degradation and influence lipid synthesis and transport. Research indicates that Golgi stress and Golgiphagy play important roles in glucose and lipid metabolism-related diseases. For example, the leucine zipper protein (LZIP) under Golgi stress conditions can promote hepatic steatosis. In mouse primary cells and human tissues, LZIP induces the expression of apolipoprotein A-IV (APOA4), which increases peripheral free fatty acid uptake, resulting in lipid accumulation in the liver and contributing to the development of fatty liver disease. This review systematically outlines the structure and function of the Golgi apparatus, the molecular regulatory mechanisms of Golgi stress and Golgiphagy, and their synergistic roles. It further elaborates on how Golgi stress and Golgiphagy participate in the regulation of glucose and lipid metabolism, discusses their clinical significance in related diseases such as diabetes, fatty liver disease, and obesity, and highlights potential novel therapeutic strategies from the perspective of Golgi-targeted medicine. Finally, this article addresses the challenges and future directions in Golgi-targeted interventions, aiming to advance the clinical translation of such strategies and foster breakthroughs in the treatment of glucose and lipid metabolism-related disorders.