Dominant Myocardial Fibrosis and Complex Immune Microenvironment Jointly Shape the Pathogenesis of Arrhythmogenic Right Ventricular Cardiomyopathy

BACKGROUND: Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a heritable life-threatening myocardial disease characterized by ventricular arrhythmias and sudden cardiac death. Few studies used RNA-sequencing (RNA-seq) technology to analyze gene expression profiles, hub genes, dominant pathogenic processes, immune microenvironment in ARVC. This study aimed to explore these questions via integrated bioinformatics analysis. METHODS: RNA-sequencing datasets of GSE107475, GSE107311, GSE107156, and GSE107125 were obtained from the Gene Expression Omnibus database, including right and left ventricular myocardium from ARVC patients and normal controls. Weighted gene co-expression network analysis identified the ARVC hub modules and genes. Functional enrichment and protein-protein interaction analysis were performed by Metascape and STRING. Single-sample gene-set enrichment analysis (ssGSEA) was applied to assess immune cell infiltration. Transcription regulator (TF) analysis was performed by TRRUST. RESULTS: Three ARVC hub modules with 25 hub genes were identified. Functional enrichment analysis of the hub genes indicated that myocardial fibrosis was the dominant pathogenic process. Higher myocardial fibrosis activity existed in ARVC than in normal controls. A complex immune microenvironment was discovered that type 2 T helper cell, type 1 T helper cell, regulatory T cell, plasmacytoid dendritic cell, neutrophil, mast cell, central memory CD4 T cell, macrophage, CD56dim natural killer cell, myeloid-derived suppressor cell, memory B cell, natural killer T cell, and activated CD8 T cell were highly infiltrated in ARVC myocardium. The immune-related hub module was enriched in immune processes and inflammatory disease pathways, with hub genes including CD74, HLA-DRA, ITGAM, CTSS, CYBB, and IRF8. A positive linear correlation existed between immune cell infiltration and fibrosis activity in ARVC. NFKB1 and RELA were the shared TFs of ARVC hub genes and immune-related hub module genes, indicating the critical role of NFκB signaling in both mechanisms. Finally, the potential lncRNA–miRNA–mRNA interaction network for ARVC hub genes was constructed. CONCLUSION: Myocardial fibrosis is the dominant pathogenic process in end-stage ARVC patients. A complex immune microenvironment exists in the diseased myocardium of ARVC, in which T cell subsets are the primary category. A tight relationship exists between myocardial fibrosis activity and immune cell infiltration. NFκB signaling pathway possibly contributes to both mechanisms.