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.

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.

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.

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.

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.

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.

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.

Zebrafish Zic Genes Mediate Developmental Signaling.

The introduction of genomics into the field of developmental biology led to a vast expansion of knowledge about developmental genes and signaling mechanisms they are involved in. Unlike mammals, the zebrafish features seven Zic genes. This provides an interesting insight into Zic gene evolution. In addition, an unprecedented bioimaging capability of semitransparent zebrafish embryos turns to be a crucial factor in medium- to large-scale analysis of the activity of potential regulatory elements. The Zic family of zinc finger proteins plays an important, relatively well-established, role in the regulation of stem cells and neural development and, in particular, during neural fate commitment and determination. At the same time, some Zic genes are expressed in mesodermal lineages, and their deficiency causes a number of developmental defects in axis formation, establishing body symmetry and cardiac morphogenesis. In stem cells, Zic genes are required to maintain pluripotency by binding to the proximal promoters of pluripotency genes (Oct4, Nanog, Sox2, etc.). During embryogenesis, the dynamic nature of Zic transcriptional regulation is manifested by the interaction of these factors with distal enhancers and other regulatory elements associated with the control of gene transcription and, in particular, with the Nodal and Wnt signaling pathways that play a role in establishing basic organization of the vertebrate body. Zic transcription factors may regulate development through acting alone as well as in combination with other transcription factors. This is achieved due to Zic binding to sites adjacent to the binding sites of other transcription factors, including Gli. This probably leads to the formation of multi-transcription factor complexes associated with enhancers.

Cytoplasmic polyadenylation-mediated translational control of maternal mRNAs directs maternal-to-zygotic transition.

In the earliest stages of animal development following fertilization, maternally deposited mRNAs direct biological processes to the point of zygotic genome activation (ZGA). These maternal mRNAs undergo cytoplasmic polyadenylation (CPA), suggesting translational control of their activation. To elucidate the biological role of CPA during embryogenesis, we performed genome-wide polysome profiling at several stages of zebrafish development. Our analysis revealed a correlation between CPA and polysome-association dynamics, demonstrating a coupling of translation to the CPA of maternal mRNAs. Pan-embryonic CPA inhibition disrupted the maternal-to-zygotic transition (MZT), causing a failure of developmental progression beyond the mid-blastula transition and changes in global gene expression that indicated a failure of ZGA and maternal mRNA clearance. Among the genes that were differentially expressed were those encoding chromatin modifiers and key transcription factors involved in ZGA, including nanog, pou5f3 and sox19b, which have distinct CPA dynamics. Our results establish the necessity of CPA for ensuring progression of the MZT. The RNA-seq data generated in this study represent a valuable zebrafish resource for the discovery of novel elements of the early embryonic transcriptome.

Changing Faces of Transcriptional Regulation Reflected by Zic3.

The advent of genomics in the study of developmental mechanisms has brought a trove of information on gene datasets and regulation during development, where the Zic family of zinc-finger proteins plays an important role. Genomic analysis of the modes of action of Zic3 in pluripotent cells demonstrated its requirement for maintenance of stem cells pluripotency upon binding to the proximal regulatory regions (promoters) of genes associated with cell pluripotency (Nanog, Sox2, Oct4, etc.) as well as cell cycle, proliferation, oncogenesis and early embryogenesis. In contrast, during gastrulation and neurulation Zic3 acts by binding the distal regulatory regions (enhancers, etc) associated with control of gene transcription in the Nodal and Wnt signaling pathways, including genes that act to break body symmetry. This illustrates a general role of Zic3 as a transcriptional regulator that acts not only alone, but in many instances in conjunction with other transcription factors. The latter is done by binding to adjacent sites in the context of multi-transcription factor complexes associated with regulatory elements.

Impaired development of neural-crest cell-derived organs and intellectual disability caused by MED13L haploinsufficiency.

MED13L is a component subunit of the Mediator complex, an important regulator of transcription that is highly conserved across eukaryotes. Here, we report MED13L disruption in a translocation t(12;19) breakpoint of a patient with Pierre-Robin syndrome, moderate intellectual disability, craniofacial anomalies, and muscular defects. The phenotype is similar to previously described patients with MED13L haploinsufficiency. Knockdown of MED13L orthologue in zebrafish, med13b, showed early defective migration of cranial neural crest cells (NCCs) that contributed to cartilage structure deformities in the later stage, recapitulating craniofacial anomalies seen in human patients. Notably, we observed abnormal distribution of developing neurons in different brain regions of med13b morphant embryos, which could be rescued upon introduction of full-length human MED13L mRNA. To compare with mammalian system, we suppressed MED13L expression by short-hairpin RNA in ES-derived human neural progenitors, and differentiated them into neurons. Transcriptome analysis revealed differential expression of components of Wnt and FGF signaling pathways in MED13L-deficient neurons. Our finding provides a novel insight into the mechanism of overlapping phenotypic outcome targeting NCCs derivatives organs in patients with MED13L haploinsufficiency, and emphasizes a clinically recognizable syndromic phenotype in these patients.

Expression of components of Wnt and Hedgehog pathways in different tissue layers during lung development in Xenopus laevis.

Although Wnt and Hedgehog (Hh) signaling pathways play important roles in mouse lung development, these have not been explored in the development of Xenopus lung. This may be due to the lack of specific molecular markers for different layers of tissue in Xenopus lung and/or insufficient knowledge on expression patterns of Wnt and Hh signaling components in Xenopus lung. In this study, we first described the early morphogenesis of Xenopus laevis lung by using surfactant protein C (sftpc) as a marker of lung epithelium and compared it with the expression patterns of several genes of Wnt and Hh pathways in Xenopus lungs. Our data showed that wnt7b was expressed in the entire lung epithelium from stage 37 to stage 45, while two other Wnt signaling components, wnt5a and wif1 (wnt inhibitory factor 1), were expressed in the mesenchyme layer of the entire lungs through stages 39-41. We also found that sonic hedgehog (shh) was expressed at stage 41 only in the anterior, but not in the posterior part of the lungs. These results show the expression of wnt5a, wnt7b, wif1 and shh in different layers of tissue of Xenopus lungs at early developmental stages, which implies different roles of these genes in the early development of Xenopus lungs. Our study for the first time defined specific molecular markers for description of early lung development in Xenopus, as well as provided information about expression of components of Wnt and Hh pathways in early Xenopus lungs, which should be useful for future functional studies.

Mercury-induced hepatotoxicity in zebrafish: in vivo mechanistic insights from transcriptome analysis, phenotype anchoring and targeted gene expression validation.

BACKGROUND: Mercury is a prominent environmental contaminant that causes detrimental effects to human health. Although the liver has been known to be a main target organ, there is limited information on in vivo molecular mechanism of mercury-induced toxicity in the liver. By using transcriptome analysis, phenotypic anchoring and validation of targeted gene expression in zebrafish, mercury-induced hepatotoxicity was investigated and a number of perturbed cellular processes were identified and compared with those captured in the in vitro human cell line studies. RESULTS: Hepato-transcriptome analysis of mercury-exposed zebrafish revealed that the earliest deregulated genes were associated with electron transport chain, mitochondrial fatty acid beta-oxidation, nuclear receptor signaling and apoptotic pathway, followed by complement system and proteasome pathway, and thereafter DNA damage, hypoxia, Wnt signaling, fatty acid synthesis, gluconeogenesis, cell cycle and motility. Comparative meta-analysis of microarray data between zebrafish liver and human HepG2 cells exposed to mercury identified some common toxicological effects of mercury-induced hepatotoxicity in both models. Histological analyses of liver from mercury-exposed fish revealed morphological changes of liver parenchyma, decreased nucleated cell count, increased lipid vesicles, glycogen and apoptotic bodies, thus providing phenotypic evidence for anchoring of the transcriptome analysis. Validation of targeted gene expression confirmed deregulated gene-pathways from enrichment analysis. Some of these genes responding to low concentrations of mercury may serve as toxicogenomic-based markers for detection and health risk assessment of environmental mercury contaminations. CONCLUSION: Mercury-induced hepatotoxicity was triggered by oxidative stresses, intrinsic apoptotic pathway, deregulation of nuclear receptor and kinase activities including Gsk3 that deregulates Wnt signaling pathway, gluconeogenesis, and adipogenesis, leading to mitochondrial dysfunction, endocrine disruption and metabolic disorders. This study provides important mechanistic insights into mercury-induced liver toxicity in a whole-animal physiology context, which will help in understanding the syndromes caused by mercury poisoning. The molecular conservation of mercury-induced hepatotoxicity between zebrafish and human cell line reveals the feasibility of using zebrafish to model molecular toxicity in human for toxicant risk assessments.

The role of vasculature and blood circulation in zebrafish swimbladder development.

BACKGROUND: Recently we have performed a detailed analysis of early development of zebrafish swimbladder, a homologous organ of tetrapod lung; however, the events of swimbladder development are still poorly characterized. Many studies have implicated the role of vascular system in development of many organs in vertebrates. As the swimbladder is lined with an intricate network of blood capillaries, it is of interest to investigate the role of the vascular system during early development of swimbladder. RESULTS: To investigate the role of endothelial cells (ECs) and blood circulation during development of the swimbladder, phenotypes of swimbladder were analysed at three different stages (approximately 2, 3 and 5 dpf [day postfertilization]) in cloche (clo) mutant and Tnnt2 morphants, in the background of transgenic lines Et(krt4:EGFP)sq33-2 and Et(krt4:EGFP)sqet3 which express EGFP in the swimbladder epithelium and outer mesothelium respectively. Analyses of the three tissue layers of the swimbladder were performed using molecular markers hb9, fgf10a, acta2, and anxa5 to distinguish epithelium, mesenchyme, and outer mesothelium. We showed that the budding stage was independent of ECs and blood flow, while early epithelial growth, mesenchymal organization and its differentiation into smooth muscle, as well as outer mesothelial organization, were dependent on ECs. Blood circulation contributed to later stage of epithelial growth, smooth muscle differentiation, and organization of the outer mesothelium. Inflation of the swimbladder was also affected as a result of absence of ECs and blood flow. CONCLUSION: Our data demonstrated that the vascular system, though not essential in swimbladder budding, plays an important role in the development of the swimbladder starting from the early growth stage, including mesenchyme organization and smooth muscle differentiation, and outer mesothelial organization, which in turn may be essential for the function of the swimbladder as reflected in its eventual inflation.

Development of zebrafish swimbladder: The requirement of Hedgehog signaling in specification and organization of the three tissue layers.

The swimbladder is a hydrostatic organ in fish postulated as a homolog of the tetrapod lung. While lung development has been well studied, the molecular mechanism of swimbladder development is essentially uncharacterized. In the present study, swimbladder development in zebrafish was analyzed by using several molecular markers: hb9 (epithelium), fgf10a and acta2 (mesenchyme), and anxa5 (mesothelium), as well as in vivo through enhancer trap transgenic lines Et(krt4:EGFP)(sq33-2) and Et(krt4:EGFP)(sqet3) that showed strong EGFP expression in the swimbladder epithelium and outer mesothelium respectively. We defined three phases of swimbladder development: epithelial budding between 36 and 48 hpf, growth with the formation of two additional mesodermal layers up to 4.5 dpf, and inflation of posterior and anterior chambers at 4.5 and 21 dpf respectively. Similar to those in early lung development, conserved expression of Hedgehog (Hh) genes, shha and ihha, in the epithelia, and Hh receptor genes, ptc1 and ptc2, as well as fgf10a in mesenchyme was observed. By analyzing several mutants affecting Hh signaling and Ihha morphants, we demonstrated an essential role of Hh signaling in swimbladder development. Furthermore, time-specific Hh inhibition by cyclopamine revealed different requirements of Hh signaling in the formation and organization of all three tissue layers of swimbladder.

Requirement of vasculogenesis and blood circulation in late stages of liver growth in zebrafish.

BACKGROUND: Early events in vertebrate liver development have been the major focus in previous studies, however, late events of liver organogenesis remain poorly understood. Liver vasculogenesis in vertebrates occurs through the interaction of endoderm-derived liver epithelium and mesoderm-derived endothelial cells (ECs). In zebrafish, although it has been found that ECs are not required for liver budding, how and when the spatio-temporal pattern of liver growth is coordinated with ECs remains to be elucidated. RESULTS: To study the process of liver development and vasculogenesis in vivo, a two-color transgenic zebrafish line Tg(lfabf:dsRed; elaA:EGFP) was generated and named LiPan for liver-specific expression of DsRed RFP and exocrine pancreas-specific expression of GFP. Using the LiPan line, we first followed the dynamic development of liver from live embryos to adult and showed the formation of three distinct yet connected liver lobes during development. The LiPan line was then crossed with Tg(fli1:EGFP)y1 and vascular development in the liver was traced in vivo. Liver vasculogenesis started at 55-58 hpf when ECs first surrounded hepatocytes from the liver bud surface and then invaded the liver to form sinusoids and later the vascular network. Using a novel non-invasive and label-free fluorescence correction spectroscopy, we detected blood circulation in the liver starting at approximately 72 hpf. To analyze the roles of ECs and blood circulation in liver development, both cloche mutants (lacking ECs) and Tnnt2 morphants (no blood circulation) were employed. We found that until 70 hpf liver growth and morphogenesis depended on ECs and nascent sinusoids. After 72 hpf, a functional sinusoidal network was essential for continued liver growth. An absence of blood circulation in Tnnt2 morphants caused defects in liver vasculature and small liver. CONCLUSION: There are two phases of liver development in zebrafish, budding and growth. In the growth phase, there are three distinct stages: avascular growth between 50-55 hpf, where ECs are not required; endothelium-dependent growth, where ECs or sinusoids are required for liver growth between 55-72 hpf before blood circulation in liver sinusoids; and circulation-dependent growth, where the circulation is essential to maintain vascular network and to support continued liver growth after 72 hpf.

Transcriptome kinetics of arsenic-induced adaptive response in zebrafish liver.

Arsenic is a prominent environmental toxicant and carcinogen; however, its molecular mechanism of toxicity and carcinogenicity remains poorly understood. In this study, we performed microarray-based expression profiling on liver of zebrafish exposed to 15 parts/million (ppm) arsenic [As(V)] for 8-96 h to identify global transcriptional changes and biological networks involved in arsenic-induced adaptive responses in vivo. We found that there was an increase of transcriptional activity associated with metabolism, especially for biosyntheses, membrane transporter activities, cytoplasm, and endoplasmic reticulum in the 96 h of arsenic treatment, while transcriptional programs for proteins in catabolism, energy derivation, and stress response remained active throughout the arsenic treatment. Many differentially expressed genes encoding proteins involved in heat shock proteins, DNA damage/repair, antioxidant activity, hypoxia induction, iron homeostasis, arsenic metabolism, and ubiquitin-dependent protein degradation were identified, suggesting strongly that DNA and protein damage as a result of arsenic metabolism and oxidative stress caused major cellular injury. These findings were comparable with those reported in mammalian systems, suggesting that the zebrafish liver coupled with the available microarray technology present an excellent in vivo toxicogenomic model for investigating arsenic toxicity. We proposed an in vivo, acute arsenic-induced adaptive response model of the zebrafish liver illustrating the relevance of many transcriptional activities that provide both global and specific information of a coordinated adaptive response to arsenic in the liver.