Single Nucleus Transcriptomics: Apical Resection in Newborn Pigs Extends the Time Window of Cardiomyocyte Proliferation and Myocardial Regeneration
Yuji Nakada,Yang Zhou,Wuming Gong,Eric Y. Zhang,Erik Skie,Thanh Nguyen,Yuhua Wei,Meng Zhao,Wangping Chen,Jiacheng Sun,S. Naqi Raza,Jake Y. Chen,Gregory P. Walcott,Daniel J. Garry,Jianyi (Jay) Zhang
DOI: https://doi.org/10.1161/circulationaha.121.056995
IF: 37.8
2022-01-01
Circulation
Abstract:HomeCirculationVol. 145, No. 23Single Nucleus Transcriptomics: Apical Resection in Newborn Pigs Extends the Time Window of Cardiomyocyte Proliferation and Myocardial Regeneration Free AccessLetterPDF/EPUBAboutView PDFView EPUBSections ToolsAdd to favoritesDownload citationsTrack citationsPermissions ShareShare onFacebookTwitterLinked InMendeleyReddit Jump toFree AccessLetterPDF/EPUBSingle Nucleus Transcriptomics: Apical Resection in Newborn Pigs Extends the Time Window of Cardiomyocyte Proliferation and Myocardial Regeneration Yuji Nakada, PhD, Yang Zhou, PhD, Wuming Gong, PhD, Eric Y. Zhang, PhD, Erik Skie, BS, Thanh Nguyen, PhD, Yuhua Wei, MS, Meng Zhao, MD, PhD, Wangping Chen, MD, Jiacheng Sun, MD, S. Naqi Raza, BS, Jake Y. Chen, PhD, Gregory P. Walcott, MD, Daniel J. Garry, MD, PhD and Jianyi (Jay) Zhang, MD, PhD Yuji NakadaYuji Nakada https://orcid.org/0000-0001-8229-4775 Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , Yang ZhouYang Zhou https://orcid.org/0000-0002-9217-5208 Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , Wuming GongWuming Gong Department of Medicine, School of Medicine, University of Minnesota, Minneapolis (W.G., E.S., D.J.G.). , Eric Y. ZhangEric Y. Zhang Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , Erik SkieErik Skie Department of Medicine, School of Medicine, University of Minnesota, Minneapolis (W.G., E.S., D.J.G.). , Thanh NguyenThanh Nguyen https://orcid.org/0000-0002-8440-1594 Informatics Institute (T.N., J.Y.C.), University of Alabama at Birmingham. , Yuhua WeiYuhua Wei https://orcid.org/0000-0003-0140-873X Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , Meng ZhaoMeng Zhao Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , Wangping ChenWangping Chen Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , Jiacheng SunJiacheng Sun Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , S. Naqi RazaS. Naqi Raza Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. , Jake Y. ChenJake Y. Chen Informatics Institute (T.N., J.Y.C.), University of Alabama at Birmingham. , Gregory P. WalcottGregory P. Walcott Department of Medicine, Cardiovascular Disease (G.P.W., J.Z.), University of Alabama at Birmingham. , Daniel J. GarryDaniel J. Garry Department of Medicine, School of Medicine, University of Minnesota, Minneapolis (W.G., E.S., D.J.G.). and Jianyi (Jay) ZhangJianyi (Jay) Zhang Correspondence to: Jianyi (Jay) Zhang, MD, PhD, Department of Biomedical Engineering, School of Medicine and School of Engineering, University of Alabama at Birmingham, 1670 University Blvd, Volker Hall G094J, Birmingham, AL 35233. Email E-mail Address: [email protected] Department of Biomedical Engineering (Y.N., Y.Z., E.Y.Z., Y.W., M.Z., W.C., J.S., S.N.R., J.Z.), University of Alabama at Birmingham. Department of Medicine, Cardiovascular Disease (G.P.W., J.Z.), University of Alabama at Birmingham. Originally published6 Jun 2022https://doi.org/10.1161/CIRCULATIONAHA.121.056995Circulation. 2022;145:1744–1747We have shown that when permanent occlusion of left anterior descending coronary artery surgery is performed in newborn large mammals (pigs) on postnatal day 1 (P1),1 the animals can recover completely from a myocardial infarction (MI) that occurs on postnatal day 28 (P28) with no evidence of scarring or decline in contractile performance.2 Our results also suggested that cardiac apical resection (AR) on P1 preserved cardiomyocyte cell cycle activity as the animals aged; thus, we analyzed single-nucleus RNA sequencing data sets to compare the transcriptomes of cardiomyocytes from fetal pigs (embryonic day 80) with those from pigs that underwent AR surgery on P1 (ARP1), AR surgery on P1 followed by left anterior descending coronary artery occlusion MI induction on P28 (ARP1MIP28), MI on P28 without previous AR (MIP28), or neither surgical procedure (controls). Hearts from ARP1MIP28 animals were explanted from postnatal day 30 (P30) through postnatal day 56 (P56) and hearts from all other groups were explanted on P1, P28, or P56 (Figure [A–D]). All experimental protocols were approved by the Institutional Animal Care and Use Committee of the University of Alabama at Birmingham and performed under the National Institutes of Health Guide for the Care and Use of Laboratory Animals. University of Alabama at Birmingham is licensed as an animal research facility (64-R-0004) by the US Department of Agriculture and has an Animal Welfare Assurance (A3255-01) on file with the Office of Laboratory Animal Welfare. Data were processed using the Seurat toolkit, latent representations of single nuclei were obtained with scvi-tools, and dimension reduction was performed using uniform manifold approximation and projection.Download figureDownload PowerPointFigure. Single-nucleus RNA sequencing analysis identified 6 cardiomyocyte subpopulations in injured and uninjured pig hearts. A, Timeline of experimental procedures. B, Coronary angiographic image with identification of the ligation site (arrow). C and D, Images of (C) the resected apex and (D) left ventricular sections; the infarcted region is circled. E and F, Cluster analysis for single-nucleus RNA sequencing data from (E) all cells and (F) cardiomyocytes; results are presented according to (Ei and Fi) experimental group/time point, (Eii) cell type (defined by marker expression), (Eiii) abundance of each of 4 cardiomyocyte-specific genes, and (Fii) cardiomyocyte cluster. G, Distribution of cardiomyocytes from each experimental group/time point across cardiomyocyte clusters. H, Heat map of expression levels for genes that were explicitly associated with each cardiomyocyte cluster. I, The expression of genes associated with each cell cycle phase was analyzed and used to calculate a G2M phase probability score for each cardiomyocyte cluster. Ji, Representative images from each section: fetal heart (FH)–embryo, ARP1MIP28–P30, and MIP28–P30 were stained with PKM2 (pyruvate kinase M2) antibody, cardiac troponin T, and DAPI. Jii, Quantification of PKM2 intensity by Image J (n=3 each). Square dots are for high magnification. Scale bar, 50 um. **P<0.01. K, The expression of glycolysis genes in each cardiomyocyte cluster was compared with their expression in high (FH–embryo) and low (control–P56) cardiomyocytes. L, The proportion of cardiomyocytes with high, medium, and low levels of expression for genes associated with each phase of the cell cycle. M, Violin plot of the sparse model score. N, Pathways that were significantly upregulated in ARP1MIP28–P30 and ARP1MIP28–P35 cardiomyocytes were identified through Gene Ontology Analysis. O, The proportion of cardiomyocytes with high, medium, and low levels of expression for genes associated with adrenergic signaling, Z disc development, and contractile activity. AR indicates apical resection; CM, cardiomyocyte; LAD, left anterior descending; MI, myocardial infarction; P, postnatal day; and UMAP, uniform manifold approximation and projection.All single-cell transcriptomic data, including the raw and processed files, are publicly available at Gene Expression Omnibus (accession number, GSE185289; URL: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE185289). A total of 218 945 high-quality nuclei were captured and Louvain clustering analysis for known lineage markers identified all 8 major cardiac cell types (Figure [E]). The Figure (Ei) shows that most of the noncardiomyocyte clusters contain nuclei from all heart samples regardless of the treatment or age, indicating limited variation between different batches of sample collections. Complete gene nucleus data were available for 94 844 cardiomyocytes and variations between the anterior apical zone and the remote zone within each group were negligible when compared with variations between groups. Cardiomyocytes were distributed among 6 clusters (cardiomyocyte 1–cardiomyocyte 6; Figure [F]). Three of the clusters were almost entirely composed of cardiomyocytes from a single experimental group (cardiomyocyte 2 [96.2% control–P56], cardiomyocyte 3 [95.7% control–P1], and cardiomyocyte 6 [92.4% fetal heart–embryo]); cardiomyocyte 4 and 5 clusters contained primarily ARP1MIP28–P35 cardiomyocytes (68.6% and 97.8%, respectively). The cardiomyocyte 1 cluster included cardiomyocytes from all other animal groups and time points as well as a small number of ARP1MIP28–P35 cardiomyocytes but few (<1%) fetal or control–P56 cardiomyocytes (Figure [G]).We found that cardiomyocyte clusters 3 through 6 have a significantly higher G2M score, representing more proliferative capability (Figure [I]). Cardiomyocyte cluster 4 had elevated expression of NTRK2 (neurotrophic receptor tyrosine kinase 2), which is involved in the Hippo signaling pathway (Figure [H]). Cardiomyocyte cluster 5 was characterized by highly expressed PKM2 (pyruvate kinase M2; Figure [H]). Transient overexpression of Pkm2 has been linked to increases in cardiomyocyte proliferation and improvements in cardiac function after ischemic injury in mice. Immunofluorescence analyses confirmed that ARP1MIP28–P30 cardiomyocytes expressed PKM2 protein (Figure [Ji and Jii]) with ample evidence of PKM2 in ARP1MIP28–P30 to ARP1MIP28–P42 cardiomyocytes. Furthermore, glycolytic molecular programs were activated in cardiomyocyte cluster 5, almost similar to cardiomyocyte clusters 3 and 6 (Figure [K]).To evaluate the proliferative capacity of cardiomyocytes from each group and time point, we analyzed our single-nucleus RNA sequencing data sets using Support Vector Machine (an algorithm for dividing a highly dimensional data set3 into 2 categories, such as proliferative and nonproliferative) and the linear sparse technique.4,5 Our results suggested that the proportion of control cardiomyocytes in each phase of the cardiomyocyte cell cycle declined from P1 to P56, whereas the proportions for control–P1 cardiomyocytes and cardiomyocytes from all other injury groups and at all other time points were broadly similar (Figure [L]). We also found 506 genes that were positively correlated with highly proliferative (ARP1MIP28–P30 and ARP1MIP28–P35) cardiomyocytes (Figure [M]) and Gene Ontology analysis indicated that these genes participated in pathways that regulated heart development, cell proliferation, and cardiomyocyte proliferation (Figure [N]). Furthermore, assessments of adrenergic signaling, Z disc development, and contractile activity (Figure [O]) increased from P1 to P56 in control cardiomyocytes, whereas measurements in cardiomyocytes from all other injury groups and time points were ≈40% to 60% of the measurements in control–P56 cardiomyocytes.The results from our single-nucleus RNA sequencing analyses identified 6 distinct cardiomyocyte subpopulations. Two of the subpopulations (cardiomyocytes 4 and 5) were cardiomyocytes with increased cell cycle activity and proliferation and very active by combined ARP1MIP28 injury, which contributed to the remuscularization of the injured left ventricle. The interpretation of our results is limited by the multinuclear nature of cardiomyocytes and the lack of data for cytoplasmic RNA transcripts. It will be interesting to examine these 2 distinct cell populations from ARP1MIP28 injury in future studies. These observations provide a foundation for future investigations of the molecular mechanisms and signaling molecules that regulate the injury-induced preservation of cardiomyocyte cell cycle activity in newborn large mammals.Article InformationAcknowledgmentsThe authors thank the University of Alabama at Birmingham Flow Cytometry and Single Cell Core and Genomics Core. Drs Nakada, Zhou, Gong, E.Y. Zhang, Nguyen, Zhao, W. Chen, Sun, J.Y. Chen, and Walcott and E. Skie, Y. Wei, and S.N. Raza conducted the experiments and analyzed the data. Drs Nakada, Zhou, Gong, E.Y. Zhang, Garry, and J. Zhang designed the study and wrote the manuscript.Sources of FundingThis work was supported by National Institutes of Health grants R01 HL114120, R01 HL 131017, R01 HL149137, and U01 HL134764.Nonstandard Abbreviations and AcronymsARapical resectionMImyocardial infarctionNTRK2neurotrophic receptor tyrosine kinase 2P1postnatal day 1P28postnatal day 28P30postnatal day 30P56postnatal day 56PKM2pyruvate kinase M2Disclosures None.FootnotesFor Sources of Funding and Disclosures, see page 1747.https://www.ahajournals.org/journal/circCorrespondence to: Jianyi (Jay) Zhang, MD, PhD, Department of Biomedical Engineering, School of Medicine and School of Engineering, University of Alabama at Birmingham, 1670 University Blvd, Volker Hall G094J, Birmingham, AL 35233. Email [email protected]eduReferences1. Zhu W, Zhang E, Zhao M, Chong Z, Fan C, Tang Y, Hunter JD, Borovjagin AV, Walcott GP, Chen JY, et al. Regenerative potential of neonatal porcine hearts.Circulation. 2018; 138:2809–2816. doi: 10.1161/CIRCULATIONAHA.118.034886LinkGoogle Scholar2. Zhao M, Zhang E, Wei Y, Zhou Y, Walcott GP, Zhang J. Apical resection prolongs the cell cycle activity and promotes myocardial regeneration after left ventricular injury in neonatal pig.Circulation. 2020; 142:913–916. doi: 10.1161/CIRCULATIONAHA.119.044619LinkGoogle Scholar3. Haghverdi L, Buettner F, Theis FJ. Diffusion maps for high-dimensional single-cell analysis of differentiation data.Bioinformatics. 2015; 31:2989–2998. doi: 10.1093/bioinformatics/btv325CrossrefMedlineGoogle Scholar4. Bi J, Bennett K, Embrechts M, Breneman C, Song M. Dimensionality reduction via sparse support vector machines.J Mach Learn Res. 2003; 3:1229–1243. doi: 10.1162/153244303322753643Google Scholar5. Chkifa A, Cohen A, Schwab C. Breaking the curse of dimensionality in sparse polynomial approximation of parametric PDEs.J Maths Pures Appls. 2015; 103:400–428. doi: 10.1016/j.matpur.2014.04.009CrossrefGoogle Scholar Previous Back to top Next FiguresReferencesRelatedDetails June 7, 2022Vol 145, Issue 23 Advertisement Article InformationMetrics © 2022 American Heart Association, Inc.https://doi.org/10.1161/CIRCULATIONAHA.121.056995PMID: 35666813 Originally publishedJune 6, 2022 Keywordsmyocytes, cardiaccoronary vesselsregenerationswinePDF download Advertisement SubjectsMechanismsMyocardial InfarctionMyocardial Regeneration