Lineage-specific regulatory changes in the pathological cardiac remodeling of hypertrophy cardiomyopathy unraveled by single-nucleus RNA-seq and spatial transcriptomics

BACKGROUNDHypertrophy cardiomyopathy (HCM) is the most common cardiac genetic disorder with the histopathological features of cardiomyocyte hypertrophy and cardiac fibrosis. The pathological remodeling that occurs in the myocardium of HCM patients may ultimately progress to heart failure and death. A thorough understanding of the cell type-specific changes in the pathological cardiac remodeling of HCM is crucial for developing successful medical therapies to prevent or mitigate the progression of this disease. METHODSWe performed single-nucleus RNA-seq of the cardiac tissues from 10 HCM patients and 2 healthy donors, and conducted spatial transcriptomic assays of 4 cardiac tissue sections from 3 HCM patients. Comparative analyses were performed to explore the lineage-specific changes in expression profile, subpopulation composition and intercellular communication in the cardiac tissues of HCM patients. Based on the results of independent analyses including pseudotime ordering, differential expression analysis, and differential regulatory network analysis, we prioritized candidate therapeutic targets for mitigating the progression to heart failure or attenuating the cardiac fibrosis in HCM. Using the spatial transcriptomic data, we examined the spatial activity patterns of the key candidate genes, pathways and subpopulations. RESULTSUnbiased clustering of 55,122 nuclei from HCM and healthy conditions revealed 9 cell lineages and 28 clusters. Significant expansion of vascular-related lineages and contraction of cardiomyocytes, fibroblasts and myeloid cells in HCM were observed. The transcriptomic dynamics during the transition towards the failing state of cardiomyocytes in HCM were uncovered. Candidate target genes for mitigating the progression to heart failure in HCM were obtained such as FGF12, IL31RA, BDNF, S100A1, CRYAB and PROS1. The transcriptomic dynamics underlying the fibroblast activation were also uncovered, and candidate targets for attenuating the cardiac fibrosis in HCM were obtained such as RUNX1, MEOX1, AEBP1, LEF1 and NRXN3. CONCLUSIONSWe provided a comprehensive analysis of the lineage-specific regulatory changes in HCM. Our analysis identified a vast array of candidate therapeutic target genes and pathways to prevent or attenuate the pathological remodeling of HCM. Our datasets constitute a valuable resource to examine the lineage-specific expression changes of HCM at single-nucleus and spatial resolution. We developed a web-based interface (http://snsthcm.fwgenetics.org/) for further exploration.