BACKGROUND:The inflammatory response of microglia is closely related to neuronal survival,regeneration,and functional recovery after spinal cord injury.Peroxiredoxin 1 is not only involved in the regulation of oxidative stress,but also has an important effect on cell proliferation,apoptosis,and inflammatory response.OBJECTIVE:To investigate the role and mechanism of peroxiredoxin 1 in the inflammatory response of microglia following spinal cord injury.METHODS:(1)Twelve female C57BL/6 mice were randomly divided into sham-operated(n=6)and spinal cord injury(n=6)groups.The sham-operated group was not modeled and acute spinal cord injury models were constructed in the spinal cord injury group using the modified Allen's method.Spinal cord tissue at the injured site was taken at 7 days after modeling and transcriptome sequencing was performed to identify differentially expressed genes.The expression of peroxiredoxin 1 in spinal cord tissues was verified using western blot and RT-qPCR.(2)Mouse microglia BV2 were divided into two groups:the control group was stimulated with lipopolysaccharide for 6 hours,and in the knockout group,lipopolysaccharide stimulation was applied for 6 hours at 24 hours after peroxiredoxin 1 was knocked down in the cells.RT-qPCR was performed to detect mRNA expression of peroxiredoxin 1,inflammatory factors(interleukin 1β,interleukin 6,inducible nitric oxide synthase,tumor necrosis factor α,C-C motif chemokine ligand 2,and C-X-C motif chemokine ligand 2),and western blot was performed to detect the expression of peroxiredoxin 1,inducible nitric oxide synthase,and reactive oxygen/mitogen-activated protein kinase signaling pathway proteins.Mouse microglia BV2 were treated in two groups:the control group was stimulated by hydrogen peroxide for 4 hours,and the knockout group was stimulated by hydrogen peroxide for 4 hours at 24 hours after knockdown of peroxiredoxin 1.The level of reactive oxygen species was detected by 2,7-dichlorodihydrofluorescein diacetate probe.RESULTS AND CONCLUSION:(1)Results from transcriptome sequencing,western blot and RT-qPCR confirmed that peroxiredoxin 1 expression levels in mouse spinal cord tissues were significantly higher in the spinal cord injury group than the sham-operated group(P＜0.05).(2)Peroxiredoxin 1 knockdown in microglial cells led to decreased expression of peroxiredoxin 1 mRNA and protein(P＜0.05),increased mRNA expression of interleukin 1β,interleukin 6,inducible nitric oxide synthase,tumor necrosis factor α,C-C motif chemokine ligand 2,and C-X-C motif chemokine ligand 2(P＜0.05),increased protein expression of inducible nitric oxide synthase,P-P38,P-JNK and P-ERK proteins(P＜0.05),and increased level of reactive oxygen species(P＜0.05).To conclude,peroxiredoxin 1 regulates microglial inflammation by targeting the reactive oxygen species/mitogen-activated protein kinase signaling pathway.

Objective: To explore the influencing factors of anxiety symptoms and to construct a predictive model based on machine learning algorithms, so as to provide support for the prevention and management of anxiety symptoms among junior high school students. Methods: From April to May 2023, a stratified random cluster sampling method was adopted to select 8 176 junior high school students from Zhengzhou and Shangqiu citys. All participants completed the Adolescent Self rating Life Events Checklist, the 10item Connor-Davidson Resilience Scale, the School Connectedness Scale, the Parent-Child Cohesion Questionnaire, and the 7 item Generalized Anxiety Disorder Scale. Logistic regression analysis identified the associated factors of anxiety symptoms among junior high school students. Predictive models were constructed using Logistic regression, Random Forest, and eXtreme Gradient Boosting (XGBoost) algorithms, with SHapley Additive exPlanations analysis explaining the optimal model. Results: The detection rate of anxiety symptoms among junior high school students was 16.3%. Logistic regression analysis showed that junior high school students who were female ( OR =1.22), in the ninth grade ( OR =1.27), living in urban areas ( OR =1.37), having a father with a college education or above ( OR =1.26), having a mother with a senior high school education ( OR =1.26), and experiencing higher levels of negative life events ( OR =1.05) reported a higher risk of anxiety symptoms(all P <0.05). In contrast, those with moderate family economic status ( OR =0.71), moderate academic burden ( OR =0.59), low academic burden ( OR =0.54), moderate sleep quality ( OR =0.46), good sleep quality ( OR =0.26), excellent sleep quality ( OR =0.15), higher levels of psychological resilience ( OR =0.96), higher levels of school connectedness ( OR =0.96), and higher levels of parent-child cohesion ( OR =0.98) reported a lower risk of anxiety symptoms (all P <0.05). Three machine learning models demonstrated good predictive performance for anxiety symptoms among junior high school students (all AUC>0.8), with the XGBoost model achieving the best predictive performance. SHAP analysis revealed that negative life events, sleep quality, school connectedness, psychological resilience and parent-child cohesion were the top five relevant factors for predicting anxiety symptoms. Conclusions The detection rate of anxiety symptoms among junior high school students is relatively high. The XGBoost model is the optimal predictive model for anxiety symptoms in the population. Negative life events, sleep quality, school connectedness, psychological resilience, and parent-child cohesion are significant correlates of anxiety symptoms among junior high school students.

ObjectiveTo construct an automatic measurement model for palmar tilt and radial inclination suitable for traditional Chinese medicine （TCM） clinical scenarios， and to validate its accuracy and efficiency in TCM manipulative reduction settings. MethodsData on anteroposterior （AP） and lateral X-rays of distal radius fractures were collected from patients admitted to 18 TCM/ integrated TCM and western medicine hospitals in Jiangsu province between September 1st， 2023， and September 1st， 2024， via the Jiangsu Diagnosis and Treatment Big Data Platform for TCM Dominant Diseases. A medical image segmentation framework based on multi-scale feature fusion and edge-awareness was employed， combined with anatomical knowledge specific to TCM orthopedics， to optimize the feature extraction strategy of an artificial intelligence （AI） model. This framework enabled automatic segmentation of fracture regions and measurement of distal radius palmar tilt and radial inclination. The accuracy of the AI model in measuring radial inclination and volar tilt was validated， and the measurement time and average time gain rate of the AI model were compared to those of manual measurement. ResultsA total of 15,444 AP and lateral X-ray images of distal radius fractures were collected， and were divided into a training set （11,144 images， 5066 AP and 6078 lateral）， a validation set （3700 images， 1840 AP and 1860 lateral）， and an independent test set （600 images， 300 AP and 300 lateral） after preprocessing. In the measurement of 300 AP X-rays in the independent test set for radial inclination， when the degree error between AI measurement and manual measurement was ＜3° and ＜5°， AI measurement accuracy was 83% and 93%， respectively. In 300 lateral X-rays in the test set for palmar tilt， when AI measurements had an error of ＜3° and ＜5° compared to manual measurements， corresponding accuracy rate was 78% and 90%， respectively. For 50 X-ray images， AI measurement time was （1.37±0.05） min for radial inclination while manual measurement time was （22.57±2.52） min （P＜0.001）； in terms of palmar tilt， the AI measurement time was （1.33±0.14） min， shorter than （23.70±2.80） min for manual measurement time （P＜0.001）. Average time gain rates for manual and AI measurements were 93.93% and 94.39% respectively. ConclusionAn automatic measurement model for palmar tilt and radial inclination in distal radius fractures has been established， enabling more accurate and efficient assessment as well as providing a tool to support the quantitative evaluation of the efficacy of TCM manipulative reduction and large-sample clinical research.

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.

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.

ObjectiveTo investigate the expression of prefrontal cortical neurons， methyltransferase-like 3 （METTL3）， fat mass and obesity-associated gene （FTO）， and AlkB homolog 5 （ALKBH5） proteins in a mouse model of post-traumatic stress disorder （PTSD）. MethodsA PTSD mouse model was established using a single prolonged stress and foot shock stimulation （SPSS） method. The despair， anxiety， and learning and memory functions of PTSD mice were assessed through the open field test， Y-maze test， and forced swimming test. Neuronal damage was detected via HE and Nissl staining. The expression levels of METTL3， FTO， ALKBH5， and neuronal nuclear protein （NEUN） were assessed by Western blot and immunofluorescence staining. ResultsCompared to control group， PTSD mice subjected to SPSS exhibited signs of despair， anxiety， and impaired learning and memory. HE and Nissl staining results showed neuronal damage in the prefrontal cortex of PTSD mice. Western blot and immunofluorescence staining results showed that the expression of the m6A-related proteins METTL3 and FTO decreased， while the expression of ALKBH5 increased in the prefrontal cortex. Additionally， NEUN protein levels showed a declining trend. ConclusionThe pathogenesis of PTSD may be associated with neuronal damage in the prefrontal cortex and alterations in m6A methylation proteins.

SMYD5 is a ribosomal methyltransferase with SET and MYND structural domains， which is a member of the SMYD family and is expressed in a variety of tissues， including ovary and testis. This enzyme participates in biological processes such as gene expression regulation， cell development and differentiation， and maintenance of genomic stability through ribosomal protein methylation modification. In recent years， research on SMYD5 has increased in cancers including hepatocellular carcinoma， gastric adenocarcinoma， and lung cancer. Studies have revealed that SMYD5 exhibits high expression levels in various diseases including hepatocellular carcinoma， gastric adenocarcinoma， lung cancer， and inflammatory bowel disease， influencing the progression of these conditions. This review summarizes the role of SMYD5 in hepatocellular carcinoma， inflammatory bowel disease， and other biological functions， aiming to provide a reference for related disease research.

Objective To investigate the temporal changes in pathological damage and macrophage polarization in liver and spleen tissues of mice infected with Angiostrongylus cantonensis, and to preliminarily unravel the peripheral immune responses during the early stage of A. cantonensis infection. Methods Forty female BALB/c mice at ages of 6 to 8 weeks were randomly divided into four groups, including the control group and 7-, 14-, and 21-day infection groups, with 10 mice in each group. Each mouse in the infection groups was inoculated with 30 third-stage (L3) larvae of A. cantonensis by oral gavage, and five mice were randomly selected from each infection group on days 7, 14, and 21 post-infection, while mice in the control group were given the same volume of physiological saline and five mice were randomly selected from the control group on the day of oral gavage. Mouse liver and spleen tissues were sampled. The histopathological changes of mouse liver and spleen tissues were observed using hematoxylin and eosin (HE) staining, and the percentage of positive staining area and the co-localization positive rates of the macrophage surface antigens F4/80, CD86, and CD206 were quantified in mouse liver and spleen tissues using immunohistochemical and immunofluorescence staining. In addition, five mice were collected from each infection group on days 7, 14, and 21 post-infection, and five mice were collected from the control group on the day of oral gavage. Mouse liver and spleen tissues were sampled for detection of macrophage markers CD86 and CD206 and macrophage phenotyping using flow cytometry, and the expression of M1 macrophage markers, including inducible nitric oxide synthase (Nos2), tumor necrosis factor-α (TNF-α), interleukin-1β (IL-1β) and M2 markers, including arginase 1 (Arg1), mannose receptor C-type 1 (Mrc1) and chitinase-like protein 3 (Chil3) was quantified in mouse liver and spleen tissues using real-time quantitative PCR (RT-qPCR) assay. Results Proliferative lesions of the hepatocyte were observed in mouse liver tissues and the follicular structures of the mouse spleen white pulp were disrupted 21 days post-infection with A. cantonensis. Immunohistochemical staining showed that there were significant differences in the percentages of F4/80, CD86 and CD206 positive staining areas in the liver and spleen tissues among the four groups of mice (F = 242.40, 197.14, 183.19, 157.65, 242.35 and 146.24; all P values < 0.001), and the percentages of positive staining in the liver and spleen tissues of mice in the 14-day infection group [(4.45 ± 0.51)%, (3.74 ± 0.67)%, (8.32 ± 0.72)%, (16.56 ± 1.14)%, (11.62 ± 0.52)%, and (8.29 ± 0.72)%, respectively] and the 21-day infection group [(3.70 ± 0.11)%, (3.22 ± 0.43)%, (11.53 ± 1.03)%, (12.59 ± 1.05)%, (9.02 ± 0.83)%, and (11.67 ± 1.10)%, respectively] were higher than in the control group [(0.35 ± 0.16)%, (0.40 ± 0.02)%, (0.93 ± 0.05)%, (2.78 ± 0.26)%, (2.33 ± 0.20)%, and (1.85 ± 0.20)%, respectively] (all P values < 0.05). Immunofluorescence staining showed significant differences in the positive rates of F4/80 co-localization with CD86 and CD206 in mouse liver and spleen tissues among the four groups (F = 24.42, 25.28, 54.51 and 130.55; all P values < 0.001). Flow cytometry detected significant differences in the proportions of CD86⁺ and CD206⁺ macrophages in mouse liver and spleen tissues among the four groups (F = 67.98, 18.41, 29.77, 172.80; all P values < 0.001), and the proportions of CD206⁺ macrophages in the liver and spleen of the 21-day infection group were significantly higher than those in the control group [(9.25 ± 2.55)% vs (3.83 ± 0.72)%, and (4.22 ± 0.56)% vs (0.47 ± 0.18)%, respectively] (both P values < 0.05). In addition, RT-qPCR assay quantified significant differences in the relative mRNA expression of M1 macrophage markers (IL-1β, TNF-α and Nos2) and M2 macrophage markers (Arg1, Chil3 and Mrc1) in mouse liver and spleen tissues among the four groups (F = 41.30, 31.82, 199.33, 19.96, 62.01, 119.76, 23.67, 95.90, 72.27, 82.59, 123.41 and 29.75; all P values < 0.05). Conclusions A. cantonensis infection may cause progressive pathological damage in mouse liver and spleen tissues, accompanied by dynamic temporal changes in macrophage polarization. M1 macrophage polarization predominates at the early stage of A. cantonensis infection and shifts towards M2 polarization at the later stages, suggesting that M2 polarization may participate in immune regulation at late stages of A. cantonensis infection by suppressing excessive inflammatory responses and promoting tissue repair.

BACKGROUND:The inflammatory response of microglia is closely related to neuronal survival,regeneration,and functional recovery after spinal cord injury.Peroxiredoxin 1 is not only involved in the regulation of oxidative stress,but also has an important effect on cell proliferation,apoptosis,and inflammatory response.OBJECTIVE:To investigate the role and mechanism of peroxiredoxin 1 in the inflammatory response of microglia following spinal cord injury.METHODS:(1)Twelve female C57BL/6 mice were randomly divided into sham-operated(n=6)and spinal cord injury(n=6)groups.The sham-operated group was not modeled and acute spinal cord injury models were constructed in the spinal cord injury group using the modified Allen's method.Spinal cord tissue at the injured site was taken at 7 days after modeling and transcriptome sequencing was performed to identify differentially expressed genes.The expression of peroxiredoxin 1 in spinal cord tissues was verified using western blot and RT-qPCR.(2)Mouse microglia BV2 were divided into two groups:the control group was stimulated with lipopolysaccharide for 6 hours,and in the knockout group,lipopolysaccharide stimulation was applied for 6 hours at 24 hours after peroxiredoxin 1 was knocked down in the cells.RT-qPCR was performed to detect mRNA expression of peroxiredoxin 1,inflammatory factors(interleukin 1β,interleukin 6,inducible nitric oxide synthase,tumor necrosis factor α,C-C motif chemokine ligand 2,and C-X-C motif chemokine ligand 2),and western blot was performed to detect the expression of peroxiredoxin 1,inducible nitric oxide synthase,and reactive oxygen/mitogen-activated protein kinase signaling pathway proteins.Mouse microglia BV2 were treated in two groups:the control group was stimulated by hydrogen peroxide for 4 hours,and the knockout group was stimulated by hydrogen peroxide for 4 hours at 24 hours after knockdown of peroxiredoxin 1.The level of reactive oxygen species was detected by 2,7-dichlorodihydrofluorescein diacetate probe.RESULTS AND CONCLUSION:(1)Results from transcriptome sequencing,western blot and RT-qPCR confirmed that peroxiredoxin 1 expression levels in mouse spinal cord tissues were significantly higher in the spinal cord injury group than the sham-operated group(P＜0.05).(2)Peroxiredoxin 1 knockdown in microglial cells led to decreased expression of peroxiredoxin 1 mRNA and protein(P＜0.05),increased mRNA expression of interleukin 1β,interleukin 6,inducible nitric oxide synthase,tumor necrosis factor α,C-C motif chemokine ligand 2,and C-X-C motif chemokine ligand 2(P＜0.05),increased protein expression of inducible nitric oxide synthase,P-P38,P-JNK and P-ERK proteins(P＜0.05),and increased level of reactive oxygen species(P＜0.05).To conclude,peroxiredoxin 1 regulates microglial inflammation by targeting the reactive oxygen species/mitogen-activated protein kinase signaling pathway.

Objective To systematically analyze the spatiotemporal distribution characteristics and epidemiological trends of tracheal, bronchus, and lung cancer (TBL) disease burden attributed to air pollution globally and in China and the United States from 1990 to 2021, and to assess the patterns of disease burden changes from 2022 to 2031 based on predictive models, providing a scientific basis for formulating targeted TBL prevention and control strategies. Methods Based on the Global Burden of Disease (GBD) 2021 database, we analyzed the disease burden data of TBL attributed to air pollution globally and in China and the United States from 1990 to 2021. R Studio 4.3.2 software was used to analyze the corresponding trends and the Bayesian age-period-cohort (BAPC) prediction model was used to predict the status of the disease burden of TBL attributed to air pollution in the world and in China and the United States from 2022 to 2031. Results In 2021, China had the highest number of deaths and disability-adjusted life years attributed to air pollution (211 400 patients and 4.8947 million person-years), followed by the United States (6 000 patients and 124 300 person-years). The age-standardized mortality rate (ASMR) and age-standardized disability-adjusted life years rate (ASDR) of TBL due to air pollution in the world and in China and the United States showed a decreasing trend. From 1990 to 2021, the ASMR and ASDR of TBL in China due to air pollution were much higher than those in the United States and the global average. In terms of gender, from 1990 to 2021, the disease burden of male patients with TBL attributed to air pollution was much higher than that of female patients. The BAPC prediction model showed that from 2022 to 2031, the ASMR and ASDR of TBL attributed to air pollution showed an upward trend globally, while they showed a downward trend in China and the United States. Conclusion Over the past 30 years, the air pollution-related TBL disease burden in the world and in China and the United States has continued to decline, but China's disease burden is still significantly higher than the global average. The disease burden in men far exceeds that in women, with men and the population aged ≥50 years being high-risk groups. In the future, the global disease trend may reverse and rise, while China and the United States are expected to continuously decline. However, precise prevention and control for high-risk groups remains a key challenge.