ObjectiveThis paper aims to observe the syndrome improvement and negative antigen conversion rate of Qingxuan Daozhi formula in the treatment of influenza in children （heat accumulation in the lung and stomach syndrome）. MethodsThrough a multi-center randomized controlled methodology design，confirmed influenza cases were collected from October 2022 to April 2023 in the pediatrics department of eight hospitals，such as Dongfang Hospital of Beijing University of Chinese Medicine. A total of 180 children with influenza and heat accumulation in the lung and stomach syndrome conforming to the standard were recruited through the clinic. The sick children meeting the inclusion criteria were randomly divided into groups by a block-randomized method. The children in the experimental group were treated with Qingxuan Daozhi formula for five days，and those in the control group were treated with Oseltamivir Phosphate Granules for five days. The primary efficacy indicator was the negative conversion rate of influenza antigen detection. Secondary efficacy indicators were the Canadian acute respiratory illness and flu scale （CARIFS） and the incidence of complications，severe cases， and critical cases. Follow-up observation was conducted on the day of enrollment，48 hours after medication，72 hours after medication， and （6+1） d after medication. ResultsOne hundred and eighty participants were randomly assigned to the experimental group （90 cases） or the control group （90 cases）. All participants were followed up during the study. Comparison of influenza antigen detection results in the primary efficacy indicators showed that the average time of negative influenza antigen conversion in the experimental group was （5.29±1.25） d，and that in the control group was （5.40±1.68） d，without a statistically significant difference. After five days of intervention，52 cases in the experimental group and 51 cases in the control group converted to negative，without a statistically significant difference. CARIFS score results in the secondary efficacy indicators showed that during 72 hours after intervention，there were statistically significant differences between the experimental group and the control group in three dimensions， including headache，muscle soreness， and the need for extra care （P<0.05）. On the （6+1） days after the intervention，the differences in both the experimental group and the control group were statistically significant in 10 dimensions， including sore throat，bad sleep，uncomfortable feeling，poor spirit and fatigue，crying more than usual，the need for extra care，symptom，function，influence on parents，and total score （P<0.05）. The comparison results within the group in the dimensional scores of symptom， function， and influence on parents，as well as the CARIFS total score showed that with the delay of follow-up time，scores of both groups decreased significantly，with a statistically significant difference （P<0.01）. Inter-group comparison results showed that the mean score of the experimental group was higher than that of the control group at the time of enrollment. With the progress of intervention，the score of the experimental group was significantly decreased compared with that of the control group. At the end of follow-up，the mean score of the experimental group was lower than that of the control group，with no statistically significant difference. In terms of the incidence of complications，severe cases， and critical cases， there were no complications，severe cases， and critical cases in the two groups，without a statistically significant difference. ConclusionThe symptom improvement effect and negative antigen conversion rate of Qingxuan Daozhi formula in the treatment of influenza in children （heat accumulation in the lung and stomach syndrome） are not inferior to Oseltamivir Phosphate granules， and children's acceptance is better. It can be more widely used in clinical treatment of influenza in children （heat accumulation in the lung and stomach syndrome）.

OBJECTIVE To confirm trigger items for adverse drug reaction （ADR） induced by bevacizumab， to identify and analyze the occurrence of related ADR， and to establish a predictive model for severe adverse reaction （SAR） caused by this drug. METHODS Based on the global trigger tool （GTT） theory， and referencing the GTT White Paper， drug package inserts and relevant literature， trigger items for bevacizumab-related ADR were confirmed using a single-round Delphi method. Utilizing these established items， electronic medical records of relevant patients at Guilin People’s Hospital from January 2020 to September 2024 were actively screened via the China Hospital Pharmacovigilance System. Pharmacists then identified and tallied the occurrence of bevacizumab-induced ADR. Data from patients with any positive trigger item served as the study subjects （divided into training and test sets at a ratio of 7∶3）， candidate feature variables were selected from 39 related variables using the Boruta algorithm， and the multivariable Logistic regression analysis was performed with the occurrence of SAR as the dependent variable. Based on these candidate features， Logistic Regression， Extreme Gradient Boosting， Light Gradient Boosting Machine， Random Forest， and Categorical Boosting models were constructed. Model performance was evaluated using metrics including the area under the curve （AUC） of receiver operating characteristic curve and recall rate. The Shapley Additive exPlanations （SHAP） method was applied to analyze and interpret the contribution of each variable. A nomogram was constructed based on the optimal model. RESULTS A total of 38 trigger items for active monitoring of bevacizumab-related ADR were determined， comprising 17 laboratory indicators， 13 clinical manifestations， and 8 intervention measures. In total， 483 patients with positive trigger items were included， and 318 patients with bevacizumab-induced ADR were identified， including 83 SARs. The positive predictive values for the trigger items and cases were 43.57% （708/1 625） and 63.84% （318/483）， respectively. Bevacizumab-induced ADR involved 7 systems/organs， with the hematological system being the most frequently involved （64.15%）. The Boruta algorithm selected 7 vari ables： serum potassium， hematocrit， albumin-to-globulin ratio， prealbumin， hypertension history， age and red blood cell count. Multivariable Logistic regression showed that elevated serum potassium levels were associated with a decreased risk of bevacizumab-induced SAR （OR＝0.234， P ＝0.002）， while a history of hypertension （OR＝2.642， P ＝0.006） and increased age （OR＝1.040， P ＝0.025） were associated with an increased risk. The Logistic Regression model demonstrated superior performance with higher AUC， F1 score and recall rate （0.761， 0.447， 0.607）， compared to other models. SHAP evaluation results indicated that variables such as serum potassium， hematocrit， and age ranked highest in importance. CONCLUSIONS Totally 38 trigger entries have been successfully identified for active screening of bevacizumab-related ADR. Elevated serum potassium levels are a protective factor against bevacizumab-induced SAR， whereas the hypertension history and increased age are risk factors. The Logistic Regression model is the optimal predictive model.

Objective To analyze the clinical and therapeutic characteristics of Epstein-Barr virus (EBV)-negative posttransplant lymphoproliferative disease (PTLD) with diffuse large B-cell lymphoma (DLBCL) in the context of specific cases and literature. Methods A case of EBV-negative DLBCL (GCB type) after kidney transplantation is reported. The patient was a 45-year-old male who underwent living-related kidney transplantation in 2016 and has been receiving triple immunosuppressive therapy with tacrolimus, mycophenolate mofetil and methylprednisolone since then. In 2024, the patient presented with intermittent fever, night sweats and gastrointestinal symptoms. The diagnosis was confirmed by endoscopic pathology, immunohistochemical staining and positron emission tomography/computed tomography. The R-CDOP regimen (rituximab + cyclophosphamide + liposomal doxorubicin + vincristine + dexamethasone) was used for treatment. Results The patient was diagnosed with EBV-negative DLBCL (GCB type, Ann Arbor stage Ⅳ B). After 4 cycles of R-CDOP chemotherapy, the efficacy assessment was partial remission, and the transplant kidney function remained stable. Conclusions For EBV-negative PTLD after kidney transplantation, it is necessary to break through the "virus-dependent" diagnostic thinking. In clinical practice, the focus should be on protecting the transplant kidney, and individualized treatment plans should be developed for patients.

ObjectiveTo explore the potential mechanism of Xiaoqinglongtang（XQL） in the prevention and treatment of high altitude pulmonary edema（HAPE） by network pharmacology， molecular docking， and molecular dynamics simulation， and to verify it by in vivo animal model. MethodsIn this study， the active ingredients， drug targets， and HAPE-related targets of XQL were collected from BATMAN-TCM， GeneCards， and Online Mendelian Inheritance in Man（OMIM） databases. The protein-protein interaction（PPI） network was constructed by using intersection targets， and the core targets were screened and visualized by Cytoscape software. Functional annotation and pathway analysis of the intersection targets were performed by gene ontology （GO） and Kyoto encyclopedia of genes and genomes （KEGG） functional enrichment. AutoDock and GROMACS were used to evaluate the binding ability of active ingredients to key targets. In the experimental verification part， a mouse model of HAPE induced by hypobaric hypoxia（simulated 6 000 m altitude for 48 h） was established. The control effect was evaluated by hematoxylin-eosin（HE） staining， lung tissue water content， lung tissue wet/dry weight ratio， real-time quantitative polymerase chain reaction（Real-time PCR） detection of gene expression levels， and immunohistochemistry and Western blot detection of key protein expression. ResultsA total of 355 active ingredients of XQL， 2 142 targets， 716 HAPE-related targets， and 236 intersection targets were obtained by network pharmacology analysis. Key core targets such as interleukin （IL）-6， tumor necrosis factor （TNF）， protein kinase B1 （Akt1）， and hypoxia-inducible factor-1α （HIF-1α） were screened. The results of GO analysis of common targets involved 738 biological processes（BP）， 72 cellular components（CC）， and 135 molecular functions（MF）. KEGG analysis effectively enriched two important signaling pathways： Phosphoinositol 3-kinase （PI3K）/Akt and HIF-1α. The results of molecular docking and molecular dynamics simulation showed that the screened active ingredients had good binding ability with key targets. In the HAPE model induced by hypobaric hypoxia（6 000 m， 48 h）， the lung tissue water content， lung tissue wet/dry weight ratio， and pathological injury score of the model group were significantly increased（P<0.01）， accompanied by exudation of a large number of red blood cells in the alveoli and alveolar interstitium， a significant increase in inflammatory cells， a significant widening of the alveolar septum， and mutual fusion between the alveoli. The XQL administration group significantly improved the above pathological changes（P<0.01）. The results of inflammatory factor expression showed that compared with the control group， the model group showed significantly up-regulated expression of TNF-α， IL-6， and IL-1β in the lung tissue（P<0.01）. Compared with the model group， the XQL administration group had significantly decreased expression of inflammatory factors（P<0.05， P<0.01）. The mRNA expression of key pathway related genes PI3K， Akt1， mammalian target of rapamycin（mTOR）， and HIF-1α was significantly increased in the model group（P<0.01）， and decreased in a concentration-dependent manner after XQL administration（P<0.05， P<0.01）. The expression levels of key proteins PI3K， phosphorylation（p）-PI3K， Akt1， p-Akt1， mTOR， p-mTOR， and HIF-1α in lung tissue were analyzed by immunohistochemistry and Western blot. Compared with the blank group， the model group showed increased expression of key proteins（P<0.05， P<0.01）. Compared with the model group， the XQL administration group exhibited decreased expression of key proteins（P<0.05， P<0.01）. ConclusionXQL can reduce lung inflammation and improve HAPE. The mechanism may be related to the regulation of PI3K/Akt/mTOR and HIF-1α pathways. This study provides a new idea and a theoretical basis for the treatment of HAPE with XQL.

ObjectiveTo explore the potential mechanism of Xiaoqinglongtang（XQL） in the prevention and treatment of high altitude pulmonary edema（HAPE） by network pharmacology， molecular docking， and molecular dynamics simulation， and to verify it by in vivo animal model. MethodsIn this study， the active ingredients， drug targets， and HAPE-related targets of XQL were collected from BATMAN-TCM， GeneCards， and Online Mendelian Inheritance in Man（OMIM） databases. The protein-protein interaction（PPI） network was constructed by using intersection targets， and the core targets were screened and visualized by Cytoscape software. Functional annotation and pathway analysis of the intersection targets were performed by gene ontology （GO） and Kyoto encyclopedia of genes and genomes （KEGG） functional enrichment. AutoDock and GROMACS were used to evaluate the binding ability of active ingredients to key targets. In the experimental verification part， a mouse model of HAPE induced by hypobaric hypoxia（simulated 6 000 m altitude for 48 h） was established. The control effect was evaluated by hematoxylin-eosin（HE） staining， lung tissue water content， lung tissue wet/dry weight ratio， real-time quantitative polymerase chain reaction（Real-time PCR） detection of gene expression levels， and immunohistochemistry and Western blot detection of key protein expression. ResultsA total of 355 active ingredients of XQL， 2 142 targets， 716 HAPE-related targets， and 236 intersection targets were obtained by network pharmacology analysis. Key core targets such as interleukin （IL）-6， tumor necrosis factor （TNF）， protein kinase B1 （Akt1）， and hypoxia-inducible factor-1α （HIF-1α） were screened. The results of GO analysis of common targets involved 738 biological processes（BP）， 72 cellular components（CC）， and 135 molecular functions（MF）. KEGG analysis effectively enriched two important signaling pathways： Phosphoinositol 3-kinase （PI3K）/Akt and HIF-1α. The results of molecular docking and molecular dynamics simulation showed that the screened active ingredients had good binding ability with key targets. In the HAPE model induced by hypobaric hypoxia（6 000 m， 48 h）， the lung tissue water content， lung tissue wet/dry weight ratio， and pathological injury score of the model group were significantly increased（P<0.01）， accompanied by exudation of a large number of red blood cells in the alveoli and alveolar interstitium， a significant increase in inflammatory cells， a significant widening of the alveolar septum， and mutual fusion between the alveoli. The XQL administration group significantly improved the above pathological changes（P<0.01）. The results of inflammatory factor expression showed that compared with the control group， the model group showed significantly up-regulated expression of TNF-α， IL-6， and IL-1β in the lung tissue（P<0.01）. Compared with the model group， the XQL administration group had significantly decreased expression of inflammatory factors（P<0.05， P<0.01）. The mRNA expression of key pathway related genes PI3K， Akt1， mammalian target of rapamycin（mTOR）， and HIF-1α was significantly increased in the model group（P<0.01）， and decreased in a concentration-dependent manner after XQL administration（P<0.05， P<0.01）. The expression levels of key proteins PI3K， phosphorylation（p）-PI3K， Akt1， p-Akt1， mTOR， p-mTOR， and HIF-1α in lung tissue were analyzed by immunohistochemistry and Western blot. Compared with the blank group， the model group showed increased expression of key proteins（P<0.05， P<0.01）. Compared with the model group， the XQL administration group exhibited decreased expression of key proteins（P<0.05， P<0.01）. ConclusionXQL can reduce lung inflammation and improve HAPE. The mechanism may be related to the regulation of PI3K/Akt/mTOR and HIF-1α pathways. This study provides a new idea and a theoretical basis for the treatment of HAPE with XQL.

ObjectiveTo investigate the effect of different visual inputs and task conditions on balance function and lower limb biomechanical characteristics in children with spastic hemiplegia. MethodsFrom March to July, 2025, 30 children aged six to nine years old with spastic hemiplegia (hemiplegia group) and 30 healthy children (control group) were selected. A 2×2×2 mixed experimental design was employed, involving groups (hemiplegia vs. healthy), tasks (single-task vs. dual-task), and vision (eyes open vs. eyes closed). One week before test, they were evaluated with Fugl-Meyer Assessment-Lower Extremities (FMA-LE), Wee Function Independent Measurement (WeeFIM) and Wechsler Intelligence Scale for Children, Fourth Edition (WISC-Ⅳ). A 3D gait analysis system and a plantar pressure testing system were used to collect spatio-temporal parameters of walking speed, stride length, cadence and step width, kinematic parameters of hip/knee/ankle joint angles, kinetic parameters of forefoot loading ratio, center of pressure (CoP) displacement, while dual-task cost (DTC) was caculated. ResultsSpatio-temporal parameters showed that under dual-task and eyes-closed conditions, walking speed and stride length decreased while step width increased in the hemiplegia group. Significant interaction effects among group, task, and vision were observed for speed, stride length, and step width (F > 4.886, P < 0.05). Kinematic parameters indicated that during dual-tasks, the hemiplegia group exhibited increased hip flexion and decreased ankle dorsiflexion; under eyes-closed conditions, knee flexion increased. The interaction of the three factors significantly affected all joint angles (F > 4.876, P < 0.05). Kinetic parameters showed that under dual-task and eyes-closed conditions, the forefoot loading ratio and anteroposterior CoP displacement decreased, while mediolateral CoP displacement increased. The interaction of the three factors significantly affected CoP displacement (F > 4.355, P < 0.05). All the DTC was significantly higher in the hemiplegia group than in the control group, except DTC of the cadence (|t| > 14.393, P < 0.001). Correlation analysis revealed that the score of FMA-LE was strongly negatively correlated with DTC (|r| > 0.731, P < 0.01). The Functional Independence Measure for Children and Working Memory Index showed moderate negative correlations with the DTC of walking speed and cadence (|r| > 0.462, P < 0.05). ConclusionThe gait and balance of children with spastic hemiplegia are concurrently influenced by dual-tasking and visual input. The superposition of visual deprivation and dual-tasks significantly exacerbates gait abnormalities. Furthermore, is strongly correlated with motor function and working memory.

ObjectiveTo investigate the effect of different visual inputs and task conditions on balance function and lower limb biomechanical characteristics in children with spastic hemiplegia. MethodsFrom March to July, 2025, 30 children aged six to nine years old with spastic hemiplegia (hemiplegia group) and 30 healthy children (control group) were selected. A 2×2×2 mixed experimental design was employed, involving groups (hemiplegia vs. healthy), tasks (single-task vs. dual-task), and vision (eyes open vs. eyes closed). One week before test, they were evaluated with Fugl-Meyer Assessment-Lower Extremities (FMA-LE), Wee Function Independent Measurement (WeeFIM) and Wechsler Intelligence Scale for Children, Fourth Edition (WISC-Ⅳ). A 3D gait analysis system and a plantar pressure testing system were used to collect spatio-temporal parameters of walking speed, stride length, cadence and step width, kinematic parameters of hip/knee/ankle joint angles, kinetic parameters of forefoot loading ratio, center of pressure (CoP) displacement, while dual-task cost (DTC) was caculated. ResultsSpatio-temporal parameters showed that under dual-task and eyes-closed conditions, walking speed and stride length decreased while step width increased in the hemiplegia group. Significant interaction effects among group, task, and vision were observed for speed, stride length, and step width (F > 4.886, P < 0.05). Kinematic parameters indicated that during dual-tasks, the hemiplegia group exhibited increased hip flexion and decreased ankle dorsiflexion; under eyes-closed conditions, knee flexion increased. The interaction of the three factors significantly affected all joint angles (F > 4.876, P < 0.05). Kinetic parameters showed that under dual-task and eyes-closed conditions, the forefoot loading ratio and anteroposterior CoP displacement decreased, while mediolateral CoP displacement increased. The interaction of the three factors significantly affected CoP displacement (F > 4.355, P < 0.05). All the DTC was significantly higher in the hemiplegia group than in the control group, except DTC of the cadence (|t| > 14.393, P < 0.001). Correlation analysis revealed that the score of FMA-LE was strongly negatively correlated with DTC (|r| > 0.731, P < 0.01). The Functional Independence Measure for Children and Working Memory Index showed moderate negative correlations with the DTC of walking speed and cadence (|r| > 0.462, P < 0.05). ConclusionThe gait and balance of children with spastic hemiplegia are concurrently influenced by dual-tasking and visual input. The superposition of visual deprivation and dual-tasks significantly exacerbates gait abnormalities. Furthermore, is strongly correlated with motor function and working memory.

ObjectiveTo investigate the effect of different visual inputs and task conditions on balance function and lower limb biomechanical characteristics in children with spastic hemiplegia. MethodsFrom March to July, 2025, 30 children aged six to nine years old with spastic hemiplegia (hemiplegia group) and 30 healthy children (control group) were selected. A 2×2×2 mixed experimental design was employed, involving groups (hemiplegia vs. healthy), tasks (single-task vs. dual-task), and vision (eyes open vs. eyes closed). One week before test, they were evaluated with Fugl-Meyer Assessment-Lower Extremities (FMA-LE), Wee Function Independent Measurement (WeeFIM) and Wechsler Intelligence Scale for Children, Fourth Edition (WISC-Ⅳ). A 3D gait analysis system and a plantar pressure testing system were used to collect spatio-temporal parameters of walking speed, stride length, cadence and step width, kinematic parameters of hip/knee/ankle joint angles, kinetic parameters of forefoot loading ratio, center of pressure (CoP) displacement, while dual-task cost (DTC) was caculated. ResultsSpatio-temporal parameters showed that under dual-task and eyes-closed conditions, walking speed and stride length decreased while step width increased in the hemiplegia group. Significant interaction effects among group, task, and vision were observed for speed, stride length, and step width (F > 4.886, P < 0.05). Kinematic parameters indicated that during dual-tasks, the hemiplegia group exhibited increased hip flexion and decreased ankle dorsiflexion; under eyes-closed conditions, knee flexion increased. The interaction of the three factors significantly affected all joint angles (F > 4.876, P < 0.05). Kinetic parameters showed that under dual-task and eyes-closed conditions, the forefoot loading ratio and anteroposterior CoP displacement decreased, while mediolateral CoP displacement increased. The interaction of the three factors significantly affected CoP displacement (F > 4.355, P < 0.05). All the DTC was significantly higher in the hemiplegia group than in the control group, except DTC of the cadence (|t| > 14.393, P < 0.001). Correlation analysis revealed that the score of FMA-LE was strongly negatively correlated with DTC (|r| > 0.731, P < 0.01). The Functional Independence Measure for Children and Working Memory Index showed moderate negative correlations with the DTC of walking speed and cadence (|r| > 0.462, P < 0.05). ConclusionThe gait and balance of children with spastic hemiplegia are concurrently influenced by dual-tasking and visual input. The superposition of visual deprivation and dual-tasks significantly exacerbates gait abnormalities. Furthermore, is strongly correlated with motor function and working memory.

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.