Multifaceted Spatial and Functional Zonation of Cardiac Cells in Adult Human Heart
Liang Chen,Kui Hua,Ningning Zhang,Jingyu Wang,Jian Meng,Zhenliang Hu,Haoxiang Gao,Fanhong Li,Yang Chen,Jie Ren,Liyang Liu,Qiming Zhou,Jin Gu,Jiangping Song,Xuegong Zhang,Shengshou Hu
DOI: https://doi.org/10.1161/circulationaha.121.055690
IF: 37.8
2022-01-01
Circulation
Abstract:HomeCirculationVol. 145, No. 4Multifaceted Spatial and Functional Zonation of Cardiac Cells in Adult Human Heart Free AccessLetterPDF/EPUBAboutView PDFView EPUBSections ToolsAdd to favoritesDownload citationsTrack citationsPermissions ShareShare onFacebookTwitterLinked InMendeleyRedditDiggEmail Jump toFree AccessLetterPDF/EPUBMultifaceted Spatial and Functional Zonation of Cardiac Cells in Adult Human Heart Liang Chen, MD, PhD, Kui Hua, PhD, Ningning Zhang, BS, Jingyu Wang, PhD, Jian Meng, BS, Zhenliang Hu, MD, Haoxiang Gao, BS, Fanhong Li, MS, Yang Chen, PhD, Jie Ren, MD, Liyang Liu, MS, Qiming Zhou, PhD, Jin Gu, PhD, Jiangping Song, MD, PhD, Xuegong Zhang, PhD and Shengshou Hu, MD Liang ChenLiang Chen https://orcid.org/0000-0001-6374-112X State Key Laboratory of Cardiovascular Disease, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China (L.C., N.Z., J.M., Z.H., J.R., J.S., S.H.). , Kui HuaKui Hua MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. , Ningning ZhangNingning Zhang State Key Laboratory of Cardiovascular Disease, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China (L.C., N.Z., J.M., Z.H., J.R., J.S., S.H.). , Jingyu WangJingyu Wang MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. Center for Synthetic and Systems Biology and School of Medicine (J.W., X.Z.), Tsinghua University, Beijing, China. , Jian MengJian Meng State Key Laboratory of Cardiovascular Disease, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China (L.C., N.Z., J.M., Z.H., J.R., J.S., S.H.). , Zhenliang HuZhenliang Hu State Key Laboratory of Cardiovascular Disease, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China (L.C., N.Z., J.M., Z.H., J.R., J.S., S.H.). , Haoxiang GaoHaoxiang Gao MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. , Fanhong LiFanhong Li https://orcid.org/0000-0002-7646-9759 MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. , Yang ChenYang Chen MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. , Jie RenJie Ren State Key Laboratory of Cardiovascular Disease, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China (L.C., N.Z., J.M., Z.H., J.R., J.S., S.H.). , Liyang LiuLiyang Liu MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. , Qiming ZhouQiming Zhou MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. , Jin GuJin Gu MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. , Jiangping SongJiangping Song Correspondence to: Jiangping Song, 167A Beilishi Rd, Beijing, 100037, China. Email E-mail Address: [email protected] https://orcid.org/0000-0002-9488-7823 State Key Laboratory of Cardiovascular Disease, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China (L.C., N.Z., J.M., Z.H., J.R., J.S., S.H.). , Xuegong ZhangXuegong Zhang Correspondence to: Xuegong Zhang, 30 Shuangqing Rd, Beijing 100084, China. Email E-mail Address: [email protected] https://orcid.org/0000-0002-9684-5643 MOE Key Laboratory of Bioinformatics, Bioinformatics Division, Beijing National Research Center for Information Science and Technology and Department of Automation (K.H., J.W., H.G., F.L., Y.C., L.L., Q.Z., J.G., X.Z.), Tsinghua University, Beijing, China. Center for Synthetic and Systems Biology and School of Medicine (J.W., X.Z.), Tsinghua University, Beijing, China. School of Life Sciences (X.Z.), Tsinghua University, Beijing, China. and Shengshou HuShengshou Hu Correspondence to: Shengshou Hu, 167A Beilishi Rd, Beijing, 100037, China. Email E-mail Address: [email protected] https://orcid.org/0000-0001-8136-0224 State Key Laboratory of Cardiovascular Disease, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China (L.C., N.Z., J.M., Z.H., J.R., J.S., S.H.). Originally published24 Jan 2022https://doi.org/10.1161/CIRCULATIONAHA.121.055690Circulation. 2022;145:315–318The heart is a well-architected organ, and study of cellular signatures across the structures of human heart is in an early stage. Our group and others have previously revealed cardiac cell remodeling under diseased and developmental scenarios.1 In addition, 2 studies reported the cellular transcriptome across 6 cardiac regions, outlining comprehensive landscapes at an unprecedented depth.2,3 A distinct heterogeneity between atria and ventricle cardiomyocytes was identified as expected, whereas the spatial heterogeneity within atria or ventricle was not well elaborated. Given that myocardial diseases preferentially affect specific regions of the heart, such as noncompaction caused by trabecular dysplasia and atrial fibrillation derived from left atria, we aim to dissect the heterogeneity of cardiac cells at more detailed and multifaceted spatial dimensions, including the previously unreported transmural axis of the ventricular wall.2,3All experiment protocols were approved by the Ethics Committee of Fuwai Hospital. The donor hearts were selected from those initially considered and processed for heart transplantation but abandoned because of mismatched heart sizes or blood types. Written informed consent for research use was obtained from families of 2 male donors who were brain-dead at the ages of 42 and 50 years. The donor hearts showed normal function in vivo and were perfused with cardioplegia followed by refrigerated transferred to the laboratory within 4 hours. In this study, we performed single nucleus RNA sequencing and obtained 160 775 nuclei from 12 cardiac sites (Figure [A]). All of the major cardiac cell types were identified with robust molecular signatures (Figure [B]). We developed a spatial heterogeneity index on the basis of the Hierarchical Graph-Based Clustering algorithm4 to quantify the heterogeneities of cells at 3 analytic levels (Figure [C]).Download figureDownload PowerPointFigure. Spatial zonation of cardiomyocytes and functional zonation of vascular cells in adult human heart. A, Schematic overview of the experimental design. An illustration of 12 anatomic sites, including left atrium (LA), right atrium (RA), interatrial septum (IAS), left ventricle (LV), right ventricle (RV), and interventricular septum (IVS). The ventricular specimens were vertically fine cut into multiple strata including 4 layers in LV (trabecular [LV-tra], subendocardial [LV-endo], middle [LV-mid] and subepicardial [LV-epi] layers), 2 layers in RV (trabecular [RV-tra] and compact [RV-com] layers), and 3 layers in IVS from left to right side (IVS-L, IVS-M, and IVS-R). B, t-SNE visualization of 12 major heart cell types from the human heart. Computation methods: log-normalization pipeline in Seurat v3 was adopted for standard analysis. Donor effect was corrected by “Harmony.” C, Spatial Heterogeneity Index (SHI) of different cell types at the 3 levels. Cardiomyocytes exhibited most dramatic heterogeneities at levels 1 and 2, and at the transmural dimension of LV (level 3), but not in IVS and RV. Mural cells and immune cells presented minimal differences at all 3 levels. Computation methods: SHI is defined as the consistency between unsupervised clustering and the anatomic regions on the basis of the Adjusted Rand Index (ARI) with the Hierarchical Graph-Based Clustering (HGC) algorithm4 cell clustering. Di, t-SNE map of cardiomyocytes colored by anatomic regions showing chambers heterogeneities. Especially atrial cardiomyocytes were identified as 2 distinct clusters, and IAS cardiomyocytes were clustered into LA and RA clusters, respectively. ii, Sankey plot showing the individual consistency of spatial heterogeneities of CMs across heart chambers. CMs from Donor 1 (reference heart) were used to train a classifier in predicting the spatial regions of CMs from Donor 2 (query heart) with a high accuracy of 0.86. Ei, Immunohistochemistry verifying chamber specific proteins in LA (MYOT) and RA (FHL2) and their layered expression in IAS. ii, Immunohistochemistry and (iii) Western blot verifying the zonated genes expression from trabecular to epicardial layers. F, Transmural zonation of CMs within LV. i, Diffusion map showing the pseudo order of LV CMs from trabecular to epicardial layer using the R package “destiny.” ii, Correlation dot plot colored by the density of dots (cells) showing the favorable individual consistency of transmural zonation of LV CMs. Computational methods: we took the pseudo-order of Donor 1 (reference heart) cells as labels and built a KNN (k=25) regressor for predicting the pseudo-order of Donor 2 (query heart) cells. G, Zonated gene expression patterns of CMs along transmural LV axis identified by R package “tradeSeq” (P adjusted<0.01). i, Heatmap of top 1000 zonated genes, clustered into 8 modules using R package “cluster,” (ii) and their enriched pathways. iii, Western blot validation of the zonated genes involved in fatty acid catabolic pathways. H, Enriched pathways in the comparative analysis of CMs between left ventricular noncompaction (LVNC) and normal hearts by Gene Set Enrichment Analysis (GSEA). I, Heatmap showing the zonated expression of genes involved in regional specific pathways (shown in H) in normal LV CMs. J, Vasculature zonation of ECs. i, Diffusion map showing the artery-capillary-vein vasculature zonation of vascular ECs. ii, Heatmap of the top 1000 zonated genes, which were clustered into 6 modules with distinct enriched pathways. K, Schematic diagram of vascular cells distribution. L, Vasculature zonation of mural cells. i, Diffusion map showing the vasculature zonation of pericytes (artery-capillary/vein) and SMCs (smooth muscle cells; artery-arteriole-vein) individually. The zonations of 2 mural cell types presented head-to-head directions with overlap at the venous site. ii, Heatmap of the top 1000 zonated genes of mural cells along the vasculature order. Novel markers in arteriole (PLD5, HMCN2) and capillary (venule) pericytes (PDZD2, CPM) were identified. M, Cozonated genes along the vasculature between pericytes and ECs. The arteriole of both ECs and pericytes coexpressed genes were involved in ECM organization (ELN), whereas capillary-derived ECs and pericytes coexpressed lipid binding genes (LPL, CD36) and VEGF receptor (FLT1). N, Multiplex fluorescence immunohistochemistry validated the distribution of GJA5/PSCK5 (in arteriole), CD36 (in capillary), and ACKR1 (in venule) in vascular ECs, the specific expression of potassium channels including KCNN3 in arteriole, KCNB1 in capillary, and KCNIP4 in venule ECs, and coexpression genes of ECs and pericytes at arteriole (ELN) and capillary (SH3RF3 and FLT1) vessels. PECAM1 is a pan-marker for ECs, and PDGFRB for pericytes. CM indicates cardiomyocyte; CM-A/V, CM-atria/ventricle; EC, endothelial cell; EndoEC, endocardial cell; IVS/IAS, interventricular/atrial septum; LEC, lymphatic EC; LV, left ventricle; LVNC, left ventricular noncompaction; RV, right ventricle; and SMC, smooth muscle cell.Cardiomyocytes exhibited dramatic heterogeneities not only at levels 1 and 2 but also in transmural dimension of the left ventricle (level 3), whereas mural cells and immune cells presented minimal spatial differences. Subclustering of cardiomyocytes identified high agreement between clusters and anatomic chambers, explaining the high spatial heterogeneity index at levels 1 and 2 (Figure [Di]). The interindividual consistency of chamber specificity was confirmed by a high accuracy rate of 0.86 when training a classifier on 1 reference heart to predict cell locations of the other (Figure [Dii]). Beyond previous studies,2,3 atrial cardiomyocytes were identified as 2 distinct clusters dominated by the left atria and right atria, respectively. Left atria cardiomyocytes were functionally enriched in heart contraction, whereas right atria were enriched in cardiac rhythm. Immunohistochemistry validated distinct gene expression (MYOT, FHL2) between left atria and right atria, and layered patterns in the atrial septum (Figure [Ei]). To illuminate the high spatial heterogeneity index of the left ventricle cardiomyocytes at level 3, we built a trajectory using a diffusion map and demonstrated a zonation pattern along the trabecular-to-epicardial axis (Figure [Fi]), which was conserved between individuals (Figure Fii). Along the transmural axis, 7168 out of the 17 089 (42%, P adjusted<0.01) genes presented dynamic alterations, of which the top 1000 variable genes were clustered into 8 expression modules with distinct biological functions (Figure [Gi, Gii, Eii, and Eiii]). Module 2 showed a linearly decreasing trend of genes enriched in myofibril. Modules 4 and 5 showed gradually increasing patterns and were enriched in circadian rhythm and fatty acid use pathways (ie, CPT1A, CD36, ACOX1) (Figure [Giii]).We further involved 2 diseased hearts with left ventricular noncompaction and performed single nucleus RNA sequencing for trabecular and compact layers of the left ventricle, respectively (Figure [H]). Compared with the normal heart, the compact layer of left ventricular noncompaction was enriched in cardiomyocyte development pathways (ie, HEY2, the marker of compact component), whereas the trabecular layer was enriched in myofibril (including left ventricular noncompaction pathogenic genes MYH7 and TNNT2) and RNA splicing process (ie, RBM20, a novel left ventricular noncompaction pathogenic gene through regulating splicing of myofibril).5 It is interesting that 150 (84%) out of the 178 enriched genes were zonated in the normal heart, suggesting that zonated genes were essential in maintaining heart structure and function (Figure [I]).The vascular endothelial cells were aligned by trajectory using a diffusion map (Figure [Ji and K]). The consecutive expression of markers and enriched pathways confirmed an artery-capillary-vein vasculature zonation along the trajectory (Figure [Jii and Jiii]). Novel vasculature genes were identified corresponding to the functions of blood vessels (Figure [Jii and N]). It is interesting that vascular endothelial cells expressed distinct potassium channels along the vasculature axis (ie, KCNN3, KCNIP4, and KCNB1) (Figure [Jii and N]). Smooth muscle cells and pericytes, the perivascular mural cell types, also showed zonation along the vasculature separately, and presented head-to-head directions with overlap at the venous site (Figure [K, Li, and Lii]). Of note, we identified a set of cozonated genes shared by endothelial cells and pericytes, confirming their synergistic functions in maintaining vascular function (Figure [M and N]).We built an interactive database (http://xglab.tech/hahca) to visualize cellular gene expression and zonation. The data and methods used in the analysis are available to any researchers for purposes of reproducing the results or replicating the procedure.In conclusion, our high-resolution cellular landscape provided a fundamental reference for understanding adult human hearts. Our results pinpointed the spatial and functional zonations of cardiomyocytes and vascular cells, which should be widely recognized for cardiovascular studies, because these profound discoveries will pave the path for in-depth investigation into myocardial diseases such as atrial cardiomyopathy and noncompaction, as well as regulation of angiogenesis and inflammatory infiltration during myocardial injury.Article InformationAcknowledgmentsThe authors thank Dr. Hairong Lyu for help in early coordination of this project and Hao Sun for help with the webpage and database. The authors thank Dr Peng Zhang for all the discussion about this work.Sources of FundingThis work was supported in part by the National Natural Science Foundation of China Projects (82100377, 61721003, 61922047, 62050178), Beijing Nova Program (2021038), National Key R&D Program of China (2021YFF1200900, 2016YFC1300902), and Tsinghua-Fuzhou Institute of Data Technologies (TFIDT2021003).Disclosures None.Footnotes*L. Chen and K. Hua contributed equally.For Sources of Funding and Disclosures, see page 318.Correspondence to: Shengshou Hu, 167A Beilishi Rd, Beijing, 100037, China. Email [email protected]orgCorrespondence to: Xuegong Zhang, 30 Shuangqing Rd, Beijing 100084, China. Email [email protected]edu.cnCorrespondence to: Jiangping Song, 167A Beilishi Rd, Beijing, 100037, China. Email [email protected]comReferences1. Asp M, Giacomello S, Larsson L, Wu C, Fürth D, Qian X, Wärdell E, Custodio J, Reimegård J, Salmén F, et al.. A spatiotemporal organ-wide gene expression and cell atlas of the developing human heart.Cell. 2019; 179:1647–1660.e19. doi: 10.1016/j.cell.2019.11.025CrossrefMedlineGoogle Scholar2. Tucker NR, Chaffin M, Fleming SJ, Hall AW, Parsons VA, Bedi KC, Akkad AD, Herndon CN, Arduini A, Papangeli I, et al.. Transcriptional and cellular diversity of the human heart.Circulation. 2020; 142:466–482. doi: 10.1161/CIRCULATIONAHA.119.045401LinkGoogle Scholar3. Litviňuková M, Talavera-López C, Maatz H, Reichart D, Worth CL, Lindberg EL, Kanda M, Polanski K, Heinig M, Lee M, et al.. Cells of the adult human heart.Nature. 2020; 588:466–472. doi: 10.1038/s41586-020-2797-4CrossrefMedlineGoogle Scholar4. Zou Z, Hua K, Zhang X. HGC: fast hierarchical clustering for large-scale single-cell data.Bioinformatics. 2021; 37:3964–3965. doi: 10.1093/bioinformatics/btab420CrossrefGoogle Scholar5. Sedaghat-Hamedani F, Haas J, Zhu F, Geier C, Kayvanpour E, Liss M, Lai A, Frese K, Pribe-Wolferts R, Amr A, et al.. Clinical genetics and outcome of left ventricular non-compaction cardiomyopathy.Eur Heart J. 2017; 38:3449–3460. doi: 10.1093/eurheartj/ehx545CrossrefMedlineGoogle Scholar Previous Back to top Next FiguresReferencesRelatedDetailsCited By Chen S, Luo Y, Gao H, Li F, Chen Y, Li J, You R, Hao M, Bian H, Xi X, Li W, Li W, Ye M, Meng Q, Zou Z, Li C, Li H, Zhang Y, Cui Y, Wei L, Chen F, Wang X, Lv H, Hua K, Jiang R and Zhang X (2022) hECA: The cell-centric assembly of a cell atlas, iScience, 10.1016/j.isci.2022.104318, 25:5, (104318), Online publication date: 1-May-2022. January 25, 2022Vol 145, Issue 4Article InformationMetrics © 2022 American Heart Association, Inc.https://doi.org/10.1161/CIRCULATIONAHA.121.055690PMID: 35073179 Originally publishedJanuary 24, 2022 Keywordsendothelial cellshearttranscriptomecardiomyocytesPDF download Advertisement SubjectsBasic Science ResearchCell Biology/Structural BiologyGene Expression and Regulation