Cardiomyocyte heterogeneity during zebrafish development and regeneration

ZebraReg-a novel platform for discovering regulators of cardiac regeneration using zebrafish

Cardiovascular disease is the leading cause of death worldwide with myocardial infarction being the most prevalent. Currently, no cure is available to either prevent or revert the massive death of cardiomyocytes that occurs after a myocardial infarction. Adult mammalian hearts display a limited regeneration capacity, but it is insufficient to allow complete myocardial recovery. In contrast, the injured zebrafish heart muscle regenerates efficiently through robust proliferation of pre-existing myocardial cells. Thus, zebrafish allows its exploitation for studying the genetic programs behind cardiac regeneration, which may be present, albeit dormant, in the adult human heart. To this end, we have established ZebraReg, a novel and versatile automated platform for studying heart regeneration kinetics after the specific ablation of cardiomyocytes in zebrafish larvae. In combination with automated heart imaging, the platform can be integrated with genetic or pharmacological approaches and used for medium-throughput screening of presumed modulators of heart regeneration. We demonstrate the versatility of the platform by identifying both anti- and pro-regenerative effects of genes and drugs. In conclusion, we present a tool which may be utilised to streamline the process of target validation of novel gene regulators of regeneration, and the discovery of new drug therapies to regenerate the heart after myocardial infarction.

Distinct features of the regenerating heart uncovered through comparative single-cell profiling

Adult humans respond to heart injury by forming a permanent scar, yet other vertebrates are capable of robust and complete cardiac regeneration. Despite progress towards characterizing the mechanisms of cardiac regeneration in fish and amphibians, the large evolutionary gulf between mammals and regenerating vertebrates complicates deciphering which cellular and molecular features truly enable regeneration. To better define these features, we compared cardiac injury responses in zebrafish and medaka, two fish species that share similar heart anatomy and common teleost ancestry but differ in regenerative capability. We used single-cell transcriptional profiling to create a time-resolved comparative cell atlas of injury responses in all major cardiac cell types across both species. With this approach, we identified several key features that distinguish cardiac injury response in the non-regenerating medaka heart. By comparing immune responses to injury, we found altered cell recruitment and a distinct pro-inflammatory gene program in medaka leukocytes, and an absence of the injury-induced interferon response seen in zebrafish. In addition, we found a lack of pro-regenerative signals, including nrg1 and retinoic acid, from medaka endothelial and epicardial cells. Finally, we identified alterations in the myocardial structure in medaka, where they lack primordial layer cardiomyocytes and fail to employ a cardioprotective gene program shared by regenerating vertebrates. Our findings reveal notable variation in injury response across nearly all major cardiac cell types in zebrafish and medaka, demonstrating how evolutionary divergence influences the hidden cellular features underpinning regenerative potential in these seemingly similar vertebrates.

Distinct features of the regenerating heart uncovered through comparative single-cell profiling

Adult humans respond to heart injury by forming a permanent scar, yet other vertebrates are capable of robust and complete cardiac regeneration. Despite progress towards characterizing the mechanisms of cardiac regeneration in fish and amphibians, the large evolutionary gulf between mammals and regenerating vertebrates complicates deciphering which cellular and molecular features truly enable regeneration. To better define these features, we compared cardiac injury responses in zebrafish and medaka, two fish species that share similar heart anatomy and common teleost ancestry but differ in regenerative capability. We used single-cell transcriptional profiling to create a time-resolved comparative cell atlas of injury responses in all major cardiac cell types across both species. With this approach, we identified several key features that distinguish cardiac injury response in the non-regenerating medaka heart. By comparing immune responses to injury, we found altered cell recruitment and a distinct pro-inflammatory gene program in medaka leukocytes, and an absence of the injury-induced interferon response seen in zebrafish. In addition, we found a lack of pro-regenerative signals, including nrg1 and retinoic acid, from medaka endothelial and epicardial cells. Finally, we identified alterations in the myocardial structure in medaka, where they lack embryonic-like primordial layer cardiomyocytes, and fail to employ a cardioprotective gene program shared by regenerating vertebrates. Our findings reveal notable variation in injury response across nearly all major cardiac cell types in zebrafish and medaka, demonstrating how evolutionary divergence influences the hidden cellular features underpinning regenerative potential in these seemingly similar vertebrates.

Analysis of cardiac single-cell RNA-sequencing data can be improved by the use of artificial-intelligence-based tools

Single-cell RNA sequencing (scRNAseq) enables researchers to identify and characterize populations and subpopulations of different cell types in hearts recovering from myocardial infarction (MI) by characterizing the transcriptomes in thousands of individual cells. However, the effectiveness of the currently available tools for processing and interpreting these immense datasets is limited. We incorporated three Artificial Intelligence (AI) techniques into a toolkit for evaluating scRNAseq data: AI Autoencoding separates data from different cell types and subpopulations of cell types (cluster analysis); AI Sparse Modeling identifies genes and signaling mechanisms that are differentially activated between subpopulations (pathway/gene set enrichment analysis), and AI Semisupervised Learning tracks the transformation of cells from one subpopulation into another (trajectory analysis). Autoencoding was often used in data denoising; yet, in our pipeline, Autoencoding was exclusively used for cell embedding and clustering. The performance of our AI scRNAseq toolkit and other highly cited non-AI tools was evaluated with three scRNAseq datasets obtained from the Gene Expression Omnibus database. Autoencoder was the only tool to identify differences between the cardiomyocyte subpopulations found in mice that underwent MI or sham-MI surgery on postnatal day (P) 1. Statistically significant differences between cardiomyocytes from P1-MI mice and mice that underwent MI on P8 were identified for six cell-cycle phases and five signaling pathways when the data were analyzed via Sparse Modeling, compared to just one cell-cycle phase and one pathway when the data were analyzed with non-AI techniques. Only Semisupervised Learning detected trajectories between the predominant cardiomyocyte clusters in hearts collected on P28 from pigs that underwent apical resection (AR) on P1, and on P30 from pigs that underwent AR on P1 and MI on P28. In another dataset, the pig scRNAseq data were collected after the injection of CCND2-overexpression Human-induced Pluripotent Stem Cell-derived cardiomyocytes (^CCND2hiPSC) into injured P28 pig heart; only the AI-based technique could demonstrate that the host cardiomyocytes increase proliferating by through the HIPPO/YAP and MAPK signaling pathways. For the cluster, pathway/gene set enrichment, and trajectory analysis of scRNAseq datasets generated from studies of myocardial regeneration in mice and pigs, our AI-based toolkit identified results that non-AI techniques did not discover. These different results were validated and were important in explaining myocardial regeneration.

Long-term culture of patient-derived cardiac organoids recapitulated Duchenne muscular dystrophy cardiomyopathy and disease progression

Duchenne Muscular Dystrophy (DMD) is an X-linked neuromuscular disease which to date is incurable. The major cause of death is dilated cardiomyopathy however, its pathogenesis is unclear as existing cellular and animal models do not fully recapitulate the human disease phenotypes. In this study, we generated cardiac organoids from patient-derived induced pluripotent stem cells (DMD-COs) and isogenic-corrected controls (DMD-Iso-COs) and studied if DMD-related cardiomyopathy and disease progression occur in the organoids upon long-term culture (up to 93 days). Histological analysis showed that DMD-COs lack initial proliferative capacity, displayed a progressive loss of sarcoglycan localization and high stress in endoplasmic reticulum. Additionally, cardiomyocyte deterioration, fibrosis and aberrant adipogenesis were observed in DMD-COs over time. RNA sequencing analysis confirmed a distinct transcriptomic profile in DMD-COs which was associated with functional enrichment in hypertrophy/dilated cardiomyopathy, arrhythmia, adipogenesis and fibrosis pathways. Moreover, five miRNAs were identified to be crucial in this dysregulated gene network. In conclusion, we generated patient-derived cardiac organoid model that displayed DMD-related cardiomyopathy and disease progression phenotypes in long-term culture. We envision the feasibility to develop a more complex, realistic and reliable in vitro 3D human cardiac-mimics to study DMD-related cardiomyopathies.

Cardiomyocyte Cell-Cycle Regulation in Neonatal Large Mammals: Single Nucleus RNA-Sequencing Data Analysis via an Artificial-Intelligence–Based Pipeline

Adult mammalian cardiomyocytes have very limited capacity to proliferate and repair the myocardial infarction. However, when apical resection (AR) was performed in pig hearts on postnatal day (P) 1 (AR_P1) and acute myocardial infarction (MI) was induced on P28 (MI_P28), the animals recovered with no evidence of myocardial scarring or decline in contractile performance. Furthermore, the repair process appeared to be driven by cardiomyocyte proliferation, but the regulatory molecules that govern the AR_P1-induced enhancement of myocardial recovery remain unclear. Single-nucleus RNA sequencing (snRNA-seq) data collected from fetal pig hearts and the hearts of pigs that underwent AR_P1, MI_P28, both AR_P1 and MI, or neither myocardial injury were evaluated via autoencoder, cluster analysis, sparse learning, and semisupervised learning. Ten clusters of cardiomyocytes (CM1–CM10) were identified across all experimental groups and time points. CM1 was only observed in AR_P1 hearts on P28 and was enriched for the expression of T-box transcription factors 5 and 20 (TBX5 and TBX20, respectively), Erb-B2 receptor tyrosine kinase 4 (ERBB4), and G Protein-Coupled Receptor Kinase 5 (GRK5), as well as genes associated with the proliferation and growth of cardiac muscle. CM1 cardiomyocytes also highly expressed genes for glycolysis while lowly expressed genes for adrenergic signaling, which suggested that CM1 were immature cardiomyocytes. Thus, we have identified a cluster of cardiomyocytes, CM1, in neonatal pig hearts that appeared to be generated in response to AR injury on P1 and may have been primed for activation of CM cell-cycle activation and proliferation by the upregulation of TBX5, TBX20, ERBB4, and GRK5.

A Vegfc-Emilin2a-Cxcl8a Signaling Axis Required for Zebrafish Cardiac Regeneration

BACKGROUND

Ischemic heart disease following the obstruction of coronary vessels leads to the death of cardiac tissue and the formation of a fibrotic scar. In contrast to adult mammals, zebrafish can regenerate their heart after injury, enabling the study of the underlying mechanisms. One of the earliest responses following cardiac injury in adult zebrafish is coronary revascularization. Defects in this process lead to impaired cardiomyocyte repopulation and scarring. Hence, identifying and investigating factors that promote coronary revascularization holds great therapeutic potential.

METHODS

We used wholemount imaging, immunohistochemistry and histology to assess various aspects of zebrafish cardiac regeneration. Deep transcriptomic analysis allowed us to identify targets and potential effectors of Vegfc (vascular endothelial growth factor C) signaling. We used newly generated loss- and gain-of-function genetic tools to investigate the role of Emilin2a (elastin microfibril interfacer 2a) and Cxcl8a (chemokine (C-X-C) motif ligand 8a)-Cxcr1 (chemokine (C-X-C) motif receptor 1) signaling in cardiac regeneration.

RESULTS

We first show that regenerating coronary endothelial cells upregulate vegfc upon cardiac injury in adult zebrafish and that Vegfc signaling is required for their proliferation during regeneration. Notably, blocking Vegfc signaling also significantly reduces cardiomyocyte dedifferentiation and proliferation. Using transcriptomic analyses, we identified emilin2a as a target of Vegfc signaling and found that manipulation of emilin2a expression can modulate coronary revascularization as well as cardiomyocyte proliferation. Mechanistically, Emilin2a induces the expression of the chemokine gene cxcl8a in epicardium-derived cells, while cxcr1, the Cxcl8a receptor gene, is expressed in coronary endothelial cells. We further show that Cxcl8a-Cxcr1 signaling is also required for coronary endothelial cell proliferation during cardiac regeneration.

CONCLUSIONS

These data show that after cardiac injury, coronary endothelial cells upregulate vegfc to promote coronary network reestablishment and cardiac regeneration. Mechanistically, Vegfc signaling upregulates epicardial emilin2a and cxcl8a expression to promote cardiac regeneration. These studies aid in understanding the mechanisms underlying coronary revascularization in zebrafish, with potential therapeutic implications to enhance revascularization and regeneration in injured human hearts.

The stress responsive gene ankrd1a is dynamically regulated during skeletal muscle development and upregulated following cardiac injury in border zone cardiomyocytes in adult zebrafish

Ankyrin repeat domain 1 (ANKRD1) is a functionally pleiotropic protein found in the nuclei and sarcomeres of cardiac and skeletal muscles, with a proposed role in linking myofibrilar stress and transcriptional regulation. Rapid upregulation of its expression in response to both physiological and pathological stress supports the involvement of ANKRD1 in muscle tissue adaptation and remodeling. However, the exact role of ANKRD1 remains poorly understood. To begin to investigate its function at higher resolution, we have generated and characterized a TgBAC(ankrd1a:EGFP) zebrafish line. This reporter line displays transgene expression in slow skeletal muscle fibers during development and exercise responsiveness in adult cardiac muscle. To better understand the role of Ankrd1a in pathological conditions in adult zebrafish, we assessed ankrd1a expression after cardiac ventricle cryoinjury and observed localized upregulation in cardiomyocytes in the border zone. We show that this expression in injured hearts is recapitulated by the TgBAC(ankrd1a:EGFP) reporter. Our results identify novel expression domains of ankrd1a and suggest an important role for Ankrd1a in the early stress response and regeneration of cardiac tissue. This new reporter line will help decipher the role of Ankrd1a in striated muscle stress response, including after cardiac injury.