Objective To investigate the distribution characteristics of Oncomelania hupensis snails in Yunnan Province fol-lowing interruption of schistosomiasis transmission, so as to provide the evidence for assessing the risk of schistosomiasis transmission and scientifically formulating the schistosomiasis surveillance program. Methods According to the requirements of the National Schistosomiasis Surveillance Scheme (2020 Edition), O. hupensis snail surveillance data were collected from 18 schistosomiasis-endemic counties (cities, districts) in Yunnan Province from 2020 to 2024, including area of snail survey, area of snail habitats, area of re-emerging snail habitats, number of frames surveyed, number of frames with O. hupensis snails, number of O. hupensis snails captured, and number of living snails, and the occurrence of frames with snails and mean density of living snails were calculated. Changes in snail status over the 5-year period from 2020 to 2024 and the differences in snail distributions specified by epidemic intensity, environmental type, and vegetation type were analyzed. Results The areas of snail survey increased from 1 727.96 hm² in 2020 to 3 894.45 hm² in 2024 (peak) across 18 schistosomiasis-endemic counties (cities, districts) in Yunnan Province during the period from 2020 through 2024. The areas of snail habitats increased from 70.36 hm² in 2020 to a peak in 2023 (172.04 hm²), followed by a reduction to 132.36 hm² in 2024, and the areas of re-emerging snail habitats increased from 42.71 hm² in 2020 to a peak in 2022 (78.43 hm²), followed by a reduction to 40.21 hm² in 2024. The occurrence of frames with snails and mean density of living snails increased from 1.24% (3 025/244 404) and (0.033 2 ± 0.038 7) snails/0.1 m² in 2020 to peaks at 2.03% (6 231/307 563) and (0.066 9 ± 0.068 4) snails/0.1 m² in 2023, followed by reductions to 1.04% (5 829/559 941) and (0.032 6 ± 0.057 7) snails/0.1 m² in 2024, respectively. There was a significant difference in the occurrence of frames with snails over the 5-year study period (χ² = 1 962.95, P < 0.05), and the occurrence of frames with snails reduced by 48.71% in 2024 relative to in 2023 (χ² = 1 411.05, P < 0.005); however, there was no significant difference in the mean density of living snails over the 5 years (H = 5.310, P > 0.05). There were significant differences in the occurrence of frames with snails (χ² = 481.27, P < 0.05) and mean density of living snails (H = 6.872, P < 0.05) in schistosomiasis-endemic areas with different epidemic intensities. The occurrence of frames with snails (χ² = 25.32 and 38.70, both P values < 0.017) and mean density of living snails (Z = 28.55 and 49.96, both P values < 0.017) were higher in schistosomiasis transmission-interrupted and eliminated areas with snails than in schistosomiasis-eliminated areas without snails, and the occurrence of frames with snails (χ² = 453.54, P < 0.017) and mean density of living snails (Z = −56.97, P < 0.017) were higher in schistosomiasis-eliminated areas with snails than in schistosomiasis transmission-interrupted areas with snails. O. hupensis snails were mainly distributed in paddy fields, dry farmlands and ditches; however, the occurrence of frames with snails (13.40%, 424/3 164) and mean density of living snails [(0.252 8 ± 0.158 7) snails/0.1 m²] were higher in ponds/weirs than in other types of environments (both P values < 0.05). Rice, dry farmland crops and weeds were main vegetations in which O. hupensis snails were distributed, and the occurrence of frames with snails (2.29%, 7 111/310 140) and mean density of living snails [(0.072 3 ± 0.018 9) snails/0.1 m²] were higher in weeds than in other types of environments (both P values < 0.05). Conclusions O. hupensis snails have been effectively controlled in Yunnan Province following implementation of integrated schistosomiasis control measures; however, there are still risk factors for schistosomiasis transmission, including reduced attention to schistosomiasis control and snail re-emergence. Improved control efforts and surveillance system construction and timely identification of risk factors of snail status and timely management are recommended to ensure the achievement of the target of schistosomiasis elimination as scheduled.

Cationic carriers have demonstrated broad application prospects in drug delivery due to their excellent drug-loading capacity and delivery performance. However, their high-density positive surface charge often leads to systemic toxicity and nonspecific uptake, posing significant barriers to clinical translation. In recent years, the emergence of charge shielding and stimuli-responsive strategies has provided effective avenues for modulating biocompatibility and targeting specificity. This review systematically summarizes the applications of chemical modification, natural polymer coating, and biomimetic membrane strategies in charge shielding. Furthermore, it explores the roles of endogenous stimuli such as pH, enzymes, and reactive oxygen species, as well as exogenous triggers like light and ultrasound, in achieving precise activation and controlled release. With the integration of multi-functional modules and the development of intelligent delivery platforms, cationic carriers are progressively advancing from laboratory research toward clinical translation. This review also discusses the translational potential and critical technical bottlenecks of related delivery systems, aiming to provide a theoretical framework and some reference for the design of next-generation smart delivery systems.

Objective To construct a multivariable prediction model for assessing the risk of cardiopulmonary complication after thoracoscopic lobectomy in patients with stage Ⅰ-Ⅱ non-small cell lung cancer(NSCLC).Methods Clinical data of 600 patients with stage Ⅰ-Ⅱ NSCLC who underwent thoracoscopic lobectomy were retrospectively analyzed.Patients were divided into the complication group(84 cases)and the non-complication group(516 cases)based on the occurrence of postoperative complication within 7 days,including atelectasis,pulmonary embolism,respiratory failure,chylothorax,massive pleural effusion,hypoxemia and atrial fibrillation.Demographic characteristics,preoperative pulmonary function assessment,pathological features and perioperative indicators were collected.Multivariate Logistic regression analysis was used to identify independent risk factors influencing postoperative cardiopulmonary complication in patients and construct a predictive model.Internal validation was performed using the Bootstrap resampling method(1 000 iterations)to evaluate the discrimination,calibration and clinical decision-making value of the model.Results Multivariate Logistic regression analysis identified the following independent risk factors of postoperative cardiopulmonary complication,including age(OR=1.832,95%CI:1.537-2.183),history of chronic obstructive pulmonary disease(COPD)(OR=6.782,95%CI:2.685-17.130),Karnofsky performance status(KPS)score(OR=0.926,95%CI:0.888-0.965),the percentage of forced expiratory volume in the first second to the predictive value(FEV1%pred)(OR=0.906,95%CI:0.845-0.972),the percentage of diffusing capacity for carbon monoxide to the expected value(DLCO%pred)(OR=0.901,95%CI:0.832-0.975),intraoperative blood loss(OR=1.025,95%CI:1.014-1.036)and one-lung ventilation time(OR=1.057,95%CI:1.034-1.080).The area under the curve(AUC)of the combined diagnosis was 0.977(95%CI:0.965-0.989),with 96.4%sensitivity and 87.6%specificity.The Hosmer-Lemeshow test indicated excellent calibration(χ2=1.285,P=0.994).Decision curve analysis demonstrated significant clinical net benefit when the risk threshold probability ranged between 20%and 98%.Conclusion The multivariable prediction model for cardiopulmonary complication after thoracoscopic lobectomy in stage Ⅰ-Ⅱ NSCLC patients exhibits strong predictive performance.

OBJECTIVE: In Pakistan, traditional medicines are an important component of the medical system, with numerous varieties and great demands. However, due to the scattered resources and the lack of systematic collection and collation, adulteration of traditional Pakistani medicine (TPM) is common, which severely affects the safety of their medicinal use and the import and export trades. Therefore, it is urgent to systematically organize and unify the management of TPM and establish a set of standards and operable methods for the identification of TPM. METHODS: We collected and organized the information on 128 TPMs with regard to their medicinal parts, efficacy, usage, and genetic material, based on Pakistan Hamdard Pharmacopoeia of Eastern Medicine: Pharmaceutical Codex. The genetic information of TPM is summarized from national center for biotechnology information (NCBI) and global pharmacopoeia genome database (GPGD). Furthermore, we utilized bioinformatics technology to supplement the chloroplast genome (cp-genome) data of 12 TPMs. To build the web server, we used the Linux + Apache + MySQL + PHP (LAMP) system and constructed the webpage on a PHP: Hypertext Preprocessor (PHP) model view controller (MVC) framework. RESULTS: We constructed a new genomic database, the traditional Pakistani medicine genomic database (TPMGD). This database comprises five entries, namely homepage, medicinal species, species identification, basic local alignment search tool (BLAST), and download. Currently, TPMGD contains basic profiles of 128 TPMs and genetic information of 102 TPMs, including 140 cytochrome c oxidase subunit I (COI) sequences and 119 mitochondrial genome sequences from Bombyx mori, 1 396 internal transcribed spacer 2 (ITS2) sequences and 1 074 intergenic region (psbA-trnH) sequences specific to 92 and 83 plant species, respectively. Additionally, TPMGD includes 199 cp-genome sequences of 82 TPMs. CONCLUSION TPMGD is a multifunctional database that integrates species description, functional information inquiry, genetic information storage, molecular identification of TPM, etc. The database not only provides convenience for TPM information queries but also establishes the scientific basis for the medication safety, species identification, and resource protection of TPM.

Objective:To investigate the feasibility and therapeutic efficacy of robot-assisted radical resection for hilar cholangiocarcinoma.Methods:This is a retrospective case series study. The clinical data of 29 patients who underwent robot-assisted radical resection for hilar cholangiocarcinoma at the Department of Minimally Invasive Hepatic Surgery,the First Affiliated Hospital of Harbin Medical University from July 2021 to February 2025 were retrospectively collected. There were 16 males and 13 females, aged ( M(IQR)) 68.0 (10.0) years (range:36 to 78 years), and body mass index (24.0±2.9) kg/m 2 (range:17.5 to 29.1 kg/m 2). Bismuth-Corlette classification: 12 cases type Ⅰ, 4 cases type Ⅱ, 6 cases type Ⅲb, and 7 cases type Ⅳ. Preoperative CA19-9 was 161.7(320.9) U/ml (range:7.1 to 1 000.0 U/ml), and carcinoembryonic antigen was 2.8(2.1)μg/L (range:0.3 to 203.1 μg/L). Preoperative total bilirubin was 134.2 (348.9) μmol/L (range:10.4 to 557.9 μmol/L), direct bilirubin was 90.8 (264.1) μmol/L (range:2.5 to 418.7 μmol/L), ALT was 136.4 (134.8) U/L (range:13.0 to 569.9 U/L), AST was 122.2 (119.9) U/L (range:16.0 to 384.0 U/L), and albumin was (34.5±6.3) g/L (range:21.7 to 41.3 g/L). Comparison of quantitative data at different time points using paired t-test or Mann-Whitney U test. Cox univariate analysis was performed for the relevant variables, and Cox multivariate regression analysis was used to screen the independent prognostic factors of patients after robot-assisted radical resection for hilar cholangiocarcinoma. Results:All the 29 patients successfully underwent robot-assisted radical resection of hilar cholangiocarcinoma, and the R0 resection rate was 93.1% (27/29) without conversion to laparotomy. The operation time was 295.0 (87.5) minutes (range:195 to 590 minutes), the intraoperative blood loss was 100.0 (150.0) ml (range:20 to 1 000 ml), the intraoperative blood transfusion rate was 20.1% (6/29), the number of lymph nodes dissected was 10.0 (7.0) pieces (range: 6 to 18 pieces), the first postoperative deflatus time was 3.0 (1.0) days (range:2 to 4 days), The oral feeding time was 5.0 (1.0) days (range: 4 to 7 days), the drainage tube removal time was 8.0 (2.0) days (range: 6 to 26 days), and the postoperative hospital stay time was 10.0 (6.0) days (range:7 to 27 days). The incidence of complications above grade Ⅱ of the Clavien-Dindo complication grading system was 24.1% (7/29), including 3 cases of gastrointestinal bleeding with recurrent high fever, 1 case of delayed gastric emptying, 1 case of bile leakage, and 5 cases of hypoalbuminemia. The total bilirubin was 42.8 (66.8) μmol/L (range:6.8 to 195.9 μmol/L), direct bilirubin was 28.1 (38.5) μmol/L (range:4.3 to 88.6 μmol/L), ALT was 55.8 (56.0) U/L (range:9.9 to 207.1 U/L), AST was 33.9 (17.9) U/L (range:10.6 to 122.7 U/L), and albumin was (32.1±3.8) g/L (range:22.8 to 37.7 g/L), the levels of transaminase and bilirubin in the postoperative liver function indexes were significantly improved compared with those before operation, and the differences were statistically significant (all P<0.05). The mortality rate of patients without perioperative death was 3.4% (1/29) at 90 days after surgery. The results of Cox multivariate regression analysis showed that R0 resection was an independent prognostic factor for survival at 1 year after surgery ( P<0.05). The follow-up time was 15.0 (12.0) months (range:6 to 30 months), 1 of the 29 patients died of intra-abdominal infection 1 week after discharge, and the remaining 28 patients were completely followed up, of which 20 patients had no recurrence and metastasis during the follow-up period, and the tumor-free survival was 15.0 (12.0) months (range:6 to 30 months), the tumor-free survival rate was 65.5% (19/29), the overall survival rate was 68.9% (20/29), and 8 patients with postoperative recurrence and metastasis. One patient with liver metastasis survived after reoperation, and one patient underwent postoperative chemoradiotherapy and died due to recurrence. There were 8 deaths during the follow-up, of which 7 died due to tumor recurrence and metastasis, and 1 died due to previous underlying diseases. Conclusion:Robot-assisted radical resection for hilar cholangiocarcinoma is feasible and effective.

Objective To construct a multivariable prediction model for assessing the risk of cardiopulmonary complication after thoracoscopic lobectomy in patients with stage Ⅰ-Ⅱ non-small cell lung cancer(NSCLC).Methods Clinical data of 600 patients with stage Ⅰ-Ⅱ NSCLC who underwent thoracoscopic lobectomy were retrospectively analyzed.Patients were divided into the complication group(84 cases)and the non-complication group(516 cases)based on the occurrence of postoperative complication within 7 days,including atelectasis,pulmonary embolism,respiratory failure,chylothorax,massive pleural effusion,hypoxemia and atrial fibrillation.Demographic characteristics,preoperative pulmonary function assessment,pathological features and perioperative indicators were collected.Multivariate Logistic regression analysis was used to identify independent risk factors influencing postoperative cardiopulmonary complication in patients and construct a predictive model.Internal validation was performed using the Bootstrap resampling method(1 000 iterations)to evaluate the discrimination,calibration and clinical decision-making value of the model.Results Multivariate Logistic regression analysis identified the following independent risk factors of postoperative cardiopulmonary complication,including age(OR=1.832,95%CI:1.537-2.183),history of chronic obstructive pulmonary disease(COPD)(OR=6.782,95%CI:2.685-17.130),Karnofsky performance status(KPS)score(OR=0.926,95%CI:0.888-0.965),the percentage of forced expiratory volume in the first second to the predictive value(FEV1%pred)(OR=0.906,95%CI:0.845-0.972),the percentage of diffusing capacity for carbon monoxide to the expected value(DLCO%pred)(OR=0.901,95%CI:0.832-0.975),intraoperative blood loss(OR=1.025,95%CI:1.014-1.036)and one-lung ventilation time(OR=1.057,95%CI:1.034-1.080).The area under the curve(AUC)of the combined diagnosis was 0.977(95%CI:0.965-0.989),with 96.4%sensitivity and 87.6%specificity.The Hosmer-Lemeshow test indicated excellent calibration(χ2=1.285,P=0.994).Decision curve analysis demonstrated significant clinical net benefit when the risk threshold probability ranged between 20%and 98%.Conclusion The multivariable prediction model for cardiopulmonary complication after thoracoscopic lobectomy in stage Ⅰ-Ⅱ NSCLC patients exhibits strong predictive performance.

Objective:To investigate the feasibility and therapeutic efficacy of robot-assisted radical resection for hilar cholangiocarcinoma.Methods:This is a retrospective case series study. The clinical data of 29 patients who underwent robot-assisted radical resection for hilar cholangiocarcinoma at the Department of Minimally Invasive Hepatic Surgery,the First Affiliated Hospital of Harbin Medical University from July 2021 to February 2025 were retrospectively collected. There were 16 males and 13 females, aged ( M(IQR)) 68.0 (10.0) years (range:36 to 78 years), and body mass index (24.0±2.9) kg/m 2 (range:17.5 to 29.1 kg/m 2). Bismuth-Corlette classification: 12 cases type Ⅰ, 4 cases type Ⅱ, 6 cases type Ⅲb, and 7 cases type Ⅳ. Preoperative CA19-9 was 161.7(320.9) U/ml (range:7.1 to 1 000.0 U/ml), and carcinoembryonic antigen was 2.8(2.1)μg/L (range:0.3 to 203.1 μg/L). Preoperative total bilirubin was 134.2 (348.9) μmol/L (range:10.4 to 557.9 μmol/L), direct bilirubin was 90.8 (264.1) μmol/L (range:2.5 to 418.7 μmol/L), ALT was 136.4 (134.8) U/L (range:13.0 to 569.9 U/L), AST was 122.2 (119.9) U/L (range:16.0 to 384.0 U/L), and albumin was (34.5±6.3) g/L (range:21.7 to 41.3 g/L). Comparison of quantitative data at different time points using paired t-test or Mann-Whitney U test. Cox univariate analysis was performed for the relevant variables, and Cox multivariate regression analysis was used to screen the independent prognostic factors of patients after robot-assisted radical resection for hilar cholangiocarcinoma. Results:All the 29 patients successfully underwent robot-assisted radical resection of hilar cholangiocarcinoma, and the R0 resection rate was 93.1% (27/29) without conversion to laparotomy. The operation time was 295.0 (87.5) minutes (range:195 to 590 minutes), the intraoperative blood loss was 100.0 (150.0) ml (range:20 to 1 000 ml), the intraoperative blood transfusion rate was 20.1% (6/29), the number of lymph nodes dissected was 10.0 (7.0) pieces (range: 6 to 18 pieces), the first postoperative deflatus time was 3.0 (1.0) days (range:2 to 4 days), The oral feeding time was 5.0 (1.0) days (range: 4 to 7 days), the drainage tube removal time was 8.0 (2.0) days (range: 6 to 26 days), and the postoperative hospital stay time was 10.0 (6.0) days (range:7 to 27 days). The incidence of complications above grade Ⅱ of the Clavien-Dindo complication grading system was 24.1% (7/29), including 3 cases of gastrointestinal bleeding with recurrent high fever, 1 case of delayed gastric emptying, 1 case of bile leakage, and 5 cases of hypoalbuminemia. The total bilirubin was 42.8 (66.8) μmol/L (range:6.8 to 195.9 μmol/L), direct bilirubin was 28.1 (38.5) μmol/L (range:4.3 to 88.6 μmol/L), ALT was 55.8 (56.0) U/L (range:9.9 to 207.1 U/L), AST was 33.9 (17.9) U/L (range:10.6 to 122.7 U/L), and albumin was (32.1±3.8) g/L (range:22.8 to 37.7 g/L), the levels of transaminase and bilirubin in the postoperative liver function indexes were significantly improved compared with those before operation, and the differences were statistically significant (all P<0.05). The mortality rate of patients without perioperative death was 3.4% (1/29) at 90 days after surgery. The results of Cox multivariate regression analysis showed that R0 resection was an independent prognostic factor for survival at 1 year after surgery ( P<0.05). The follow-up time was 15.0 (12.0) months (range:6 to 30 months), 1 of the 29 patients died of intra-abdominal infection 1 week after discharge, and the remaining 28 patients were completely followed up, of which 20 patients had no recurrence and metastasis during the follow-up period, and the tumor-free survival was 15.0 (12.0) months (range:6 to 30 months), the tumor-free survival rate was 65.5% (19/29), the overall survival rate was 68.9% (20/29), and 8 patients with postoperative recurrence and metastasis. One patient with liver metastasis survived after reoperation, and one patient underwent postoperative chemoradiotherapy and died due to recurrence. There were 8 deaths during the follow-up, of which 7 died due to tumor recurrence and metastasis, and 1 died due to previous underlying diseases. Conclusion:Robot-assisted radical resection for hilar cholangiocarcinoma is feasible and effective.

BackgroundSibling relationships play a critical role in shaping anxiety symptoms in firstborn children. Anxiety symptoms often originate in early childhood and can persist into adolescence and adulthood. However， there is insufficient research on anxiety symptoms in preschool children， especially firstborn preschool children. ObjectiveTo explore the influencing factors of anxiety symptoms among firstborn preschool children， so as to provide references for the intervention of anxiety symptom for children in families with multiple children. MethodsFrom October to December 2021， a total of 8 449 children from 234 kindergartens in Longhua District of Shenzhen were included using a cluster sampling method. Sibling Inventory of Behavior （SIB） and Spence Preschool Anxiety Scale （SPAS） were used to investigate. Logistic regression analysis was used to identify influencing factors of anxiety symptoms in firstborn preschool children. ResultsA total of 8 419 （99.64%） valid questionnaires were collected. Anxiety symptoms were detected in 344（4.09%） firstborn preschool children. Statistically significant differences were observed between anxiety group and non-anxiety group in terms of household registration， monthly family income， maternal age， maternal education level， paternal education level， family living conditions and whether they are left-behind children （χ²/t=9.906， 33.490， 5.136， 13.485， 9.690， 17.332， 21.975， P<0.05 or 0.01）. Compared with non-anxiety group， children in the anxiety group scored higher on the SIB dimensions of rivalry， aggression and avoidance （t=165.322， 74.471， 286.419， P<0.01）， and lower on companionship， empathy and teaching （t=59.133， 42.417， 39.112， P<0.01）. Risk factors for anxiety symptoms in firstborn preschool children included left-behind children， as well as negative sibling relationships characterized by rivalry and avoidance （OR=1.195， 1.143， 1.260， P<0.05 or 0.01）. ConclusionFirstborn preschool children who are left-behind are more susceptible to anxiety symptoms. Negative sibling relationships， characterized by competition and avoidance， may also contribute to the emergence of anxiety symptoms in firstborn preschool children.

Objective To predict the potential geographic distribution of Oncomelania hupensis in Yunnan Province using random forest (RF) and maximum entropy (MaxEnt) models, so as to provide insights into O. hupensis surveillance and control in Yunnan Province. Methods The O. hupensis snail survey data in Yunnan Province from 2015 to 2016 were collected and converted into O. hupensis snail distribution site data. Data of 22 environmental variables in Yunnan Province were collected, including twelve climate variables (annual potential evapotranspiration, annual mean ground surface temperature, annual precipitation, annual mean air pressure, annual mean relative humidity, annual sunshine duration, annual mean air temperature, annual mean wind speed, ≥ 0 ℃ annual accumulated temperature, ≥ 10 ℃ annual accumulated temperature, aridity and index of moisture), eight geographical variables (normalized difference vegetation index, landform type, land use type, altitude, soil type, soil textureclay content, soil texture-sand content and soil texture-silt content) and two population and economic variables (gross domestic product and population). Variables were screened with Pearson correlation test and variance inflation factor (VIF) test. The RF and MaxEnt models and the ensemble model were created using the biomod2 package of the software R 4.2.1, and the potential distribution of O. hupensis snails after 2016 was predicted in Yunnan Province. The predictive effects of models were evaluated through cross-validation and independent tests, and the area under the receiver operating characteristic curve (AUC), true skill statistics (TSS) and Kappa statistics were used for model evaluation. In addition, the importance of environmental variables was analyzed, the contribution of environmental variables output by the models with AUC values of > 0.950 and TSS values of > 0.850 were selected for normalization processing, and the importance percentage of environmental variables was obtained to analyze the importance of environmental variables. Results Data of 148 O. hupensis snail distribution sites and 15 environmental variables were included in training sets of RF and MaxEnt models, and both RF and MaxEnt models had high predictive performance, with both mean AUC values of > 0.900 and all mean TSS values and Kappa values of > 0.800, and significant differences in the AUC (t = 19.862, P < 0.05), TSS (t = 10.140, P < 0.05) and Kappa values (t = 10.237, P < 0.05) between two models. The AUC, TSS and Kappa values of the ensemble model were 0.996, 0.954 and 0.920, respectively. Independent data verification showed that the AUC, TSS and Kappa values of the RF model and the ensemble model were all 1, which still showed high performance in unknown data modeling, and the MaxEnt model showed poor performance, with TSS and Kappa values of 0 for 24%(24/100) of the modeling results. The modeling results of 79 RF models, 38 MaxEnt models and their ensemble models with AUC values of > 0.950 and TSS values of > 0.850 were included in the evaluation of importance of environmental variables. The importance of annual sunshine duration (SSD) was 32.989%, 37.847% and 46.315% in the RF model, the MaxEnt model and their ensemble model, while the importance of annual mean relative humidity (RHU) was 30.947%, 15.921% and 28.121%, respectively. Important environment variables were concentrated in modeling results of the RF model, dispersed in modeling results of the MaxEnt model, and most concentrated in modeling results of the ensemble model. The potential distribution of O. hupensis snails after 2016 was predicted to be relatively concentrated in Yunnan Province by the RF model and relatively large by the MaxEnt model, and the distribution of O. hupensis snails predicted by the ensemble model was mostly the joint distribution of O. hupensis snails predicted by RF and MaxEnt models. Conclusions Both RF and MaxEnt models are effective to predict the potential distribution of O. hupensis snails in Yunnan Province, which facilitates targeted O. hupensis snail control.

Schistosomiasis was once hyper-endemic in Yunnan Province. Following concerted efforts for over 70 years, remarkable achievements have been made for schistosomiasis control in the province. In 2004, the Mid- and Long-term Plan for Schistosomiasis Prevention and Control in Yunnan Province was initiated in Yunnan Province, and the target for transmission control of schistosomiasis was achieved in the province in 2009. Following the subsequent implementation of the Outline for Key Projects in Integrated Schistosomiasis Control Program (2009—2015) and the 13th Five - year Plan for Schistosomiasis Control in Yunnan Province, no acute schistosomiasis had been identified in Yunnan Province for successive 12 years, and no local Schistosoma japonicum infections had been detected in humans, animals or Oncomelania hupensis snails for successive 6 years in the province by the end of 2020. The transmission of schistosomiasis was interrupted in Yunnan Province in 2020. This review summarizes the history of schistosomiasis, changes in schistosomiasis prevalence and progress of schistosomiasis control in Yunnan Province, and proposes the future priorities for schistosomiasis control in the province.