Heart failure (HF) with preserved ejection fraction (HFpEF) constitutes half of all HF but lacks effective therapy. Understanding of its myocardial biology remains limited because of a paucity of heart tissue molecular analysis.

We performed RNA sequencing on right ventricular septal endomyocardial biopsies prospectively obtained from patients meeting consensus criteria for HFpEF (n=41) contrasted with right ventricular septal tissue from patients with HF with reduced ejection fraction (HFrEF, n=30) and donor controls (n=24). Principal component analysis and hierarchical clustering tested for transcriptomic distinctiveness between groups, effect of comorbidities, and differential gene expression with pathway enrichment contrasted HF groups and donor controls. Within HFpEF, non–negative matrix factorization and weighted gene coexpression analysis identified molecular subgroups, and the resulting clusters were correlated with hemodynamic and clinical data.

Patients with HFpEF were more often women (59%), African American (68%), obese (median body mass index 41), and hypertensive (98%), with clinical HF characterized by 65% New York Heart Association Class III or IV, nearly all on a loop diuretic, and 70% with a HF hospitalization in the previous year. Principal component analysis separated HFpEF from HFrEF and donor controls with minimal overlap, and this persisted after adjusting for primary comorbidities: body mass index, sex, age, diabetes, and renal function. Hierarchical clustering confirmed group separation. Nearly half the significantly altered genes in HFpEF versus donor controls (1882 up, 2593 down) changed in the same direction in HFrEF; however, 5745 genes were uniquely altered between HF groups. Compared with controls, uniquely upregulated genes in HFpEF were enriched in mitochondrial adenosine triphosphate synthesis/electron transport, pathways downregulated in HFrEF. HFpEF-specific downregulated genes engaged endoplasmic reticulum stress, autophagy, and angiogenesis. Body mass index differences largely accounted for HFpEF upregulated genes, whereas neither this nor broader comorbidity adjustment altered pathways enriched in downregulated genes. Non–negative matrix factorization identified 3 HFpEF transcriptomic subgroups with distinctive pathways and clinical correlates, including a group closest to HFrEF with higher mortality, and a mostly female group with smaller hearts and proinflammatory signaling. These groupings remained after sex adjustment. Weighted gene coexpression analysis yielded analogous gene clusters and clinical groupings.

HFpEF exhibits distinctive broad transcriptomic signatures and molecular subgroupings with particular clinical features and outcomes. The data reveal new signaling targets to consider for precision therapeutics.