Identification of Potential Differentially-Methylated/Expressed Genes in Chronic Obstructive Pulmonary Disease

Transcriptomic Insights into Different Stimulation Intensity of Electroacupuncture in Treating COPD in Rat Models

Background

Electroacupuncture (EA), with varying stimulation intensities, has demonstrated therapeutic potentials in both animal and clinical studies for the treatment of chronic obstructive pulmonary disease (COPD). However, a comprehensive investigation of the intensity-related effects, particularly 1mA and 3mA of EA, and the underlying mechanisms remains lacking.

Methods

A COPD rat model was established by prolonged exposure to cigarette smoke and intermittent intratracheal instillation of lipopolysaccharide. EA treatment was administered at acupoints BL13 (Feishu) and ST36 (Zusanli), 20 minutes daily for 2 weeks, with intensities of 1mA and 3mA. EA effectiveness was evaluated by pulmonary function, histopathological change, serum level of inflammatory cytokines, and level of oxidative stress markers in serum and lung tissues. Transcriptome profiling and weighted gene co-expression network analysis (WGCNA) were performed to reveal gene expression patterns and identify hub genes. Real-time quantitative PCR (RT-qPCR) and Western blot (WB) were performed to detect the mRNA and protein expression levels, respectively.

Results

EA at both 1mA and 3mA exerted differing therapeutic effects by improving lung function and reducing inflammation and oxidative stress in COPD rats. Transcriptome analysis revealed distinct expression patterns between the two groups, functionally corresponding to shared and intensity-specific (1mA and 3mA) enriched pathways. Eight candidate genes were identified, including Aqp9, Trem1, Mrc1, and Gpnmb that were downregulated by EA and upregulated in COPD. Notably, Msr1 and Slc26a4 exclusively downregulated in EA-1mA, while Pde3a and Bmp6 upregulated solely in EA-3mA. WGCNA constructed 5 key modules and elucidated the module-trait relationship, with the aforementioned 8 genes being highlighted. Additionally, their mRNA and protein levels were validated by RT-qPCR and WB.

Conclusion

Our results demonstrated that 1mA and 3mA intensities induce distinct gene expression patterns at the transcriptional level, associated with shared and 1mA vs 3mA-specific enriched pathways. Genes Mrc1, Gpnmb, Trem1, and Aqp9 emerge as promising targets, and further studies are needed to elucidate their functional consequences in COPD.

Different Transcriptome Features of Peripheral Blood Mononuclear Cells in Non-Emphysematous Chronic Obstructive Pulmonary Disease

Non-emphysematous chronic obstructive pulmonary disease (COPD), which is defined based on chest computed tomography findings, presented different transcriptome features of peripheral blood mononuclear cells (PBMCs) compared with emphysematous COPD. Enrichment analysis of transcriptomic data in COPD demonstrated that the "Hematopoietic cell lineage" pathway in Kyoto Encyclopedia of Genes and Genomes pathway analysis was highly upregulated, suggesting that cellular dynamic dysregulation in COPD lungs is affected by pathologically modified PBMCs. The differentially expressed genes (DEGs) upregulated in PBMCs reflected the disease state of non-emphysematous COPD. Upregulated DEGs such as XCL1, PRKCZ, TMEM102, CD200R1, and AQP1 activate T lymphocytes and eosinophils. Upregulating keratan sulfate biosynthesis and metabolic processes is associated with protection against the destruction of the distal airways. ITGA3 upregulation augments interactions with extracellular matrix proteins, and COL6A1 augments the profibrotic mast cell phenotype during alveolar collagen VI deposition. Upregulating HSPG2, PDGFRB, and PAK4 contributes to the thickening of the airway wall, and upregulating SERPINF1 expression explains the better-preserved vascular bed. Therefore, gene expression and pathway analysis in PBMCs in patients with non-emphysematous COPD represented type 2 immune responses and airway remodeling features. Therefore, these patients have asthmatic potential despite no clinical signs of asthma, in contrast to those with emphysematous COPD.