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 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 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.

Objective To construct an index system for assessment of schistosomiasis transmission risk following natural disasters such as rainstorms, floods, earthquakes, mudslides, and landslides, so as to provide insights into rapid identification of schistosomiasis transmission risk post-disasters and formulation of targeted schistosomiasis control strategies. Methods An initial framework for the index system for assessment of schistosomiasis transmission risk following natural disasters was drafted through literature review, brainstorming, and focus group discussions. Two rounds of expert correspondence consultations were conducted using the Delphi method to refine and finalize the system, and the degrees of expert activeness, authority and endorse ment, and consensus were evaluated. In addition, the weights of each index were calculated using the analytic hierarchy process. Results A total of 18 experts participated in the consultation. The expert positive coefficients were 100.00% and 94.44% for two rounds of consultations, with authority coefficients of 0.92 and 0.94, respectively. The coefficients of coordination on the index importance, rationality and operability were 0.209, 0.185, 0.222 and 0.407, 0.214, 0.257 for two rounds of consultations, respectively, and all consistency tests were statistically significant (χ² = 246.771 to 505.278, all P values < 0.001). Following two rounds of expert consultations, an index system consisting of 6 first-level indicators, 15 second-level indicators, and 49 third-level indicators was ultimately constructed. In terms of first-level indicators, “disaster situation”, “previous epidemics”, “healthcare guarantee”, “response capacity” and “emergency recovery” had the highest weights, each at 18.18%. Regarding second-level indicators, “Schistosoma japonicum infections in animals”, “S. japonicum infections in snails” and “medical treatment” had the highest weights, each at 7.35%. In terms of third-level indicators, ten items had the highest weights, including “identification of schistosomiasis cases”, “detection of S. japonicum infections in wild feces”, “detection of S. japonicum infections in snails”, “reserves of schistosomiasis diagnostic/testing reagents and consumables”, “reserves of chemotherapy agents for human and animal schistosomiasis”, “reserves of cercariacides”, “periodical surveillance on schistosomiasis”, “identification of schistosomiasis transmission risk and timely response”, “normal provision of diagnosis and treatment services” and “post-disaster schistosomiasis surveillance”, each at 2.40%. Conclusion A scientific, systematic, and practical index system has been constructed for assessment of schistosomiasis transmission risk following natural disasters, which may provide insights into rapid post-disaster identification of schistosomiasis transmission risk, formulation of targeted schistosomiasis control strategies and optimization of resource allocation.

Medical parasitology, as a course bridging basic medical sciences and clinical medicine, has an important disciplinary value in the medical education system. This study investigated the composition of parasitology teachers from multiple medical colleges and universities across China. The results showed that there was a significant difference in the proportion of teachers with clinical medicine background knowledge, and there was common dilemma that there were insufficient clinical medicine knowledge reserves among teachers in some medical colleges and universities, who encountered severe teaching challenges. Based on this issue, this study constructed a basic-clinical medicine collaborative problem-based learning (PBL) teaching model. This model integrated theoretical teaching, case analyses, and experimental operations, and combined transdisciplinary team building and multidimensional teacher training, which significantly improved the clinical teaching capability among parasitology teachers, and effectively compensated the impact of insufficient clinical medicine knowledge reserves on teaching. Following teaching reform, students' scores significantly improved, and their case analysis capability enhanced. This study provides a practical path to address the shortage of clinical medicine background knowledge among parasitology teachers, which facilitates the progress of educational reform of medical parasitology and improvement of teaching quality.

Background Co-exposure to noise and microwave radiation occurs frequently. The central nervous system has been identified as a sensitive target organ for both noise and microwave exposure individually, and the underlying mechanisms remain poorly understood. The specific biological effects resulting from co-exposure to these two factors have yet to be fully elucidated. Objective To clarify the effects of co-exposure to noise and microwave on neurobehavior and hippocampal tissue structure, and to explore the underlying mechanism through the assessment of serum cytokines. Methods C57BL/6N mice were selected and randomly assigned to a blank control group, a noise group, a microwave group, and a combined noise & microwave exposure group. To establish the exposure models, the noise group was subjected to broadband noise at 100 dB for 2 h, while the microwave group received radiation at a central frequency of 9.375 GHz with an average power density of 12 mW·cm⁻² and a specific absorption rate of 2.58 W·kg⁻¹ for 15 min. Open field and tail suspension tests assessed anxiety-like emotional behaviour; novel object recognition and Y-maze tests evaluated cognitive function. Histological changes in hippocampal tissue were examined using haematoxylin and eosin (HE) staining, and Nissl staining under light microscopy. Serum cytokine levels were measured using radioimmunoassay and enzyme-linked immunosorbent assay (ELISA). Results After 3 d of exposure, the noise, microwave, and combined exposure groups showed significant reductions in exploration frequency, duration, and distance within the central zone of the open field test compared to the control group (P < 0.01); the combined exposure group exhibited increased ratios of peripheral-to-central exploration time and distance (P < 0.05). After 7 d of exposure, compared with the control group, the noise group maintained a decrease in central zone exploration time (P < 0.01), while the combined exposure group showed persistent decline across all central zone metrics (P < 0.05) and elevated peripheral-to-central ratios (P < 0.05); compared to the microwave group, the combined exposure group showed significant less time in the central zone (P < 0.05) and higher peripheral-to-central ratios (P < 0.05). Regarding behaviour and cognition, compared with the control group, the combined exposure group showed increased immobility time in the tail suspension test after 3 d of exposure (P < 0.01). At this interval, all exposure groups demonstrated reduced frequency and duration of novel object recognition (P < 0.05), with the combined exposure group showing a marked decrease in novel arm exploration time (P < 0.01). After 7 d of exposure, compared with the control group, the noise group showed reduced novel object recognition frequency (P < 0.05), and both the noise and microwave groups exhibited decreased novel arm exploration time (P < 0.05). Pathological alterations including an increased number of hyperchromatic nuclei and depleted Nissl bodies were observed in the CA3 and DG regions across all exposure groups with the most severe lesions observed in the combined exposure group. Serum levels of central nervous system-specific protein β (S-100β), glial fibrillary acidic protein (GFAP), and corticosterone (CORT) were significantly elevated in all exposure groups compared with the control group (P < 0.05). Aquaporin-4 (AQP4) levels increased in the combined exposure group (P < 0.05), while CXC chemokine ligand 10 (CXCL10) levels rose in both the noise and combined groups compared with the control group (P < 0.05). Specifically, S-100β and CXCL10 levels in the combined exposure group were higher than those in the microwave group (P < 0.05); moreover, levels of S-100β, GFAP, CORT, AQP4, and CXCL10 in the combined exposure group were significantly higher than those in the noise group (P < 0.05). Conclusion Combined exposure to noise and microwave radiation induces pathological changes in the hippocampus of mice, increases levels of serum stress hormones and neuro-specific biomarkers. These impairments are more severe than those observed following single-factor exposure. The underlaying mechanism may be related to systemic stress response, neuronal damage, astrocyte activation, and changes in blood-brain barrier permeability, leading to emotional behavioral abnormalities and cognitive decline.

Background In 2021, chronic obstructive pulmonary disease (COPD) emerged as the forth leading cause of death in the world. However, the impact of air pollutants on COPD is still inconsistent across current studies. Objective To analyze the relationship between ambient sulfur dioxide (SO₂) exposure and hospital admissions for COPD in Lanzhou, and to examine the modified effects of SO₂ across different genders, age groups, and seasons. Methods A total of 20718 hospital admission records for COPD were collected from Lanzhou between 2016 and 2020, alongside synchronized data on SO₂ concentration and meteorological variables. A generalized additive model integrated with a distributed lag nonlinear model was used to assess the relationship and the lag effects between SO₂ exposure and COPD admissions. Potential confounders, including meteorological factors, day-of-the-week effect, and holidays, were controlled for, with a maximum lag period set at 7 d. Two-pollutant models and sensitivity analyses were conducted to ensure robustness. Results are expressed as relative risks (RR) with 95% confidence intervals (95%CI). Results During the study period, the average daily number of COPD admissions was 11.34±7.03. The average mass concentration of SO₂ was (23.72±12.35) μg·m⁻³. SO₂ exposure was positively correlated with COPD admissions; specifically, each 10 μg·m⁻³ increase in SO₂ concentration was associated with an RR of 1.440 (95%CI: 1.317, 1.573). Stratified analyses found that at a cumulative lag of 7 d, the RRs were higher among females, individuals aged ≥65 years, and during the cold season. Conversely, no significant increase in COPD admission risk related to SO₂ was observed during the warm season. For every 10 μg·m⁻³ increase in SO₂ concentration, the RR was 1.483 (95%CI: 1.311, 1.678) for females, 1.441 (95%CI: 1.311, 1.583) for those aged > 65 years， and 1.595 (95%CI: 1.313, 1.937) during the cold season. Conclusion Short-term exposure to ambient SO₂ exposure in Lanzhou is associated with an increased risk of COPD admission. The association is more pronounced among females, the elderly (aged> 65 years) and during the cold season.

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.

Background Sunshine duration is closely associated with population mental health and emotional states, although its relationship with mental and behavioral disorders (MBD) remains insufficiently studied. Objective To analyze the effect of sunshine duration on hospital admissions for MBD in Zigong City, Sichuan Province. Methods Hospital admission records for MBD from 10 medical institutions, meteorological data, and ambient air pollutant concentrations were collected in Zigong City from January 1, 2019 to December 31, 2024. A distributed lag non-linear model (DLNM) was employed to calculate single-day and cumulative lag effects of different sunshine duration exposures—0 h (P₀, P₅, P₂₅), 6 h (P₇₅), and 10.4 h (P₉₅)—on hospitalization risks for MBD, stratified by diagnostic category, sex, and age groups. Results This study analyzed 36597 hospital admissions for MBD in Zigong City. The exposure-response curve for sunshine duration at lag 0 d exhibited an inverted N-shaped pattern, with minimum relative risk (RR) of hospitalizations observed at 14 h sunshine duration. The immediate effects of different levels of sunshine duration on the hospitalization risk for MBD, schizophrenia, schizotypal and delusional disorders (SSDD), and neurotic, stress-related and somatoform disorders (NSSD) were significant, with the cumulative effect peaking at lags 12–13 days. For extremely short sunshine duration 0 h, the maximum relative risk (RR) values (with 95% confidence intervals, 95%CI) of hospitalization risk were 1.209 (1.001, 1.459) at lag 1 d for SSDD patients and 1.398 (1.133, 1.724) at lag 0 d for NSSD patients, respectively. The cumulative effect on hospitalization risk induced by short sunshine duration 6 h was higher than that induced by extremely short sunshine duration 0 h and long sunshine duration 10.4 h. Among populations stratified by gender and age, male patients and those aged over 70 years had higher hospitalization risks. For short sunshine duration 6 h, the cumulative RR (CRR) value (95%CI) of hospitalization risk for male MBD patients was 3.373 (1.362, 8.358) at cumulative lag 0–12 d, which was higher than that for female patients [2.943 (1.471, 5.889) at cumulative lag 0–13 d]. The population aged over 70 years was more sensitive to single-day and cumulative lag effects on hospitalization risk than other age groups. For short sunshine duration 6 h, the RR and CRR values (95%CI) of single-day and cumulative lag effects on hospitalization risk were 1.336 (1.002, 1.783) at lag 0 day and 5.085 (1.680, 15.392) at cumulative lag 0–12 d, respectively. At lag 13 d, long sunshine duration 10.4 h exposure was associated with depression hospitalizations in patients <20 years with RR=1.127 (95%CI: 1.002, 1.270). Conclusion When the daily sunshine duration is less than 14 h, it will increase the hospitalization risk of MBD in the population to varying degrees, with both single-day and cumulative lag effects. It is recommended that the public, especially vulnerable groups such as people with MBD aged 70 years and above and those with depression under 20 years, strengthen health protection to reduce the health risks caused by changes in sunshine.