Development of the cardiac conduction system in zebrafish

Fibroblast Growth Factor (FGF) 13

Fibroblast Growth Factor (FGF) 13, also referred to as FGF homologous factor (FHF) 2, is a member of the FGF11 subfamily that is characterized as having sequence similarities to classical FGF receptor (FGFR)-binding FGFs, but functionally do not bind FGFRs. In this primer mini-review, we summarize current knowledge regarding FGF13 expression, mutant analyses, and gene and protein structure. Similar to other FHFs, FGF13 has been considered a non-secreted protein that lacks an amino signal and is prominently expressed in developing and mature neurons of the central and peripheral nervous systems, as well as the heart. The expression of FGF13 is not limited to early embryonic stages and has been shown to persist in adult tissues. As well, FGF13 is known to localize subcellularly, both within the cytoplasm and the nucleus. FGF13 is extremely adaptable, as it interacts with MAPK scaffolding protein islet brain 2 (IB2), stabilizes microtubules, or binds to voltage-gated sodium channels. Fgf13 mutant mouse lines display various neurological pathologies. Through sequence mapping, FGF13 is considered a candidate causative gene that is mutated in multiple human X-linked neurological diseases.

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.

POPDC1 Variants Cause Atrioventricular Node Dysfunction and Arrhythmogenic Changes in Cardiac Electrophysiology and Intracellular Calcium Handling in Zebrafish

Popeye domain-containing (POPDC) proteins selectively bind cAMP and mediate cellular responses to sympathetic nervous system (SNS) stimulation. The first discovered human genetic variant (POPDC1S201F) is associated with atrioventricular (AV) block, which is exacerbated by increased SNS activity. Zebrafish carrying the homologous mutation (popdc1S191F) display a similar phenotype to humans. To investigate the impact of POPDC1 dysfunction on cardiac electrophysiology and intracellular calcium handling, homozygous popdc1S191F and popdc1 knock-out (popdc1KO) zebrafish larvae and adult isolated popdc1S191F hearts were studied by functional fluorescent analysis. It was found that in popdc1S191F and popdc1KO larvae, heart rate (HR), AV delay, action potential (AP) and calcium transient (CaT) upstroke speed, and AP duration were less than in wild-type larvae, whereas CaT duration was greater. SNS stress by β-adrenergic receptor stimulation with isoproterenol increased HR, lengthened AV delay, slowed AP and CaT upstroke speed, and shortened AP and CaT duration, yet did not result in arrhythmias. In adult popdc1S191F zebrafish hearts, there was a higher incidence of AV block, slower AP upstroke speed, and longer AP duration compared to wild-type hearts, with no differences in CaT. SNS stress increased AV delay and led to further AV block in popdc1S191F hearts while decreasing AP and CaT duration. Overall, we have revealed that arrhythmogenic effects of POPDC1 dysfunction on cardiac electrophysiology and intracellular calcium handling in zebrafish are varied, but already present in early development, and that AV node dysfunction may underlie SNS-induced arrhythmogenesis associated with popdc1 mutation in adults.

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.

A phenotype-based forward genetic screen identifies Dnajb6 as a sick sinus syndrome gene

Previously we showed the generation of a protein trap library made with the gene-break transposon (GBT) in zebrafish (Danio rerio) that could be used to facilitate novel functional genome annotation towards understanding molecular underpinnings of human diseases (Ichino et al, 2020). Here, we report a significant application of this library for discovering essential genes for heart rhythm disorders such as sick sinus syndrome (SSS). SSS is a group of heart rhythm disorders caused by malfunction of the sinus node, the heart's primary pacemaker. Partially owing to its aging-associated phenotypic manifestation and low expressivity, molecular mechanisms of SSS remain difficult to decipher. From 609 GBT lines screened, we generated a collection of 35 zebrafish insertional cardiac (ZIC) mutants in which each mutant traps a gene with cardiac expression. We further employed electrocardiographic measurements to screen these 35 ZIC lines and identified three GBT mutants with SSS-like phenotypes. More detailed functional studies on one of the arrhythmogenic mutants, GBT411, in both zebrafish and mouse models unveiled Dnajb6 as a novel SSS causative gene with a unique expression pattern within the subpopulation of sinus node pacemaker cells that partially overlaps with the expression of hyperpolarization activated cyclic nucleotide gated channel 4 (HCN4), supporting heterogeneity of the cardiac pacemaker cells.

Genetics of Atavism

Abstract

Atavisms have attracted people’s attention for a long time. First, atavisms excited their imagination and created fertile ground for myths and superstitions. With the development of science, atavisms became the subject of investigation, which soon provided evidence to support evolutionary theory. However, at the molecular level, the formation of atavisms remained insufficiently understood. Recent progress in comparative genomics and molecular developmental biology has helped in understanding the processes underlying the formation of one of the human atavisms: the vestigial tail.

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.

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.

Ablation of Mto1 in zebrafish exhibited hypertrophic cardiomyopathy manifested by mitochondrion RNA maturation deficiency

Deficient maturations of mitochondrial transcripts are linked to clinical abnormalities but their pathophysiology remains elusive. Previous investigations showed that pathogenic variants in MTO1 for the biosynthesis of τm5U of tRNAGlu, tRNAGln, tRNALys, tRNATrp and tRNALeu(UUR) were associated with hypertrophic cardiomyopathy (HCM). Using mto1 knock-out(KO) zebrafish generated by CRISPR/Cas9 system, we demonstrated the pleiotropic effects of Mto1 deficiency on mitochondrial RNA maturations. The perturbed structure and stability of tRNAs caused by mto1 deletion were evidenced by conformation changes and sensitivity to S1-mediated digestion of tRNAGln, tRNALys, tRNATrp and tRNALeu(UUR). Notably, mto1KO zebrafish exhibited the global decreases in the aminoacylation of mitochondrial tRNAs with the taurine modification. Strikingly, ablated mto1 mediated the expression of MTPAP and caused the altered polyadenylation of cox1, cox3, and nd1 mRNAs. Immunoprecipitation assay indicated the interaction of MTO1 with MTPAP related to mRNA polyadenylation. These alterations impaired mitochondrial translation and reduced activities of oxidative phosphorylation complexes. These mitochondria dysfunctions caused heart development defects and hypertrophy of cardiomyocytes and myocardial fiber disarray in ventricles. These cardiac defects in the mto1KO zebrafish recapitulated the clinical phenotypes in HCM patients carrying the MTO1 mutation(s). Our findings highlighted the critical role of MTO1 in mitochondrial transcript maturation and their pathological consequences in hypertrophic cardiomyopathy.

Utility of Zebrafish Models of Acquired and Inherited Long QT Syndrome

Long-QT Syndrome (LQTS) is a cardiac electrical disorder, distinguished by irregular heart rates and sudden death. Accounting for ∼40% of cases, LQTS Type 2 (LQTS2), is caused by defects in the Kv11.1 (hERG) potassium channel that is critical for cardiac repolarization. Drug block of hERG channels or dysfunctional channel variants can result in acquired or inherited LQTS2, respectively, which are typified by delayed repolarization and predisposition to lethal arrhythmia. As such, there is significant interest in clear identification of drugs and channel variants that produce clinically meaningful perturbation of hERG channel function. While toxicological screening of hERG channels, and phenotypic assessment of inherited channel variants in heterologous systems is now commonplace, affordable, efficient, and insightful whole organ models for acquired and inherited LQTS2 are lacking. Recent work has shown that zebrafish provide a viable in vivo or whole organ model of cardiac electrophysiology. Characterization of cardiac ion currents and toxicological screening work in intact embryos, as well as adult whole hearts, has demonstrated the utility of the zebrafish model to contribute to the development of therapeutics that lack hERG-blocking off-target effects. Moreover, forward and reverse genetic approaches show zebrafish as a tractable model in which LQTS2 can be studied. With the development of new tools and technologies, zebrafish lines carrying precise channel variants associated with LQTS2 have recently begun to be generated and explored. In this review, we discuss the present knowledge and questions raised related to the use of zebrafish as models of acquired and inherited LQTS2. We focus discussion, in particular, on developments in precise gene-editing approaches in zebrafish to create whole heart inherited LQTS2 models and evidence that zebrafish hearts can be used to study arrhythmogenicity and to identify potential anti-arrhythmic compounds.

The Intrinsic Cardiac Nervous System and Its Role in Cardiac Pacemaking and Conduction

The cardiac autonomic nervous system (CANS) plays a key role for the regulation of cardiac activity with its dysregulation being involved in various heart diseases, such as cardiac arrhythmias. The CANS comprises the extrinsic and intrinsic innervation of the heart. The intrinsic cardiac nervous system (ICNS) includes the network of the intracardiac ganglia and interconnecting neurons. The cardiac ganglia contribute to the tight modulation of cardiac electrophysiology, working as a local hub integrating the inputs of the extrinsic innervation and the ICNS. A better understanding of the role of the ICNS for the modulation of the cardiac conduction system will be crucial for targeted therapies of various arrhythmias. We describe the embryonic development, anatomy, and physiology of the ICNS. By correlating the topography of the intracardiac neurons with what is known regarding their biophysical and neurochemical properties, we outline their physiological role in the control of pacemaker activity of the sinoatrial and atrioventricular nodes. We conclude by highlighting cardiac disorders with a putative involvement of the ICNS and outline open questions that need to be addressed in order to better understand the physiology and pathophysiology of the ICNS.

Deletion of Gtpbp3 in zebrafish revealed the hypertrophic cardiomyopathy manifested by aberrant mitochondrial tRNA metabolism

GTPBP3 is a highly conserved tRNA modifying enzyme for the biosynthesis of τm⁵U at the wobble position of mitochondrial tRNA^Glu, tRNA^Gln, tRNA^Lys, tRNA^Trp and tRNA^Leu(UUR). The previous investigations showed that GTPBP3 mutations were associated with hypertrophic cardiomyopathy (HCM). However, the pathophysiology of GTPBP3 deficiency remains elusively. Using the gtpbp3 knockout zebrafish generated by CRISPR/Cas9 system, we demonstrated the aberrant mitochondrial tRNA metabolism in gtpbp3 knock-out zebrafish. The deletion of gtpbp3 may alter functional folding of tRNA, indicated by conformation changes and sensitivity to S1-mediated digestion of tRNA^Glu, tRNA^Lys, tRNA^Trp and tRNA^Leu(UUR). Strikingly, gtpbp3 knock-out zebrafish displayed the global increases in the aminoacylated efficiencies of mitochondrial tRNAs. The aberrant mitochondrial tRNA metabolisms impaired mitochondrial translation, produced proteostasis stress and altered activities of respiratory chain complexes. These mitochondria dysfunctions caused the alterations in the embryonic heart development and reduced fractional shortening of ventricles in mutant zebrafish. Notably, the gtpbp3 knock-out zebrafish exhibited hypertrophy of cardiomyocytes and myocardial fiber disarray in ventricles. These cardiac defects in the gtpbp3 knock-out zebrafish recapitulated the clinical phenotypes in HCM patients carrying the GTPBP3 mutation(s). Our findings highlight the fundamental role of defective nucleotide modifications of tRNAs in mitochondrial biogenesis and their pathological consequences in hypertrophic cardiomyopathy.

Toxic effects of flufenoxuron on development and vascular formation during zebrafish embryogenesis

Flufenoxuron, a chitin synthesis inhibitor that is widely used in developed countries as an insecticide, is rarely degraded in the environment. In addition to that in insects, flufenoxuron-mediated non-targeted death in organisms such as lizards and bees has been reported. However, the toxic effects of this compound on vascular development during embryogenesis, as well as the underlying mechanism, have not yet been elucidated. In the present study, we assessed abnormal development and cardiovascular damage induced by flufenoxuron in zebrafish embryos. Exposed zebrafish had malformed eyes and pathological characteristics such as heart and yolk sac edema. In accordance with developmental inhibition, cell cycle regulatory genes were dysregulated in zebrafish embryos upon exposure to flufenoxuron. We also discovered that this agent can disrupt vascular formation by interfering with angiogenesis-associated genes including the genes encoding vascular endothelial growth factor Aa (vegfaa), vegfc, fms-related tyrosine kinase 1 (flt1), and flt4 in zebrafish embryos. These anti-angiogenic effects of flufenoxuron were further verified using a well-known angiogenesis model, namely human umbilical vein endothelial cells. In conclusion, our results suggest that flufenoxuron inhibits overall development and angiogenesis during embryogenesis.

Muscarinic receptors promote pacemaker fate at the expense of secondary conduction system tissue in zebrafish

Deterioration or inborn malformations of the cardiac conduction system (CCS) interfere with proper impulse propagation in the heart and may lead to sudden cardiac death or heart failure. Patients afflicted with arrhythmia depend on antiarrhythmic medication or invasive therapy, such as pacemaker implantation. An ideal way to treat these patients would be CCS tissue restoration. This, however, requires precise knowledge regarding the molecular mechanisms underlying CCS development. Here, we aimed to identify regulators of CCS development. We performed a compound screen in zebrafish embryos and identified tolterodine, a muscarinic receptor antagonist, as a modifier of CCS development. Tolterodine provoked a lower heart rate, pericardiac edema, and arrhythmia. Blockade of muscarinic M3, but not M2, receptors induced transcriptional changes leading to amplification of sinoatrial cells and loss of atrioventricular identity. Transcriptome data from an engineered human heart muscle model provided additional evidence for the contribution of muscarinic M3 receptors during cardiac progenitor specification and differentiation. Taken together, we found that muscarinic M3 receptors control the CCS already before the heart becomes innervated. Our data indicate that muscarinic receptors maintain a delicate balance between the developing sinoatrial node and the atrioventricular canal, which is probably required to prevent the development of arrhythmia.

A novel conserved enhancer at zebrafish zic3 and zic6 loci drives neural expression

Background

Identifying enhancers and deciphering their putative roles represent a major step to better understand the mechanism of metazoan gene regulation, development, and the role of regulatory elements in disease. Comparative genomics and transgenic assays have been used with some success to identify critical regions that are involved in regulating the spatiotemporal expression of genes during embryogenesis.

Results

We identified two novel tetrapod‐teleost conserved noncoding elements within the vicinity of the zic3 and zic6 loci in the zebrafish genome and demonstrated their ability to drive tissue‐specific expression in a transgenic zebrafish assay. The syntenic analysis and robust green fluorescent expression in the developing habenula in the stable transgenic line were correlated with known sites of endogenous zic3 and zic6 expression.

Conclusion

This transgenic line that expresses green fluorescent protein in the habenula is a valuable resource for studying a specific population of cells in the zebrafish central nervous system. Our observations indicate that a genomic sequence that is conserved between humans and zebrafish acts as an enhancer that likely controls zic3 and zic6 expression.

Identified a novel enhancer near zebrafish zic3/zic6 locus.

The novel enhancer drives tissue‐specific expression in the habenula.

Zebrafish transgenic line generated in this study can be a useful resource for studying development of habenula.

Intrinsic regulation of sinoatrial node function and the zebrafish as a model of stretch effects on pacemaking

Excitation of the heart occurs in a specialised region known as the sinoatrial node (SAN). Tight regulation of SAN function is essential for the maintenance of normal heart rhythm and the response to (patho-)physiological changes. The SAN is regulated by extrinsic (central nervous system) and intrinsic (neurons, peptides, mechanics) factors. The positive chronotropic response to stretch in particular is essential for beat-by-beat adaptation to changes in hemodynamic load. Yet, the mechanism of this stretch response is unknown, due in part to the lack of an appropriate experimental model for targeted investigations. We have been investigating the zebrafish as a model for the study of intrinsic regulation of SAN function. In this paper, we first briefly review current knowledge of the principal components of extrinsic and intrinsic SAN regulation, derived primarily from experiments in mammals, followed by a description of the zebrafish as a novel experimental model for studies of intrinsic SAN regulation. This mini-review is followed by an original investigation of the response of the zebrafish isolated SAN to controlled stretch. Stretch causes an immediate and continuous increase in beating rate in the zebrafish isolated SAN. This increase reaches a maximum part way through a period of sustained stretch, with the total change dependent on the magnitude and direction of stretch. This is comparable to what occurs in isolated SAN from most mammals (including human), suggesting that the zebrafish is a novel experimental model for the study of mechanisms involved in the intrinsic regulation of SAN function by mechanical effects.