Identification of the Immune Status of Hypertrophic Cardiomyopathy by Integrated Analysis of Bulk- and Single-Cell RNA Sequencing Data
Authors: Zhao Wei , Wu Tianyu , Zhan Jian , Dong Zishuang 
Published: October 2022

Objectives. Hypertrophic cardiomyopathy (HCM) is the most common hereditary cardiomyopathy and immune infiltration is considered an indispensable factor involved in its pathogenesis. In this study, we attempted to combine bulk sequencing and single-cell sequencing to map the immune infiltration-related genes in hypertrophic cardiomyopathy. Methods. The GSE36961, GSE160997, and GSE122930 datasets were obtained from the Gene Expression Omnibus database. The compositional patterns of the 18 types of immune cell fraction and pathway enrichment score in control and HCM patients were estimated based on the GSE36961 cohort using xCell algorithm. The Weighted Gene Coexpression Network Analysis (WGCNA) was performed to identify genes associated with immune infiltration for hypertrophic cardiomyopathy. The area under the curve (AUC) value was obtained and used to evaluate the discriminatory ability of common immune-related DEGs. “NetworkAnalyst” platform was used to identify TF-gene and TF-miRNA interaction with identified common genes. Heat map was used to determine the association between common DEGs and various immune cells. Results. Immune infiltration analysis by the xCell algorithm showed a higher level of CD8+ naive T cells, CD8+ T cells, as well as a lower level of activated dendritic cells (aDC), dendritic cells (DC), immature dendritic cells (iDC), conventional dendritic cells (cDC), macrophages, M1 macrophages, monocytes, and NKT cell in HCM compared with the control group in GSE36961 dataset. aDC, macrophages, and M1 macrophages were the top three discriminators between HCM and control groups with the area under the curve (AUC) of 0.907, 0.867, and 0.941. WGCNA analysis showed that 1258 immune-related genes were included in four different modules. Of these modules, the turquoise module showed a pivotal correlation with HCM. 13 common immune-related DEGs were found by intersecting common DEGs in GSE36961 and GSE160997 datasets with genes from the genes in turquoise module. 5 hub immune-related genes (S100A9, TYROBP, FCER1G, CD14, and S100A8) were identified by protein interaction network. Through analysis of single-cell sequencing data, S100a9, TYROBP, FCER1G, and S100a8 were mainly expressed by infiltrated M1 proinflammatory cells, especially Ccr2-M1 proinflammatory macrophage cells in the heart immune microenvironment while Cd14 was expressed by infiltrated M1 proinflammatory macrophage cells and M2 macrophages in transverse aortic constriction (TAC) mice at 1 week. Higher M2 macrophage and M1 proinflammatory macrophage infiltration as well as lower Ccr2-M1 proinflammatory macrophage and dendritic cells were shown in TAC 1week mice compared with sham mice. Conclusions. There was a difference in immune infiltration between HCM patients and normal groups. aDC, macrophages, and M1 macrophages were the top three discriminator immune cell subsets between HCM and control groups. S100A9, TYROBP, FCER1G, CD14, and S100A8 were identified as potential biomarkers to discriminate HCM from the control group. S100a9, TYROBP, FCER1G, and S100a8 were mainly expressed by infiltrated M1 proinflammatory cells, especially Ccr2-M1 proinflammatory cells in the heart immune microenvironment while Cd14 was expressed by M2 macrophages in transverse aortic constriction (TAC) mice at 1 week.