Identification of Potential Key Genes in the Pathogenesis of Chronic Obstructive Pulmonary Disease Through Bioinformatics Analysis

Exosome microRNA-125a-5p derived from epithelium promotes M1 macrophage polarization by targeting IL1RN in chronic obstructive pulmonary disease

BACKGROUND

The interplay between airway epithelium and macrophages plays a pivotal role in Chronic Obstructive Pulmonary Disease (COPD) pathogenesis. Exosomes, which transport miRNA cargo, have emerged as novel mediators of intercellular communication. MicroRNA-125a-5p (miR-125a-5p) has been implicated in macrophage polarization.This study aims to investigate the role of exosomal miR-125a-5p in the dysfunctional epithelium-macrophage cross-talk in cigarette smoke (CS)-induced COPD.

METHODS

In cell models, THP-1 monocytic cells were differentiated into macrophages (M0). Human bronchial epithelial cells treated with CS extract (CSE) were co-cultured with M0. Exosomes were isolated from culture media using commercial kits and characterized using nanoparticle tracking analysis (NTA) and transmission electron microscopy (TEM). Exosomes labeled with PKH26 red fluorescent cell linker kits were incubated with macrophages. Luciferase reporter assay was used to confirm the target gene of miR-125a-5p. In mouse experiments, inhibiting miR-125a-5p was utilized to examine its role in macrophage polarization. Furthermore, the underlying mechanism was explored.

RESULTS

In vitro results indicated that CSE treatment led to upregulation of miR-125a-5p in HBE cells, and exosomes contained miR-125a-5p. PKH26-labeled exosomes were internalized by macrophages. Co-culture experiments between bronchial epithelial cells and miR-125a-5p mimic resulted in significant increase in M1 macrophage markers (TNF-α, iNOS-2, IL-1β) and decrease in M2 markers (IL-10 and Arg-1). In COPD mouse models, miR-125a-5p inhibitor reduced levels of TNF-α, IL-1β, and IL-6. Luciferase assays revealed that miR-125a-5p inhibitors enhanced the relative luciferase activity of IL1RN. Mechanistic experiments demonstrated that HBE-derived exosomes transfected with miR-125a-5p mimics promoted upregulation of MyD88, TRAF6, p65, iNOS-2, and downregulation of Arg-1.

CONCLUSION

This study suggests that exosomal miR-125a-5p may act as a mediator in the cross-talk between airway epithelium and macrophage polarization in COPD. Exosomal miR-125a-5p targeting IL1RN may promote M1 macrophage polarization via the MyD88/NF-κB pathway.

Combined analysis of bulk RNA and single-cell RNA sequencing to identify pyroptosis-related markers and the role of dendritic cells in chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is characterized by dyspnea caused by airflow limitation. Further development may lead to decreased lung function and other lung diseases. Pyroptosis is a type of programmed cell death that involves multiple pathways. For example, the pathway induced by the NLR family pyrin domain containing 3 (NLRP3) inflammasome is closely associated with COPD exacerbation. Therefore, in this study, various machine learning algorithms were applied to screen for diagnostically relevant pyroptosis-related genes from the GEO dataset, and the results were verified using external datasets. The results showed that deep neural networks and logistic regression algorithms had the highest AUC of 0.91 and 0.74 in the internal and external test sets, respectively. Here, we explored the immune landscape of COPD using diagnosis-related genes. We found that the infiltrating abundance of dendritic cells significantly differed between the COPD and control groups. Finally, the communication patterns of each cell type were explored based on scRNA-seq data. The critical role of significant pathways involved in communication between DCS and other cell populations in the occurrence and progression of COPD was identified.

Identification of genes and key pathways underlying the pathophysiological association between sarcopenia and chronic obstructive pulmonary disease

PURPOSE

Chronic obstructive pulmonary disease (COPD) patients are likely to develop sarcopenia, while the exact mechanism underlying the association between sarcopenia and COPD is still not clear. This cohort study aims to explore the genes, signaling pathways, and transcription factors (TFs) that are related to the molecular pathogenesis of sarcopenia and COPD.

METHODS

According to the strict inclusion criteria, two gene sets (GSE8479 for sarcopenia and GSE76925 for COPD) were obtained from the Gene Expression Omnibus (GEO) platform. Overlapping differentially expressed genes (DEGs) in sarcopenia and COPD were detected, and comprehensive bioinformatics analysis was conducted, including functional annotation, enrichment analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), construction of a protein-protein interaction (PPI) network, co-expression analysis, identification and validation of hub genes, and TFs prediction and verification.

RESULTS

In total, 118 downregulated and 92 upregulated common DEGs were detected. Functional analysis revealed that potential pathogenesis involves oxidoreductase activity and ferroptosis. Thirty hub genes were detected, and ATP metabolic process and oxidative phosphorylation were identified to be closely related to the hub genes. Validation analysis revealed that SAA1, C3, and ACSS2 were significantly upregulated, whereas ATF4, PPARGC1A, and MCTS1 were markedly downregulated in both sarcopenia and COPD. In addition, six TFs (NFKB1, RELA, IRF7, SP1, MYC, and JUN) were identified to regulate the expression of these genes, and SAA1 was found to be coregulated by NFKB1 and RELA.

CONCLUSION

This study uncovers potential common mechanisms of COPD complicated by sarcopenia. The hub gene SAA1 and the NF-κB signaling pathway could be involved, and oxidative phosphorylation and ferroptosis might be important contributors to this comorbidity.

The Effects of Ellagic Acid on Experimental Corrosive Esophageal Burn Injury

This study aimed to investigate the antioxidant effect of Ellagic acid (EA) on wound healing in sodium hydroxide (NaOH)-induced corrosive esophageal burn injury. The interaction networks and functional annotations were conducted using Cytoscape software. A total of 24 Wistar albino rats were divided into control, corrosive esophageal burn (CEB) and CEB + EA groups. Burn injury was created by 20% NaOH and 30 mg/kg EA was per oral administered to rats. At the end of the 28-day experimental period, Malondialdehyde (MDA) content was measured. Esophageal tissue samples were processed for histological staining. The EA-target interaction network was revealed to be involved in regulating crucial cellular mechanisms for burn wound healing, with epidermal growth factor (EGF) identified as a central mediator. An increase in animal weight in the CEB + EA group was observed in the EA-treated group after CEB injury. Burn injury increased MDA content, but EA treatment decreased its level after CEB injury. Stenosis index, collagen degeneration, inflammation, fibrosis and necrosis levels were increased after CEB injury. EA treatment improved histopathology in the CEB + EA group compared to the CEB group. The expression of EGF was decreased in the CEB group but upregulated in the EA-treated group, suggesting a potential involvement of EA in cellular processes and tissue regeneration. EA, through its antioxidative and tissue regenerative properties, significantly contributes to alleviating the adverse effects of CEB injury, promoting wound healing.

Comprehensive Analysis of a Competing Endogenous RNA Co-Expression Network in Chronic Obstructive Pulmonary Disease

Purpose

Chronic obstructive pulmonary disease (COPD) is the main cause of mortality world widely. Non-coding RNAs (lncRNAs) and associated competitive endogenous RNAs (ceRNAs) networks were recently proved to lead to mRNA gene expression downregulation but were still unclear in COPD. This study aims to investigate and elucidate the mechanisms underlying the involvement of ceRNA co-expression networks in COPD pathogenesis.

Methods

Obtained expression signature of data from the Gene Expression Omnibus database and compared the differentially expression of mRNAs and miRNAs between COPD patients and healthy smokers. Predicted the miRNA-lncRNA and miRNA-mRNA interaction using online library and employed CIBERSORT to measure the proportions of the 22 immune cells in the COPD and control groups.

Results

Established a ceRNA-network comprising 11 lncRNAs, 5 miRNAs, and 16 mRNAs. Using the weighted correlation network analysis method, we identified hub genes and hub miRNAs and obtained one core sub-network, XIST, FGD5-AS1, KCNQ1OT1, HOXA11-AS, LINC00667, H19, PRKCQ-AS1, NUTM2A-AS1/has-mir-454-3p/ZNF678, PRRG4. COPD patients had different proportions of immune cells than controls, and these variations were associated with the magnitude of pulmonary function parameters.

Conclusion

The ceRNA-network, particularly the core sub-network, may be a putative goal for COPD, in which specific immune cells were involved.

Integrated bioinformatics analysis for the identification of idiopathic pulmonary fibrosis-related genes and potential therapeutic drugs

OBJECTIVE

The pathogenesis of idiopathic pulmonary fibrosis (IPF) remains unclear. We sought to identify IPF-related genes that may participate in the pathogenesis and predict potential targeted traditional Chinese medicines (TCMs).

METHODS

Using IPF gene-expression data, Wilcoxon rank-sum tests were performed to identify differentially expressed genes (DEGs). Protein-protein interaction (PPI) networks, hub genes, and competitive endogenous RNA (ceRNA) networks were constructed or identified by Cytoscape. Quantitative polymerase chain reaction (qPCR) experiments in TGF-β1-induced human fetal lung (HFL) fibroblast cells and a pulmonary fibrosis mouse model verified gene reliability. The SymMap database predicted potential TCMs targeting IPF. The reliability of TCMs was verified in TGF-β1-induced MRC-5 cells.

MATERIALS

Multiple gene-expression profile data of normal lung and IPF tissues were downloaded from the Gene Expression Omnibus database. HFL fibroblast cells and MRC-5 cells were purchased from Wuhan Procell Life Science and Technology Co., Ltd. (Wuhan, China). C57BL/12 mice were purchased from Beijing Vital River Laboratory Animal Technology Co., Ltd. (Beijing, China).

RESULTS

In datasets GSE134692 and GSE15197, DEGs were identified using Wilcoxon rank-sum tests (both p < 0.05). Among them, 1885 DEGs were commonly identified, and 87% (1640 genes) had identical dysregulation directions (binomial test, p < 1.00E-16). A PPI network with 1623 nodes and 8159 edges was constructed, and 18 hub genes were identified using the Analyze Network plugin in Cytoscape. Of 18 genes, CAV1, PECAM1, BMP4, VEGFA, FYN, SPP1, and COL1A1 were further validated in the GeneCards database and independent dataset GSE24206. ceRNA networks of VEGFA, SPP1, and COL1A1 were constructed. The genes were verified by qPCR in samples of TGF-β1-induced HFL fibroblast cells and pulmonary fibrosis mice. Finally, Sea Buckthorn and Gnaphalium Affine were predicted as potential TCMs for IPF. The TCMs were verified by qPCR in TGF-β1-induced MRC-5 cells.

CONCLUSION

This analysis strategy may be useful for elucidating novel mechanisms underlying IPF at the transcriptome level. The identified hub genes may play key roles in IPF pathogenesis and therapy.

Germline variants of the genes involved in NF-kB activation are associated with the risk of COPD and lung cancer development

Chronic obstructive pulmonary disease (COPD) and lung cancer (LC) are closely related diseases associated with smoking history and dysregulated immune response. However, not all smokers develop the disease, indicating that genetic susceptibility could be important. Therefore, the aim of this study was to search for the potential overlapping genetic biomarkers, with a focus on single nucleotide polymorphisms (SNPs) located in the regulatory regions of immune-related genes. Additionally, the aim was to see if an identified SNP has potentially an effect on proinflamma-tory cytokine concentration in the serum of COPD patients. We extracted summary data of variants in 1511 immune-related genes from COPD and LC genome-wide association studies (GWAS) from the UK Biobank. The LC data had 203 cases, patients diagnosed with LC, and 360 938 controls, while COPD data had 1 897 cases and 359 297 controls. Assuming 1 association/gene, SNPs with a p-value < 3.3 × 10^-5 were considered statistically significantly associated with the disease. We identified seven SNPs located in different genes (BAG6, BTNL2, TNF, HCP5, MICB, NCR3, ABCF1, TCF7L1) to be associated with the COPD risk and two with the LC risk (HLA-C, HLA-B), with statistical significance. We also identified two SNPs located in the IL2RA gene associated with LC (rs2386841; p = 1.86 × 10^-4) and COPD (rs11256442; p = 9.79 × 10^-3) but with lower significance. Functional studies conducted on COPD patients showed that RNA expression of IL2RA, IFNγ and related proinflammatory cytokines in blood serum did not correlate with a specific genotype. Although results presented in this study do not fully support our hypothesis, it is worth to mention that the identified genes/SNPs that were associated with either COPD or LC risk, all were involved in the activation of the NF-κB transcription factor which is closely related to the regulation of the inflammatory response, a condition associated with both pathologies.

Identification of candidate aberrant differentially methylated/expressed genes in asthma

BACKGROUND

Asthma is an important non-communicable disease worldwide. DNA methylation is associated with the occurrence and development of asthma. We are aimed at assuring differential expressed genes (DEGs) modified by aberrantly methylated genes (DMGs) and pathways related to asthma by integrating bioinformatics analysis.

METHODS

One mRNA dataset (GSE64913) and one gene methylation dataset (GSE137716) were selected from the Gene Expression Omnibus (GEO) database. Functional enrichment analysis was performed using GeneCodies 4.0 database. All gene expression matrices were analyzed by Gene set enrichment analysis (GSEA) software. STRING was applied to construct a protein-protein interaction (PPI) network to find the hub genes. Then, electronic validation was performed to verify the hub genes, followed by the evaluation of diagnostic value. Eventually, quantitative real-time polymerase chain reaction (qRT-PCR) was utilized to detect the expression of hub genes.

RESULTS

In total, 14 hypomethylated/high-expression genes and 10 hypermethylated/low-expression genes were obtained in asthma. Among them, 10 hub genes were identified in the PPI network. Functional analysis demonstrated that the differentially methylated/expressed genes were primarily associated with the lung development, cytosol and protein binding. Notably, HLA-DOA was enriched in asthma. FKBP5, WNT5A, TM4SF1, PDK4, EPAS1 and GMPR had potential diagnostic value for asthma.

CONCLUSION

The project explored the pathogenesis of asthma, which may provide a research basis for the prediction and the drug development of asthma.