Immune infiltration landscape and immune-marker molecular typing of pulmonary fibrosis with pulmonary hypertension

Bioinformatics analysis of the immune cell infiltration characteristics and correlation with crucial diagnostic markers in pulmonary arterial hypertension

BACKGROUND

Pulmonary arterial hypertension (PAH) is a pathophysiological syndrome, characterized by pulmonary vascular remodeling. Immunity and inflammation are progressively recognized properties of PAH, which are crucial for the initiation and maintenance of pulmonary vascular remodeling. This study explored immune cell infiltration characteristics and potential biomarkers of PAH using comprehensive bioinformatics analysis.

METHODS

Microarray data of GSE117261, GSE113439 and GSE53408 datasets were downloaded from Gene Expression Omnibus database. The differentially expressed genes (DEGs) were identified in GSE117261 dataset. The proportions of infiltrated immune cells were evaluated by CIBERSORT algorithm. Feature genes of PAH were selected by least absolute shrinkage and selection operator (LASSO) regression analysis and validated by fivefold cross-validation, random forest and logistic regression. The GSE113439 and GSE53408 datasets were used as validation sets and logistic regression and receiver operating characteristic (ROC) curve analysis were performed to evaluate the prediction value of PAH. The PAH-associated module was identified by weighted gene association network analysis (WGCNA). The intersection of genes in the modules screened and DEGs was used to construct protein-protein interaction (PPI) network and the core genes were selected. After the intersection of feature genes and core genes, the hub genes were identified. The correlation between hub genes and immune cell infiltration was analyzed by Pearson correlation analysis. The expression level of LTBP1 in the lungs of monocrotaline-induced PAH rats was determined by Western blotting. The localization of LTBP1 and CD4 in lungs of PAH was assayed by immunofluorescence.

RESULTS

A total of 419 DEGs were identified, including 223 upregulated genes and 196 downregulated genes. Functional enrichment analysis revealed that a significant enrichment in inflammation, immune response, and transforming growth factor β (TGFβ) signaling pathway. CIBERSORT analysis showed that ten significantly different types of immune cells were identified between PAH and control. Resting memory CD4+ T cells, CD8+ T cells, γδ T cells, M1 macrophages, and resting mast cells in the lungs of PAH patients were significantly higher than control. Seventeen feature genes were identified by LASSO regression for PAH prediction. WGCNA identified 15 co-expression modules. PPI network was constructed and 100 core genes were obtained. Complement C3b/C4b receptor 1 (CR1), thioredoxin reductase 1 (TXNRD1), latent TGFβ binding protein 1 (LTBP1), and toll-like receptor 1 (TLR1) were identified as hub genes and LTBP1 has the highest diagnostic efficacy for PAH (AUC = 0.968). Pearson correlation analysis showed that LTBP1 was positively correlated with resting memory CD4+ T cells, but negatively correlated with monocytes and neutrophils. Western blotting showed that the protein level of LTBP1 was increased in the lungs of monocrotaline-induced PAH rats. Immunofluorescence of lung tissues from rats with PAH showed increased expression of LTBP1 in pulmonary arteries as compared to control and LTBP1 was partly colocalized with CD4+ cells in the lungs.

CONCLUSION

LTBP1 was correlated with immune cell infiltration and identified as the critical diagnostic maker for PAH.

Inhibition of Th17 cells by donepezil ameliorates experimental lung fibrosis and pulmonary hypertension

Rationale: Pulmonary hypertension (PH) secondary to lung fibrosis belongs to WHO Group III, one of the most common subgroups of PH; however, it lacks effective treatment options. Cholinesterase inhibitor donepezil (DON) has been shown to effectively improve Group I PH. However, its effects on Group III PH are unknown. Methods: A lung fibrosis-induced PH mouse model was constructed using a single intratracheal instillation of bleomycin (BLM), after which DON was administered daily. Pulmonary artery and right ventricle (RV) remodeling were evaluated at the end of the study. Lung tissue in each group was analyzed using RNA sequencing, and the results were further verified with datasets from patients with PH. The mechanisms underlying DON-induced effects on PH were verified both in vivo and in vitro. Results: DON effectively improved pulmonary artery and RV remodeling in the BLM-induced mouse model. Transcriptomic profiles of lung tissue indicated that the expression of inflammatory and fibrotic genes was significantly changed in this process. In the animal model and patients with PH, T helper 17 lymphocytes (Th17) were the most common inflammatory cells infiltrating the lung tissue. DON significantly inhibited lung fibroblast activation; thus, preventing lung fibrosis and reducing the inflammatory response and Th17 cell infiltration in the BLM-induced lung tissue. In addition, Th17 cells could activate lung fibroblasts by secreting IL17A, and DON-mediated inhibition of Th17 cell differentiation was found to depend on the α7nAchR-JAK2-STAT3 pathway. Conclusion: DON can alleviate lung fibrosis and PH in an experimental mouse model. It inhibited pro-inflammatory Th17 cell differentiation, which is dependent on a cholinergic receptor pathway, thereby regulating fibroblast activation.

Weighted gene co-expression network analysis identifies dysregulated B-cell receptor signaling pathway and novel genes in pulmonary arterial hypertension

Background

Pulmonary arterial hypertension (PAH) is a devastating cardio-pulmonary vascular disease in which chronic elevated pulmonary arterial pressure and pulmonary vascular remodeling lead to right ventricular failure and premature death. However, the exact molecular mechanism causing PAH remains unclear.

Methods

RNA sequencing was used to analyze the transcriptional profiling of controls and rats treated with monocrotaline (MCT) for 1, 2, 3, and 4 weeks. Weighted gene co-expression network analysis (WGCNA) was employed to identify the key modules associated with the severity of PAH. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to explore the potential biological processes and pathways of key modules. Real-time PCR and western blot analysis were used to validate the gene expression. The hub genes were validated by an independent dataset obtained from the Gene Expression Omnibus database.

Results

A total of 26 gene modules were identified by WGCNA. Of these modules, two modules showed the highest correlation with the severity of PAH and were recognized as the key modules. GO analysis of key modules showed the dysregulated inflammation and immunity, particularly B-cell-mediated humoral immunity in MCT-induced PAH. KEGG pathway analysis showed the significant enrichment of the B-cell receptor signaling pathway in the key modules. Pathview analysis revealed the dysregulation of the B-cell receptor signaling pathway in detail. Moreover, a series of humoral immune response-associated genes, such as BTK, BAFFR, and TNFSF4, were found to be differentially expressed in PAH. Additionally, five genes, including BANK1, FOXF1, TLE1, CLEC4A1, and CLEC4A3, were identified and validated as the hub genes.

Conclusion

This study identified the dysregulated B-cell receptor signaling pathway, as well as novel genes associated with humoral immune response in MCT-induced PAH, thereby providing a novel insight into the molecular mechanisms underlying inflammation and immunity and therapeutic targets for PAH.