Spatially resolved RNA-sequencing of the embryonic heart identifies a role for Wnt/β-catenin signaling in autonomic control of heart rate

scRNA-seq reveals the diversity of the developing cardiac cell lineage and molecular players in heart rhythm regulation

We utilized scRNA-seq to delineate the diversity of cell types in the zebrafish heart. Transcriptome profiling of over 50,000 cells at 48 and 72 hpf defined at least 18 discrete cell lineages of the developing heart. Utilizing well-established gene signatures, we identified a population of cells likely to be the primary pacemaker and characterized the transcriptome profile defining this critical cell type. Two previously uncharacterized genes, atp1b3b and colec10, were found to be enriched in the sinoatrial cardiomyocytes. CRISPR/Cas9-mediated knockout of these two genes significantly reduced heart rate, implicating their role in cardiac development and conduction. Additionally, we describe other cardiac cell lineages, including the endothelial and neural cells, providing their expression profiles as a resource. Our results established a detailed atlas of the developing heart, providing valuable insights into cellular and molecular mechanisms, and pinpointed potential new players in heart rhythm regulation.

Single cell RNA-sequencing and RNA-tomography of the avian embryo extending body axis

Introduction: Vertebrate body axis formation initiates during gastrulation and continues within the tail bud at the posterior end of the embryo. Major structures in the trunk are paired somites, which generate the musculoskeletal system, the spinal cord-forming part of the central nervous system, and the notochord, with important patterning functions. The specification of these different cell lineages by key signalling pathways and transcription factors is essential, however, a global map of cell types and expressed genes in the avian trunk is missing. Methods: Here we use high-throughput sequencing approaches to generate a molecular map of the emerging trunk and tailbud in the chick embryo. Results and Discussion: Single cell RNA-sequencing (scRNA-seq) identifies discrete cell lineages including somites, neural tube, neural crest, lateral plate mesoderm, ectoderm, endothelial and blood progenitors. In addition, RNA-seq of sequential tissue sections (RNA-tomography) provides a spatially resolved, genome-wide expression dataset for the avian tailbud and emerging body, comparable to other model systems. Combining the single cell and RNA-tomography datasets, we identify spatially restricted genes, focusing on somites and early myoblasts. Thus, this high-resolution transcriptome map incorporating cell types in the embryonic trunk can expose molecular pathways involved in body axis development.

Development of the Cardiac Conduction System

The electrical impulses that coordinate the sequential, rhythmic contractions of the atria and ventricles are initiated and tightly regulated by the specialized tissues of the cardiac conduction system. In the mature heart, these impulses are generated by the pacemaker cardiomyocytes of the sinoatrial node, propagated through the atria to the atrioventricular node where they are delayed and then rapidly propagated to the atrioventricular bundle, right and left bundle branches, and finally, the peripheral ventricular conduction system. Each of these specialized components arise by complex patterning events during embryonic development. This chapter addresses the origins and transcriptional networks and signaling pathways that drive the development and maintain the function of the cardiac conduction system.

Single-cell profiling of the developing embryonic heart in Drosophila

Drosophila is an important model for studying heart development and disease. Yet, single-cell transcriptomic data of its developing heart have not been performed. Here, we report single-cell profiling of the entire fly heart using ∼3000 Hand-GFP embryos collected at five consecutive developmental stages, ranging from bilateral migrating rows of cardiac progenitors to a fused heart tube. The data revealed six distinct cardiac cell types in the embryonic fly heart: cardioblasts, both Svp+ and Tin+ subtypes; and five types of pericardial cell (PC) that can be distinguished by four key transcription factors (Eve, Odd, Ct and Tin) and include the newly described end of the line PC. Notably, the embryonic fly heart combines transcriptional signatures of the mammalian first and second heart fields. Using unique markers for each heart cell type, we defined their number and location during heart development to build a comprehensive 3D cell map. These data provide a resource to track the expression of any gene in the developing fly heart, which can serve as a reference to study genetic perturbations and cardiac diseases.

Cfdp1 Is Essential for Cardiac Development and Function

Cardiovascular diseases (CVDs) are the prevalent cause of mortality worldwide. A combination of environmental and genetic effectors modulates the risk of developing them. Thus, it is vital to identify candidate genes and elucidate their role in the manifestation of the disease. Large-scale human studies have revealed the implication of Craniofacial Development Protein 1 (CFDP1) in Coronary Artery Disease (CAD). CFDP1 belongs to the evolutionary conserved Bucentaur (BCNT) family, and to date, its function and mechanism of action in Cardiovascular Development are still unclear. We utilized zebrafish to investigate the role of cfdp1 in the developing heart due to the high genomic homology, similarity in heart physiology, and ease of experimental manipulations. We showed that cfdp1 was expressed during development, and we tested two morpholinos and generated a cfdp1 mutant line. The cfdp1-/- embryos developed arrhythmic hearts and exhibited defective cardiac performance, which led to a lethal phenotype. Findings from both knockdown and knockout experiments showed that abrogation of cfdp1 leads to downregulation of Wnt signaling in embryonic hearts during valve development but without affecting Notch activation in this process. The cfdp1 zebrafish mutant line provides a valuable tool for unveiling the novel mechanism of regulating cardiac physiology and function. cfdp1 is essential for cardiac development, a previously unreported phenotype most likely due to early lethality in mice. The detected phenotype of bradycardia and arrhythmias is an observation with potential clinical relevance for humans carrying heterozygous CFDP1 mutations and their risk of developing CAD.

Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes of zebrafish

Maintenance of cardiomyocyte identity is vital for normal heart development and function. However, our understanding of cardiomyocyte plasticity remains incomplete. Here, we show that sustained expression of the zebrafish transcription factor Nr2f1a prevents the progressive acquisition of ventricular cardiomyocyte (VC) and pacemaker cardiomyocyte (PC) identities within distinct regions of the atrium. Transcriptomic analysis of flow-sorted atrial cardiomyocytes (ACs) from nr2f1a mutant zebrafish embryos showed increased VC marker gene expression and altered expression of core PC regulatory genes, including decreased expression of nkx2.5, a critical repressor of PC differentiation. At the arterial (outflow) pole of the atrium in nr2f1a mutants, cardiomyocytes resolve to VC identity within the expanded atrioventricular canal. However, at the venous (inflow) pole of the atrium, there is a progressive wave of AC transdifferentiation into PCs across the atrium toward the arterial pole. Restoring Nkx2.5 is sufficient to repress PC marker identity in nr2f1a mutant atria and analysis of chromatin accessibility identified an Nr2f1a-dependent nkx2.5 enhancer expressed in the atrial myocardium directly adjacent to PCs. CRISPR/Cas9-mediated deletion of the putative nkx2.5 enhancer leads to a loss of Nkx2.5-expressing ACs and expansion of a PC reporter, supporting that Nr2f1a limits PC differentiation within venous ACs via maintaining nkx2.5 expression. The Nr2f-dependent maintenance of AC identity within discrete atrial compartments may provide insights into the molecular etiology of concurrent structural congenital heart defects and associated arrhythmias.

The opportunities and challenges of using Drosophila to model human cardiac diseases

The Drosophila heart tube seems simple, yet it has notable anatomic complexity and contains highly specialized structures. In fact, the development of the fly heart tube much resembles that of the earliest stages of mammalian heart development, and the molecular-genetic mechanisms driving these processes are highly conserved between flies and humans. Combined with the fly's unmatched genetic tools and a wide variety of techniques to assay both structure and function in the living fly heart, these attributes have made Drosophila a valuable model system for studying human heart development and disease. This perspective focuses on the functional and physiological similarities between fly and human hearts. Further, it discusses current limitations in using the fly, as well as promising prospects to expand the capabilities of Drosophila as a research model for studying human cardiac diseases.

WNT/β-catenin pathway is a key regulator of cardiac function and energetic metabolism

The WNT/β-catenin pathway is a master regulator of cardiac development and growth, and its activity is low in healthy adult hearts. However, even this low activity is essential for maintaining normal heart function. Acute activation of the WNT/β-catenin signaling cascade is considered to be cardioprotective after infarction through the upregulation of prosurvival genes and reprogramming of metabolism. Chronically high WNT/β-catenin pathway activity causes profibrotic and hypertrophic effects in the adult heart. New data suggest more complex functions of β-catenin in metabolic maturation of the perinatal heart, establishing an adult pattern of glucose and fatty acid utilization. Additionally, low basal activity of the WNT/β-catenin cascade maintains oxidative metabolism in the adult heart, and this pathway is reactivated by physiological or pathological stimuli to meet the higher energy needs of the heart. This review summarizes the current state of knowledge of the organization of canonical WNT signaling and its function in cardiogenesis, heart maturation, adult heart function, and remodeling. We also discuss the role of the WNT/β-catenin pathway in cardiac glucose, lipid metabolism, and mitochondrial physiology.

Advances and Challenges in Spatial Transcriptomics for Developmental Biology

Development from single cells to multicellular tissues and organs involves more than just the exact replication of cells, which is known as differentiation. The primary focus of research into the mechanism of differentiation has been differences in gene expression profiles between individual cells. However, it has predominantly been conducted at low throughput and bulk levels, challenging the efforts to understand molecular mechanisms of differentiation during the developmental process in animals and humans. During the last decades, rapid methodological advancements in genomics facilitated the ability to study developmental processes at a genome-wide level and finer resolution. Particularly, sequencing transcriptomes at single-cell resolution, enabled by single-cell RNA-sequencing (scRNA-seq), was a breath-taking innovation, allowing scientists to gain a better understanding of differentiation and cell lineage during the developmental process. However, single-cell isolation during scRNA-seq results in the loss of the spatial information of individual cells and consequently limits our understanding of the specific functions of the cells performed by different spatial regions of tissues or organs. This greatly encourages the emergence of the spatial transcriptomic discipline and tools. Here, we summarize the recent application of scRNA-seq and spatial transcriptomic tools for developmental biology. We also discuss the limitations of current spatial transcriptomic tools and approaches, as well as possible solutions and future prospects.

Full-Length Spatial Transcriptomics Reveals the Unexplored Isoform Diversity of the Myocardium Post-MI

We introduce Single-cell Nanopore Spatial Transcriptomics (scNaST), a software suite to facilitate the analysis of spatial gene expression from second- and third-generation sequencing, allowing to generate a full-length near-single-cell transcriptional landscape of the tissue microenvironment. Taking advantage of the Visium Spatial platform, we adapted a strategy recently developed to assign barcodes to long-read single-cell sequencing data for spatial capture technology. Here, we demonstrate our workflow using four short axis sections of the mouse heart following myocardial infarction. We constructed a de novo transcriptome using long-read data, and successfully assigned 19,794 transcript isoforms in total, including clinically-relevant, but yet uncharacterized modes of transcription, such as intron retention or antisense overlapping transcription. We showed a higher transcriptome complexity in the healthy regions, and identified intron retention as a mode of transcription associated with the infarct area. Our data revealed a clear regional isoform switching among differentially used transcripts for genes involved in cardiac muscle contraction and tissue morphogenesis. Molecular signatures involved in cardiac remodeling integrated with morphological context may support the development of new therapeutics towards the treatment of heart failure and the reduction of cardiac complications.

Stx4 is required to regulate cardiomyocyte Ca2+ handling during vertebrate cardiac development

Requirements for vesicle fusion within the heart remain poorly understood, despite the multitude of processes that necessitate proper intracellular trafficking within cardiomyocytes. Here, we show that Syntaxin 4 (STX4), a target-Soluble N-ethylmaleimide sensitive factor attachment receptor (t-SNARE) protein, is required for normal vertebrate cardiac conduction and vesicular transport. Two patients were identified with damaging variants in STX4. A patient with a homozygous R240W missense variant displayed biventricular dilated cardiomyopathy, ectopy, and runs of non-sustained ventricular tachycardia, sensorineural hearing loss, global developmental delay, and hypotonia, while a second patient displayed severe pleiotropic abnormalities and perinatal lethality. CRISPR/Cas9-generated stx4 mutant zebrafish exhibited defects reminiscent of these patients' clinical presentations, including linearized hearts, bradycardia, otic vesicle dysgenesis, neuronal atrophy, and touch insensitivity by 3 days post fertilization. Imaging of Vamp2+ vesicles within stx4 mutant zebrafish hearts showed reduced docking to the cardiomyocyte sarcolemma. Optical mapping of the embryonic hearts coupled with pharmacological modulation of Ca²⁺ handling together support that zebrafish stx4 mutants have a reduction in L-type Ca²⁺ channel modulation. Transgenic overexpression of zebrafish Stx4^R241W, analogous to the first patient's STX4^R240W variant, indicated that the variant is hypomorphic. Thus, these data show an in vivo requirement for SNAREs in regulating normal embryonic cardiac function and that variants in STX4 are associated with pleiotropic human disease, including cardiomyopathy.

Transcriptome-scale methods for uncovering subcellular RNA localization mechanisms

Across a variety of systems, thousands of RNAs are localized to specific subcellular locations. However, for the vast majority of these RNAs, the mechanisms that underlie their transport are unknown. Historically, these mechanisms were uncovered for a single transcript at a time by laboriously testing the ability of RNA fragments to direct transcript localization. Recently developed methods profile the content of subcellular transcriptomes using high-throughput sequencing, allowing the analysis of the localization of thousands of transcripts at once. By identifying commonalities shared among multiple localized transcripts, these methods have the potential to rapidly expand our understanding of RNA localization mechanisms.

A dual SHOX2:GFP; MYH6:mCherry knockin hESC reporter line for derivation of human SAN-like cells

The sinoatrial node (SAN) is the primary pacemaker of the heart. The human SAN is poorly understood due to limited primary tissue access and limitations in robust in vitro derivation methods. We developed a dual SHOX2:GFP; MYH6:mCherry knockin human embryonic stem cell (hESC) reporter line, which allows the identification and purification of SAN-like cells. Using this line, we performed several rounds of chemical screens and developed an efficient strategy to generate and purify hESC-derived SAN-like cells (hESC-SAN). The derived hESC-SAN cells display molecular and electrophysiological characteristics of bona fide nodal cells, which allowed exploration of their transcriptional profile at single-cell level. In sum, our dual reporter system facilitated an effective strategy for deriving human SAN-like cells, which can potentially be used for future disease modeling and drug discovery.

Asymmetric Hapln1a drives regionalized cardiac ECM expansion and promotes heart morphogenesis in zebrafish development

AIMS

Vertebrate heart development requires the complex morphogenesis of a linear tube to form the mature organ, a process essential for correct cardiac form and function, requiring coordination of embryonic laterality, cardiac growth, and regionalized cellular changes. While previous studies have demonstrated broad requirements for extracellular matrix (ECM) components in cardiac morphogenesis, we hypothesized that ECM regionalization may fine tune cardiac shape during heart development.

METHODS AND RESULTS

Using live in vivo light sheet imaging of zebrafish embryos, we describe a left-sided expansion of the ECM between the myocardium and endocardium prior to the onset of heart looping and chamber ballooning. Analysis using an ECM sensor revealed the cardiac ECM is further regionalized along the atrioventricular axis. Spatial transcriptomic analysis of gene expression in the heart tube identified candidate genes that may drive ECM expansion. This approach identified regionalized expression of hapln1a, encoding an ECM cross-linking protein. Validation of transcriptomic data by in situ hybridization confirmed regionalized hapln1a expression in the heart, with highest levels of expression in the future atrium and on the left side of the tube, overlapping with the observed ECM expansion. Analysis of CRISPR-Cas9-generated hapln1a mutants revealed a reduction in atrial size and reduced chamber ballooning. Loss-of-function analysis demonstrated that ECM expansion is dependent upon Hapln1a, together supporting a role for Hapln1a in regionalized ECM modulation and cardiac morphogenesis. Analysis of hapln1a expression in zebrafish mutants with randomized or absent embryonic left-right asymmetry revealed that laterality cues position hapln1a-expressing cells asymmetrically in the left side of the heart tube.

CONCLUSION

We identify a regionalized ECM expansion in the heart tube which promotes correct heart development, and propose a novel model whereby embryonic laterality cues orient the axis of ECM asymmetry in the heart, suggesting these two pathways interact to promote robust cardiac morphogenesis.

Modeling Human Cardiac Arrhythmias: Insights from Zebrafish

Cardiac arrhythmia, or irregular heart rhythm, is associated with morbidity and mortality and is described as one of the most important future public health challenges. Therefore, developing new models of cardiac arrhythmia is critical for understanding disease mechanisms, determining genetic underpinnings, and developing new therapeutic strategies. In the last few decades, the zebrafish has emerged as an attractive model to reproduce in vivo human cardiac pathologies, including arrhythmias. Here, we highlight the contribution of zebrafish to the field and discuss the available cardiac arrhythmia models. Further, we outline techniques to assess potential heart rhythm defects in larval and adult zebrafish. As genetic tools in zebrafish continue to bloom, this model will be crucial for functional genomics studies and to develop personalized anti-arrhythmic therapies.

When Bigger Is Better: 3D RNA Profiling of the Developing Head in the Catshark Scyliorhinus canicula

We report the adaptation of RNA tomography, a technique allowing spatially resolved, genome-wide expression profiling, to a species occupying a key phylogenetic position in gnathostomes, the catshark Scyliorhinus canicula. We focused analysis on head explants at an embryonic stage, shortly following neural tube closure and of interest for a number of developmental processes, including early brain patterning, placode specification or the establishment of epithalamic asymmetry. As described in the zebrafish, we have sequenced RNAs extracted from serial sections along transverse, horizontal and sagittal planes, mapped the data onto a gene reference taking advantage of the high continuity genome recently released in the catshark, and projected read counts onto a digital model of the head obtained by confocal microscopy. This results in the generation of a genome-wide 3D atlas, containing expression data for most protein-coding genes in a digital model of the embryonic head. The digital profiles obtained for candidate forebrain regional markers along antero-posterior, dorso-ventral and left-right axes reproduce those obtained by in situ hybridization (ISH), with expected relative organizations. We also use spatial autocorrelation and correlation as measures to analyze these data and show that they provide adequate statistical tools to extract novel expression information from the model. These data and tools allow exhaustive searches of genes exhibiting any predefined expression characteristic, such a restriction to a territory of interest, thus providing a reference for comparative analyses across gnathostomes. This methodology appears best suited to species endowed with large embryo or organ sizes and opens novel perspectives to a wide range of evo-devo model organisms, traditionally counter-selected on size criterion.

Identification of astroglia-like cardiac nexus glia that are critical regulators of cardiac development and function

Glial cells are essential for functionality of the nervous system. Growing evidence underscores the importance of astrocytes; however, analogous astroglia in peripheral organs are poorly understood. Using confocal time-lapse imaging, fate mapping, and mutant genesis in a zebrafish model, we identify a neural crest-derived glial cell, termed nexus glia, which utilizes Meteorin signaling via Jak/Stat3 to drive differentiation and regulate heart rate and rhythm. Nexus glia are labeled with gfap, glast, and glutamine synthetase, markers that typically denote astroglia cells. Further, analysis of single-cell sequencing datasets of human and murine hearts across ages reveals astrocyte-like cells, which we confirm through a multispecies approach. We show that cardiac nexus glia at the outflow tract are critical regulators of both the sympathetic and parasympathetic system. These data establish the crucial role of glia on cardiac homeostasis and provide a description of nexus glia in the PNS.

Mechanosensitive Notch-Dll4 and Klf2-Wnt9 signaling pathways intersect in guiding valvulogenesis in zebrafish

In the zebrafish embryo, the onset of blood flow generates fluid shear stress on endocardial cells, which are specialized endothelial cells that line the interior of the heart. High levels of fluid shear stress activate both Notch and Klf2 signaling, which play crucial roles in atrioventricular valvulogenesis. However, it remains unclear why only individual endocardial cells ingress into the cardiac jelly and initiate valvulogenesis. Here, we show that lateral inhibition between endocardial cells, mediated by Notch, singles out Delta-like-4-positive endocardial cells. These cells ingress into the cardiac jelly, where they form an abluminal cell population. Delta-like-4-positive cells ingress in response to Wnt9a, which is produced in parallel through an Erk5-Klf2-Wnt9a signaling cascade also activated by blood flow. Hence, mechanical stimulation activates parallel mechanosensitive signaling pathways that produce binary effects by driving endocardial cells toward either luminal or abluminal fates. Ultimately, these cell fate decisions sculpt cardiac valve leaflets.

3D-cardiomics: A spatial transcriptional atlas of the mammalian heart

Understanding the spatial gene expression and regulation in the heart is key to uncovering its developmental and physiological processes, during homeostasis and disease. Numerous techniques exist to gain gene expression and regulation information in organs such as the heart, but few utilize intuitive true-to-life three-dimensional representations to analyze and visualise results. Here we combined transcriptomics with 3D-modelling to interrogate spatial gene expression in the mammalian heart. For this, we microdissected and sequenced transcriptome-wide 18 anatomical sections of the adult mouse heart. Our study has unveiled known and novel genes that display complex spatial expression in the heart sub-compartments. We have also created 3D-cardiomics, an interface for spatial transcriptome analysis and visualization that allows the easy exploration of these data in a 3D model of the heart. 3D-cardiomics is accessible from http://3d-cardiomics.erc.monash.edu/.

Transcriptome profile of the sinoatrial ring reveals conserved and novel genetic programs of the zebrafish pacemaker

BACKGROUND

Sinoatrial Node (SAN) is part of the cardiac conduction system, which controls the rhythmic contraction of the vertebrate heart. The SAN consists of a specialized pacemaker cell population that has the potential to generate electrical impulses. Although the SAN pacemaker has been extensively studied in mammalian and teleost models, including the zebrafish, their molecular nature remains inadequately comprehended.

RESULTS

To characterize the molecular profile of the zebrafish sinoatrial ring (SAR) and elucidate the mechanism of pacemaker function, we utilized the transgenic line sqet33mi59BEt to isolate cells of the SAR of developing zebrafish embryos and profiled their transcriptome. Our analyses identified novel candidate genes and well-known conserved signaling pathways involved in pacemaker development. We show that, compared to the rest of the heart, the zebrafish SAR overexpresses several mammalian SAN pacemaker signature genes, which include hcn4 as well as those encoding calcium- and potassium-gated channels. Moreover, genes encoding components of the BMP and Wnt signaling pathways, as well as members of the Tbx family, which have previously been implicated in pacemaker development, were also overexpressed in the SAR. Among SAR-overexpressed genes, 24 had human homologues implicated in 104 different ClinVar phenotype entries related to various forms of congenital heart diseases, which suggest the relevance of our transcriptomics resource to studying human heart conditions. Finally, functional analyses of three SAR-overexpressed genes, pard6a, prom2, and atp1a1a.2, uncovered their novel role in heart development and physiology.

CONCLUSION

Our results established conserved aspects between zebrafish and mammalian pacemaker function and revealed novel factors implicated in maintaining cardiac rhythm. The transcriptome data generated in this study represents a unique and valuable resource for the study of pacemaker function and associated heart diseases.

Integrating single-cell and spatial transcriptomics to elucidate intercellular tissue dynamics

Single-cell RNA sequencing (scRNA-seq) identifies cell subpopulations within tissue but does not capture their spatial distribution nor reveal local networks of intercellular communication acting in situ. A suite of recently developed techniques that localize RNA within tissue, including multiplexed in situ hybridization and in situ sequencing (here defined as high-plex RNA imaging) and spatial barcoding, can help address this issue. However, no method currently provides as complete a scope of the transcriptome as does scRNA-seq, underscoring the need for approaches to integrate single-cell and spatial data. Here, we review efforts to integrate scRNA-seq with spatial transcriptomics, including emerging integrative computational methods, and propose ways to effectively combine current methodologies.

Genomic and physiological analyses of the zebrafish atrioventricular canal reveal molecular building blocks of the secondary pacemaker region

The atrioventricular canal (AVC) is the site where key structures responsible for functional division between heart regions are established, most importantly, the atrioventricular (AV) conduction system and cardiac valves. To elucidate the mechanism underlying AVC development and function, we utilized transgenic zebrafish line sqet31Et expressing EGFP in the AVC to isolate this cell population and profile its transcriptome at 48 and 72 hpf. The zebrafish AVC transcriptome exhibits hallmarks of mammalian AV node, including the expression of genes implicated in its development and those encoding connexins forming low conductance gap junctions. Transcriptome analysis uncovered protein-coding and noncoding transcripts enriched in AVC, which have not been previously associated with this structure, as well as dynamic expression of epithelial-to-mesenchymal transition markers and components of TGF-β, Notch, and Wnt signaling pathways likely reflecting ongoing AVC and valve development. Using transgenic line Tg(myl7:mermaid) encoding voltage-sensitive fluorescent protein, we show that abolishing the pacemaker-containing sinoatrial ring (SAR) through Isl1 loss of function resulted in spontaneous activation in the AVC region, suggesting that it possesses inherent automaticity although insufficient to replace the SAR. The SAR and AVC transcriptomes express partially overlapping species of ion channels and gap junction proteins, reflecting their distinct roles. Besides identifying conserved aspects between zebrafish and mammalian conduction systems, our results established molecular hallmarks of the developing AVC which underlies its role in structural and electrophysiological separation between heart chambers. This data constitutes a valuable resource for studying AVC development and function, and identification of novel candidate genes implicated in these processes.

Towards spatio-temporally resolved developmental cardiac gene regulatory networks in zebrafish

Heart formation in the zebrafish involves a rapid, complex series of morphogenetic events in three-dimensional space that spans cardiac lineage specification through to chamber formation and maturation. This process is tightly orchestrated by a cardiac gene regulatory network (GRN), which ensures the precise spatio-temporal deployment of genes critical for heart formation. Alterations of the timing or spatial localisation of gene expression can have a significant impact in cardiac ontogeny and may lead to heart malformations. Hence, a better understanding of the cellular and molecular basis of congenital heart disease relies on understanding the behaviour of cardiac GRNs with precise spatiotemporal resolution. Here, we review the recent technical advances that have expanded our capacity to interrogate the cardiac GRN in zebrafish. In particular, we focus on studies utilising high-throughput technologies to systematically dissect gene expression patterns, both temporally and spatially during heart development.

Bridging scales: From cell biology to physiology using in situ single-cell technologies

Biological organization crosses multiple spatial scales: from molecular, cellular, to tissues and organs. The proliferation of molecular profiling technologies enables increasingly detailed cataloging of the components at each scale. However, the scarcity of spatial profiling has made it challenging to bridge across these scales. Emerging technologies based on highly multiplexed in situ profiling are paving the way to study the spatial organization of cells and tissues in greater detail. These new technologies provide the data needed to cross the scale from cell biology to physiology and identify the fundamental principles that govern tissue organization. Here, we provide an overview of these key technologies and discuss the current and future insights these powerful techniques enable.

From Stripes to a Beating Heart: Early Cardiac Development in Zebrafish

The heart is the first functional organ to form during vertebrate development. Congenital heart defects are the most common type of human birth defect, many originating as anomalies in early heart development. The zebrafish model provides an accessible vertebrate system to study early heart morphogenesis and to gain new insights into the mechanisms of congenital disease. Although composed of only two chambers compared with the four-chambered mammalian heart, the zebrafish heart integrates the core processes and cellular lineages central to cardiac development across vertebrates. The rapid, translucent development of zebrafish is amenable to in vivo imaging and genetic lineage tracing techniques, providing versatile tools to study heart field migration and myocardial progenitor addition and differentiation. Combining transgenic reporters with rapid genome engineering via CRISPR-Cas9 allows for functional testing of candidate genes associated with congenital heart defects and the discovery of molecular causes leading to observed phenotypes. Here, we summarize key insights gained through zebrafish studies into the early patterning of uncommitted lateral plate mesoderm into cardiac progenitors and their regulation. We review the central genetic mechanisms, available tools, and approaches for modeling congenital heart anomalies in the zebrafish as a representative vertebrate model.

Atrial and Sinoatrial Node Development in the Zebrafish Heart

Proper development and function of the vertebrate heart is vital for embryonic and postnatal life. Many congenital heart defects in humans are associated with disruption of genes that direct the formation or maintenance of atrial and pacemaker cardiomyocytes at the venous pole of the heart. Zebrafish are an outstanding model for studying vertebrate cardiogenesis, due to the conservation of molecular mechanisms underlying early heart development, external development, and ease of genetic manipulation. Here, we discuss early developmental mechanisms that instruct appropriate formation of the venous pole in zebrafish embryos. We primarily focus on signals that determine atrial chamber size and the specialized pacemaker cells of the sinoatrial node through directing proper specification and differentiation, as well as contemporary insights into the plasticity and maintenance of cardiomyocyte identity in embryonic zebrafish hearts. Finally, we integrate how these insights into zebrafish cardiogenesis can serve as models for human atrial defects and arrhythmias.

Complementary and divergent functions of zebrafish Tango1 and Ctage5 in tissue development and homeostasis

Coat protein complex II (COPII) factors mediate cargo export from the endoplasmic reticulum (ER), but bulky collagens and lipoproteins are too large for traditional COPII vesicles. Mammalian CTAGE5 and TANGO1 have been well characterized individually as specialized cargo receptors at the ER that function with COPII coats to facilitate trafficking of bulky cargoes. Here, we present a genetic interaction study in zebrafish of deletions in ctage5, tango1, or both to investigate their distinct and complementary potential functions. We found that Ctage5 and Tango1 have different roles related to organogenesis, collagen versus lipoprotein trafficking, stress-pathway activation, and survival. While disruption of both ctage5 and tango1 compounded phenotype severity, mutation of either factor alone revealed novel tissue-specific defects in the building of heart, muscle, lens, and intestine, in addition to previously described roles in the development of neural and cartilage tissues. Together, our results demonstrate that Ctage5 and Tango1 have overlapping functions, but also suggest divergent roles in tissue development and homeostasis.

From whole-mount to single-cell spatial assessment of gene expression in 3D

Unravelling spatio-temporal patterns of gene expression is crucial to understanding core biological principles from embryogenesis to disease. Here we review emerging technologies, providing automated, high-throughput, spatially resolved quantitative gene expression data. Novel techniques expand on current benchmark protocols, expediting their incorporation into ongoing research. These approaches digitally reconstruct patterns of embryonic expression in three dimensions, and have successfully identified novel domains of expression, cell types, and tissue features. Such technologies pave the way for unbiased and exhaustive recapitulation of gene expression levels in spatial and quantitative terms, promoting understanding of the molecular origin of developmental defects, and improving medical diagnostics.

Spatial modeling of prostate cancer metabolic gene expression reveals extensive heterogeneity and selective vulnerabilities

Spatial heterogeneity is a fundamental feature of the tumor microenvironment (TME), and tackling spatial heterogeneity in neoplastic metabolic aberrations is critical for tumor treatment. Genome-scale metabolic network models have been used successfully to simulate cancer metabolic networks. However, most models use bulk gene expression data of entire tumor biopsies, ignoring spatial heterogeneity in the TME. To account for spatial heterogeneity, we performed spatially-resolved metabolic network modeling of the prostate cancer microenvironment. We discovered novel malignant-cell-specific metabolic vulnerabilities targetable by small molecule compounds. We predicted that inhibiting the fatty acid desaturase SCD1 may selectively kill cancer cells based on our discovery of spatial separation of fatty acid synthesis and desaturation. We also uncovered higher prostaglandin metabolic gene expression in the tumor, relative to the surrounding tissue. Therefore, we predicted that inhibiting the prostaglandin transporter SLCO2A1 may selectively kill cancer cells. Importantly, SCD1 and SLCO2A1 have been previously shown to be potently and selectively inhibited by compounds such as CAY10566 and suramin, respectively. We also uncovered cancer-selective metabolic liabilities in central carbon, amino acid, and lipid metabolism. Our novel cancer-specific predictions provide new opportunities to develop selective drug targets for prostate cancer and other cancers where spatial transcriptomics datasets are available.

On Zebrafish Disease Models and Matters of the Heart

Coronary artery disease (CAD) is the leading form of cardiovascular disease (CVD), which is the primary cause of mortality worldwide. It is a complex disease with genetic and environmental risk factor contributions. Reports in human and mammalian models elucidate age-associated changes in cardiac function. The diverse mechanisms involved in cardiac diseases remain at the center of the research interest to identify novel strategies for prevention and therapy. Zebrafish (Danio rerio) have emerged as a valuable vertebrate model to study cardiovascular development over the last few decades. The facile genetic manipulation via forward and reverse genetic approaches combined with noninvasive, high-resolution imaging and phenotype-based screening has provided new insights to molecular pathways that orchestrate cardiac development. Zebrafish can recapitulate human cardiac pathophysiology due to gene and regulatory pathways conservation, similar heart rate and cardiac morphology and function. Thus, generations of zebrafish models utilize the functional analysis of genes involved in CAD, which are derived from large-scale human population analysis. Here, we highlight recent studies conducted on cardiovascular research focusing on the benefits of the combination of genome-wide association studies (GWAS) with functional genomic analysis in zebrafish. We further summarize the knowledge obtained from zebrafish studies that have demonstrated the architecture of the fundamental mechanisms underlying heart development, homeostasis and regeneration at the cellular and molecular levels.

Optogenetic sensors in the zebrafish heart: a novel in vivo electrophysiological tool to study cardiac arrhythmogenesis

Cardiac arrhythmias are among the most challenging human disorders to diagnose and treat due to their complex underlying pathophysiology. Suitable experimental animal models are needed to study the mechanisms causative for cardiac arrhythmogenesis. To enable in vivo analysis of cardiac cellular electrophysiology with a high spatial and temporal resolution, we generated and carefully validated two zebrafish models, one expressing an optogenetic voltage indicator (chimeric VSFP-butterfly CY) and the other a genetically encoded calcium indicator (GCaMP6f) in the heart.

Methods: High-speed epifluorescence microscopy was used to image chimeric VSFP-butterfly CY and GCaMP6f in the embryonic zebrafish heart, providing information about the spatiotemporal patterning of electrical activation, action potential configuration and intracellular Ca²⁺ dynamics. Plotting VSFP or GCaMP6f signals on a line along the myocardial wall over time facilitated the visualization and analysis of electrical impulse propagation throughout the heart. Administration of drugs targeting the sympathetic nervous system or cardiac ion channels was used to validate sensitivity and kinetics of both zebrafish sensor lines. Using the same microscope setup, we imaged transparent juvenile casper fish expressing GCaMP6f, demonstrating the feasibility of imaging cardiac optogenetic sensors at later stages of development.

Results: Isoproterenol slightly increased heart rate, diastolic Ca²⁺ levels and Ca²⁺ transient amplitudes, whereas propranolol caused a profound decrease in heart rate and Ca²⁺ transient parameters in VSFP-Butterfly and GCaMP6f embryonic fish. Ik_r blocker E-4031 decreased heart rate and increased action potential duration in VSFP-Butterfly fish. I_Ca,L blocker nifedipine caused total blockade of Ca²⁺ transients in GCaMP6f fish and a reduced heart rate, altered ventricular action potential duration and disrupted atrial-ventricular electrical conduction in VSFP-Butterfly fish. Imaging of juvenile animals demonstrated the possibility of employing an older zebrafish model for in vivo cardiac electrophysiology studies. We observed differences in atrial and ventricular Ca²⁺ recovery dynamics between 3 dpf and 14 dpf casper fish, but not in Ca²⁺ upstroke dynamics.

Conclusion: By introducing the optogenetic sensors chimeric VSFP-butterfly CY and GCaMP6f into the zebrafish we successfully generated an in vivo cellular electrophysiological readout tool for the zebrafish heart. Complementary use of both sensor lines demonstrated the ability to study heart rate, cardiac action potential configuration, spatiotemporal patterning of electrical activation and intracellular Ca²⁺ homeostasis in embryonic zebrafish. In addition, we demonstrated the first successful use of an optogenetic sensor to study cardiac function in older zebrafish. These models present a promising new research tool to study the underlying mechanisms of cardiac arrhythmogenesis.

Decoding the Heart through Next Generation Sequencing Approaches

: Vertebrate organs develop through a complex process which involves interaction between multiple signaling pathways at the molecular, cell, and tissue levels. Heart development is an example of such complex process which, when disrupted, results in congenital heart disease (CHD). This complexity necessitates a holistic approach which allows the visualization of genome-wide interaction networks, as opposed to assessment of limited subsets of factors. Genomics offers a powerful solution to address the problem of biological complexity by enabling the observation of molecular processes at a genome-wide scale. The emergence of next generation sequencing (NGS) technology has facilitated the expansion of genomics, increasing its output capacity and applicability in various biological disciplines. The application of NGS in various aspects of heart biology has resulted in new discoveries, generating novel insights into this field of study. Here we review the contributions of NGS technology into the understanding of heart development and its disruption reflected in CHD and discuss how emerging NGS based methodologies can contribute to the further understanding of heart repair.