Objective: To explore the nonlinear dynamic effects of social capital, adverse childhood experiences (ACEs) and depressive symptoms on suicidal behavior among vocational high school students, so as to provide theoretical basis and practical references for formulating suicide prevention strategies. Methods: A convenience sampling method was employed to include 668 students from a vocational high school from Wuhan in March 2023. Social capital was used as the asymmetry variable, while ACEs and depressive symptoms were used as bifurcation variables, a cusp catastrophe model was constructed to analyze the nonlinear changes in suicidal behavior among vocational high school students, and its fit was compared with linear and Logistic regression models. Results: Among students in the health vocational high school in Wuhan, only suicidal ideation accounted for 8.5%, only suicide attempt for 18.6%, neither accounted for 31.9%, and both for 41.0%. Gender, left behind experience, family economic status, parental parenting styles, depressive symptoms, social capital, and ACEs were all related factors influencing suicidal behavior among vocational high school students ( χ 2/H=19.03, 13.33, 21.11, 46.70, 144.38, 24.61, 118.77, all P <0.05). Violin plots showed a bimodal distribution of suicidal behavior, indicating nonlinear variation characteristics. The cusp catastrophe model results showed that social capital was negatively correlated with suicidal behavior, but the relationship was bifurcated by ACEs ( α social capital = -0.006 , β ACEs =0.075) and depressive symptoms ( α social capital =-0.013, β depressive =0.028) (all P <0.05). When both ACEs and depressive symptoms coexisted, the impact of ACEs was stronger ( β ACEs =0.077, β depressive =0.014) (both P <0.05). The cusp catastrophe model fitted ( R 2=0.886, 0.881, 0.882) better than the linear ( R 2=0.258, 0.219, 0.258) and Logistic regression models ( R 2= 0.242, 0.211 , 0.176). Gender stratified analysis results showed that bifurcation effect of ACEs was stronger in males than in females( β boys =0.224, β girls =0.086); in females, both ACEs and depressive symptoms had a bifurcation effect, with the former showing a stronger effect ( β ACEs =0.062, β depressive =0.015) (all P <0.05). Conclusions Suicidal behavior among vocational high school students exhibits nonlinear characteristics. Improving social capital to reducing ACEs and depressive symptoms may contribute to decreasing adolescent suicidal behaviors.

This paper introduces a real-world cohort research model for community-acquired pneumonia （CAP） based on the Jiangsu Traditional Chinese Medicine （TCM） Dominant Diseases Diagnosis and Treatment Data Platform. Firstly， data cleaning is performed by standardizing diagnosis， symptoms， treatment and imaging， intelligently extracting unstructured information， and cleaning and constructing a standardized database. Secondly， for cohort establishment， CAP patients across the province are screened in accordance with CAP diagnostic criteria to build a high-quality disease-specific cohort. Lastly， in terms of protocol design， the characteristics of TCM research and the CAP disease profile are considered to determine appropriate inclusion and exclusion criteria， estimate sample size， define interventions， outcomes and economic evaluations， providing a reference for real-world TCM research on CAP.

This paper introduces a real-world cohort research model for community-acquired pneumonia （CAP） based on the Jiangsu Traditional Chinese Medicine （TCM） Dominant Diseases Diagnosis and Treatment Data Platform. Firstly， data cleaning is performed by standardizing diagnosis， symptoms， treatment and imaging， intelligently extracting unstructured information， and cleaning and constructing a standardized database. Secondly， for cohort establishment， CAP patients across the province are screened in accordance with CAP diagnostic criteria to build a high-quality disease-specific cohort. Lastly， in terms of protocol design， the characteristics of TCM research and the CAP disease profile are considered to determine appropriate inclusion and exclusion criteria， estimate sample size， define interventions， outcomes and economic evaluations， providing a reference for real-world TCM research on CAP.

ObjectiveTo identify potential influencing factors of urate crystal deposition at ankle/foot in patients with hyperuricemia （HUA）， and to analyze the predictive value of inflammatory indicators for urate crystal deposition in patients with different traditional Chinese medicine （TCM） syndromes， so as to provide potential reference for clinical risk assessment and individualized TCM intervention. MethodsA retrospective study was carried out with the enrollment of 231 HUA patients from The Third Affiliated People's Hospital of Fujian University of Traditional Chinese Medicine between January 2021 and December 2024. The enrolled patients were further divided into a crystal deposition-positive group （143 cases） and a crystal deposition-negative group （88 cases） according to the results of dual-energy computed tomography （CT）. Sociodemographic data， living habits， serum uric acid levels， and inflammatory indicators of the enrolled patients were collcted， and TCM syndrome differentiation was performed. Furthermore， univariate analysis was used to compare inter-group differences in clinical characteristics. MMultivariate Logistic regression was applied to identify the influencing factors of urate crystal deposition. In addition， the receiver operating characteristic （ROC） curves were plotted to evaluate the predictive efficacy of inflammatory indicators for crystal deposition across different TCM syndromes. ResultsThere were statistically significant inter-group differences in the proportion of males， age， body mass index， proportion of mental labor， rate of low water intake， and rate of high-sugar beverage consumption （P<0.05），whereas no significant difference in low exercise intensity was found between the two groups. Furthermore， compared with the negative group， the positive group had higher serum uric acid level， neutrophil-to-lymphocyte ratio （NLR）， and platelet-to-lymphocyte ratio （PLR）， but lower systemic immune-inflammation index （SIRI） （P<0.05）. Regarding the distribution of TCM syndromes， the positive group was dominated by the dampness-heat accumulation syndrome （55/143，38.46%）， while the negative group was mainly characterized by the phlegm-turbidity obstruction syndrome （44/88，50.00%）. Multivariate Logistic regression analysis revealed that high-sugar beverage consumption， elevated NLR， and elevated PLR were risk factors for urate crystal deposition ［odd ratio （OR） = 8.002， 5.377， 1.034， respectively； 95% CI 1.572-40.732， 2.179-13.270， 1.013-1.054，all P<0.05］， while SIRI was a protective factor （OR = 0.869， 95% CI 0.778-0.971， P<0.05）. In the positive group， patients with the dampness-heat accumulation syndrome exhibited the highest NLR， while the lowest PLR and SIRI， showing statistically significant differences with those of other syndromes （all P<0.05）. In addition， ROC curve analysis indicated that for the dampness-heat accumulation syndrome， the combined "NLR + PLR" model had an area under the curve （AUC） of 0.901 （95% CI 0.850-0.951， P<0.01）， with a sensitivity of 89.1% and a specificity of 79.5%； for the blood stasis-heat obstruction syndrome， the combined "NLR + PLR" model had an AUC of 0.880 （95% CI 0.825-0.934， P<0.01）， with a sensitivity of 100.0% and a specificity of 67.3%； for the liver-kidney Yin-deficiency syndrome， the single PLR model had an AUC of 0.842 （95% CI 0.731-0.952， P<0.01）， with a sensitivity of 83.3% and a specificity of 84.0%. ConclusionUrate crystal deposition in HUA patients exhibits intimate associations with high-sugar beverage consumption as well as elevated NLR and PLR levels. Meanwhile， TCM syndrome differentiation has potential correlation with inflammatory characteristics. The inflammatory indicator-based prediction model constructed based on TCM syndromes exhibits good predictive value.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

ObjectiveTranscranial magneto-acoustic stimulation (TMAS) is an emerging non-invasive neuromodulation technique that may provide a novel non-pharmacological intervention strategy for Parkinson's disease (PD). PD is characterized by the progressive degeneration of dopaminergic neurons in the substantia nigra pars compacta (SNc), leading to motor impairments such as bradykinesia, tremor, and rigidity. Increasing evidence indicates that mitochondrial dysfunction and impaired mitochondrial quality control are central mechanisms underlying dopaminergic neuronal loss. In particular, abnormalities in mitophagy and mitochondrial fission-fusion balance contribute substantially to oxidative stress, energy metabolic failure, and neuronal injury. At present, most clinical treatments for PD mainly alleviate symptoms but do not effectively halt disease progression. Therefore, exploring new interventions targeting the core pathological mechanisms is of considerable significance. This study aims to investigate whether TMAS can improve neural damage and motor dysfunction in PD mice by regulating mitophagy and the fission/fusion dynamic balance, thereby providing theoretical and experimental support for its application in PD treatment. MethodsMale C57BL/6 mice were used in this study. A PD model was established by intraperitoneal injection of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) for 7 consecutive days. After model induction, mice in the intervention group received TMAS once daily for 14 consecutive days, whereas the corresponding control group received sham stimulation. The stimulation target was positioned over the primary motor cortex (M1). Motor performance was evaluated using the pole test and the open-field test. To verify the activation effect of TMAS on the target cortical region, c-Fos immunohistochemistry was performed in the M1. To assess nigral dopaminergic neuronal injury, tyrosine hydroxylase (TH) immunohistochemistry was used to quantify TH-positive neurons in the SNc. Mitochondrial function was evaluated by measuring reactive oxygen species (ROS) levels and adenosine triphosphate (ATP) content in the SNc. Western blot was further performed to determine the expression of mitophagy-related proteins, including PINK1, Parkin, LC3-II, and p62, as well as mitochondrial dynamics-related proteins, including Drp1 and Opa1. ResultsTMAS significantly increased the number of c-Fos-positive cells in M1 (P<0.000 1), indicating effective activation of neurons in the targeted cortical region. Compared with the control group, MPTP-treated mice exhibited marked motor dysfunction, including a significant reduction in total distance traveled in the open-field test (P<0.000 1) and mean speed (P=0.000 1), as well as significant prolongation of turn time and total climbing time in the pole test (P<0.000 1). These behavioral impairments were accompanied by a substantial loss of TH-positive dopaminergic neurons in the SNc, whereas TMAS significantly increased TH-positive neuron survival (P<0.000 1). In parallel, MPTP induced a pronounced increase in ROS levels and a significant reduction in ATP content, indicating severe mitochondrial dysfunction and energy metabolism impairment (P<0.01). TMAS treatment significantly improved motor performance, as reflected by the reversal of MPTP-induced impairment in the open-field and pole tests, and significantly reduced ROS accumulation (P<0.01) while restoring ATP production (P<0.001). At the molecular level, MPTP markedly downregulated PINK1 and Parkin, decreased p62 expression, increased LC3-II accumulation, elevated Drp1 expression, and reduced Opa1 expression, whereas TMAS significantly reversed these abnormalities, suggesting restoration of mitophagy-related mitochondrial quality control and re-establishment of mitochondrial fission-fusion balance. Collectively, these findings indicate that TMAS ameliorates MPTP-induced neurotoxicity and restores mitochondrial homeostasis and energy metabolism. ConclusionTMAS effectively attenuates neural damage and improves motor dysfunction in MPTP-induced PD mice. Its neuroprotective effects are closely associated with multidimensional regulation of the mitochondrial quality control system, including restoration of PINK1/Parkin-mediated mitophagy and rebalancing of Drp1/Opa1-related mitochondrial dynamics. Rather than acting only as a symptomatic neuromodulatory intervention, TMAS may influence a key pathological axis of PD by improving mitochondrial homeostasis in SNc and protecting nigral dopaminergic neurons. These findings provide experimental evidence supporting TMAS as a promising non-invasive physical intervention for PD.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

ObjectiveTranscranial magneto-acoustic stimulation (TMAS) is an emerging non-invasive neuromodulation technique that may provide a novel non-pharmacological intervention strategy for Parkinson's disease (PD). PD is characterized by the progressive degeneration of dopaminergic neurons in the substantia nigra pars compacta (SNc), leading to motor impairments such as bradykinesia, tremor, and rigidity. Increasing evidence indicates that mitochondrial dysfunction and impaired mitochondrial quality control are central mechanisms underlying dopaminergic neuronal loss. In particular, abnormalities in mitophagy and mitochondrial fission-fusion balance contribute substantially to oxidative stress, energy metabolic failure, and neuronal injury. At present, most clinical treatments for PD mainly alleviate symptoms but do not effectively halt disease progression. Therefore, exploring new interventions targeting the core pathological mechanisms is of considerable significance. This study aims to investigate whether TMAS can improve neural damage and motor dysfunction in PD mice by regulating mitophagy and the fission/fusion dynamic balance, thereby providing theoretical and experimental support for its application in PD treatment. MethodsMale C57BL/6 mice were used in this study. A PD model was established by intraperitoneal injection of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) for 7 consecutive days. After model induction, mice in the intervention group received TMAS once daily for 14 consecutive days, whereas the corresponding control group received sham stimulation. The stimulation target was positioned over the primary motor cortex (M1). Motor performance was evaluated using the pole test and the open-field test. To verify the activation effect of TMAS on the target cortical region, c-Fos immunohistochemistry was performed in the M1. To assess nigral dopaminergic neuronal injury, tyrosine hydroxylase (TH) immunohistochemistry was used to quantify TH-positive neurons in the SNc. Mitochondrial function was evaluated by measuring reactive oxygen species (ROS) levels and adenosine triphosphate (ATP) content in the SNc. Western blot was further performed to determine the expression of mitophagy-related proteins, including PINK1, Parkin, LC3-II, and p62, as well as mitochondrial dynamics-related proteins, including Drp1 and Opa1. ResultsTMAS significantly increased the number of c-Fos-positive cells in M1 (P<0.000 1), indicating effective activation of neurons in the targeted cortical region. Compared with the control group, MPTP-treated mice exhibited marked motor dysfunction, including a significant reduction in total distance traveled in the open-field test (P<0.000 1) and mean speed (P=0.000 1), as well as significant prolongation of turn time and total climbing time in the pole test (P<0.000 1). These behavioral impairments were accompanied by a substantial loss of TH-positive dopaminergic neurons in the SNc, whereas TMAS significantly increased TH-positive neuron survival (P<0.000 1). In parallel, MPTP induced a pronounced increase in ROS levels and a significant reduction in ATP content, indicating severe mitochondrial dysfunction and energy metabolism impairment (P<0.01). TMAS treatment significantly improved motor performance, as reflected by the reversal of MPTP-induced impairment in the open-field and pole tests, and significantly reduced ROS accumulation (P<0.01) while restoring ATP production (P<0.001). At the molecular level, MPTP markedly downregulated PINK1 and Parkin, decreased p62 expression, increased LC3-II accumulation, elevated Drp1 expression, and reduced Opa1 expression, whereas TMAS significantly reversed these abnormalities, suggesting restoration of mitophagy-related mitochondrial quality control and re-establishment of mitochondrial fission-fusion balance. Collectively, these findings indicate that TMAS ameliorates MPTP-induced neurotoxicity and restores mitochondrial homeostasis and energy metabolism. ConclusionTMAS effectively attenuates neural damage and improves motor dysfunction in MPTP-induced PD mice. Its neuroprotective effects are closely associated with multidimensional regulation of the mitochondrial quality control system, including restoration of PINK1/Parkin-mediated mitophagy and rebalancing of Drp1/Opa1-related mitochondrial dynamics. Rather than acting only as a symptomatic neuromodulatory intervention, TMAS may influence a key pathological axis of PD by improving mitochondrial homeostasis in SNc and protecting nigral dopaminergic neurons. These findings provide experimental evidence supporting TMAS as a promising non-invasive physical intervention for PD.

Background As climate change intensifies and extreme temperature events become more frequent, non-optimal temperature has emerged as a significant contributor to the global disease burden, representing a pressing public health challenge. Objective To analyze the disease burden, temporal trends, and health inequalities attributable to non-optimal, high, and low temperatures in China from 1990 to 2021, and to compare these findings with global levels to provide a scientific basis for targeted prevention strategies. Methods Using data from the Global Burden of Disease 2021 (GBD 2021), we extracted mortality rates and disability-adjusted life year (DALY) rates, and other indicators attributable to non-optimal, high, and low temperatures by sex, age, region, and cause. Joinpoint regression was applied to examine temporal trends. Decomposition analysis identified driving factors of change, while the slope index of inequality (SII) and concentration index (CI) quantified disparities across socio-demographic index (SDI) levels. Results From 1990 to 2021, the age-standardized mortality rates (ASMR) and age-standardized DALY rates (ASDR) attributable to non-optimal temperature in China exhibited a downward trend, decreasing from 66.48 (95%UI: 58.09, 76.56) to 32.70 (95%UI: 27.26, 39.26) per 100000 population, and from 1219.59 (95%UI: 1056.28, 1418.37) to 493.22 (95%UI: 403.88, 609.32) per 100000 population, respectively. Burdens attributable to non-optimal temperature and low temperature were higher than the global average, whereas the high temperature burden was lower. Males consistently experienced higher ASMR and ASDR attributable to non-optimal temperature than females. Cardiovascular diseases, chronic respiratory diseases, and respiratory infections and tuberculosis were the top three causes of non-optimal temperature-attributable burdens. Decomposition analysis revealed that population aging and growth were the primary drivers of increased burden, while epidemiological changes primarily drove the decline. Health inequalities were most predominant between extreme SDI regions but narrowed over time. Conclusion Despite the overall decline in burden attributable to non-optimal temperature in China, significant challenges remain, including high risks from cold exposure, gender disparities, and the compounding effects of an aging population with cardiovascular or respiratory diseases. Policy makers should prioritize climate change adaptation, focusing on elderly health and regional equity while strengthening the public health workforce.

Background As climate change intensifies and extreme temperature events become more frequent, non-optimal temperature has emerged as a significant contributor to the global disease burden, representing a pressing public health challenge. Objective To analyze the disease burden, temporal trends, and health inequalities attributable to non-optimal, high, and low temperatures in China from 1990 to 2021, and to compare these findings with global levels to provide a scientific basis for targeted prevention strategies. Methods Using data from the Global Burden of Disease 2021 (GBD 2021), we extracted mortality rates and disability-adjusted life year (DALY) rates, and other indicators attributable to non-optimal, high, and low temperatures by sex, age, region, and cause. Joinpoint regression was applied to examine temporal trends. Decomposition analysis identified driving factors of change, while the slope index of inequality (SII) and concentration index (CI) quantified disparities across socio-demographic index (SDI) levels. Results From 1990 to 2021, the age-standardized mortality rates (ASMR) and age-standardized DALY rates (ASDR) attributable to non-optimal temperature in China exhibited a downward trend, decreasing from 66.48 (95%UI: 58.09, 76.56) to 32.70 (95%UI: 27.26, 39.26) per 100000 population, and from 1219.59 (95%UI: 1056.28, 1418.37) to 493.22 (95%UI: 403.88, 609.32) per 100000 population, respectively. Burdens attributable to non-optimal temperature and low temperature were higher than the global average, whereas the high temperature burden was lower. Males consistently experienced higher ASMR and ASDR attributable to non-optimal temperature than females. Cardiovascular diseases, chronic respiratory diseases, and respiratory infections and tuberculosis were the top three causes of non-optimal temperature-attributable burdens. Decomposition analysis revealed that population aging and growth were the primary drivers of increased burden, while epidemiological changes primarily drove the decline. Health inequalities were most predominant between extreme SDI regions but narrowed over time. Conclusion Despite the overall decline in burden attributable to non-optimal temperature in China, significant challenges remain, including high risks from cold exposure, gender disparities, and the compounding effects of an aging population with cardiovascular or respiratory diseases. Policy makers should prioritize climate change adaptation, focusing on elderly health and regional equity while strengthening the public health workforce.