Comprehensive analysis of scRNA-seq and bulk RNA-seq reveals the non-cardiomyocytes heterogeneity and novel cell populations in dilated cardiomyopathy.

Background: Dilated cardiomyopathy (DCM) is one of the most common causes of heart failure. Infiltration and alterations in non-cardiomyocytes of the human heart involve crucially in the occurrence of DCM and associated immunotherapeutic approaches.

Methods: We constructed a single-cell transcriptional atlas of DCM and normal patients. Then, the xCell algorithm, EPIC algorithm, MCP counter algorithm, and CIBERSORT method were applied to identify DCM-related cell types with a high degree of precision and specificity using RNA-seq datasets. We further analyzed the heterogeneity among cell types, performed trajectory analysis, examined transcription factor regulatory networks, investigated metabolic heterogeneity, and conducted intercellular communication analysis. Finally, we used bulk RNA-seq data to confirm the roles of M2-like2 subpopulations and GAS6 in DCM.

Results: We integrated and analyzed Single-cell sequencing (scRNA-seq) data from 7 DCM samples and 3 normal heart tissue samples, totaling 70,958 single-cell data points. Based on gene-specific expression and prior marker genes, we identified 9 distinct subtypes, including fibroblasts, endothelial cells, myeloid cells, pericytes, T/NK cells, smooth muscle cells, neuronal cells, B cells, and cardiomyocytes. Using machine learning methods to quantify bulk RNA-seq data, we found significant differences in fibroblasts, T cells, and macrophages between DCM and normal samples. Further analysis revealed high heterogeneity in tissue preference, gene expression, functional enrichment, immunodynamics, transcriptional regulatory factors, metabolic changes, and communication patterns in fibroblasts and myeloid cells. Among fibroblast subpopulations, proliferative F3 cells were implicated in the fibroblast transition process in DCM, while myofibroblast F6 cells promoted the fibroblast transition to a late cell state in DCM. Additionally, two subpopulations of M2 macrophages, M2-like1 and M2-like2, were identified with distinct features. The M2-like2 cell subpopulation, which was enriched in glycolysis and fatty acid metabolism, involved in inflammation inhibition and fibrosis promotion. Cell‒cell communication analysis indicated the GAS6-MERTK axis might exhibit interaction between M2 macrophage and M2-like1 macrophage. Furthermore, deconvolution analysis for bulk RNA-seq data revealed a significant increase in M2-like2 subpopulations in DCM, suggesting a more important role for this cell population in DCM.

Conclusions: We revealed the heterogeneity of non-cardiomyocytes in DCM and identified subpopulations of myofibroblast and macrophages engaged in DCM, which suggested a potential significance of non-cardiomyocytes in treatment of DCM.