Pla2g12b drives expansion of triglyceride-rich lipoproteins

Vertebrates transport hydrophobic triglycerides through the circulatory system by packaging them within amphipathic particles called Triglyceride-Rich Lipoproteins. Yet, it remains largely unknown how triglycerides are loaded onto these particles. Mutations in Phospholipase A2 group 12B (PLA2G12B) are known to disrupt lipoprotein homeostasis, but its mechanistic role in this process remains unclear. Here we report that PLA2G12B channels lipids within the lumen of the endoplasmic reticulum into nascent lipoproteins. This activity promotes efficient lipid secretion while preventing excess accumulation of intracellular lipids. We characterize the functional domains, subcellular localization, and interacting partners of PLA2G12B, demonstrating that PLA2G12B is calcium-dependent and tightly associated with the membrane of the endoplasmic reticulum. We also detect profound resistance to atherosclerosis in PLA2G12B mutant mice, suggesting an evolutionary tradeoff between triglyceride transport and cardiovascular disease risk. Here we identify PLA2G12B as a key driver of triglyceride incorporation into vertebrate lipoproteins.

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.

Transcriptional profiling of zebrafish identifies host factors controlling susceptibility to Shigella flexneri

Shigella flexneri is a human-adapted pathovar of Escherichia coli that can invade the intestinal epithelium, causing inflammation and bacillary dysentery. Although an important human pathogen, the host response to S. flexneri has not been fully described. Zebrafish larvae represent a valuable model for studying human infections in vivo. Here, we use a Shigella-zebrafish infection model to generate mRNA expression profiles of host response to Shigella infection at the whole-animal level. Immune response-related processes dominate the signature of early Shigella infection (6 h post-infection). Consistent with its clearance from the host, the signature of late Shigella infection (24 h post-infection) is significantly changed, and only a small set of immune-related genes remain differentially expressed, including acod1 and gpr84. Using mutant lines generated by ENU, CRISPR mutagenesis and F0 crispants, we show that acod1- and gpr84-deficient larvae are more susceptible to Shigella infection. Together, these results highlight the power of zebrafish to model infection by bacterial pathogens and reveal the mRNA expression of the early (acutely infected) and late (clearing) host response to Shigella infection.

Selenoprotein deficiency disorder predisposes to aortic aneurysm formation

Aortic aneurysms, which may dissect or rupture acutely and be lethal, can be a part of multisystem disorders that have a heritable basis. We report four patients with deficiency of selenocysteine-containing proteins due to selenocysteine Insertion Sequence Binding Protein 2 (SECISBP2) mutations who show early-onset, progressive, aneurysmal dilatation of the ascending aorta due to cystic medial necrosis. Zebrafish and male mice with global or vascular smooth muscle cell (VSMC)-targeted disruption of Secisbp2 respectively show similar aortopathy. Aortas from patients and animal models exhibit raised cellular reactive oxygen species, oxidative DNA damage and VSMC apoptosis. Antioxidant exposure or chelation of iron prevents oxidative damage in patient's cells and aortopathy in the zebrafish model. Our observations suggest a key role for oxidative stress and cell death, including via ferroptosis, in mediating aortic degeneration.

Gaucher disease protects against tuberculosis

Biallelic mutations in the glucocerebrosidase (GBA1) gene cause Gaucher disease, characterized by lysosomal accumulation of glucosylceramide and glucosylsphingosine in macrophages. Gaucher and other lysosomal diseases occur with high frequency in Ashkenazi Jews. It has been proposed that the underlying mutations confer a selective advantage, in particular conferring protection against tuberculosis. Here, using a zebrafish Gaucher disease model, we find that the mutation GBA1 N370S, predominant among Ashkenazi Jews, increases resistance to tuberculosis through the microbicidal activity of glucosylsphingosine in macrophage lysosomes. Consistent with lysosomal accumulation occurring only in homozygotes, heterozygotes remain susceptible to tuberculosis. Thus, our findings reveal a mechanistic basis for protection against tuberculosis by GBA1 N370S and provide biological plausibility for its selection if the relatively mild deleterious effects in homozygotes were offset by significant protection against tuberculosis, a rampant killer of the young in Europe through the Middle Ages into the 19th century.

P38 Mediates Tumor Suppression through Reduced Autophagy and Actin Cytoskeleton Changes in NRAS-Mutant Melanoma

Hotspot mutations in the NRAS gene are causative genetic events associated with the development of melanoma. Currently, there are no FDA-approved drugs directly targeting NRAS mutations. Previously, we showed that p38 acts as a tumor suppressor in vitro and in vivo with respect to NRAS-mutant melanoma. We observed that because of p38 activation through treatment with the protein synthesis inhibitor, anisomycin leads to a transient upregulation of several targets of the cAMP pathway, representing a stressed cancer cell state that is often observed by therapeutic doses of MAPK inhibitors in melanoma patients. Meanwhile, genetically induced p38 or its stable transduction leads to a distinct cellular transcriptional state. Contrary to previous work showing an association of invasiveness with high p38 levels in BRAF-mutated melanoma, there was no correlation of p38 expression with NRAS-mutant melanoma invasion, highlighting the difference in BRAF and NRAS-driven melanomas. Although the role of p38 has been reported to be that of both tumor suppressor and oncogene, we show here that p38 specifically plays the role of a tumor suppressor in NRAS-mutant melanoma. Both the transient and stable activation of p38 elicits phosphorylation of mTOR, reported to be a master switch in regulating autophagy. Indeed, we observed a correlation between elevated levels of phosphorylated mTOR and a reduction in LC3 conversion (LCII/LCI), indicative of suppressed autophagy. Furthermore, a reduction in actin intensity in p38-high cells strongly suggests a role of mTOR in regulating actin and a remodeling in the NRAS-mutant melanoma cells. Therefore, p38 plays a tumor suppressive role in NRAS-mutant melanomas at least partially through the mechanism of mTOR upregulation, suppressed autophagy, and reduced actin polymerization. One or more combinations of MEK inhibitors with either anisomycin, rapamycin, chloroquine/bafilomycin, and cytochalasin modulate p38 activation, mTOR phosphorylation, autophagy, and actin polymerization, respectively, and they may provide an alternate route to targeting NRAS-mutant melanoma.

Conserved roles for Hnf4 family transcription factors in zebrafish development and intestinal function

Transcription factors play important roles in the development of the intestinal epithelium and its ability to respond to endocrine, nutritional, and microbial signals. Hepatocyte nuclear factor 4 family nuclear receptors are liganded transcription factors that are critical for the development and function of multiple digestive organs in vertebrates, including the intestinal epithelium. Zebrafish have 3 hepatocyte nuclear factor 4 homologs, of which, hnf4a was previously shown to mediate intestinal responses to microbiota in zebrafish larvae. To discern the functions of other hepatocyte nuclear factor 4 family members in zebrafish development and intestinal function, we created and characterized mutations in hnf4g and hnf4b. We addressed the possibility of genetic redundancy amongst these factors by creating double and triple mutants which showed different rates of survival, including apparent early lethality in hnf4a; hnf4b double mutants and triple mutants. RNA sequencing performed on digestive tracts from single and double mutant larvae revealed extensive changes in intestinal gene expression in hnf4a mutants that were amplified in hnf4a; hnf4g mutants, but limited in hnf4g mutants. Changes in hnf4a and hnf4a; hnf4g mutants were reminiscent of those seen in mice including decreased expression of genes involved in intestinal function and increased expression of cell proliferation genes, and were validated using transgenic reporters and EdU labeling in the intestinal epithelium. Gnotobiotics combined with RNA sequencing also showed hnf4g has subtler roles than hnf4a in host responses to microbiota. Overall, phenotypic changes in hnf4a single mutants were strongly enhanced in hnf4a; hnf4g double mutants, suggesting a conserved partial genetic redundancy between hnf4a and hnf4g in the vertebrate intestine.

mTOR-regulated mitochondrial metabolism limits mycobacterium-induced cytotoxicity

Necrosis of macrophages in the granuloma, the hallmark immunological structure of tuberculosis, is a major pathogenic event that increases host susceptibility. Through a zebrafish forward genetic screen, we identified the mTOR kinase, a master regulator of metabolism, as an early host resistance factor in tuberculosis. We found that mTOR complex 1 protects macrophages from mycobacterium-induced death by enabling infection-induced increases in mitochondrial energy metabolism fueled by glycolysis. These metabolic adaptations are required to prevent mitochondrial damage and death caused by the secreted mycobacterial virulence determinant ESAT-6. Thus, the host can effectively counter this early critical mycobacterial virulence mechanism simply by regulating energy metabolism, thereby allowing pathogen-specific immune mechanisms time to develop. Our findings may explain why Mycobacterium tuberculosis, albeit humanity's most lethal pathogen, is successful in only a minority of infected individuals.

Loss of slc39a14 causes simultaneous manganese hypersensitivity and deficiency in zebrafish

Manganese neurotoxicity is a hallmark of hypermanganesemia with dystonia 2, an inherited manganese transporter defect caused by mutations in SLC39A14. To identify novel potential targets of manganese neurotoxicity, we performed transcriptome analysis of slc39a14-/- mutant zebrafish that were exposed to MnCl2. Differentially expressed genes mapped to the central nervous system and eye, and pathway analysis suggested that Ca2+ dyshomeostasis and activation of the unfolded protein response are key features of manganese neurotoxicity. Consistent with this interpretation, MnCl2 exposure led to decreased whole-animal Ca2+ levels, locomotor defects and changes in neuronal activity within the telencephalon and optic tectum. In accordance with reduced tectal activity, slc39a14-/- zebrafish showed changes in visual phototransduction gene expression, absence of visual background adaptation and a diminished optokinetic reflex. Finally, numerous differentially expressed genes in mutant larvae normalised upon MnCl2 treatment indicating that, in addition to neurotoxicity, manganese deficiency is present either subcellularly or in specific cells or tissues. Overall, we assembled a comprehensive set of genes that mediate manganese-systemic responses and found a highly correlated and modulated network associated with Ca2+ dyshomeostasis and cellular stress. This article has an associated First Person interview with the first author of the paper.

Translational relevance of forward genetic screens in animal models for the study of psychiatric disease

Psychiatric disorders represent a significant burden in our societies. Despite the convincing evidence pointing at gene and gene-environment interaction contributions, the role of genetics in the etiology of psychiatric disease is still poorly understood. Forward genetic screens in animal models have helped elucidate causal links. Here we discuss the application of mutagenesis-based forward genetic approaches in common animal model species: two invertebrates, nematodes (Caenorhabditis elegans) and fruit flies (Drosophila sp.); and two vertebrates, zebrafish (Danio rerio) and mice (Mus musculus), in relation to psychiatric disease. We also discuss the use of large scale genomic studies in human populations. Despite the advances using data from human populations, animal models coupled with next-generation sequencing strategies are still needed. Although with its own limitations, zebrafish possess characteristics that make them especially well-suited to forward genetic studies exploring the etiology of psychiatric disorders.

Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants

In model organisms, RNA-sequencing (RNA-seq) is frequently used to assess the effect of genetic mutations on cellular and developmental processes. Typically, animals heterozygous for a mutation are crossed to produce offspring with different genotypes. Resultant embryos are grouped by genotype to compare homozygous mutant embryos to heterozygous and wild-type siblings. Genes that are differentially expressed between the groups are assumed to reveal insights into the pathways affected by the mutation. Here we show that in zebrafish, differentially expressed genes are often over-represented on the same chromosome as the mutation due to different levels of expression of alleles from different genetic backgrounds. Using an incross of haplotype-resolved wild-type fish, we found evidence of widespread allele-specific expression, which appears as differential expression when comparing embryos homozygous for a region of the genome to their siblings. When analysing mutant transcriptomes, this means that the differential expression of genes on the same chromosome as a mutation of interest may not be caused by that mutation. Typically, the genomic location of a differentially expressed gene is not considered when interpreting its importance with respect to the phenotype. This could lead to pathways being erroneously implicated or overlooked due to the noise of spurious differentially expressed genes on the same chromosome as the mutation. These observations have implications for the interpretation of RNA-seq experiments involving outbred animals and non-inbred model organisms.

Ankk1 Loss of Function Disrupts Dopaminergic Pathways in Zebrafish

Ankyrin repeat and kinase domain containing 1 (ANKK1) is a member of the receptor-interacting protein serine/threonine kinase family, known to be involved in cell proliferation, differentiation and activation of transcription factors. Genetic variation within the ANKK1 locus is suggested to play a role in vulnerability to addictions. However, ANKK1 mechanism of action is still poorly understood. It has been suggested that ANKK1 may affect the development and/or functioning of dopaminergic pathways. To test this hypothesis, we generated a CRISPR-Cas9 loss of function ankk1 zebrafish line causing a 27 bp insertion that disrupts the ankk1 sequence introducing an early stop codon. We found that ankk1 transcript levels were significantly lower in ankk1 mutant (ankk127ins) fish compared to their wild type (ankk1+/+) siblings. In ankk1+/+ adult zebrafish brain, ankk1 protein was detected in isocortex, hippocampus, basolateral amygdala, mesencephalon, and cerebellum, resembling the mammalian distribution pattern. In contrast, ankk1 protein was reduced in the brain of ankk127ins/27ins fish. Quantitative polymerase chain reaction analysis revealed an increase in expression of drd2b mRNA in ankk127ins at both larval and adult stages. In ankk1+/+ adult zebrafish brain, drd2 protein was detected in cerebral cortex, cerebellum, hippocampus, and caudate homolog regions, resembling the pattern in humans. In contrast, drd2 expression was reduced in cortical regions of ankk127ins/27ins being predominantly found in the hindbrain. No differences in the number of cell bodies or axonal projections detected by anti-tyrosine hydroxylase immunostaining on 3 days post fertilization (dpf) larvae were found. Behavioral analysis revealed altered sensitivity to effects of both amisulpride and apomorphine on locomotion and startle habituation, consistent with a broad loss of both pre and post synaptic receptors. Ankk127ins mutants showed reduced sensitivity to the effect of the selective dopamine receptor antagonist amisulpride on locomotor responses to acoustic startle and were differentially sensitive to the effects of the non-selective dopamine agonist apomorphine on both locomotion and habituation. Taken together, our findings strengthen the hypothesis of a functional relationship between ANKK1 and DRD2, supporting a role for ANKK1 in the maintenance and/or functioning of dopaminergic pathways. Further work is needed to disentangle ANKK1’s role at different developmental stages.

Behavioral and Gene Regulatory Responses to Developmental Drug Exposures in Zebrafish

Developmental consequences of prenatal drug exposure have been reported in many human cohorts and animal studies. The long-lasting impact on the offspring-including motor and cognitive impairments, cranial and cardiac anomalies and increased prevalence of ADHD-is a socioeconomic burden worldwide. Identifying the molecular changes leading to developmental consequences could help ameliorate the deficits and limit the impact. In this study, we have used zebrafish, a well-established behavioral and genetic model with conserved drug response and reward pathways, to identify changes in behavior and cellular pathways in response to developmental exposure to amphetamine, nicotine or oxycodone. In the presence of the drug, exposed animals showed altered behavior, consistent with effects seen in mammalian systems, including impaired locomotion and altered habituation to acoustic startle. Differences in responses seen following acute and chronic exposure suggest adaptation to the presence of the drug. Transcriptomic analysis of exposed larvae revealed differential expression of numerous genes and alterations in many pathways, including those related to cell death, immunity and circadian rhythm regulation. Differential expression of circadian rhythm genes did not correlate with behavioral changes in the larvae, however, two of the circadian genes, arntl2 and per2, were also differentially expressed at later stages of development, suggesting a long-lasting impact of developmental exposures on circadian gene expression. The immediate-early genes, egr1, egr4, fosab, and junbb, which are associated with synaptic plasticity, were downregulated by all three drugs and in situ hybridization showed that the expression for all four genes was reduced across all neuroanatomical regions, including brain regions implicated in reward processing, addiction and other psychiatric conditions. We anticipate that these early changes in gene expression in response to drug exposure are likely to contribute to the consequences of prenatal exposure and their discovery might pave the way to therapeutic intervention to ameliorate the long-lasting deficits.

Total Nucleic Acid Extraction from Single Zebrafish Embryos for Genotyping and RNA-seq

RNA sequencing allows for the quantification of the transcriptome of embryos to investigate transcriptional responses to various perturbations (e.g., mutations, infections, drug treatments). Previous protocols either lack the option to genotype individual samples, or are laborious and therefore difficult to do at a large scale. We have developed a protocol to extract total nucleic acid from individual zebrafish embryos. Individual embryos are lysed in 96-well plates and nucleic acid is extracted using SPRI beads. The total nucleic acid can be genotyped and then DNase I treated to produce RNA samples for sequencing. This protocol allows for processing large numbers of individual samples, with the ability to genotype each sample, which makes it possible to undertake transcriptomic analysis on mutants at timepoints before the phenotype is visible. Graphic abstract: Extraction of total nucleic acid from individual zebrafish embryos for genotyping and RNA-seq.

NRASQ61K melanoma tumor formation is reduced by p38-MAPK14 activation in zebrafish models and NRAS-mutated human melanoma cells

Oncogenic BRAF and NRAS mutations drive human melanoma initiation. We used transgenic zebrafish to model NRAS-mutant melanoma, and the rapid tumor onset allowed us to study candidate tumor suppressors. We identified P38α-MAPK14 as a potential tumor suppressor in The Cancer Genome Atlas melanoma cohort of NRAS-mutant melanomas, and overexpression significantly increased the time to tumor onset in transgenic zebrafish with NRAS-driven melanoma. Pharmacological activation of P38α-MAPK14 using anisomycin reduced in vitro viability of melanoma cultures, which we confirmed by stable overexpression of p38α. We observed that the viability of MEK inhibitor resistant melanoma cells could be reduced by combined treatment of anisomycin and MEK inhibition. Our study demonstrates that activating the p38α-MAPK14 pathway in the presence of oncogenic NRAS abrogates melanoma in vitro and in vivo. SIGNIFICANCE: The significance of our study is in the accountability of NRAS mutations in melanoma. We demonstrate here that activation of p38α-MAPK14 pathway can abrogate NRAS-mutant melanoma which is contrary to the previously published role of p38α-MAPK14 pathway in BRAF mutant melanoma. These results implicate that BRAF and NRAS-mutant melanoma may not be identical biologically. We also demonstrate the translational benefit of our study by using a small molecule compound-anisomycin (already in use for other diseases in clinical trials) to activate p38α-MAPK14 pathway.

PRL3-DDX21 Transcriptional Control of Endolysosomal Genes Restricts Melanocyte Stem Cell Differentiation

Melanocytes, replenished throughout life by melanocyte stem cells (MSCs), play a critical role in pigmentation and melanoma. Here, we reveal a function for the metastasis-associated phosphatase of regenerating liver 3 (PRL3) in MSC regeneration. We show that PRL3 binds to the RNA helicase DDX21, thereby restricting productive transcription by RNAPII at master transcription factor (MITF)-regulated endolysosomal vesicle genes. In zebrafish, this mechanism controls premature melanoblast expansion and differentiation from MSCs. In melanoma patients, restricted transcription of this endolysosomal vesicle pathway is a hallmark of PRL3-high melanomas. Our work presents the conceptual advance that PRL3-mediated control of transcriptional elongation is a differentiation checkpoint mechanism for activated MSCs and has clinical relevance for the activity of PRL3 in regenerating tissue and cancer.

A point mutation decouples the lipid transfer activities of microsomal triglyceride transfer protein

Apolipoprotein B-containing lipoproteins (B-lps) are essential for the transport of hydrophobic dietary and endogenous lipids through the circulation in vertebrates. Zebrafish embryos produce large numbers of B-lps in the yolk syncytial layer (YSL) to move lipids from yolk to growing tissues. Disruptions in B-lp production perturb yolk morphology, readily allowing for visual identification of mutants with altered B-lp metabolism. Here we report the discovery of a missense mutation in microsomal triglyceride transfer protein (Mtp), a protein that is essential for B-lp production. This mutation of a conserved glycine residue to valine (zebrafish G863V, human G865V) reduces B-lp production and results in yolk opacity due to aberrant accumulation of cytoplasmic lipid droplets in the YSL. However, this phenotype is milder than that of the previously reported L475P stalactite (stl) mutation. MTP transfers lipids, including triglycerides and phospholipids, to apolipoprotein B in the ER for B-lp assembly. In vitro lipid transfer assays reveal that while both MTP mutations eliminate triglyceride transfer activity, the G863V mutant protein unexpectedly retains ~80% of phospholipid transfer activity. This residual phospholipid transfer activity of the G863V mttp mutant protein is sufficient to support the secretion of small B-lps, which prevents intestinal fat malabsorption and growth defects observed in the mttpstl/stl mutant zebrafish. Modeling based on the recent crystal structure of the heterodimeric human MTP complex suggests the G865V mutation may block triglyceride entry into the lipid-binding cavity. Together, these data argue that selective inhibition of MTP triglyceride transfer activity may be a feasible therapeutic approach to treat dyslipidemia and provide structural insight for drug design. These data also highlight the power of yolk transport studies to identify proteins critical for B-lp biology.

Identification of slit3 as a locus affecting nicotine preference in zebrafish and human smoking behaviour

To facilitate smoking genetics research we determined whether a screen of mutagenized zebrafish for nicotine preference could predict loci affecting smoking behaviour. From 30 screened F3 sibling groups, where each was derived from an individual ethyl-nitrosurea mutagenized F0 fish, two showed increased or decreased nicotine preference. Out of 25 inactivating mutations carried by the F3 fish, one in the slit3 gene segregated with increased nicotine preference in heterozygous individuals. Focussed SNP analysis of the human SLIT3 locus in cohorts from UK (n=863) and Finland (n=1715) identified two variants associated with cigarette consumption and likelihood of cessation. Characterisation of slit3 mutant larvae and adult fish revealed decreased sensitivity to the dopaminergic and serotonergic antagonist amisulpride, known to affect startle reflex that is correlated with addiction in humans, and increased htr1aa mRNA expression in mutant larvae. No effect on neuronal pathfinding was detected. These findings reveal a role for SLIT3 in development of pathways affecting responses to nicotine in zebrafish and smoking in humans.

Common and distinct transcriptional signatures of mammalian embryonic lethality

The Deciphering the Mechanisms of Developmental Disorders programme has analysed the morphological and molecular phenotypes of embryonic and perinatal lethal mouse mutant lines in order to investigate the causes of embryonic lethality. Here we show that individual whole-embryo RNA-seq of 73 mouse mutant lines (>1000 transcriptomes) identifies transcriptional events underlying embryonic lethality and associates previously uncharacterised genes with specific pathways and tissues. For example, our data suggest that Hmgxb3 is involved in DNA-damage repair and cell-cycle regulation. Further, we separate embryonic delay signatures from mutant line-specific transcriptional changes by developing a baseline mRNA expression catalogue of wild-type mice during early embryogenesis (4-36 somites). Analysis of transcription outside coding sequence identifies deregulation of repetitive elements in Morc2a mutants and a gene involved in gene-specific splicing. Collectively, this work provides a large scale resource to further our understanding of early embryonic developmental disorders.

The gene regulatory basis of genetic compensation during neural crest induction

The neural crest (NC) is a vertebrate-specific cell type that contributes to a wide range of different tissues across all three germ layers. The gene regulatory network (GRN) responsible for the formation of neural crest is conserved across vertebrates. Central to the induction of the NC GRN are AP-2 and SoxE transcription factors. NC induction robustness is ensured through the ability of some of these transcription factors to compensate loss of function of gene family members. However the gene regulatory events underlying compensation are poorly understood. We have used gene knockout and RNA sequencing strategies to dissect NC induction and compensation in zebrafish. We genetically ablate the NC using double mutants of tfap2a;tfap2c or remove specific subsets of the NC with sox10 and mitfa knockouts and characterise genome-wide gene expression levels across multiple time points. We find that compensation through a single wild-type allele of tfap2c is capable of maintaining early NC induction and differentiation in the absence of tfap2a function, but many target genes have abnormal expression levels and therefore show sensitivity to the reduced tfap2 dosage. This separation of morphological and molecular phenotypes identifies a core set of genes required for early NC development. We also identify the 15 somites stage as the peak of the molecular phenotype which strongly diminishes at 24 hpf even as the morphological phenotype becomes more apparent. Using gene knockouts, we associate previously uncharacterised genes with pigment cell development and establish a role for maternal Hippo signalling in melanocyte differentiation. This work extends and refines the NC GRN while also uncovering the transcriptional basis of genetic compensation via paralogues.

Embryonic development is an intricate process that requires genes to be active at the right time and place. Organisms have evolved mechanisms that ensure faithful execution of developmental programmes even if genes fail to function. For example, in a process called genetic compensation, one or more genes become activated in response to loss of function of another. In this work we use the zebrafish model to investigate how two related genes, tfap2a and tfap2c, interact to ensure establishment of the neural crest, a vertebrate-specific cell type that contributes to many different tissues. Losing tfap2a activity causes mild morphological defects and losing tfap2c has no visible effect. Yet when both are inactive, embryos are severely abnormal due to lack of neural crest-derived tissues. Here we show that loss of tfap2a triggers upregulation of tfap2c which prevents the loss of neural crest tissue. However, the genes normally regulated by tfap2a respond differently to tfap2c allowing us to identify the first tier of the Ap2 network and new players in neural crest biology. Our work demonstrates that the expression signature of partial, but morphologically sufficient, genetic compensation provides an opportunity to dissect gene regulatory networks.

Dicer1 is required for pigment cell and craniofacial development in zebrafish

The multidomain RNase III endoribonuclease DICER is required for the generation of most functional microRNAs (miRNAs). Loss of Dicer affects developmental processes at different levels. Here, we characterized the zebrafish Dicer1 mutant, dicer1^sa9205, which has a single point mutation induced by N-ethyl-N-nitrosourea mutagenesis. Heterozygous dicer1^sa9205 developed normally, being phenotypically indistinguishable from wild-type siblings. Homozygous dicer1^sa9205 mutants display smaller eyes, abnormal craniofacial development and aberrant pigmentation. Reduced numbers of both iridophores and melanocytes were observed in the head and ventral trunk of dicer1^sa9205 homozygotes; the effect on melanocytes was stronger and detectable earlier in development. The expression of microphthalmia-associated transcription factor a (mitfa), the master gene for melanocytes differentiation, was enhanced in dicer1-depleted fish. Similarly, the expression of SRY-box containing gene 10 (sox10), required for mitfa activation, was higher in mutants than in wild types. In silico and in vivo analyses of either sox10 or mitfa 3'UTRs revealed conserved potential miRNA binding sites likely involved in the post-transcriptional regulation of both genes. Based on these findings, we propose that dicer1 participates in the gene regulatory network governing zebrafish melanocyte differentiation by controlling the expression of mitfa and sox10.

Compensatory growth renders Tcf7l1a dispensable for eye formation despite its requirement in eye field specification

The vertebrate eye originates from the eye field, a domain of cells specified by a small number of transcription factors. In this study, we show that Tcf7l1a is one such transcription factor that acts cell-autonomously to specify the eye field in zebrafish. Despite the much-reduced eye field in tcf7l1a mutants, these fish develop normal eyes revealing a striking ability of the eye to recover from a severe early phenotype. This robustness is not mediated through genetic compensation at neural plate stage; instead, the smaller optic vesicle of tcf7l1a mutants shows delayed neurogenesis and continues to grow until it achieves approximately normal size. Although the developing eye is robust to the lack of Tcf7l1a function, it is sensitised to the effects of additional mutations. In support of this, a forward genetic screen identified mutations in hesx1, cct5 and gdf6a, which give synthetically enhanced eye specification or growth phenotypes when in combination with the tcf7l1a mutation.

Left and right eyes develop independently, yet they consistently grow to roughly the same size in humans and other creatures. How they do this remains a mystery, though scientists have learned that both eyes originate from a single group of cells in the developing nervous system called the eye field. As development progresses, the eye field splits in two, and buds into the two separate compartments from which each eye forms. As the eyes grow, the cells in each compartment specialize, or ‘differentiate’, to make working left and right eyes. Scientists often study eye development in zebrafish embryos because it is easy to see each step in the process.

Now, Young at al. show that zebrafish with a mutation that causes the eye field to be half its normal size go on to form normal-sized eyes. Somehow these developing embryos overcome this deleterious mutation. It turns out that the eyes of zebrafish with this mutation grow for a longer period of time than typical zebrafish eyes. This change allows the mutant fish’s eyes to catch up and reach normal size. When Young et al. removed some cells from one of the forming eyes of normal zebrafish embryos they found that same thing happened. The smaller eye developed for a longer time and delayed its differentiation until both eyes were the same size. Conversely, when eyes developed from a larger than normal eye field, growth stopped prematurely and differentiation began early preventing the eyes from ending up oversized.

Though the fish were able to overcome the effects of one mutation to develop normal-sized eyes, adding a second mutation that affected eye development led to unusual sized eyes or absence of eyes. Together the experiments identify genes and mechanisms essential for the formation and size of the eyes. Given that the processes underlying eye formation are very similar in many animals, this new information should help scientists to better understand eye abnormalities in humans.

scRNA-Seq reveals distinct stem cell populations that drive hair cell regeneration after loss of Fgf and Notch signaling

Loss of sensory hair cells leads to deafness and balance deficiencies. In contrast to mammalian hair cells, zebrafish ear and lateral line hair cells regenerate from poorly characterized support cells. Equally ill-defined is the gene regulatory network underlying the progression of support cells to differentiated hair cells. scRNA-Seq of lateral line organs uncovered five different support cell types, including quiescent and activated stem cells. Ordering of support cells along a developmental trajectory identified self-renewing cells and genes required for hair cell differentiation. scRNA-Seq analyses of fgf3 mutants, in which hair cell regeneration is increased, demonstrates that Fgf and Notch signaling inhibit proliferation of support cells in parallel by inhibiting Wnt signaling. Our scRNA-Seq analyses set the foundation for mechanistic studies of sensory organ regeneration and is crucial for identifying factors to trigger hair cell production in mammals. The data is searchable and publicly accessible via a web-based interface.

Deep phenotyping in zebrafish reveals genetic and diet-induced adiposity changes that may inform disease risk

The regional distribution of adipose tissues is implicated in a wide range of diseases. For example, proportional increases in visceral adipose tissue increase the risk for insulin resistance, diabetes, and CVD. Zebrafish offer a tractable model system by which to obtain unbiased and quantitative phenotypic information on regional adiposity, and deep phenotyping can explore complex disease-related adiposity traits. To facilitate deep phenotyping of zebrafish adiposity traits, we used pairwise correlations between 67 adiposity traits to generate stage-specific adiposity profiles that describe changing adiposity patterns and relationships during growth. Linear discriminant analysis classified individual fish according to an adiposity profile with 87.5% accuracy. Deep phenotyping of eight previously uncharacterized zebrafish mutants identified neuropilin 2b as a novel gene that alters adipose distribution. When we applied deep phenotyping to identify changes in adiposity during diet manipulations, zebrafish that underwent food restriction and refeeding had widespread adiposity changes when compared with continuously fed, equivalently sized control animals. In particular, internal adipose tissues (e.g., visceral adipose) exhibited a reduced capacity to replenish lipid following food restriction. Together, these results in zebrafish establish a new deep phenotyping technique as an unbiased and quantitative method to help uncover new relationships between genotype, diet, and adiposity.

The age of heterozygous telomerase mutant parents influences the adult phenotype of their offspring irrespective of genotype in zebrafish

Background: Mutations in proteins involved in telomere maintenance lead to a range of human diseases, including dyskeratosis congenita, idiopathic pulmonary fibrosis and cancer. Telomerase functions to add telomeric repeats back onto the ends of chromosomes, however non-canonical roles of components of telomerase have recently been suggested.

Methods: Here we use a zebrafish telomerase mutant which harbours a nonsense mutation in tert to investigate the adult phenotypes of fish derived from heterozygous parents of different ages. Furthermore we use whole genome sequencing data to estimate average telomere lengths.

Results: We show that homozygous offspring from older heterozygotes exhibit signs of body wasting at a younger age than those of younger parents, and that offspring of older heterozygous parents weigh less irrespective of genotype. We also demonstrate that tert homozygous mutant fish have a male sex bias, and that clutches from older parents also have a male sex bias in the heterozygous and wild-type populations. Telomere length analysis reveals that the telomeres of younger heterozygous parents are shorter than those of older heterozygous parents.

Conclusions: These data indicate that the phenotypes observed in offspring from older parents cannot be explained by telomere length. Instead we propose that Tert functions outside of telomere length maintenance in an age-dependent manner to influence the adult phenotypes of the next generation.

mRNA processing in mutant zebrafish lines generated by chemical and CRISPR-mediated mutagenesis produces unexpected transcripts that escape nonsense-mediated decay

As model organism-based research shifts from forward to reverse genetics approaches, largely due to the ease of genome editing technology, a low frequency of abnormal phenotypes is being observed in lines with mutations predicted to lead to deleterious effects on the encoded protein. In zebrafish, this low frequency is in part explained by compensation by genes of redundant or similar function, often resulting from the additional round of teleost-specific whole genome duplication within vertebrates. Here we offer additional explanations for the low frequency of mutant phenotypes. We analyzed mRNA processing in seven zebrafish lines with mutations expected to disrupt gene function, generated by CRISPR/Cas9 or ENU mutagenesis methods. Five of the seven lines showed evidence of altered mRNA processing: one through a skipped exon that did not lead to a frame shift, one through nonsense-associated splicing that did not lead to a frame shift, and three through the use of cryptic splice sites. These results highlight the need for a methodical analysis of the mRNA produced in mutant lines before making conclusions or embarking on studies that assume loss of function as a result of a given genomic change. Furthermore, recognition of the types of adaptations that can occur may inform the strategies of mutant generation.

The recent rise of reverse genetic, gene targeting methods has allowed researchers to readily generate mutations in any gene of interest with relative ease. Should these mutations have the predicted effect on the mRNA and encoded protein, we would expect many more abnormal phenotypes than are typically being seen in reverse genetic screens. Here we set out to explore some of the reasons for this discrepancy by studying seven separate mutations in zebrafish. We present evidence that thorough cDNA sequence analysis is a key step in assessing the likelihood that a given mutation will produce hypomorphic or null alleles. This study reveals that mRNA processing in the mutant background often produces transcripts that escape nonsense-mediated decay, thereby potentially preserving gene function. By understanding the ways that cells avoid the deleterious consequences of mutations, researchers can better design reverse genetic strategies to increase the likelihood of gene disruption.

A high-resolution mRNA expression time course of embryonic development in zebrafish

We have produced an mRNA expression time course of zebrafish development across 18 time points from 1 cell to 5 days post-fertilisation sampling individual and pools of embryos. Using poly(A) pulldown stranded RNA-seq and a 3′ end transcript counting method we characterise temporal expression profiles of 23,642 genes. We identify temporal and functional transcript co-variance that associates 5024 unnamed genes with distinct developmental time points. Specifically, a class of over 100 previously uncharacterised zinc finger domain containing genes, located on the long arm of chromosome 4, is expressed in a sharp peak during zygotic genome activation. In addition, the data reveal new genes and transcripts, differential use of exons and previously unidentified 3′ ends across development, new primary microRNAs and temporal divergence of gene paralogues generated in the teleost genome duplication. To make this dataset a useful baseline reference, the data can be browsed and downloaded at Expression Atlas and Ensembl.

KDM2A integrates DNA and histone modification signals through a CXXC/PHD module and direct interaction with HP1

Functional genomic elements are marked by characteristic DNA and histone modification signatures. How combinatorial chromatin modification states are recognized by epigenetic reader proteins and how this is linked to their biological function is largely unknown. Here we provide a detailed molecular analysis of chromatin recognition by the lysine demethylase KDM2A. Using biochemical approaches we identify a nucleosome interaction module within KDM2A consisting of a CXXC type zinc finger, a PHD domain and a newly identified Heterochromatin Protein 1 (HP1) interaction motif that mediates direct binding between KDM2A and HP1. This nucleosome interaction module enables KDM2A to decode nucleosomal H3K9me3 modification in addition to CpG methylation signals. The multivalent engagement with DNA and HP1 results in a nucleosome binding circuit in which KDM2A can be recruited to H3K9me3-modified chromatin through HP1, and HP1 can be recruited to unmodified chromatin by KDM2A. A KDM2A mutant deficient in HP1-binding is inactive in an in vivo overexpression assay in zebrafish embryos demonstrating that the HP1 interaction is essential for KDM2A function. Our results reveal a complex regulation of chromatin binding for both KDM2A and HP1 that is modulated by DNA- and H3K9-methylation, and suggest a direct role for KDM2A in chromatin silencing.

The age of heterozygous telomerase mutant parents influences the adult phenotype of their offspring irrespective of genotype in zebrafish

Background: Mutations in proteins involved in telomere maintenance lead to a range of human diseases, including dyskeratosis congenita, idiopathic pulmonary fibrosis and cancer. Telomerase functions to add telomeric repeats back onto the ends of chromosomes, however non-canonical roles of components of telomerase have recently been suggested.Methods: Here we use a zebrafish telomerase mutant which harbours a nonsense mutation intertto investigate the adult phenotypes of fish derived from heterozygous parents of different ages. Furthermore we use whole genome sequencing data to estimate average telomere lengths.Results: We show that homozygous offspring from older heterozygotes exhibit signs of body wasting at a younger age than those of younger parents, and that offspring of older heterozygous parents weigh less irrespective of genotype. We also demonstrate thatterthomozygous mutant fish have a male sex bias, and that clutches from older parents also have a male sex bias in the heterozygous and wild-type populations. Telomere length analysis reveals that the telomeres of younger heterozygous parents are shorter than those of older heterozygous parents.Conclusions: These data indicate that the phenotypes observed in offspring from older parents cannot be explained by telomere length. Instead we propose that Tert functions outside of telomere length maintenance in an age-dependent manner to influence the adult phenotypes of the next generation.

Report of a Meeting on Contemporary Topics in Zebrafish Husbandry and Care

A meeting on Contemporary Topics in Zebrafish Husbandry and Care was held in the United Kingdom in 2014, with the aim of providing a discussion forum for researchers, animal technologists, and veterinarians from academia and industry to share good practice and exchange ideas. Presentation topics included protocols for optimal larval rearing, implementing the 3Rs (replacement, reduction, and refinement) in large-scale colony management, and environmental enrichment. The audience also participated in a survey of current practice relating to practical husbandry, cryopreservation, and the provision of enrichment.

Efficient identification of CRISPR/Cas9-induced insertions/deletions by direct germline screening in zebrafish

BACKGROUND

The CRISPR/Cas9 system is a prokaryotic immune system that infers resistance to foreign genetic material and is a sort of 'adaptive immunity'. It has been adapted to enable high throughput genome editing and has revolutionised the generation of targeted mutations.

RESULTS

We have developed a scalable analysis pipeline to identify CRISPR/Cas9 induced mutations in hundreds of samples using next generation sequencing (NGS) of amplicons. We have used this system to investigate the best way to screen mosaic Zebrafish founder individuals for germline transmission of induced mutations. Screening sperm samples from potential founders provides much better information on germline transmission rates and crucially the sequence of the particular insertions/deletions (indels) that will be transmitted. This enables us to combine screening with archiving to create a library of cryopreserved samples carrying known mutations. It also allows us to design efficient genotyping assays, making identifying F1 carriers straightforward.

CONCLUSIONS

The methods described will streamline the production of large numbers of knockout alleles in selected genes for phenotypic analysis, complementing existing efforts using random mutagenesis.

Zebrafish models for nemaline myopathy reveal a spectrum of nemaline bodies contributing to reduced muscle function

Nemaline myopathy is characterized by muscle weakness and the presence of rod-like (nemaline) bodies. The genetic etiology of nemaline myopathy is becoming increasingly understood with mutations in ten genes now known to cause the disease. Despite this, the mechanism by which skeletal muscle weakness occurs remains elusive, with previous studies showing no correlation between the frequency of nemaline bodies and disease severity. To investigate the formation of nemaline bodies and their role in pathogenesis, we generated overexpression and loss-of-function zebrafish models for skeletal muscle α-actin (ACTA1) and nebulin (NEB). We identify three distinct types of nemaline bodies and visualize their formation in vivo, demonstrating these nemaline bodies not only exhibit different subcellular origins, but also have distinct pathological consequences within the skeletal muscle. One subtype is highly dynamic and upon breakdown leads to the accumulation of cytoplasmic actin contributing to muscle weakness. Examination of a Neb-deficient model suggests this mechanism may be common in nemaline myopathy. Another subtype results from a reduction of actin and forms a more stable cytoplasmic body. In contrast, the final type originates at the Z-disk and is associated with myofibrillar disorganization. Analysis of zebrafish and muscle biopsies from ACTA1 nemaline myopathy patients demonstrates that nemaline bodies also possess a different protein signature. In addition, we show that the ACTA1^D286G mutation causes impaired actin incorporation and localization in the sarcomere. Together these data provide a novel examination of nemaline body origins and dynamics in vivo and identifies pathological changes that correlate with muscle weakness.

High-throughput and quantitative genome-wide messenger RNA sequencing for molecular phenotyping

BACKGROUND

We present a genome-wide messenger RNA (mRNA) sequencing technique that converts small amounts of RNA from many samples into molecular phenotypes. It encompasses all steps from sample preparation to sequence analysis and is applicable to baseline profiling or perturbation measurements.

RESULTS

Multiplex sequencing of transcript 3' ends identifies differential transcript abundance independent of gene annotation. We show that increasing biological replicate number while maintaining the total amount of sequencing identifies more differentially abundant transcripts.

CONCLUSIONS

This method can be implemented on polyadenylated RNA from any organism with an annotated reference genome and in any laboratory with access to Illumina sequencing.

A systematic genome-wide analysis of zebrafish protein-coding gene function

Since the publication of the human reference genome, the identities of specific genes associated with human diseases are being discovered at a rapid rate. A central problem is that the biological activity of these genes is often unclear. Detailed investigations in model vertebrate organisms, typically mice, have been essential for understanding the activities of many orthologues of these disease-associated genes. Although gene-targeting approaches and phenotype analysis have led to a detailed understanding of nearly 6,000 protein-coding genes, this number falls considerably short of the more than 22,000 mouse protein-coding genes. Similarly, in zebrafish genetics, one-by-one gene studies using positional cloning, insertional mutagenesis, antisense morpholino oligonucleotides, targeted re-sequencing, and zinc finger and TAL endonucleases have made substantial contributions to our understanding of the biological activity of vertebrate genes, but again the number of genes studied falls well short of the more than 26,000 zebrafish protein-coding genes. Importantly, for both mice and zebrafish, none of these strategies are particularly suited to the rapid generation of knockouts in thousands of genes and the assessment of their biological activity. Here we describe an active project that aims to identify and phenotype the disruptive mutations in every zebrafish protein-coding gene, using a well-annotated zebrafish reference genome sequence, high-throughput sequencing and efficient chemical mutagenesis. So far we have identified potentially disruptive mutations in more than 38% of all known zebrafish protein-coding genes. We have developed a multi-allelic phenotyping scheme to efficiently assess the effects of each allele during embryogenesis and have analysed the phenotypic consequences of over 1,000 alleles. All mutant alleles and data are available to the community and our phenotyping scheme is adaptable to phenotypic analysis beyond embryogenesis.