Single cell transcriptomic analyses of human heart failure with preserved ejection fraction

Background: Heart failure with preserved ejection fraction (HFpEF) is a poorly understood, multi-system disease with high morbidity and mortality. To improve our understanding of its underlying biology, we used single-nucleus RNA sequencing (snRNA-seq) to characterize cell-specific gene expression patterns in human HFpEF myocardium. Methods: Septal myocardial biopsies (2-3 mg) from 30 HFpEF patients and 29 non-failing donor controls were analyzed using the 10X Genomics platform, with nuclei isolated from combined samples (6 patients/pool). Genotype-based demultiplexing was performed with souporcell, and gene expression quantified with CellRanger and CellBender. After quality control, nuclei were clustered and annotated by cell types based on specific marker genes. Differential expression (DE) by cell-type in HFpEF vs controls was performed using limma-voom and functional analysis performed using Gene Set Enrichment Analysis. Data were compared to dilated cardiomyopathy (DCM) using prior snRNA-seq in DCM vs respective controls. Results: We successfully demultiplexed pooled myocardial biopsies, assigning >75% of droplets to individual patients. From eight pooled samples (19 HFpEF, 24 controls), we recovered 48,886 nuclei and identified 14 cell types. Cardiomyocytes (5159 differentially expressed [DE] genes, 36%) and fibroblasts (5905 DE genes, 49%) showed the most DE genes, while endothelial cells (2143), pericytes (1812), and macrophages (1405) had fewer. Enriched pathways common to multiple cell types included transcription/translation, immune activation, metabolism, and protein quality control. Of 7848 DE genes identified via pseudo-bulk snRNA-seq, 51% were DE in fibroblasts and 47% in cardiomyocytes, compared to <20% in other cell types. Unlike dilated cardiomyopathy (DCM), sub-clustering fibroblasts did not reveal an activated fibroblast population in HFpEF. Comparative analysis between HFpEF and DCM identified transcriptional differences primarily in cardiomyocytes. Conclusions: This study demonstrates the power of genotype-based demultiplexing for single-cell transcriptomic analyses of small endomyocardial biopsies and identifies cardiomyocytes as the principal cell type with distinct transcriptional changes in HFpEF versus DCM. These findings, coupled with differential gene expression and functional pathway analyses, illuminate HFpEF pathways and may nominate compelling targets for future mechanistic studies and therapeutic efforts for HFpEF.