We have generated several hundred lines of zebrafish (Danio rerio), each heterozygous for a recessive embryonic lethal mutation. Since many tumor suppressor genes are recessive lethals, we screened our colony for lines that display early mortality and/or gross evidence of tumors. We identified 12 lines with elevated cancer incidence. Fish from these lines develop malignant peripheral nerve sheath tumors, and in some cases also other tumor types, with moderate to very high frequencies. Surprisingly, 11 of the 12 lines were each heterozygous for a mutation in a different ribosomal protein (RP) gene, while one line was heterozygous for a mutation in a zebrafish paralog of the human and mouse tumor suppressor gene, neurofibromatosis type 2. Our findings suggest that many RP genes may act as haploinsufficient tumor suppressors in fish. Many RP genes might also be cancer genes in humans, where their role in tumorigenesis could easily have escaped detection up to now.

Zebrafish embryos are an exceptional system for studying vertebrate development. Historically, studies using zebrafish to uncover key players in developmentally regulated gene expression have entailed detailed analysis of transcription factors. It is now apparent that epigenetic modifications of both DNA and histone tails are equally important in the regulation of gene expression during development. As such, blocking the function of key epigenetic modifiers impairs development, albeit with surprising tissue specificity. For instance, DNA methylation is an important epigenetic mark that is depleted in embryos lacking dnmt1 and uhrf1. These embryos display developmental defects in the eye, liver, pancreas, and larval lethality. Interestingly, human tumors derived from these same organs have aberrant changes in DNA methylation and altered expression of genes that are thought to contribute to formation of these cancers. These observations have provided a mechanistic basis for treating cancer with drugs that block the enzymes that facilitate DNA and histone modifications. Thus, it is important to understand the consequences of targeting these factors in a whole animal. We review the use of zebrafish for probing the genetic, cellular, and physiological response to alterations in the epigenome and highlight exciting data illustrating that epigenetic studies using zebrafish can inform and impact cancer biology.

In addition to having a Cx43 ortholog, the zebrafish genome also contains a Cx43-like gene, Cx40.8. Here, we investigate the expression of cx40.8 in zebrafish fins and the function of Cx40.8 in HeLa cells. We find that cx40.8 is present in the same population of dividing cells as cx43. Unlike Cx43, dye coupling assays suggest that Cx40.8 only inefficiently forms functional gap junction channels. However, co-transfection reveals that Cx40.8 can co-localize with Cx43 in gap junction plaques, and that the resulting plaques contain functional gap junction channels. Together, these data suggest the possibility that Cx40.8 may functionally interact with Cx43 to regulate cell proliferation in vivo.

Emerging evidence indicates that paracrine signals from endothelial cells play a role in tissue differentiation and organ formation [1-3]. Here, we identify a novel role for endothelial cells in modulating hepatocyte polarization during liver organogenesis. We find that in zebrafish, the apical domain of the hepatocytes predicts the location of the intrahepatic biliary network. The refinement of hepatocyte polarization coincides with the invasion of endothelial cells into the liver, and these endothelial cells migrate along the maturing basal surface of the hepatocytes. Using genetic, pharmacological, and transplantation experiments, we provide evidence that endothelial cells influence the polarization of the adjacent hepatocytes. This influence of endothelial cells on hepatocytes is mediated at least in part by the cell-surface protein Heart of glass and contributes to the establishment of coordinately aligned hepatocyte apical membranes and evenly spaced intrahepatic conduits.

The unfolded protein response (UPR) is a complex network of sensors and target genes that ensure efficient folding of secretory proteins in the endoplasmic reticulum (ER). UPR activation is mediated by three main sensors, which regulate the expression of hundreds of targets. UPR activation can result in outcomes ranging from enhanced cellular function to cell dysfunction and cell death. How this pathway causes such different outcomes is unknown. Fatty liver disease (steatosis) is associated with markers of UPR activation and robust UPR induction can cause steatosis; however, in other cases, UPR activation can protect against this disease. By assessing the magnitude of activation of UPR sensors and target genes in the liver of zebrafish larvae exposed to three commonly used ER stressors (tunicamycin, thapsigargin and Brefeldin A), we have identified distinct combinations of UPR sensors and targets (i.e. subclasses) activated by each stressor. We found that only the UPR subclass characterized by maximal induction of UPR target genes, which we term a stressed-UPR, induced steatosis. Principal component analysis demonstrated a significant positive association between UPR target gene induction and steatosis. The same principal component analysis showed significant correlation with steatosis in samples from patients with fatty liver disease. We demonstrate that an adaptive UPR induced by a short exposure to thapsigargin prior to challenging with tunicamycin reduced both the induction of a stressed UPR and steatosis incidence. We conclude that a stressed UPR causes steatosis and an adaptive UPR prevents it, demonstrating that this pathway plays dichotomous roles in fatty liver disease.

Activation of the unfolded protein response (UPR) can be either adaptive or pathological. We term the pathological UPR that causes fatty liver disease a "stressed UPR." Here we investigate the mechanism of stressed UPR activation in zebrafish bearing a mutation in thetrappc11gene, which encodes a component of the transport protein particle (TRAPP) complex.trappc11mutants are characterized by secretory pathway defects, reflecting disruption of the TRAPP complex. In addition, we uncover a defect in protein glycosylation intrappc11mutants that is associated with reduced levels of lipid-linked oligosaccharides (LLOs) and compensatory up-regulation of genes in the terpenoid biosynthetic pathway that produces the LLO anchor dolichol. Treating wild-type larvae with terpenoid or LLO synthesis inhibitors phenocopies the stressed UPR seen intrappc11mutants and is synthetically lethal withtrappc11mutation. We propose that reduced LLO level causing hypoglycosylation is a mechanism of stressed UPR induction intrappc11mutants. Of importance, in human cells, depletion of TRAPPC11, but not other TRAPP components, causes protein hypoglycosylation, and lipid droplets accumulate in fibroblasts from patients with theTRAPPC11mutation. These data point to a previously unanticipated and conserved role for TRAPPC11 in LLO biosynthesis and protein glycosylation in addition to its established function in vesicle trafficking.

Cell-cell communication, facilitating the exchange of small metabolites, ions and second messengers, takes place via aqueous proteinaceous channels called gap junctions. Connexins (cx) are the subunits of a gap junction channel. Mutations in zebrafish cx43 produces the short fin (sof (b123) ) phenotype and is characterized by short fins due to reduced segment length of the bony fin rays and reduced cell proliferation. Previously established results from our lab demonstrate that Cx43 plays a dual role regulating both cell proliferation (growth) and joint formation (patterning) during the process of skeletal morphogenesis. In this study, we show that Hapln1a (Hyaluronan and Proteoglycan Link Protein 1a) functions downstream of cx43. Hapln1a belongs to the family of link proteins that play an important role in stabilizing the ECM by linking the aggregates of hyaluronan and proteoglycans. We validated that hapln1a is expressed downstream of cx43 by in situ hybridization and quantitative RT-PCR methods. Moreover, in situ hybridization at different time points revealed that hapln1a expression peaks at 3 days post amputation. Expression of hapln1a is located in the medial mesenchyme and the in the lateral skeletal precursor cells. Furthermore, morpholino mediated knock-down of hapln1a resulted in reduced fin regenerate length, reduced bony segment length and reduced cell proliferation, recapitulating all the phenotypes of cx43 knock-down. Moreover, Hyaluronic Acid (HA) levels are dramatically reduced in hapln1a knock-down fins, attesting the importance of Hapln1a in stabilizing the ECM. Attempts to place hapln1a in our previously defined cx43-sema3d pathway suggest that hapln1a functions in a parallel genetic pathway. Collectively, our data suggest that Cx43 mediates independent Sema3d and Hapln1a pathways in order to coordinate skeletal growth and patterning.

Gap junctions are proteinaceous channels that reside at the plasma membrane and permit the exchange of ions, metabolites, and second messengers between neighboring cells. Connexin proteins are the subunits of gap junction channels. Mutations in zebrafish cx43 cause the short fin (sof(b123)) phenotype which is characterized by short fins due to defects in length of the bony fin rays. Previous findings from our lab demonstrate that Cx43 is required for both cell proliferation and joint formation during fin regeneration. Here we demonstrate that semaphorin3d (sema3d) functions downstream of Cx43. Semas are secreted signaling molecules that have been implicated in diverse cellular functions such as axon guidance, cell migration, cell proliferation, and gene expression. We suggest that Sema3d mediates the Cx43-dependent functions on cell proliferation and joint formation. Using both in situ hybridization and quantitative RT-PCR, we validated that sema3d expression depends on Cx43 activity. Next, we found that knockdown of Sema3d recapitulates all of the sof(b123) and cx43-knockdown phenotypes, providing functional evidence that Sema3d acts downstream of Cx43. To identify the potential Sema3d receptor(s), we evaluated gene expression of neuropilins and plexins. Of these, nrp2a, plxna1, and plxna3 are expressed in the regenerating fin. Morpholino-mediated knockdown of plxna1 did not cause cx43-specific defects, suggesting that PlexinA1 does not function in this pathway. In contrast, morpholino-mediated knockdown of nrp2a caused fin overgrowth and increased cell proliferation, but did not influence joint formation. Moreover, morpholino-mediated knockdown of plxna3 caused short segments, influencing joint formation, but did not alter cell proliferation. Together, our findings reveal that Sema3d functions in a common molecular pathway with Cx43. Furthermore, functional evaluation of putative Sema3d receptors suggests that Cx43-dependent cell proliferation and joint formation utilize independent membrane-bound receptors to mediate downstream cellular phenotypes.

With the emergence of zebrafish as an important model organism, a concerted effort has been made to study its transcriptome. This effort is limited, however, by gaps in zebrafish annotation, which are especially pronounced concerning transcripts dynamically expressed during zygotic genome activation (ZGA). To date, short-read sequencing has been the principal technology for zebrafish transcriptome annotation. In part because these sequence reads are too short for assembly methods to resolve the full complexity of the transcriptome, the current annotation is rudimentary. By providing direct observation of full-length transcripts, recently refined long-read sequencing platforms can dramatically improve annotation coverage and accuracy. Here, we leveraged the SMRT platform to study the transcriptome of zebrafish embryos before and after ZGA. Our analysis revealed additional novelty and complexity in the zebrafish transcriptome, identifying 2539 high-confidence novel transcripts that originated from previously unannotated loci and 1835 high-confidence new isoforms in previously annotated genes. We validated these findings using a suite of computational approaches including structural prediction, sequence homology, and functional conservation analyses, as well as by confirmatory transcript quantification with short-read sequencing data. Our analyses provided insight into new homologs and paralogs of functionally important proteins and noncoding RNAs, isoform switching occurrences, and different classes of novel splicing events. Several novel isoforms representing distinct splicing events were validated through PCR experiments, including the discovery and validation of a novel 8-kb transcript spanning multiple mir-430 elements, an important driver of early development. Our study provides a significantly improved zebrafish transcriptome annotation resource.

Acquisition of cellular fate during development is initiated and maintained by well-coordinated patterns of gene expression that are dictated by the epigenetic landscape and genome organization in the nucleus. While the epigenetic marks that mediate developmental gene expression patterns during organogenesis have been well studied, less is known about how epigenetic marks influence nuclear organization during development. This study examines the relationship between nuclear structure, chromatin accessibility, DNA methylation, and gene expression during hepatic outgrowth in zebrafish larvae. We investigate the relationship between these features using mutants that lack DNA methylation. Hepatocyte nuclear morphology was established coincident with hepatocyte differentiation at 80 h post-fertilization (hpf), and nuclear shape and size continued to change until the conclusion of outgrowth and morphogenesis at 120 hpf. Integrating ATAC-Seq analysis with DNA methylation profiling of zebrafish livers at 120 hpf showed that closed and highly methylated chromatin occupies most transposable elements and that open chromatin correlated with gene expression. DNA hypomethylation, due to mutation of genes encoding ubiquitin-like, containing PHD and RING Finger Domains 1 (uhrf1) and DNA methyltransferase (dnmt1), did not block hepatocyte differentiation, but had dramatic effects on nuclear organization. Hepatocytes in uhrf1 mutants have large, deformed nuclei with multiple nucleoli, downregulation of nucleolar genes, and a complete lack of the nuclear lamina. Loss of lamin B2 staining was phenocopied by dnmt1 mutation. Together, these data show that hepatocyte nuclear morphogenesis coincides with organ morphogenesis and outgrowth, and that DNA methylation directs chromatin organization, and, in turn, hepatocyte nuclear shape and size during liver development.

Extracellular matrix plays a dynamic role during the process of wound healing, embryogenesis and tissue regeneration. Caudal fin regeneration in zebrafish is an excellent model to study tissue and skeletal regeneration. We have analyzed the expression pattern of some of the well characterized ECM proteins during the process of caudal fin regeneration in zebrafish. Our results show that a transitional matrix analogous to the one formed during newt skeletal and heart muscle regeneration is synthesized during fin regeneration. Here we demonstrate that a provisional matrix rich in hyaluronic acid, tenascin C, and fibronectin is synthesized following amputation. Additionally, we observed that the link protein Hapln1a dependent ECM, consisting of Hapln1a, hyaluronan and proteoglycan aggrecan, is upregulated during fin regeneration. Laminin, the protein characteristic of differentiated tissues, showed only modest change in the expression pattern. Our findings on zebrafish fin regeneration implicates that changes in the extracellular milieu represent an evolutionarily conserved mechanism that proceeds during tissue regeneration, yet with distinct players depending on the type of tissue that is involved.

Ubiquitin-like, containing PHD and RING finger domains 1 (uhrf1) is regulated at the transcriptional level during the cell cycle and in developing zebrafish embryos. We identify phosphorylation as a novel means of regulating UHRF1 and demonstrate that Uhrf1 phosphorylation is required for gastrulation in zebrafish. Human UHRF1 contains a conserved cyclin-dependent kinase 2 (CDK2) phosphorylation site at Ser-661 that is phosphorylated in vitro by CDK2 partnered with cyclin A2 (CCNA2), but not cyclin E. An antibody specific for phospho-Ser-661 recognizes UHRF1 in both mammalian cancer cells and in nontransformed zebrafish cells, but not in zebrafish bearing a mutation in ccna2. Depleting Uhrf1 from zebrafish embryos by morpholino injection causes arrest before gastrulation and early embryonic death. This phenotype is rescued by wild-type UHRF1, but not by UHRF1 in which the phospho-acceptor site is mutated, demonstrating that UHRF1 phosphorylation is essential for embryogenesis. UHRF1 was detected in the nucleus and cytoplasm, whereas nonphosphorylatable UHRF1 is unable to localize to the cytoplasm, suggesting the importance of localization in UHRF1 function. Together, these data point to an essential role for UHRF1 phosphorylation by CDK/CCNA2 during early vertebrate development.

Fatty liver disease (FLD) is characterized by lipid accumulation in hepatocytes and is accompanied by secretory pathway dysfunction, resulting in induction of the unfolded protein response (UPR). Activating transcription factor 6 (ATF6), one of three main UPR sensors, functions to both promote FLD during acute stress and reduce FLD during chronic stress. There is little mechanistic understanding of how ATF6, or any other UPR factor, regulates hepatic lipid metabolism to cause disease. We addressed this using zebrafish genetics and biochemical analyses and demonstrate that Atf6 is necessary and sufficient for FLD. atf6 transcription is significantly upregulated in the liver of zebrafish with alcoholic FLD and morpholino-mediated atf6 depletion significantly reduced steatosis incidence caused by alcohol. Moreover, overexpression of active, nuclear Atf6 (nAtf6) in hepatocytes caused FLD in the absence of stress. mRNA-Seq and qPCR analyses of livers from five day old nAtf6 transgenic larvae revealed upregulation of genes promoting glyceroneogenesis and fatty acid elongation, including fatty acid synthase (fasn), and nAtf6 overexpression in both zebrafish larvae and human hepatoma cells increased the incorporation of 14C-acetate into lipids. Srebp transcription factors are key regulators of lipogenic enzymes, but reducing Srebp activation by scap morpholino injection neither prevented FLD in nAtf6 transgenics nor synergized with atf6 knockdown to reduce alcohol-induced FLD. In contrast, fasn morpholino injection reduced FLD in nAtf6 transgenic larvae and synergistically interacted with atf6 to reduce alcoholic FLD. Thus, our data demonstrate that Atf6 is required for alcoholic FLD and epistatically interacts with fasn to cause this disease, suggesting triglyceride biogenesis as the mechanism of UPR induced FLD.

Individuals with congenital disorders of glycosylation (CDG) have recessive mutations in genes required for protein N-glycosylation, resulting in multi-systemic disease. Despite the well-characterized biochemical consequences in these individuals, the underlying cellular defects that contribute to CDG are not well understood. Synthesis of the lipid-linked oligosaccharide (LLO), which serves as the sugar donor for the N-glycosylation of secretory proteins, requires conversion of fructose-6-phosphate to mannose-6-phosphate via the phosphomannose isomerase (MPI) enzyme. Individuals who are deficient in MPI present with bleeding, diarrhea, edema, gastrointestinal bleeding and liver fibrosis. MPI-CDG patients can be treated with oral mannose supplements, which is converted to mannose-6-phosphate through a minor complementary metabolic pathway, restoring protein glycosylation and ameliorating most symptoms, although liver disease continues to progress. Because Mpi deletion in mice causes early embryonic lethality and thus is difficult to study, we used zebrafish to establish a model of MPI-CDG. We used a morpholino to block mpi mRNA translation and established a concentration that consistently yielded 13% residual Mpi enzyme activity at 4 days post-fertilization (dpf), which is within the range of MPI activity detected in fibroblasts from MPI-CDG patients. Fluorophore-assisted carbohydrate electrophoresis detected decreased LLO and N-glycans in mpi morphants. These deficiencies resulted in 50% embryonic lethality by 4 dpf. Multi-systemic abnormalities, including small eyes, dysmorphic jaws, pericardial edema, a small liver and curled tails, occurred in 82% of the surviving larvae. Importantly, these phenotypes could be rescued with mannose supplementation. Thus, parallel processes in fish and humans contribute to the phenotypes caused by Mpi depletion. Interestingly, mannose was only effective if provided prior to 24 hpf. These data provide insight into treatment efficacy and the broader molecular and developmental abnormalities that contribute to disorders associated with defective protein glycosylation.

Fin length in the zebrafish is achieved by the distal addition of bony segments of the correct length. Genetic and molecular data provided evidence that segment growth uses a single pulse of growth, followed by a period of stasis. Examination of cell proliferation during segment growth was predicted to expose a graphical model consistent with a single burst of cell division (e.g., constant, parabolic, or exponential decay) during the lengthening of the distal-most segment. Cell proliferation was detected either by labeling animals with bromodeoxyuridine (during S-phase) or monitoring histone3-phosphate (mitosis). Results from both methods revealed that the number of proliferating cells fluctuates in apparent pulses as a segment grows (i.e., during the growth phase). Thus, rather than segment size being the result of a single burst of proliferation, it appears that segment growth is the result of several pulses of cell division that occur approximately every 60 microns (average segment length approximately 250 microns). These results indicate that segment lengthening requires multiple pulses of cell proliferation.

How joints are correctly positioned in the vertebrate skeleton remains poorly understood. From our studies on the regenerating fin, we have evidence that the gap junction protein Cx43 suppresses joint formation by suppressing the expression of the evx1 transcription factor. Joint morphogenesis proceeds through at least two discrete stages. First, cells that will produce the joint condense in a single row on the bone matrix ("initiation"). Second, these cells separate coincident with articulation of the bone matrix. We propose that Cx43 activity is transiently reduced prior to joint initiation.

Transmembrane protein 184A (TMEM184A) was recently identified as the heparin receptor in vascular cells. Heparin binds specifically to TMEM184A and induces anti-proliferative signaling in vitro. Though it is highly conserved, the physiological function of TMEM184A remains unknown. The objective of this study was to investigate the expression and effects on vascular regeneration of TMEM184A using the adult zebrafish regenerating caudal fin as an in vivo model. Here, we show that Tmem184a is expressed in vascular endothelial cells (ECs) of mature and regenerating zebrafish fins. Transient morpholino (MO)-mediated knockdown of Tmem184a using two validated MOs results in tangled regenerating vessels that do not grow outward and limit normal overall fin regeneration. A significant increase in EC proliferation is observed. Consistent with in vitro work with tissue culture vascular cells, heparin has the opposite effect and decreases EC proliferation which also hinders overall fin regeneration. Collectively, our study suggests that Tmem184a is a novel regulator of angiogenesis.

Angiogenesis, the outgrowth of new blood vessels from existing vasculature, is critical during development, tissue formation, and wound healing. In response to vascular endothelial growth factors (VEGFs), endothelial cells are activated to proliferate and move towards the signal, extending the vessel. These events are directed by VEGF-VEGF receptor (Vegfr2) signal transduction, which in turn is modulated by heparan sulfate proteoglycans (HSPGs). HSPGs are glycoproteins covalently attached to HS glycosaminoglycan chains. Transmembrane protein 184a (Tmem184a) has been recently identified as a heparin receptor, which is believed to bind heparan sulfate chains in vivo. Therefore, Tmem184a has the potential to fine-tune interactions between VEGF and HS, modulating Vegfr2-dependent angiogenesis. The function of Tmem184a has been investigated in the regenerating zebrafish caudal fin, but its role has yet to be evaluated during developmental angiogenesis. Here we provide insights into how Tmem184a contributes to the proper formation of the vasculature in zebrafish embryos. First, we find that knockdown of Tmem184a causes a reduction in the number of intact intersegmental vessels (ISVs) in the zebrafish embryo. This phenotype mimics that of vegfr2b knockout mutants, which have previously been shown to exhibit severe defects in ISV development. We then test the importance of HS interactions by removing the binding domain within the Tmem184a protein, which has a negative effect on angiogenesis. Tmem184a is found to act synergistically with Vegfr2b, indicating that the two gene products function in a common pathway to modulate angiogenesis. Moreover, we find that knockdown of Tmem184a leads to an increase in endothelial cell proliferation but a decrease in the amount of VE-cadherin present. Together, these findings suggest that Tmem184a is necessary for ISVs to organize into mature, complete vessels.

Plexins (plxns) are transmembrane (TM) receptors involved in the guidance of vascular, lymphatic vessel, and neuron growth as well as cancer metastasis. Plxn signaling results in cytosolic GTPase-activating protein activity, and previous research implicates dimerization as important for activation of plxn signaling. Purified, soluble plxn extracellular and cytosolic domains exhibit only weak homomeric interactions, suggesting a role for the plxn TM and juxtamembrane regions in homooligomerization. In this study, we consider a heptad repeat in the Danio rerio PlxnA3 cytosolic juxtamembrane domain (JM) for its ability to influence PlxnA3 homooligomerization in TM-domain containing constructs. Site-directed mutagenesis in conjunction with the AraTM assay and bioluminescent energy transfer (BRET²) suggest an interface involving a JM heptad repeat, in particular residue M1281, regulates PlxnA3 homomeric interactions when examined in constructs containing an ectodomain, TM and JM domain. In the presence of a neuropilin-2a co-receptor and semaphorin 3F ligand, disruption to PlxnA3 homodimerization caused by an M1281F mutation is eliminated, suggesting destabilization of the PlxnA3 homodimer in the JM is not sufficient to disrupt co-receptor complex formation. In contrast, enhanced homodimerization of PlxnA3 caused by mutation M1281L remains even in the presence of ligand semaphorin 3F and co-receptor neuropilin-2a. Consistent with this pattern of PlxnA3 dimerization in the presence of ligand and co-receptor, destabilizing mutations to PlxnA3 homodimerization (M1281F) are able to rescue motor patterning defects in sidetracked zebrafish embryos, whereas mutations that enhance PlxnA3 homodimerization (M1281L) are not. Collectively, our results indicate the JM heptad repeat, in particular residue M1281, forms a switchable interface that modulates both PlxnA3 homomeric interactions and signal transduction.

During fin regeneration, osteoblasts must continually differentiate for outgrowth of the bony fin rays. Bone maturity increases in a distal-proximal manner, and osteoblast maturation can be detected similarly when following gene expression. We find that early markers for osteoblast differentiation are expressed in a discrete domain at the distal end of the fin, just proximal to the adjacent germinal compartment of dividing cells. Matrix genes, required at later stages developmentally, are expressed in a population of cells proximally to the early genes. A marker for mature osteoblasts is expressed in cells further proximal. These domains of gene expression are partially overlapping, perhaps revealing additional levels of osteoblast maturity. We suggest a model for growth where new cells are continually added to the distal-most osteoblast compartment, while osteoblasts in more proximal locations differentiate, thus translating developmental time to location on the proximal-distal axis.