Spatial Transcriptomic Analysis of Focal and Normal Areas of Myocyte Disarray in Human Hypertrophic Cardiomyopathy

Abstract

Hypertrophic Cardiomyopathy (HCM) is a common inherited disorder that can lead to heart failure and sudden cardiac death, characterized at the histological level by focal areas of myocyte disarray, hypertrophy and fibrosis, and few disease-targeted therapies exist. To identify, focal, spatially restricted alterations in transcriptional pathways and reveal novel therapeutic targets, we performed a spatial transcriptomic analysis of areas of focal myocyte disarray compared to areas of normal tissue, using a commercially available platform (GeoMx, nanoString). We analyzed surgical myectomy tissue from four patients with HCM and control interventricular septum tissue from two unused organ donor hearts that were free of cardiovascular disease. Histological sections were reviewed by an expert pathologist and 72 focal areas with varying degrees of myocyte disarray (normal, mild, moderate, severe) were chosen for analysis. Areas of interest were interrogated with the Human Cancer Transcriptome Atlas designed to profile 1800 transcripts. Differential expression analysis revealed significant changes in gene expression between HCM and Control tissue, and functional enrichment analysis indicated these genes were primarily involved in interferon production and mitochondrial energetics. Within HCM tissue, differentially expressed genes between areas of mild and moderate disarray were enriched for genes related to mitochondrial energetics (moderate disarray) and response to oxygen/cytokine levels (mild disarray). The comparison between areas of moderate and severe disarray were enriched for genes related to the c-Jun N-terminal kinase (JNK) cascade in severe disarray. Analysis of ligand-receptor pair gene expression revealed that HCM tissue exhibited downregulation of platelet-derived growth factor (PDGF), NOTCH, junctional adhesion molecule, and CD46 signaling, while showing upregulation of fibronectin, CD99, cadherin, and amyloid precursor protein signaling. A deconvolution analysis utilizing the matched single nuclei RNA-sequencing (snRNA-seq) data to determine cell type composition in areas of interest revealed significant differences in fibroblast and vascular cell composition in areas of severe disarray when compared to normal areas in HCM samples. Cell composition in normal areas from control tissue was also divergent from normal areas in HCM samples, which was consistent with the differential expression results. Overall, our data identify novel and potential disease-modifying targets for therapy in HCM.