Objective To evaluate the effectiveness of schistosomiasis surveillance in Jiangsu Province during the stage moving from transmission control to transmission interruption, and to analyze the current risk and challenges, so as to provide the evidence for achieving the target of schistosomiasis elimination. Methods Schistosomiasis surveillance data were collected from Jiangsu Province from 2012 to 2024, and the endemic areas, Schistosoma japonicum infections in humans and livestock, Oncomelania hupensis snail distribution and implementation of integrated interventions were descriptively analyzed. In addition, the trends in areas with snails, seroprevalence of human S. japonicum infections and numbers of advanced schistosomiasis cases were assessed using a Joinpoint regression model. Results The endemic areas of schistosomiasis continued to shrink in Jiangsu Province from 2012 to 2024, with the number of schistosomiasis-eliminated counties (cities, districts) increasing from 53 (75.71%) to 63 (96.92%), and interruption of schistosomiasis transmission was achieved across the province. A total of 4 600 300 person-times were tested for serum antibodies against S. japonicum, with 28 719 person-times positive detected; and 616 500 person-times were tested S. japonicum infections among local residents in Jiangsu Province from 2012 to 2024, with only 3 egg-positives detected, and no egg-positives found since 2017. A total of 187 600 herd-times were tested for schistosomiasis in livestock, and no S. japonicum infections were found. O. hupensis snail survey was performed covering 1 018 408.97 hm², and a total of 35 556.35 hm² was found with snail-infested habitats, including 174.40 hm2 of emerging snail-infested habitats. A total of 1 102 800 O. hupensis snails were identified for S. japonicum infections, and no infections were found. The areas of snail-infested habitats appeared a tendency towards a rise in Jiangsu Province from 2019 to 2023 (APC = 23.67%, P < 0.05), and the actual areas of snail-infested habitats appeared a tendency towards a decline from 2012 to 2015 (APC = −22.77%, P < 0.05), and towards a rise from 2015 to 2023 (APC = 9.76%, P < 0.01). The seroprevalence of anti-S. japonicum antibodies appeared a tendency towards a decline among residents in Jiangsu Province from 2017 to 2023 (APC = −14.92%, P < 0.01). In addition, the number of newly diagnosed advanced schistosomiasis cases appeared a tendency towards a decline from 2012 to 2024 (APC = −12.02%, P < 0.01), and the numbers of advanced schistosomiasis patients requiring treatment showed a tendency towards a decline from 2012 to 2021 (APC = −10.56%, P < 0.01) and from 2021 to 2023 (APC = −20.06%, P < 0.01). Conclusions Great progresses had been achieved in schistosomiasis control in Jiangsu Province following transmission control, and transmission interruption had been achieved; however, there are still snail-infested habitats. High-intensity surveillance and integrated control are required to be maintained to advance the achievement of the target of schistosomiasis elimination in Jiangsu Province.

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 investigate the species of sandflies and the prevalence of Leishmania infections in sandflies from selected areas of northern and northwestern China, so as to provide insights into identification of leishmaniasis vectors and assessment of epidemiological trends of leishmaniasis in China. Methods Sandfly samples were collected from Mentougou District of Beijing Municipality, Xiangning County in Linfen City of Shanxi Province, Ejin Banner in Alxa League of Inner Mongolia Autonomous Region, and Payzawat County of Karamay District of Karamay City, Gaochang District of Turpan City in Xinjiang Uygur Autonomous Region from July 2023 to July 2024. Approximately 100 intact female sandfly samples were randomly selected from each site and the species of sandflies was identified according to morphological characteristics and molecular assays. Female sandflies originating from the same habitat were grouped into pools of 10 individuals. Leishmania infection was detected using polymerase chain reaction (PCR) assay targeting the internal transcribed spacer 1 (ITS-1) gene, and the prevalence of Leishmania infection was calculated in sandflies from different sampling sites using the minimum infection rate (MIR) method. In addition, positive amplicons were sequenced and subjected to phylogenetic analysis. Results A total of 6 155 sandflies were collected from different environments at sampling sites across the six aforementioned regions from July 2023 to July 2024. Phlebotomus chinensis (96.00%) was the dominant sandfly species in Mentougou District, Beijing Municipality, with a small proportion of Ph. sergenti (4.00%), and only Ph. chinensis was found in Xiangning County, Linfen City, Shanxi Province. Ph. wui was the only sandfly species detected in Ejin Banner, Alxa League, Inner Mongolia Autonomous Region, and Payzawat County, Kashgar City, Xinjiang Uygur Autonomous Region, and Ph. caucasicus (97.70%) was the dominant sandfly species in Karamay District, Karamay City, Xinjiang Uygur Autonomous Region, with a small proportion of Ph. wui (2.30%), while Ph. alexandri was the only species in Gaochang District, Turpan City, Xinjiang Uygur Autonomous Region. A total of 40, 60, 34, 18, 18, and 22 pools of sandfly samples were tested from Mentougou District in Beijing Municipality, Xiangning County in Linfen City of Shanxi Province, Ejin Banner in Alxa League of Inner Mongolia Autonomous Region, Payzawat County in Kashgar City, Karamay District in Karamay City, and Gaochang District in Turpan City of Xinjiang Uygur Autonomous Region, respectively. L. infantum was detected in Ph. chinensis samples from Mentougou District in Beijing Municipality, and Xiangning County of Linfen City in Shanxi Province, with MIR of 0.25% to 1.00%, and L. donovani was detected in Ph. wui from Ejin Banner in Alxa League of Inner Mongolia Autonomous Region, and Payzawat County in Kashgar City of Xinjiang Uygur Autonomous Region, with MIR of 0.56% to 0.88%; however, no Leishmania infection was detected in Ph. caucasicus from Karamay District in Karamay City or Ph. alexandri from Gaochang District in Turpan City of Xinjiang Uygur Autonomous Region. Phylogenetic analysis showed that the Leishmania ITS-1 gene sequences obtained from Mentougou District in Beijing Municipality and Xiangning County in Linfen City of Shanxi Province were clustered into the same clade with the reference sequences of L. infantum ITS-1 gene, while the Leishmania ITS-1 gene sequences obtained from Ejin Banner in Alxa League of Inner Mongolia Autonomous Region and Payzawat County in Kashgar City of Xinjiang Uygur Autonomous Region were clustered into the same clade with the reference sequences of L. donovani ITS-1 gene. Conclusions There are variations in sandfly species in selected areas of northern and northwestern China, and variations in the species of Leishmania infecting sandflies. Improved surveillance of sandfly vectors and targeted control strategies with adaptations to geographical features and leishmaniasis vectors are recommended.

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 To investigate the current distribution and predict the future suitable habitats of Haemaphysalis longicornis in Nanjing City, so as to provide insights into control and early warning of ticks and management of tick-borne diseases in Nanjing City. Methods The electronic map of Nanjing City was obtained from the National Platform for Common GeoSpatial Information Services. The distribution of H. longicornis and the longitude and latitude of distribution points from 2022 to 2024 were obtained from centers for disease control and prevention across each district in Nanjing City. Climatic and environmental variable data in Nanjing City were captured from the Worldclim database. Initially, 19 bioclimatic variables in this database were selected, including annual mean temperature, mean diurnal range, isothermality, temperature seasonality, maximum temperature of the warmest month, minimum temperature of the warmest 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. The elevation and normalized difference vegetation index were obtained from Data Sharing Platform of the Center for Resources and Environmental Sciences, Chinese Academy of Sciences. Then, the distribution points of H. longicornis, elevation, vegetation index and 19 bioclimatic variables were loaded into the software MaxEnt 3.4.4 to evaluate and screen out the variables with a contribution rate of 1% and higher. ArcGIS 10.8.1 software was used to extract the elevation, vegetation index and 19 bioclimatic variables of the distribution points of H. longicornis for a correlation analysis. If the absolute value of the correlation coefficient was 0.8 and higher, the variable with the higher contribution was retained. The 2050 dataset of the BCCCSM2-MR atmospheric circulation model in the coupled model intercomparison project phase 6 (CMIP6) were obtained from the Worldclim database as climate data for 2050. Screened H. longicornis species data and environmental and climate data were loaded into the maximum entropy (MaxEnt) model with the software MaxEnt 3.4.4 for training and validation, and then, all data generated from the model were imported into the software ArcGIS 10.8.1 to generate raster data and yield the map pertaining to the distribution of H. longicornis risk in Nanjing City. The accuracy of the model was evaluated with a receiver operating characteristic (ROC) curve, and the predictive effect of the model was assessed with area under the ROC curve (AUC). The suitable habitats of H. longicornis were classified in Nanjing City with the software ArcGIS 10.8.1, and the areas of distribution of suitable habitats in various categories were recorded to create the map of current H. longicornis suitable habitats classification in Nanjing City. The climatic and geographic information data in 2050 were employed as future environmental and climatic factors, and current environmental and climatic factors and current H. longicornis distribution data were additionally used to predict the future suitable habitats of H. longicornis in Nanjing City. In addition, the contributions of environmental and climatic factors to distribution of suitable habitats of H. longicornis was evaluated with the Jackknife method in Nanjing City. Results A total of 10 environmental and climatic variables were screened for analysis of the suitability of H. longicornis in Nanjing City based on correlation analyses and contributions of the MaxEnt model, including annual mean temperature, precipitation of the warmest quarter, vegetation index, precipitation of the wettest month, temperature annual range, annual precipitation, mean temperature of the warmest quarter, elevation, mean temperature of the wettest quarter, and maximum temperature of the warmest month, and annual mean temperature (34.8%), precipitation of the warmest quarter (17.3%), vegetation index (13.1%), and precipitation of the wettest month (10.8%) contributed relatively highly to the distribution of suitable habitats of H. longicornis in Nanjing City. The mean AUC of the ROC curve was 0.810 ± 0.055 for 10 repeated modeling results of the MaxEnt model, indicating high predictive performance of the model. The potential distribution areas of H. longicornis were predicted to be mainly located in Luhe District, Pukou District, Jiangning District, Lishui District, and Gaochun District in Nanjing City with the MaxEnt model. Under current climatic conditions, the area of potential suitable habitats of H. longicornis was 4 182.42 km² in Nanjing City, including 1 252.94 km² highly suitable habitats, which accounted for 19.00% of the total area of Nanjing City. Under the climate scenario in 2050, the area of potential suitable habitats of H. longicornis was projected to increase to 5 467.58 km² in Nanjing City, accounting for 82.95% of the total area of the city, and these habitats were mainly concentrated in Luhe District, Pukou District, Jiangning District, Lishui District, and Gaochun District. The areas of suitable habitats of H. longicornis at various categories were predicted to vary greatly in 2050, and the area of highly suitable habitats of H. longicornis was projected to increase to 2 378.82 km², accounting for 36.08% of the total area of Nanjing City. Based on jackknife tests and contributions of environmental and climatic variables, 6 dominant environmental and climatic factors were screened, including annual mean temperature (34.8% contribution), precipitation of the warmest quarter (17.3% contribution), vegetation index (13.1% contribution), precipitation of the wettest month (10.8% contribution), temperature annual range (5.4% contribution), and mean temperature of the warmest quarter (5.0% contribution), with cumulative contributions of 86.4%. Conclusion The distribution of H. longicornis is strongly associated with vegetation, temperature and precipitation in Nanjing City. Future climate change may lead to an expansion of the distribution area of H. longicornis in Nanjing City.

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.

The frequent extreme climatic events post multifaceted impacts on the distribution of Oncomelania hupensis, the intermediate host of Schistosoma japonicum in the context of global climate change. This article systematically reviews the effects of four types of extreme climatic events, including floods, droughts, heat waves, and cold waves, on the survival, reproduction, and distribution of Oncomelania hupensis. Floods may expand suitable snail habitats, and increase both emerging and re-emerging snail habitats; however, the impact of floods on O. hupensis density is characterized by a lag effect of a decline followed by a rise. Droughts may cause fragmentation of suitable O. hupensis snail habitats, reduced O. hupensis snail egg production rates, and increased O. hupensis snail mortality, and heat waves may cause an increase in O. hupensis snail mortality, a reduction in numbers of O. hupensis snail populations and shrinking of O. hupensis snail distribution, while cold waves may cause a reduction in O. hupensis snail density and a rise in O. hupensis snail mortality. Extreme climate events pose both shortand long-term effects on the distribution of O. hupensis. Intensified surveillance of O. hupensis snails is required in high-risk environments.

Background Aluminum (Al) is a lightweight metal that is widely present in the environment and the human body. It has been documented to cause various adverse health effects including the suppression of humoral immunity. Objective To investigate the role of oxidative stress in Al-induced humoral immunity suppression and to evaluate the possible protective effects of iron supplementation on this process. Methods Adult C57BL/6J mice were exposed to Al at concentrations of 0, 200, or 800 μg·mL⁻¹ via drinking water for three consecutive months. The expression of major histocompatibility complex class Ⅱ (I-A), proliferating cell markr-67 (Ki-67), and 2',7'-dichlorodihydrofluorescein diacetate (DCFH-DA) in splenic B cells was evaluated through flow cytometry. Splenic B cells from the mice treated with 800 μg·mL⁻¹ Al or the control were sorted and treated in vitro with glutathione (GSH), N-Acetyl-L-cysteine (NAC), or a control vehicle. After 24 h, the expression of I-A was evaluated; and the hydroxyl radical (·OH)-generating potential, ·OH production, malondialdehyde (MDA) production, and iron content were assessed using commercial kits. Sixteen mice treated with 800 μg·mL⁻¹ Al received an intravenous injection of either a ferric chloride solution containing 0.3 g·L⁻¹ iron or a 0.9% sodium chloride solution, while eight control mice received 0.9% sodium chloride solution; the injection volume was 0.1 mL per mouse. Two and a half days after injection, I-A and Ki-67 expressions, ·OH-generating potential, ·OH production, and MDA production in splenic B cells were measured; and the concentrations of serum immunoglobulin (Ig) M and IgG were measured through (enzyme-linked immunosorbent assay) ELISA. The splenic B cells sorted from untreated mice were exposed to 0, 12.5, 25, or 50 μg·L⁻¹ Al in vitro. The splenic B cells treated with 50 μg·L⁻¹ Al and the splenic B cells sorted from 800 μg·mL⁻¹ Al-treated mice were additionally treated with GSH and NAC in vitro. The iron supplementation groups, which included the 50 μg·L⁻¹ Al-treated group and splenic B cells sorted from 800 μg·mL⁻¹ Al-treated mice, were treated with a culture medium containing 30 μmol·L⁻¹ iron in vitro. I-A and Ki-67 expressions, ·OH-generating potential, ·OH production, and MDA production in B cells were detected after a 24-h treatment period. Results In the in vivo mouse model, exposure to 800 μg·mL⁻¹ Al significantly inhibited the I-A and Ki-67 expressions (P<0.05), increased DCFH-DA expression and ·OH-generating potential (P<0.05, P<0.01), decreased iron content (P<0.01) and ·OH and MDA production (P<0.01, P<0.001) of splenic B cells, as well as serum IgM and IgG concentrations (P<0.05, P<0.01) in the mice. Exposure to 200 μg·mL⁻¹ Al showed a tendency to decrease the I-A and Ki-67 expressions, and to increase the DCFH-DA expression in splenic B cells, but these differences were not significant. In the in vitro splenic B-cell model, Al (12.5, 25, and 50 μg·L⁻¹) inhibited I-A and Ki-67 expressions (P<0.05, P<0.01) across all concentrations; 50 μg·L⁻¹ Al increased the ·OH-generating potential (P<0.05), and decreased ·OH and MDA production (P<0.01, P<0.05) in B cells. Treatment with GSH and NAC further suppressed I-A expression (P<0.05) in B cells. Iron supplementation increased the ·OH and MDA production (P<0.05), restored I-A and Ki-67 expressions (P<0.05, P<0.01) in B cells, and elevated the serum IgM and IgG concentrations (P<0.05) in Al-treated mice. Conclusion Al suppresses humoral immunity and ·OH production in B cells. The underlying mechanism may involve the decreased iron content and the subsequent retardation of the Fenton reaction in B cells. Supplementing with iron can restore the Fenton reaction in B cells and potentially reverse Al-induced impairment of humoral immunity.

Background The prevalence of chronic diseases among the Chinese occupational population is rising steadily, with hypertension and diabetes becoming important health concerns. Occupational health examinations (OHE) provide stable population coverage, standardized protocols, and fixed follow-up intervals, offering a robust data foundation for risk assessment. However, most existing hypertension prediction studies rely on cross-sectional data and mainly focus on clinic onset, failing to capture the dynamic progression and cumulation of individual risk. Objective To construct a machine learning-based risk prediction model for blood pressure abnormality in occupational populations, providing a reference for health risk stratification and targeted health interventions. Methods Longitudinal data from 2020 to 2023 were extracted from the occupational health examination database of an institution in Shanghai. After excluding individuals with hypertension in any of the first three years, 3710 workers who participated 4 consecutive years of OHE were included. Six categories of blood pressure-related indicators from 2020 to 2022 were selected, including basic information, lifestyles, occupational exposures, cardiovascular indicators, routine blood tests, and biochemical markers, comprising 30 variables in total. For feature engineering, continuous variables were processed to reflect three-year dynamic characteristics: level (mean), variability (variance), and temporal trend (slope). Categorical variables were incorporated using 2022 values. Least absolute shrinkage and selection operator regression was applied to identify variables for inclusion in the machine learning models. Five models (decision tree, logistic regression, random forest, support vector machine, and XGBoost) were employed to predict blood pressure abnormality risk in 2023. And then sensitive analysis was conducted. The optimal model was selected based on the area under the curve (AUC) of receiver operating characteristic. Feature importance and SHAP analyses were applied to interpret the final model. Results Feature engineering transformed 17 continuous variables into 51 secondary variables, combined with 13 baseline variables, resulted in 24 variables screened out by LASSO regression. The logistic regression model achieved the best performance. Feature importance analysis indicated that the dynamic trajectory of blood pressure played a central role in the model, complemented by other biochemical and lifestyle indicators. SHAP analysis further showed that the model’s ability to not only identify risks of high-normal blood pressure and hypertension but also quantify the specific contribution of individual examination indicators to predicted risk, supporting risk-stratified population management in occupational health management. Conclusion Based on an existing occupational health examination database and routinely collected data, this study has developed a risk prediction model for blood pressure abnormalities in occupational populations. By incorporating dynamic change characteristics from longitudinal examinations, the study demonstrates the feasibility of using routine health data for early risk identification and stratified management. However, its applicability is limited by the availability of repeated examination data and the focus on specific occupational groups, requiring further validation across broader populations and diverse clinic scenarios.

Background Meteorological factors are among the key extrinsic triggers for the onset and exacerbation of cardiovascular and cerebrovascular diseases (CVD). Against the backdrop of sustained global warming, elucidating the impact of ambient temperature and atmospheric pressure on CVD, especially on pre-hospital CVD emergent events, has become imperative for evidence-based prevention and emergency preparedness. Objective To quantify the temporal trends of daily mean temperature and atmospheric pressure and their associations with pre-hospital CVD emergent events in Yantai, and to explore effect modification by demographic subgroups and geographic areas, thereby providing an empirical basis for the rational allocation of emergency medical resources. Methods Pre-hospital CVD emergency data from January 1, 2016 to December 31, 2022 were selected from the Yantai 120 Emergency Medical Command System. Synchronous meteorological factors and environmental pollutant data were obtained from the websites of the National Oceanic and Atmospheric Administration and the National Centers for Environmental Information of the United States. Time-series analysis combined with distributed lag non-linear model was used to analyze the association between daily temperature, atmospheric pressure, and pre-hospital CVD emergencies. Average annual percentage changes (AAPC) were calculated using Joinpoint (version 5.2.0.0) to reflect temporal trends. Spearman correlation analysis was employed to screen variables with low collinearity for inclusion in the multi-pollutant adjusted models. Results From 2016 to 2022, a total of 268275 pre-hospital CVD emergency cases were recorded in Yantai City, showing an increasing annual trend (AAPC=7.16%, P=0.002). Emergency volume, daily temperature, atmospheric pressure, and PM₁₀ displayed marked seasonal patterns, with summer bottoms and winter peaks for emergencies, daily temperature higher in summer and lower in winter, and the inverse for atmospheric pressure. Correlation analysis showed that both daily temperature and atmospheric pressure were positively correlated with O₃, with correlation coefficients of 0.745 and 0.723, respectively; negatively correlated with PM_2.5, with correlation coefficients of −0.246 and −0.235, respectively; and negatively correlated with PM₁₀, with correlation coefficients of −0.241 and −0.229, respectively (P < 0.05). A U-shaped correlation was revealed between daily temperature and pre-hospital CVD emergency risk. Using 24 °C as the reference, risk was elevated at both low and high daily temperatures, with the strongest cold effect at –11 ℃ (RR=1.53, 95%CI: 1.35, 1.73) and strongest heat effect at 30 ℃ (RR=1.06, 95%CI: 1.01, 1.11); cold effects were prolonged (lag0–14d) compared to heat effects. Atmospheric pressure exhibited a positive association with emergencies; using a median of 1012.2 hPa as the reference, the maximum risk occurred at 1037 hPa (RR=1.25, 95%CI: 1.03, 1.52), with a 0–5 d lag. In the subgroup analysis, the impact of low temperature was more pronounced among females (RR_{−9 ℃}=1.66), people over 60 years (RR_{−9 ℃}=1.91), and rural residents (RR_{−9 ℃}=1.76). High atmospheric pressure significantly increased the risk of pre-hospital CVD emergency among females (RR_{1035 hPa}=1.54), people under 18 years (RR_{1035 hPa}=3.62), and rural residents (RR_{1035 hPa}=1.46). The sensitivity analysis results indicated that, after excluding the impact of pandemic and the degrees of freedom for time variations, the effects of daily average temperature and atmospheric pressure on pre-hospital CVD emergency volume were consistent with the primary analysis findings. Conclusion Both non-optimal temperature and atmospheric pressure significantly increase the pre-hospital CVD emergency volume; therefore, emergency medical services should develop targeted response protocols, optimize resource allocation, and intensify public health education to mitigate the pre-hospital burden attributable to temperature and pressure fluctuations.