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.

Lack of Stim2 Affects Vision-Dependent Behavior and Sensitivity to Hypoxia.

Stromal interaction molecules (STIMs) are endoplasmic reticulum-resident proteins that regulate Ca2+ homeostasis and signaling by store-operated calcium entry (SOCE). The different properties and functions of STIM1 and STIM2 have been described mostly based on work in vitro. STIM2 knockout mice do not survive until adulthood. Therefore, we generated and characterized stim2a and stim2b double-knockout zebrafish. The (stim2a;stim2b)-/- zebrafish did not have any apparent morphological phenotype. However, RNA sequencing revealed 1424 differentially expressed genes. One of the most upregulated genes was annexin A3a, which is a marker of activated microglia. This corresponded well to an increase in Neutral Red staining in the in vivo imaging of the (stim2a;stim2b)-/- zebrafish brain. The lack of Stim2 decreased zebrafish survival under low oxygen conditions. Behavioral tests, such as the visual-motor response test and dark-light preference test, indicated that (stim2a;stim2b)-/- larvae might have problems with vision. This was consistent with the downregulation of many genes that are related to light perception. The periodic acid-Schiff staining of retina sections from adult zebrafish revealed alterations of the stratum pigmentosum, suggesting the involvement of a Stim2-dependent process in visual perception. Altogether, these data reveal new functions for Stim2 in the nervous system.

xcore: an R package for inference of gene expression regulators.

BACKGROUND: Elucidating the Transcription Factors (TFs) that drive the gene expression changes in a given experiment is a common question asked by researchers. The existing methods rely on the predicted Transcription Factor Binding Site (TFBS) to model the changes in the motif activity. Such methods only work for TFs that have a motif and assume the TF binding profile is the same in all cell types. RESULTS: Given the wealth of the ChIP-seq data available for a wide range of the TFs in various cell types, we propose that gene expression modeling can be done using ChIP-seq "signatures" directly, effectively skipping the motif finding and TFBS prediction steps. We present xcore, an R package that allows TF activity modeling based on ChIP-seq signatures and the user's gene expression data. We also provide xcoredata a companion data package that provides a collection of preprocessed ChIP-seq signatures. We demonstrate that xcore leads to biologically relevant predictions using transforming growth factor beta induced epithelial-mesenchymal transition time-courses, rinderpest infection time-courses, and embryonic stem cells differentiated to cardiomyocytes time-course profiled with Cap Analysis Gene Expression. CONCLUSIONS: xcore provides a simple analytical framework for gene expression modeling using linear models that can be easily incorporated into differential expression analysis pipelines. Taking advantage of public ChIP-seq databases, xcore can identify meaningful molecular signatures and relevant ChIP-seq experiments.

Profiling subcellular localization of nuclear-encoded mitochondrial gene products in zebrafish.

Most mitochondrial proteins are encoded by nuclear genes, synthetized in the cytosol and targeted into the organelle. To characterize the spatial organization of mitochondrial gene products in zebrafish (Danio rerio), we sequenced RNA from different cellular fractions. Our results confirmed the presence of nuclear-encoded mRNAs in the mitochondrial fraction, which in unperturbed conditions, are mainly transcripts encoding large proteins with specific properties, like transmembrane domains. To further explore the principles of mitochondrial protein compartmentalization in zebrafish, we quantified the transcriptomic changes for each subcellular fraction triggered by the chchd4a -/- mutation, causing the disorders in the mitochondrial protein import. Our results indicate that the proteostatic stress further restricts the population of transcripts on the mitochondrial surface, allowing only the largest and the most evolutionary conserved proteins to be synthetized there. We also show that many nuclear-encoded mitochondrial transcripts translated by the cytosolic ribosomes stay resistant to the global translation shutdown. Thus, vertebrates, in contrast to yeast, are not likely to use localized translation to facilitate synthesis of mitochondrial proteins under proteostatic stress conditions.

A Zebrafish/Drosophila Dual System Model for Investigating Human Microcephaly.

Microcephaly presents in neurodevelopmental disorders with multiple aetiologies, including bi-allelic mutation in TUBGCP2, a component of the biologically fundamental and conserved microtubule-nucleation complex, γ-TuRC. Elucidating underlying principles driving microcephaly requires clear phenotype recapitulation and assay reproducibility, areas where go-to experimental models fall short. We present an alternative simple vertebrate/invertebrate dual system to investigate fundamental TUBGCP2-related processes driving human microcephaly and associated developmental traits. We show that antisense morpholino knockdown (KD) of the Danio rerio homolog, tubgcp2, recapitulates human TUBGCP2-associated microcephaly. Co-injection of wild type mRNA pre-empts microcephaly in 55% of KD zebrafish larvae, confirming causality. Body shortening observed in morphants is also rescued. Mitotic marker (pH3) staining further reveals aberrantly accumulated dividing brain cells in microcephalic tubgcp2 KD morphants, indicating that tubgcp2 depletion disrupts normal mitosis and/or proliferation in zebrafish neural progenitor brain cells. Drosophila melanogaster double knockouts (KO) for TUBGCP2 homologs Grip84/cg7716 also develop microcephalic brains with general microsomia. Exacerbated Grip84/cg7716-linked developmental aberration versus single mutations strongly suggests interactive or coinciding gene functions. We infer that tubgcp2 and Grip84/cg7716 affect brain size similarly to TUBGCP2 and recapitulate both microcephaly and microcephaly-associated developmental impact, validating the zebrafish/fly research model for human microcephaly. Given the conserved cross-phyla homolog function, the data also strongly support mitotic and/or proliferative disruption linked to aberrant microtubule nucleation in progenitor brain cells as key mechanistic defects for human microcephaly.

Adar-mediated A-to-I editing is required for embryonic patterning and innate immune response regulation in zebrafish.

Adenosine deaminases (ADARs) catalyze the deamination of adenosine to inosine, also known as A-to-I editing, in RNA. Although A-to-I editing occurs widely across animals and is well studied, new biological roles are still being discovered. Here, we study the role of A-to-I editing in early zebrafish development. We demonstrate that Adar, the zebrafish orthologue of mammalian ADAR1, is essential for establishing the antero-posterior and dorso-ventral axes and patterning. Genome-wide editing discovery reveals pervasive editing in maternal and the earliest zygotic transcripts, the majority of which occurred in the 3'-UTR. Interestingly, transcripts implicated in gastrulation as well as dorso-ventral and antero-posterior patterning are found to contain multiple editing sites. Adar knockdown or overexpression affect gene expression by 12 hpf. Analysis of adar-/- zygotic mutants further reveals that the previously described role of Adar in mammals in regulating the innate immune response is conserved in zebrafish. Our study therefore establishes distinct maternal and zygotic functions of RNA editing by Adar in embryonic patterning along the zebrafish antero-posterior and dorso-ventral axes, and in the regulation of the innate immune response, respectively.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements.

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Multi-omics analyses of early liver injury reveals cell-type-specific transcriptional and epigenomic shift.

BACKGROUND: Liver fibrosis is a wound-healing response to tissue injury and inflammation hallmarked by the extracellular matrix (ECM) protein deposition in the liver parenchyma and tissue remodelling. Different cell types of the liver are known to play distinct roles in liver injury response. Hepatocytes and liver endothelial cells receive molecular signals indicating tissue injury and activate hepatic stellate cells which produce ECM proteins upon their activation. Despite the growing knowledge on the molecular mechanism underlying hepatic fibrosis in general, the cell-type-specific gene regulatory network associated with the initial response to hepatotoxic injury is still poorly characterized. RESULTS: In this study, we used thioacetamide (TAA) to induce hepatic injury in adult zebrafish. We isolated three major liver cell types - hepatocytes, endothelial cells and hepatic stellate cells - and identified cell-type-specific chromatin accessibility and transcriptional changes in an early stage of liver injury. We found that TAA induced transcriptional shifts in all three cell types hallmarked by significant alterations in the expression of genes related to fatty acid and carbohydrate metabolism, as well as immune response-associated and vascular-specific genes. Interestingly, liver endothelial cells exhibit the most pronounced response to liver injury at the transcriptome and chromatin level, hallmarked by the loss of their angiogenic phenotype. CONCLUSION: Our results uncovered cell-type-specific transcriptome and epigenome responses to early stage liver injury, which provide valuable insights into understanding the molecular mechanism implicated in the early response of the liver to pro-fibrotic signals.

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.

Cardiac-specific ß-catenin deletion dysregulates energetic metabolism and mitochondrial function in perinatal cardiomyocytes.

ß-Catenin signaling pathway regulates cardiomyocytes proliferation and differentiation, though its involvement in metabolic regulation of cardiomyocytes remains unknown. We used one-day-old mice with cardiac-specific knockout of ß-catenin and neonatal rat ventricular myocytes treated with ß-catenin inhibitor to investigate the role of ß-catenin metabolism regulation in perinatal cardiomyocytes. Transcriptomics of perinatal ß-catenin-ablated hearts revealed a dramatic shift in the expression of genes involved in metabolic processes. Further analysis indicated an inhibition of lipolysis and glycolysis in both in vitro and in vivo models. Finally, we showed that ß-catenin deficiency leads to mitochondria dysfunction via the downregulation of Sirt1/PGC-1α pathway. We conclude that cardiac-specific ß-catenin ablation disrupts the energy substrate shift that is essential for postnatal heart maturation, leading to perinatal lethality of homozygous ß-catenin knockout mice.

The development of zebrafish pancreas affected by deficiency of Hedgehog signaling.

The pancreas development depends on complex regulation of several signaling pathways, including the Hedgehog (Hh) signaling via a receptor complex component, Smoothened, which deficiency blocks the Hh signaling. Such a defect in birds and mammals results in an annular pancreas. We showed that in developing zebrafish, the mutation of Smoothened or inhibition of Hh signaling by its antagonist cyclopamine caused developmental defects of internal organs, liver, pancreas, and gut. In particular, the pancreatic primordium was duplicated. The two exocrine pancreatic primordia surround the gut. This phenomenon correlates with a significant reduction of the gut's diameter, causing the annular pancreas phenotype.

Fish-Ing for Enhancers in the Heart.

Precise control of gene expression is crucial to ensure proper development and biological functioning of an organism. Enhancers are non-coding DNA elements which play an essential role in regulating gene expression. They contain specific sequence motifs serving as binding sites for transcription factors which interact with the basal transcription machinery at their target genes. Heart development is regulated by intricate gene regulatory network ensuring precise spatiotemporal gene expression program. Mutations affecting enhancers have been shown to result in devastating forms of congenital heart defect. Therefore, identifying enhancers implicated in heart biology and understanding their mechanism is key to improve diagnosis and therapeutic options. Despite their crucial role, enhancers are poorly studied, mainly due to a lack of reliable way to identify them and determine their function. Nevertheless, recent technological advances have allowed rapid progress in enhancer discovery. Model organisms such as the zebrafish have contributed significant insights into the genetics of heart development through enabling functional analyses of genes and their regulatory elements in vivo. Here, we summarize the current state of knowledge on heart enhancers gained through studies in model organisms, discuss various approaches to discover and study their function, and finally suggest methods that could further advance research in this field.

Exploring Translational Control of Maternal mRNAs in Zebrafish.

The study of translational regulation requires reliable measurement of both mRNA levels and protein synthesis. Cytoplasmic polyadenylation is a prevalent mode of translational regulation during oogenesis and early embryogenesis. Here the length of the poly(A) tail of an mRNA is coupled to its translatability. We describe a protocol to identify translationally regulated genes and measure their translation rate in the early zebrafish embryo using genome-wide polysome profiling. This protocol relies on the isolation of mRNA by means of an rRNA depletion strategy, which avoids capture bias due to short poly(A) tail that can occur when using conventional oligo(dT)-based methods. We also present a simple PCR-based method to measure the poly(A) tail length of selected mRNAs.

Targeted RNA Knockdown by a Type III CRISPR-Cas Complex in Zebrafish.

RNA interference is a powerful experimental tool for RNA knockdown, but not all organisms are amenable. Here, we provide a proof of principle demonstration that a type III Csm effector complex can be used for programmable mRNA transcript degradation in eukaryotes. In zebrafish, Streptococcus thermophilus Csm complex (StCsm) proved effective for knockdown of maternally expressed EGFP in germ cells of Tg(ddx4:ddx4-EGFP) fish. It also led to significant, albeit less drastic, fluorescence reduction at one day postfertilization in Tg(myl7:GFP) and Tg(fli1:EGFP) fish that express EGFP zygotically. StCsm targeted against the endogenous tdgf1 elicited the characteristic one-eyed phenotype with greater than 50% penetrance, and hence with similar efficiency to morpholino-mediated knockdown. We conclude that Csm-mediated knockdown is very efficient for maternal transcripts and can also be used for mixed maternal/early zygotic and early zygotic transcripts, in some cases reaching comparable efficiency to morpholino-based knockdown without significant off-target effects.

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.

Dynamics of cardiomyocyte transcriptome and chromatin landscape demarcates key events of heart development.

Organogenesis involves dynamic regulation of gene transcription and complex multipathway interactions. Despite our knowledge of key factors regulating various steps of heart morphogenesis, considerable challenges in understanding its mechanism still exist because little is known about their downstream targets and interactive regulatory network. To better understand transcriptional regulatory mechanism driving heart development and the consequences of its disruption in vivo, we performed time-series analyses of the transcriptome and genome-wide chromatin accessibility in isolated cardiomyocytes (CMs) from wild-type zebrafish embryos at developmental stages corresponding to heart tube morphogenesis, looping, and maturation. We identified genetic regulatory modules driving crucial events of heart development that contained key cardiac TFs and are associated with open chromatin regions enriched for DNA sequence motifs belonging to the family of the corresponding TFs. Loss of function of cardiac TFs Gata5, Tbx5a, and Hand2 affected the cardiac regulatory networks and caused global changes in chromatin accessibility profile, indicating their role in heart development. Among regions with differential chromatin accessibility in mutants were highly conserved noncoding elements that represent putative enhancers driving heart development. The most prominent gene expression changes, which correlated with chromatin accessibility modifications within their proximal promoter regions, occurred between heart tube morphogenesis and looping, and were associated with metabolic shift and hematopoietic/cardiac fate switch during CM maturation. Our results revealed the dynamic regulatory landscape throughout heart development and identified interactive molecular networks driving key events of heart morphogenesis.

The translational regulation of maternal mRNAs in time and space.

Since their discovery, the study of maternal mRNAs has led to the identification of mechanisms underlying their spatiotemporal regulation within the context of oogenesis and early embryogenesis. Following synthesis in the oocyte, maternal mRNAs are translationally silenced and sequestered into storage in cytoplasmic granules. At the same time, their unique distribution patterns throughout the oocyte and embryo are tightly controlled and connected to their functions in downstream embryonic processes. At certain points in oogenesis and early embryogenesis, maternal mRNAs are translationally activated to perform their functions in a timely manner. The cytoplasmic polyadenylation machinery is responsible for the translational activation of maternal mRNAs, and its role in initiating the maternal to zygotic transition events has recently come to light. Here, we summarize the current knowledge on maternal mRNA regulation, with particular focus on cytoplasmic polyadenylation as a mechanism for translational regulation.