Lineage-specific regulatory changes in hypertrophic cardiomyopathy unraveled by single-nucleus RNA-seq and spatial transcriptomics

Hypertrophic cardiomyopathy (HCM) is the most common cardiac genetic disorder characterized by cardiomyocyte hypertrophy and cardiac fibrosis. Pathological cardiac remodeling in the myocardium of HCM patients may progress to heart failure. An in-depth elucidation of the lineage-specific changes in pathological cardiac remodeling of HCM is pivotal for the development of therapies to mitigate the progression. Here, we performed single-nucleus RNA-seq of the cardiac tissues from HCM patients or healthy donors and conducted spatial transcriptomic assays on tissue sections from patients. Unbiased clustering of 55,122 nuclei from HCM and healthy conditions revealed 9 cell lineages and 28 clusters. Lineage-specific changes in gene expression, subpopulation composition, and intercellular communication in HCM were discovered through comparative analyses. According to the results of pseudotime ordering, differential expression analysis, and differential regulatory network analysis, potential key genes during the transition towards a failing state of cardiomyocytes such as FGF12, IL31RA, and CREB5 were identified. Transcriptomic dynamics underlying cardiac fibroblast activation were also uncovered, and potential key genes involved in cardiac fibrosis were obtained such as AEBP1, RUNX1, MEOX1, LEF1, and NRXN3. Using the spatial transcriptomic data, spatial activity patterns of the candidate genes, pathways, and subpopulations were confirmed on patient tissue sections. Moreover, we showed experimental evidence that in vitro knockdown of AEBP1 could promote the activation of human cardiac fibroblasts, and overexpression of AEBP1 could attenuate the TGFβ-induced activation. Our study provided a comprehensive analysis of the lineage-specific regulatory changes in HCM, which laid the foundation for targeted drug development in HCM.