Objective: To determine the causative agent, clinical manifestations and risk factors for infection during a September 2023 outbreak of acute haemorrhagic conjunctivitis (AHC) in Pak Touch village, Ratanakiri province, Cambodia. Methods: A retrospective case-control study was conducted. Cases were age-matched to controls (1:1), who were randomly selected from the village population. Twenty-one conjunctival samples were analysed using real-time reverse transcription–polymerase chain reaction (RT–PCR). RNA sequencing was additionally performed to identify the causative agent of the outbreak. Logistic regression models were used to identify significant risk factors. Results: A total of 73 cases and 73 controls were included in the analysis. Cases had a median age of 20 years (range: 1–70, mean and standard deviation: 27.7 ± 20.0), and 46.6% (34/73) were male. The overall attack rate was 12.3% (73 cases/594 residents). Clinical presentations included conjunctival hyperaemia (100%), subconjunctival haemorrhage (82.2%, 60), pain and discharge (64.4%, 47 each), eyelid swelling (57.5%, 42) and tearing (54.8%, 40). RT–PCR identified enterovirus in 52.4% (11/21) of conjunctival swabs, with RNA sequencing confirming the coxsackievirus A24 variant as the causative agent in five swabs. Statistical analysis identified significant risk factors, including physical contact with patients with acute haemorrhagic conjunctivitis (adjusted odds ratio [aOR]: 4.42, 95% confidence interval [CI]: 1.90–10.10), frequent eye rubbing (aOR: 4.56, 95% CI: 2.00–10.37) and poor hand hygiene (aOR: 3.70, 95% CI: 1.64–8.43). Discussion: The outbreak of acute haemorrhagic conjunctivitis in Pak Touch village was primarily caused by coxsackievirus A24. Significant risk factors included physical contact with infected individuals, frequent eye rubbing and poor hand hygiene. Effective hygiene measures are crucial to prevent the spread of AHC.

ObjectiveTo investigate the effects of exogenous nitric oxide （NO） on the accumulation and quality formation mechanism of flavonoids in Scutellariae Radix. MethodsFresh roots of Scutellaria baicalensis were treated with sodium nitroprusside （SNP） solutions at concentrations of 0.0， 7.5， and 20 mmol·L^-1， respectively. Kits and supporting reaction systems were used to determine the following indicators of samples in each group， including （1） reactive oxygen species： changes in the content of nitric oxide （NO）， superoxide anion （O2-·）， and hydrogen peroxide （H₂O₂）， （2） lipid oxidation products： changes in the content of malondialdehyde （MDA）， （3） antioxidant enzymes： changes in the activity of superoxide dismutase （SOD）， catalase （CAT）， and peroxidase （POD）， and （4） metabolic key enzymes and metabolic intermediates： Changes in the activity of phenylalanine ammonia-lyase （PAL） and the content of 1，3-diphosphoglycerate （1，3-DPG）. High-performance liquid chromatography （HPLC） was used to determine the changes in the content of secondary metabolites in flavonoids such as baicalin， wogonoside， baicalein， and wogonin in each group. ResultsCompared with that on day 0， during the treatment with 20 mmol·L^-1 SNP solution， the content of NO， O2-·， H₂O₂， and MDA in fresh roots of Scutellaria baicalensis increased significantly （P<0.01）， by 138.7%， 75.7%， 292.9%， and 140.3% respectively. The activity of SOD， CAT， and POD increased by 100.7%， 41.4%， and 70.4% respectively （P<0.01）， and that of PAL increased by 111.1% （P<0.01）. The content of the secondary metabolites baicalin and wogonoside remained basically unchanged， while that of baicalein and wogonin increased by 401.3% and 457.9% respectively （P<0.01）. The content of 1，3-DPG increased by 165.0% （P<0.01）. In the 7.5 mmol·L^-1 SNP solution treatment group， the stress degree was relatively mild. Compared with that on day 0， the content of NO， O2-·， H₂O₂， and MDA increased by 76.2%， 57.7%， 221.4%， and 99.7% respectively （P<0.01）. The activity of antioxidant enzymes SOD， CAT， and POD increased by 79.5%， 55.6%， and 51.9% respectively （P<0.01）， and that of PAL increased by 69.0% （P<0.01）. The content of baicalein and wogonin increased by 204.7% and 272.2% respectively （P<0.01）， and that of 1，3-DPG increased by 128.0% （P<0.01）. The 20 mmol·L^-1 SNP group showed better performance than the 7.5 mmol·L^-1 SNP group. ConclusionSNP can simulate the physiological state of Scutellaria baicalensis under stress conditions， activate the antioxidant protection mechanism， promote the biosynthesis of flavonoids in Scutellariae Radix， and significantly improve the quality of medicinal materials in Scutellariae Radix.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.

Background Long working hours, as a common risk factor for occupational stress, is closely related to the occurrence of depressive symptoms. Understanding how long working hours affect occupational stress and depressive symptoms will inform occupational health interventions. Objective To quantify the impact of long working hours exposure on occupational stress and depressive symptoms among Internet industry employees, translate black-box outputs into actionable insights, and demonstrate the value of interpretable machine learning for early-warning occupational-health surveillance. Methods A dataset was derived from a cross-sectional survey involving 2866 internet industry employees in China. This survey was part of the project Risk Assessment Of Long Working Hour Exposure And Its Adverse Health Effects, conducted by the National Institute for Occupational Health and Poisoning Control, Chinese Center for Disease Control and Prevention, from 2021 to 2023. Working hours, occupational stress and depressive symptoms were quantified with a set of structured questionnaires including the Core Occupational Stress Scale and the Patient Health Questionnaire. Pairwise associations were screened by Mantel tests and variance-inflation factors. Key predictors identified through feature selection were fed into six machine-learning risk-prediction models. Visual interpretation was provided by feature importance, Shapley additive explanations (SHAP) and local interpretable model-agnostic explanations (LIME), while directed causal effects and intervention impacts of prolonged working hours exposure on occupational stress and depressive symptoms were dissected with causal explanation of features techniques. Results The positive rates of occupational stress and depressive symptoms among internet employees were 12.9% and 77.8% respectively. Twelve core features for occupational stress and nine for depressive symptoms were retained after selection. After these features were supplied to six predictive algorithms and evaluated on five metrics, the Light Gradient Boosting Machine (LGBM) achieved the highest accuracy—0.89 for occupational stress and 0.79 for depressive symptoms on the hold-out test set. The feature-importance rankings converged on fatigue accumulation and life satisfaction as dominant drivers for both outcomes, whereas weekly working hours and daily overtime emerged as the principal exposure-related predictors. The SHAP summary plots revealed that longer weekly hours and daily overtime systematically elevated the probability of occupational stress. The causal feature explanation further quantified that ascending one category in weekly working hours increased the probability of occupational stress by 7.04%. Conclusion Exposure to long working hours is associated with both occupational stress and depressive symptoms among internet industry employees. Interpretable machine-learning frameworks translate these associations into transparent, defensible drivers, enabling precise identification of the pivotal factors and their interplay. This evidence base equips occupational-health practitioners with actionable insights for designing targeted prevention and intervention strategies.

Buyang Huanwu Decoction(BYHWD), as one of the classic formulas in traditional Chinese medicine(TCM) for the treatment of cerebral ischemic stroke(CIS), has demonstrated definite effects in clinical practice. However, the material basis and mechanism of treatment have not been systematically elucidated. This study employed network pharmacology and molecular docking to analyze the potential targets and mechanisms of blood-and brain-penetrating active components of BYHWD in reducing cell apoptosis in CIS. Cell experiments were then carried out to validate the prediction results. In the experiments, five active components including hydroxysafflor yellow A( HSYA), tetramethylpyrazine( TMP), astragaloside Ⅳ( AS-Ⅳ), amygdalin( AMY), and paeoniflorin(PF) were selected to explore the pharmacological effects of BYHWD. HT22 cells were treated with BYHWD, and the cell counting kit-8(CCK-8) method was employed to examine the toxic and side effects of BYHWD. A cell model of oxygen-glucose deprivation/reoxygenation( OGD/R) was constructed, with apoptosis and pyroptosis as the main screening indicators. The levels of lactate dehydrogenase(LDH) and glutathione(GSH) were measured to assess the cell membrane integrity. Flow cytometry was employed to detect apoptosis, and the activities of caspase-3 and caspase-1 were measured to clarify the status of apoptosis and pyroptosis. ELISA was employed to determine the levels of interleukin(IL)-1β and IL-18 to confirm pyroptosis. HSYA and AMY were identified in this study as the active components regulating apoptosis and pyroptosis. TUNEL was employed to detect the apoptosis rate, and Western blot was employed to determine the expression levels of apoptosis-related proteins B-cell lymphoma-2(Bcl-2), Bcl-2-associated X protein(Bax), and caspase-3, which confirmed that the anti-apoptotic effect of the combined component group was superior to that of the single component groups. The molecular docking results revealed strong binding affinity of HSYA and AMY with SDF-1α and CXCR4.AMD3100, a selective antagonist of CXCR4, was then used for intervention. The results of Western blot showed alterations in the expression levels of apoptosis-associated proteins, SDF-1α, and CXCR4. In conclusion, HSYA and AMY influence cellular apoptosis by modulating the SDF-1α/CXCR4 signaling cascade.
Drugs, Chinese Herbal/chemistry* ; Apoptosis/drug effects* ; Animals ; Neurons/cytology* ; Mice ; Molecular Docking Simulation ; Cell Line ; Glucose/metabolism* ; Humans ; Neuroprotective Agents/pharmacology*

OBJECTIVE: To reveal the molecular basis of knee osteoarthritis (KOA) with Yang deficiency and blood stasis syndrome by analyzing the gene expression profiles in synovial fluid and blood of KOA patients with this syndrome. METHODS: A total of 80 KOA patients were recruited from October 2022 to June 2024, including 40 cases in the non-Yang deficiency and blood stasis group (27 males and 13 females), with an average age of (61.75±3.45) years old;and 40 cases in the Yang deficiency and blood stasis group (22 males and 18 females), with an average age of (62.00±2.76) years old. The levels of body mass index (BMI), high-density lipoprotein (HDL), low-density lipoprotein (LDL), fibrinogen, total cholesterol, and D-dimer were recorded and summarized. Blood and synovial fluid samples from patients were collected for gene expression profile microarray sequencing, and then PCR and immunohistochemistry were used for clinical verification on the patients' synovial fluid and cartilage samples. RESULTS: Logistic regression analysis showed that compared with KOA patients with non-Yang deficiency and blood stasis syndrome, those with Yang deficiency and blood stasis syndrome had increased BMI, LDL, fibrinogen, total cholesterol, and D-dimer, and decreased HDL, with a clear correlation between the two groups. There were 562 differential genes in the blood, among which 322 were up-regulated and 240 were down-regulated;755 differential genes were found in the synovial fluid, with 350 up-regulated and 405 down-regulated. KEGG signaling pathway analysis of synovial fluid revealed changes in lipid metabolism-related pathways, including cholesterol metabolism, fatty acid metabolism, and PPARG signaling pathway. Analysis of the involved differential genes identified 6 genes in synovial fluid that were closely related to lipid metabolism, namely LRP1, LPL, ACOT6, TM6SF2, DGKK, and PPARG. Subsequently, PCR and immunohistochemical verification were performed using synovial fluid and cartilage samples, and the results were consistent with those of microarray sequencing. CONCLUSION This study explores the clinical and genomic correlation between traditional Chinese medicine syndromes and knee osteoarthritis from the perspective of lipid metabolism, and proves that abnormal lipid metabolism is closely related to KOA with Yang deficiency and blood stasis syndrome from both clinical and basic aspects.
Humans ; Male ; Female ; Middle Aged ; Synovial Fluid/metabolism* ; Osteoarthritis, Knee/metabolism* ; Yang Deficiency/complications* ; Aged

Aim To explore the effect of hydroxysafflor yellow A(HSYA)on lipocalin 2(LCN2)-induced fer-roptosis in HT22 cells and the related mechanism.Methods Thirty male Sprague-Dawley(SD)rats were used to establish the middle cerebral artery occlu-sion/reperfusion(MCAO/R)model by the suture method.The rats were randomly divided into the Sham group,the MCAO/R group,and the MCAO/R+HSYA group.The infarct area was measured by TTC staining,and the degree of neurological deficit was evaluated by the Z-Longa scoring method.The expressions of LCN2 and 24P3R in brain tissues were detected by Western blot.LCN2 protein was added to HT-22 cells,and the cells were divided into the normal group,the LCN2 group,and the LCN2+HSYA group.The optimal con-centration of LCN2-induced neuronal ferroptosis was screened by LDH assay and Western blot,and the ex-pression levels of ferritin,FPN1,GPX4,SLC7A11,COX2,and 24P3R were detected.LCN2 was knocked down by siRNA transfection,and the expressions of GPX4 and ferritin were detected.The contents of glu-tathione(GSH),malondialdehyde(MDA),GPX4,and Fe2+were determined by colorimetry,and the expres-sion of GPX4 was detected by immunofluorescence.The binding force between HSYA and LCN2 was ana-lyzed by molecular docking technology.Results Ani-mal experiments showed that HSYA could reduce the cerebral infarction area and decrease the neurological function score of MCAO/R rats.Compared with the sham group,the levels of LCN2 and 24P3R increased in the MCAO/R group,while HSYA inhibited their ex-pressions.Cell experiments showed that the optimal concentration of LCN2 to induce ferroptosis in HT22 cells was 2 μmol·L-1.After knocking down LCN2 by siRNA transfection,compared with the LCN2 group,the expression levels of GPX4 and ferritin in the siLCN2 group increased significantly.Compared with the nor-mal group,the expressions of SLC7A11,GPX4,FPN1,ferritin,and GSH in the LCN2 group decreased signifi-cantly,while the concentration of Fe2+,and the expres-sions of MDA,COX2,and 24P3R increased.HSYA could increase the expressions of SLC7A11,GPX4,FPN1,ferritin,and GSH,reduce the contents of Fe2+and MDA,and inhibit the expressions of COX2 and 24P3R.Molecular docking showed that the binding en-ergy between HSYA and LCN2 was-8.0 kJ·mol-1.Conclusion HSYA can inhibit LCN2-induced ferrop-tosis in HT22 cells through the SLC7A11/GPX4 signa-ling pathway.

Objective:To investigate whether hysteroscopic transcervical resection of septum (TCRS) prior to in vitro fertilization/intracytoplasmic sperm injection (IVF/ICSI) can improve cumulative live birth rates in patients with uterine septum. Methods:A retrospective cohort study was conducted to analyze data from 244 patients with partial uterine septum who underwent IVF/ICSI at the Reproductive Medicine Center of the Third Affiliated Hospital of Zhengzhou University between January 2016 and August 2022. The patients were divided into a surgical group ( n=171) and a non-surgical group ( n=73) based on whether TCRS was performed prior to IVF/ICSI. The clinical outcomes of IVF/ICSI in the two groups were analyzed, with the primary observation indicator being the cumulative live birth rate. Cox regression analysis was employed to identify determinants. Results:The age of patients in the operated group [(31.20±3.80) years] was younger than that in the non-operated group [(32.92±5.34) years, P=0.005], and the basal antral follicle count [17.0 (11.0, 24.0)] was higher than that in the non-operated group [14.0 (8.0, 21.5), P=0.039]. There were no significant differences in other baseline data (all P>0.05). The cumulative pregnancy rate [79.53% (136/171)] and the cumulative live birth rate [60.23% (103/171)] in the operated group during the 24-month follow-up period were significantly higher than those in the non-operated group [65.75% (48/73), P=0.022; 45.21% (33/73), P=0.030]. Compared with the operated group [296.0 (260.0, 430.0) d], the duration from the start of ovarian stimulation to the first live birth was significantly prolonged in the non-operated group [379.0 (329.5, 471.5) d, P<0.001]. Adjusted Cox-regression analysis showed that whether or not surgery was performed ( HR=1.683, 95% CI: 1.116-2.539, P=0.013) and the basal antral follicle count ( HR=1.032, 95% CI: 1.000-1.065, P=0.048) were independent factors affecting cumulative live birth rate. Conclusion:Performing TCRS before IVF/ICSI can improve cumulative live birth rates of patients with uterine septum.

Background Occupational stress has emerged as a critical public health concern affecting the physical and mental well-being of workers in the manufacturing sector. However, researchers typically evaluate its core driver—cumulative fatigue—using a crude binary “present/absent” variable, thereby overlooking the high-dimensional complexity and heterogeneity inherent in fatigue characteristics. This oversimplification constrains both the precision and predictive performance of occupational stress risk assessment model. Objective Leveraging a data-driven approach, to survey data on cumulative fatigue among manufacturing employees, and then use this new classification to develop and validate an occupational stress prediction model, with an ultimate aim of enhancing the accuracy and effectiveness of occupational stress assessment. Methods A set of cross-sectional survey data on 3871 manufacturing employees in 2021 were derived from the “Long working hours exposure and its adverse health effect risk assessment” program of the National Institute for Occupational Health and Poison Control, Chinese Center for Disease Control and Prevention. Occupational stress was assessed using the Core Occupational Stress Measurement Scale, while cumulative fatigue was evaluated via the Worker’s Self-Diagnostic Questionnaire for Fatigue Accumulation. Boruta method was applied to screen core variables from 20 occupational stress influencing factors. Dimensionality reduction of cumulative fatigue features was performed using a Bayesian sparse autoencoder (BASE), followed by comparison and application of clustering methods to achieve multi-class classification of the reduced features. Six machine learning classification models were then selected including logistic regression, support vector machine, decision tree, random forest, adaptive boosting, and lightweight gradient boosting machine (LightGBM) to construct and compare two occupational stress prediction models involving cumulative fatigue combined with ten other core variables: one utilizing the original binary cumulative fatigue label as a core variable, and another employing the novel fatigue classification proposed herein as a core variable. Results The positive rate of occupational stress among the manufacturing employees was 38.9%. The Boruta method identified 11 core variables of occupational stress: depressive symptoms, cumulative fatigue, age, length of service, tenure in current position, average weekly working hours, average daily overtime hours, average monthly income, low levels of exercise, life satisfaction, and sleep quality. The BSAE reduced the original cumulative fatigue factors to 12 dimensions, while the K-means clustering grouped the reduced fatigue features into three categories: no fatigue, moderate fatigue, and sever fatigue. Six occupational stress prediction models were constructed using the three-category labels of cumulative fatigue and ten additional factors. The results indicated that the LightGBM model performed best, achieving an area under the curve (AUC) of 0.78, an accuracy of 0.77, and an F1 score of 0.72. Compared with the model incorporating the traditional binary labels for cumulative fatigue, the best prediction AUC (0.72), accuracy (0.66), and F1 score (0.65) improved by 6%, 11%, and 7%, respectively. Conclusion Cumulative fatigue is a significant predictor of occupational stress among manufacturing employees. Applying data dimensionality reduction combined with three-class cluster analysis to characterize cumulative fatigue improves the performance of occupational stress prediction models. From a data-driven perspective, machine learning methods demonstrate strengths in processing complex datasets, capturing nonlinear relationships, and generating predictions based on key variables. This study therefore confirms that refined processing of core factors substantially enhances the predictive capability of machine learning models in assessing occupational stress risk in the manufacturing workforce.