ObjectiveThis paper aims to investigate the material basis of the anti-inflammatory efficacy and mechanism of action of Bushen Tongdu prescription （BSTDP）. MethodsThe chemical components of BSTDP and its blood-absorbed components in vivo were systematically identified by using ultra-performance liquid chromatography-linear ion trap-electrostatic field orbitrap high-resolution mass spectrometry （UPLC-LIT-Orbitrap-MS）. Network pharmacology was employed to screen blood-absorbed bioactive components and potential targets of this formula. A protein-protein interaction （PPI） network of core targets was constructed to conduct enrichment analysis. Molecular docking was further utilized to verify the binding affinity between key components and targets. The inflammatory model was established and verified in vivo by using a transgenic zebrafish Tg （mpx： GFP）. At three days post-fertilization （3 dpf）， larvae of zebrafish were randomly assigned to blank group， model group， positive drug dexamethasone acetate group （75 μmol·L^-1）， and BSTDP groups with low， medium， and high doses （500， 1 000， and 2 000 mg·L^-1）. The distribution and quantity of neutrophils in the yolk sac region were observed under a fluorescence microscope. The mRNA expression levels of key genes in the toll-like receptor 4 （TLR4）/myeloid differentiation factor 88 （MyD88）/nuclear factor kappa-B （NF-κB） signaling pathway and inflammatory factors including interleukin （IL）-1β， IL-6， and tumor necrosis factor-α （TNF-α） were detected by Real-time quantitative polymerase chain reaction （Real-time PCR）. ResultsA total of 120 chemical components were identified in BSTDP， among which 26 original components were confirmed by using serum pharmacochemical methods. A total of 227 common targets linking rheumatoid arthritis （RA） and the blood-absorbed components were screened by network pharmacology. It is suggested that pseudobrucine， vomicine， sinapine， rehmannioside， cinnamyl alcohol glycoside， and methylephedrine exert anti-inflammatory effects by acting on core targets including protein kinase B1 （Akt1）， signal transducer and activator of transcription 3 （STAT3）， tumor necrosis factor （TNF）， TLR4， mitogen-activated protein kinase 14 （MAPK14）， and phosphatidylinositol-4，5-bisphosphate 3-kinase catalytic subunit α （PIK3CA）， thereby modulating multiple signaling pathways such as TLR4 and NF-κB. In vivo verification in zebrafish demonstrates that the maximum tolerable concentration of Bushen Tongdu Formula is 2 000 mg·L^-1. Compared to those in the blank group， zebrafish in the model group showed a significantly higher number of neutrophils in the yolk sac region （P<0.01） and rising mRNA levels of TLR4， MyD88， NF-κB， TNF-α， IL-6， and IL-1β （P<0.01）. Compared to that in the model group， the number of neutrophils was significantly reduced in BSTDP groups with medium and high doses， as well as the dexamethasone acetate group （P<0.05， P<0.01）. There was no statistically significant difference in the low dose group. The mRNA expression levels of TLR4， MyD88， NF-κB， TNF-α， IL-6， and IL-1β were significantly down-regulated （P<0.05， P<0.01）. ConclusionThis paper identifies the material basis of the efficacy of BSTDP， demonstrating that the formula can exert an anti-inflammatory effect through the TLR4/MyD88/NF-κB signaling pathway. The results provide scientific experimental evidence for its further clinical application.

ObjectiveThis paper aims to investigate the material basis of the anti-inflammatory efficacy and mechanism of action of Bushen Tongdu prescription （BSTDP）. MethodsThe chemical components of BSTDP and its blood-absorbed components in vivo were systematically identified by using ultra-performance liquid chromatography-linear ion trap-electrostatic field orbitrap high-resolution mass spectrometry （UPLC-LIT-Orbitrap-MS）. Network pharmacology was employed to screen blood-absorbed bioactive components and potential targets of this formula. A protein-protein interaction （PPI） network of core targets was constructed to conduct enrichment analysis. Molecular docking was further utilized to verify the binding affinity between key components and targets. The inflammatory model was established and verified in vivo by using a transgenic zebrafish Tg （mpx： GFP）. At three days post-fertilization （3 dpf）， larvae of zebrafish were randomly assigned to blank group， model group， positive drug dexamethasone acetate group （75 μmol·L^-1）， and BSTDP groups with low， medium， and high doses （500， 1 000， and 2 000 mg·L^-1）. The distribution and quantity of neutrophils in the yolk sac region were observed under a fluorescence microscope. The mRNA expression levels of key genes in the toll-like receptor 4 （TLR4）/myeloid differentiation factor 88 （MyD88）/nuclear factor kappa-B （NF-κB） signaling pathway and inflammatory factors including interleukin （IL）-1β， IL-6， and tumor necrosis factor-α （TNF-α） were detected by Real-time quantitative polymerase chain reaction （Real-time PCR）. ResultsA total of 120 chemical components were identified in BSTDP， among which 26 original components were confirmed by using serum pharmacochemical methods. A total of 227 common targets linking rheumatoid arthritis （RA） and the blood-absorbed components were screened by network pharmacology. It is suggested that pseudobrucine， vomicine， sinapine， rehmannioside， cinnamyl alcohol glycoside， and methylephedrine exert anti-inflammatory effects by acting on core targets including protein kinase B1 （Akt1）， signal transducer and activator of transcription 3 （STAT3）， tumor necrosis factor （TNF）， TLR4， mitogen-activated protein kinase 14 （MAPK14）， and phosphatidylinositol-4，5-bisphosphate 3-kinase catalytic subunit α （PIK3CA）， thereby modulating multiple signaling pathways such as TLR4 and NF-κB. In vivo verification in zebrafish demonstrates that the maximum tolerable concentration of Bushen Tongdu Formula is 2 000 mg·L^-1. Compared to those in the blank group， zebrafish in the model group showed a significantly higher number of neutrophils in the yolk sac region （P<0.01） and rising mRNA levels of TLR4， MyD88， NF-κB， TNF-α， IL-6， and IL-1β （P<0.01）. Compared to that in the model group， the number of neutrophils was significantly reduced in BSTDP groups with medium and high doses， as well as the dexamethasone acetate group （P<0.05， P<0.01）. There was no statistically significant difference in the low dose group. The mRNA expression levels of TLR4， MyD88， NF-κB， TNF-α， IL-6， and IL-1β were significantly down-regulated （P<0.05， P<0.01）. ConclusionThis paper identifies the material basis of the efficacy of BSTDP， demonstrating that the formula can exert an anti-inflammatory effect through the TLR4/MyD88/NF-κB signaling pathway. The results provide scientific experimental evidence for its further clinical application.

ObjectiveTo establish a mouse model of hypoxic pulmonary hypertension (HPH) through intermittent hypoxia induction, and to develop a comprehensive and reliable evaluation system for the HPH model, thereby providing experimental evidence for mechanistic studies and translational applications related to this disease. MethodsTwenty-four male specific pathogen-free (SPF) BALB/c mice were randomly divided into a control group and a model group, with 12 mice per group. Mice in the model group were placed in a hypobaric oxygen chamber control system to establish the HPH model, whereas mice in the control group received no intervention. After 28 days of modeling, a comprehensive evaluation of the pathophysiological characteristics of both groups was performed using a general condition scoring scale, echocardiography, hemodynamic measurements, blood gas analysis, hematological tests, organ coefficient determination, histopathological examination, and enzyme-linked immunosorbent assay (ELISA). In addition, correlation analyses were conducted among echocardiographic parameters, the contents of endothelin‑1 (ET‑1) and N‑terminal pro‑B‑type natriuretic peptide (NT‑proBNP) in lung tissue, and other measured indices. ResultsWith prolonged modeling duration time, body weight, water intake, and activity level of mice in the model group were significantly decreased compared with the control group (P<0.05). Echocardiography showed that, compared with the control group, the model group exhibited increased right ventricular dimensions (P<0.01), right ventricular anterior wall thickening (P<0.01), and widened main pulmonary artery diameter (P<0.01), whereas the peak systolic velocity across the pulmonary valve, tricuspid annular plane systolic excursion, and peak systolic velocity of the tricuspid annulus were significantly decreased (P<0.001); tricuspid regurgitation was observed in some model animals. Hemodynamic results revealed that right ventricular systolic pressure was elevated in the model group compared with the control group (P<0.001). Blood gas analysis showed that pH, partial pressure of oxygen, oxygen saturation, actual bicarbonate, and total carbon dioxide were all decreased in the model group compared with the control group (P<0.01). Hematological indices demonstrated that lymphocyte counts (P<0.05) and reticulocyte counts (P<0.001) were decreased in the model group compared with the control group. Compared with the control group, the organ coefficients of heart and lung in the model group were significantly increased (P<0.05 and P<0.001). Pathological examination revealed that the right ventricular hypertrophy index was significantly increased in the model group (P<0.001), with varying degrees of damage to the right ventricle, pulmonary artery, and pulmonary vessels; both the pulmonary artery wall thickness percentage and the pulmonary wall area percentage were significantly elevated (P<0.001). ELISA results showed that the levels of ET‑1 and NT‑proBNP in lung tissues were significantly increased in the model group compared with the control group (P<0.05 and P<0.001, respectively). Correlation analysis indicated that some echocardiographic parameters were highly correlated with multiple indices in the development of HPH (P<0.05). ConclusionEchocardiography can accurately assess a series of hemodynamic changes in HPH, including right ventricular structural, functional impairment, and pulmonary hypertension. Laboratory tests not only help verify whether the model has been successfully established, but also provide deeper insights into the pathogenesis of HPH, evaluate the effects of interventions, and offer scientific evidence for clinical outcomes. Pathological examination can further confirm the alterations in pulmonary vascular remodeling and increased right heart load. This multimodal analysis provides a reliable animal model and evaluation paradigm for both basic and translational research on HPH, and is of great significance for exploring disease mechanisms and developing precision therapeutic strategies.

Background: Alzheimer’s disease (AD) is the most prevalent neurodegenerative disorder. Cellular drug screening represents a key approach in the preclinical discovery of AD therapeutics. Moreover, selecting an appropriate cellular model is essential for exploring the mechanisms underlying AD pathogenesis. Objective: This study aimed to establish a progressive strategy integrating cell model construction, monomer evaluation, and effective component group-based drug discovery for the development of novel therapeutics targeting AD. Methods: A novel AD12 cell line stably expressing the double-transgenic genes amyloid precursor protein (APP)695/presenilin 1(PS1)-ΔE9 was first constructed and subjected to lipopolysaccharide (LPS) stimulation. The resulting LPS-induced neuroinflammatory response was characterized using metabolomic analysis. Subsequently, 16 alkaloids derived from traditional Chinese medicine(TCM) were screened with the LPS-AD12 cell model. Molecular docking was then performed to assess the binding affinity of benzoylaconitine (BAC) with APP and PS1 proteins. A targeted metabolomics approach focusing on amino acids and the tricarboxylic acid (TCA) cycle was applied to BAC-treated LPS-AD12 cells. Finally, aconite extract was administered to five-familial-AD-mutations(5xFAD) transgenic mice, and therapeutic efficacy was evaluated through phenotypic and metabolomic analyses. Results: LPS stimulation induced neuroinflammation in AD12 cells, characterized by elevated levels of interleukin-6, amyloid-β (Aβ ), and PS1. The primary metabolic pathways affected in AD12 cells included alanine, aspartate, and glutamate metabolism, as well as arginine biosynthesis. Notably, BAC significantly reduced the levels of Aβ and interleukin-6 in LPS-AD12 cells. Furthermore, molecular docking analyses revealed a strong binding affinity of BAC with APP and PS1 proteins. Targeted metabolomics analysis indicated that the pharmacological activity of BAC was primarily associated with its anti-inflammatory effects, modulation of the cholinergic signaling pathway, and regulation of energy metabolism. Behavioral experiments in 5xFAD mice, together with immunofluorescence assessments of brain tissues (amyloid-β, glial fibrillary acidic protein, and ionized calcium-binding adapter molecule 1), demonstrated that aconite extract treatment significantly alleviated AD-related symptoms. In addition, metabolomics-based mechanistic analyses revealed that serine, acetylcholine, ornithine, and gamma-aminobutyric acid exhibited consistent trends in both brain tissue and serum samples from treated mice. Conclusion: This study established a novel AD12 cell model and an effective, reliable LPS-induced cell-based screening system. Additionally, we proposed a new drug discovery strategy integrating cell model construction, monomer evaluation, and effective component group-based drug discovery. This screening system has the potential to accelerate the identification of anti-AD agents and provides a feasible framework for subsequent in-depth investigations into new drug development from TCM, thereby contributing to the modernization of TCM-derived therapeutics.

Objective To investigate the relationship between the levels of soluble FMS-like tyrosine ki-nase-1(sFlt-1),monocyte chemoatgulant protein-1(MCP-1),and soluble lectin-like oxidized low-density lipo-protein receptor 1(sLOX-1)in the serum of patients with coronary heart disease and the disease state and prognosis.Methods A total of 126 patients with coronary heart disease treated in Zhongshan Hospital,Dalian University from May 2022 to May 2024 were selected as study objects.The patients were divided into good prognosis group(n=78)and poor prognosis group(n=48)according to prognostic results.The serum sFlt-1,MCP-1 and sLOX-1 levels were detected by enzyme-linked immunosorbent assay.The degree of coronary artery stenosis in patients was measured by coronary angiography,as indicated by the Gensini score.Spearman correlation analysis was used to analyze the correlation between Gensini score and serum sFlt-1,MCP-1 and sLOX-1 levels.Multivariate Logistic regression was used to analyze the factors affecting the prognosis of pa-tients with coronary heart disease.The predictive efficacy of serum sFlt-1,MCP-1 and sLOX-1 levels for the prognosis of patients with coronary heart disease was evaluated by receiver operating characteristic(ROC)curve.Results The age in the poor prognosis group was older than that in the good prognosis group(P＜0.05),Gensini score was higher than that in the good prognosis group(P＜0.05).Compared with the poor prognosis group,the serum sFlt-1,MCP-1 and sLOX-1 levels in the good prognosis group were lower(P＜0.05).Gensini score was positively correlated with the serum sFlt-1,MCP-1 and sLOX-1 levels.Multivariate Logistic regression analysis showed that old age,high Gensini score,high sFlt-1 level,high MCP-1 level and high sLOX-1 level were all risk factors effecting the prognosis of patients with coronary heart disease(P＜0.05).ROC curve analysis showed that the areas under the curve(AUC)of serum sFlt-1,MCP-1,and sLOX-1 for predicting the prognosis of coronary heart disease patients were 0.787(95％CI:0.703-0.870),0.815(95％CI:0.741-0.890)and 0.795(95％CI:0.715-0.876),respectively.The AUC for predicting the prognosis of coronary heart disease patients using a combination of three was 0.923(95％CI:0.876-0.970).Conclusion The serum levels of sFlt-1,MCP-1 and sLOX-1 in patients with coronary heart disease are signifi-cantly increased,and are positively correlated with the degree of coronary artery stenosis.The levels of sFlt-1,MCP-1 and sLOX-1 in serum have a certain predictive effect on the prognosis of patients with coronary heart disease.

Objective To explore macrophage subpopulations in ischemic stroke (IS) by using single-cell RNA sequencing (scRNA-seq) data analysis and High-Dimensional Weighted Gene Co-Expression Network Analysis (hdWGCNA). Methods Based on single-cell sequencing data, transcriptomic information for different cell types was obtained, and macrophages were selected for subpopulation identification. hdWGCNA, cell-cell communication, and pseudotime trajectory analysis were used to explore the characteristics of macrophage subpopulations following IS. Key genes related to IS were identified using microarray data and validated for diagnostic potential through Receiver Operating Characteristic (ROC) analysis. Gene Set Enrichment Analysis (GSEA) was conducted to investigate the potential functions of these genes. Results The scRNA-seq data analysis revealed significant changes in macrophage subpopulation composition after IS. A specific macrophage subpopulation enriched in the stroke group was identified and designated as MCAO-specific macrophages (MSM). Pseudotime trajectory analysis indicated that MSM cells were in an intermediate stage of macrophage differentiation. Cell-cell communication analysis uncovered complex interactions between MSM cells and other cells, with the CCL6-CCR1 signaling axis potentially playing a crucial role in neuroinflammation. Two gene modules associated with MSM were identified via hdWGCNA, significantly enriched in pathways related to NOD-like receptors and antigen processing. By integrating differentially expressed MSM genes with conventional transcriptomic data, three IS-related hub genes were identified: Arg1, CLEC4D, and CLEC4E. Conclusion This study reveals the characteristics and functions of macrophage subpopulations following IS and identifies three hub genes with potential diagnostic value, providing novel insights into the pathological mechanisms of IS.
Macrophages/metabolism* ; Computational Biology/methods* ; Single-Cell Analysis/methods* ; Transcriptome ; Ischemic Stroke/metabolism* ; Animals ; Gene Regulatory Networks ; Gene Expression Profiling ; Humans ; Male

ObjectiveTo study the effect of Tanreqing injection （TRQ） on the pulmonary flora in the rat model of chronic obstructive pulmonary disease （COPD）. MethodWistar rats were randomized into control， model， and TRQ groups. The rats in other groups except the control group were treated by smoking combined with intratracheal instillation of lipopolysaccharide for the modeling of COPD. The TRQ group was intraperitoneally injected with TRQ （2 g·kg^-1）. At the end of the experiment， after blood collection from the abdominal aorta of the rats， the lung tissue was collected for hematoxylin-eosin and picric sirius red staining to reveal the pathological changes. The lung lavage fluid was collected， and the diversity and relative abundance of lung flora in different groups were analyzed by 16S rDNA amplicon sequencing. ResultThe lungs of the control group were normal， and those of the model group showed neutrophil infiltration， telangiectasia， lung hemorrhage and emphysema in individual cases， and thickening of collagen fibers in the trachea. Compared with the model group， the TRQ group showed significantly improved lungs and recovered collagen fibers. The MLI analysis showed that compared with the control group， the model group showcased increased alveolar space （P<0.01）， which was reduced in the TRQ group （P<0.05）. Compared with the control group， the model group showed increased wall thickness （P<0.01）， and the increase was attenuated in the TRQ group （P<0.01）. TRQ increased the Simpson index and altered the α diversity of pulmonary flora. The results of principal co-ordinate analysis showed that TRQ changed the β diversity and reduced the β diversity index of pulmonary flora. At the genus level， the model group showed increased relative abundance of g_Bacillus and g_Brevundimonas and decreased relative abundance of g_Pseudomonas， compared with the control group. After treatment with TRQ， the relative abundance of g_Stenotrophomonas increased， and that of g_Bacillus decreased. The LEfSe of differential taxa between groups showed that the modeling increased the relative abundance of g_Lachnospiraceae_NK4A136_group， and TRQ treatment increased the relative abundance of g_Rhodococcus and g_Stenotrophomonas. ConclusionTRQ can regulate the diversity of pulmonary flora and restore the balance of bacterial genera in the rat model of COPD， which may be one of the mechanisms of the prevention and treatment of COPD with TRQ.

Background It is a research hotspot to study the changes of metabolites and metabolic pathways in the process of coal worker's pneumoconiosis (CWP) by metabonomics and to explore its pathogenesis. Objective To study the change of metabolites in bronchoalveolar lavage fluid (BALF) of patients with CWP and explore the metabolic regulation mechanism of the disease. Methods Patients with CWP who met the national diagnostic criteria according to Diagnosis of occupational pneumoconiosis (GBZ 70-2015) and underwent massive whole lung lavage were selected as the case group, and patients with tracheostenosis who underwent bronchoscopy were selected as the control group. BALF samples were collected from the cases and the controls. After filtering out large particles and mucus, the supernatant was stored in a −80 ℃ refrigerator. The samples were detected and analyzed by liquid chromatography-mass spectrometry after adding extraction solution, cold bath ultrasonication, and high-speed centrifugation, and the metabolic profiles and related data of CWP patients were obtained. The differential metabolites related to the occurrence and development of CWP were screened by multiple statistical analysis; furthermore, we searched the Kyoto Encyclopedia of Genes and Genomes (KEGG) database for potential metabolic pathways involved in the progression. Results There was no significant difference in the general conditions of the subjects, such as weight, height, age, and length of service among the stage I group, the stage II group, the stage III group, and the control group (P˃0.05). When comparing the CWP stage I group with the control group, 48 differential metabolites were screened out, among which 14 were up-regulated and 34 were down-regulated. A total of 66 differential metabolites were screened out between the patients with CWP stage II and the controls, 14 up-regulated and 52 down-regulated differential metabolites. Compared with the control group, 63 differential metabolites were screened out in the patients with CWP stage III, including 11 up-regulated and 52 down-regulated differential metabolites. There were 36 differential metabolites that may be related to the occurrence of CWP, among which 11 differential metabolites were up-regulated, and 25 were down-regulated. Four significant differential metabolic pathways were identified through KEGG database query: linoleic acid metabolic pathway, alanine metabolic pathway, sphingolipid metabolic pathway, and glycerophospholipid metabolic pathway. Conclusion The metabolomic study of BALF show that there are 36 different metabolites in the occurrence and development of CWP, mainly associating with linoleic acid metabolism, alanine metabolism, sphingolipid metabolism, and glycerophospholipid metabolism pathways.

Background Pneumoconiosis is a widespread occupational disease in China at present. As a type of lung diseases, its pathological damage is mainly irreversible fibrotic changes in the lungs. Several studies have shown that the occurrence and development of lung diseases such as coal workers' pneumoconiosis are closely related to intestinal flora. Objective To observe intestinal flora of coal workers' pneumoconiosis patients based on the results of 16SrDNA high-throughput sequencing and evaluate the changes of intestinal flora after treatment with tetrandrine tablets. Methods A total of 80 patients with coal workers' pneumoconiosis attending the outpatient clinic of the Department of Occupational Diseases of the Emergency General Hospital from April to July 2022 were enrolled. All patients were treated with tetrandrine tablets for 4 weeks, with group A before the treatment of tetrandrine tablets and group B after the treatment. In the same period, 24 healthy controls (group C) were set up. Stool samples were collected before and after the treatment. Using 16SrDNA high-throughput sequencing, gene V3-V4 sequencing technology, and bioinformatic analysis platform, we evaluated the intestinal flora after treatment by groups. Results The dominant flora at the phylum level and genus level were the same across three groups. The relative abundances of phylum Bacteroidetes, Bifidobacterium, Bacteroides, and Facealibacterium in groups B and C were higher than those in group A, and the relative abundances of phy-lum Actinobacteria, genus Blautia, and genus Romboutsia in groups B and C were lower than those in group A (P＜0.05). The relative abundances of genus Clostridium, genus Megamonas, and genus Lactobacillus in group C was lower than that in groups A and B (P＜0.05). The alpha diversity analysis showed that the Chao1 index was higher in group A than in group C (P＜0.01). Compared with group A, the Shannon index was higher in group B, and the increases of Simpson index were all statistically significant in stage I patients (P＜0.05), but the differences in Chao1 index were not statistically significant (P＞0.05). The differences in the values of Chao1 index, Shannon index, and Simpson index in stage Ⅱ and stage III patients were not statistically significant (P＞0.05). The beta diversity analysis showed that the difference in flora structure between group A and group C was statistically significant (P＜0.05); the differences in flora structure before and after treatment in the same stage patients were statistically significant (P＜0.05). The partial least squares discriminant analysis (PLS-DA) showed that there were significant differences between group A and group C, and between group A and group B. The LEfSe analysis showed that the significant markers contributing to the differences were basically the same in stage I, stage Ⅱ, and stage Ⅲ after treatment, which were mainly phylum Bacteroidetes and its subordinate groups, class Negativicutes, or-der Selenomonas, and genus Facealibacterium. Conclusion There are differences in the distribution of flora between coal workers' pneumoconiosis patients and healthy individuals, and the structure and relative abundance of intestinal flora are changed and the number of beneficial flora is increased after treatment with tetrandrine tablets.

Background At present, the practice of pulmonary rehabilitation for pneumoconiosis in China is in a primary stage. The basis for formulating an individualized comprehensive pulmonary rehabilitation plan is still insufficient, which is one of the factors limiting the development of community-level rehabilitation work. Objective To formulate an exercise prescription based on maximum heart rate measured by cardiopulmonary exercise test (CPET), conduct an individualized comprehensive pulmonary rehabilitation program with the exercise prescription for patients with stable pneumoconiosis, and evaluate its role in improving exercise endurance and quality of life, thus provide a basis for the application and promotion of pulmonary rehabilitation. Methods A total of 68 patients were recruited from the Occupational Disease Prevention Hospital of Jinneng Holding Coal Industry Group Co., Ltd. from April to August 2022 , and were divided into an intervention group and a control group by random number table method, with 34 cases in each group. All the pneumoconiosis patients participated in a baseline test. The control group was given routine drug treatment, while the intervention group received multidisciplinary comprehensive pulmonary rehabilitation treatment on the basis of routine drug treatment, including health education, breathing training, exercise training, nutrition guidance, psychological intervention, and sleep management, whose exercise intensity was determined according to the maximum heart rate provided by CPET. The rehabilitation training lasted for 24 weeks. Patients were evaluated at registration and the end of study respectively. CPET was used to measure peak oxygen uptake per kilogram (pVO₂/kg), anaerobic threshold (AT), carbon dioxide equivalent of ventilation (EqCO₂), maximum metabolic equivalent (METs), and maximum work (W_max). The modified British Medical Research Council Dyspnea Questionnaire (mMRC), Self-rating Anxiety Scale (SAS), Self-rating Depression Scale (SDS), Pittsburgh Sleep Quality Index (PSQI), Chronic Obstructive Pulmonary Disease Assessment Test (CAT), and Short Form of Health Survey (SF-36) were used to evaluate the potential effect of the comprehensive pulmonary rehabilitation program. Results Among the included 68 patients, 63 patients were having complete data, then 31 cases were assigned in the control group and 32 cases in the interventional group. Before the intervention, there was no significant difference in pVO₂/kg, AT, EqCO₂, METs, or W_max between the two groups (P>0.05). At the end of the trail, the indicators like pVO₂/kg [(19.81±2.38) mL·(min·kg)⁻¹], AT [(14.48±2.33) mL·(min·kg)⁻¹], METs (5.64±0.69), and W_max [(85.25±14) W] of patients in the intervention group were all higher than those [(13.90±2.37) mL·(min·kg)⁻¹, (11.70±1.94) mL·(min kg)⁻¹, (3.97±0.70), and (61.77±14.72) W, respectively] in the control group (P<0.001); there was no significant difference in EqCO₂ between the two groups (P=0.083). Before the trial, there was no significant difference in mMRC, SAS, SDS, PSQI, or CAT scores between the two groups (P>0.05). At the end of the trail, the mMRC score (1.16±0.57), SAS score (27.93±2.12), SDS score (26.48±1.44), PSQI score (1.08±0.88), and CAT score (4.34±3.28) of patients in the intervention group were lower than those [(2.03±0.83), (35.87±6.91), (34.23±6.65), (5.37±3.03), and (13.87±7.53), respectively] in the control group (P<0.001). The SF-36 scores of bodily pain (94.13±10.72), general health (87.50±5.68), vitality (95.31±5.53), mental health (99.88±0.71), and health changes (74.22±4.42) in the intervention group were higher than those [(71.87±32.72), (65.81±15.55), (74.52±16.45), (86.97±16.56), and (29.84±13.50), respectively] in the control group (P<0.001), and no significant difference was found in social functioning and role emotional scores (P>0.05). Conclusion Comprehensive pulmonary rehabilitation can increase the oxygen intake and exercise endurance of pneumoconiosis patients, ameliorate dyspnea symptoms, elevate psychological state and sleep quality, and improve the quality of life.