Search

In zebrafish, maternally deposited yolk is the source of nutrients for embryogenesis prior to digestive system maturation. Yolk nutrients are processed and secreted to the growing organism by an extra-embryonic tissue, the yolk syncytial layer (YSL). The export of lipids from the YSL occurs through the production of triacylglycerol-rich lipoproteins. Here we report that mutations in the triacylglycerol synthesis enzyme, diacylglycerol acyltransferase-2 (Dgat2), cause yolk sac opacity due to aberrant accumulation of cytoplasmic lipid droplets in the YSL. Although triacylglycerol synthesis continues, it is not properly coupled to lipoprotein production as dgat2 mutants produce fewer, smaller, ApoB-containing lipoproteins. Unlike DGAT2-null mice, which are lipopenic and die soon after birth, zebrafish dgat2 mutants are viable, fertile, and exhibit normal mass and adiposity. Residual Dgat activity cannot be explained by the activity of other known Dgat isoenzymes, as dgat1a;dgat1b;dgat2 triple mutants continue to produce YSL lipid droplets and remain viable as adults. Further, the newly identified diacylglycerol acyltransferase, Tmem68, is also not responsible for the residual triacylglycerol synthesis activity. Unlike overexpression of Dgat1a and Dgat1b, monoacylglycerol acyltransferase-3 (Mogat3b) overexpression does not rescue yolk opacity, suggesting it does not possess Dgat activity in the YSL. However, mogat3b;dgat2 double mutants exhibit increased yolk opacity and often have structural alterations of the yolk extension. Quadruple mogat3b;dgat1a;dgat1b;dgat2 mutants either have severely reduced viability and stunted growth or do not survive past 3 days post fertilization, depending on the dgat2 mutant allele present. Our study highlights the remarkable ability of vertebrates to synthesize triacylglycerol through multiple biosynthetic pathways., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2024. Published by The Company of Biologists Ltd.)

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., (© 2023. The Author(s).)

Septins are evolutionarily conserved GTP-binding proteins known for their roles in cell division and host defence against Shigella infection. Although septin group members are viewed to function as hetero-oligomeric complexes, the role of individual septins within these complexes or in isolation is poorly understood. Decades of work using mouse models has shown that some septins (including SEPT7) are essential for animal development, while others (including SEPT6) are dispensable, suggesting that some septins may compensate for the absence of others. The zebrafish genome encodes 19 septin genes, representing the full complement of septin groups described in mice and humans. In this report, we characterise null mutants for zebrafish Sept6 (a member of the SEPT6 group) and Sept15 (a member of the SEPT7 group) and test their role in zebrafish development and host defence against Shigella infection. We show that null mutants for Sept6 and Sept15 are both viable, and that expression of other zebrafish septins are not significantly affected by their mutation. Consistent with previous reports using knockdown of Sept2, Sept7b, and Sept15, we show that Sept6 and Sept15 are required for host defence against Shigella infection. These results highlight Shigella infection of zebrafish as a powerful system to study the role of individual septins in vivo., (© 2023 The Authors. Cytoskeleton published by Wiley Periodicals LLC.)

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.

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.

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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., Competing Interests: Declaration of interests L.R. and E.L.P. are advisory board members for Cell. E.L.P. is a scientific advisory board member of ImmunoMet and a founder of Rheos Medicines. For the purpose of open access, the authors have applied for a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission. This work is licensed under a Creative Commons Attribution 4.0 International License., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2022. Published by The Company of Biologists Ltd.)

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., (Copyright © 2022. Published by Elsevier Ltd.)

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., Competing Interests: RW, EM, SW, EB No competing interests declared, (© 2022, White et al.)

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 ( ankk1 27 ins ) 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 ankk1 27 ins/ 27 ins fish. Quantitative polymerase chain reaction analysis revealed an increase in expression of drd2b mRNA in ankk1 27 ins 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 ankk1 27 ins/ 27 ins 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. Ankk1 27 ins 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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Leggieri, García-González, Torres-Perez, Havelange, Hosseinian, Mech, Keatinge, Busch-Nentwich and Brennan.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Mech, Merteroglu, Sealy, Teh, White, Havelange, Brennan and Busch-Nentwich.)

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., Competing Interests: Competing interestsThe authors declare no competing interests., (Copyright © The Authors; exclusive licensee Bio-protocol LLC.)

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., (© 2020 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

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., Competing Interests: Declaration of Interests L.I.Z. is a founder and stockholder of Fate Therapeutics, CAMP4 Therapeutics, and Scholar Rock and a consultant for Celularity., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: The authors have declared that no competing interests exist.

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 F 3 sibling groups, where each was derived from an individual ethyl-nitrosurea mutagenized F 0 fish, two showed increased or decreased nicotine preference. Out of 25 inactivating mutations carried by the F 3 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., Competing Interests: JG, AB, MP, RR, DJ, AS, MT, EB, DS, AM, JK, TP, VK, RW, CB No competing interests declared, (© 2020, García-González et al.)

Immune cells congregate at specific loci to fight infections during inflammatory responses, a process that must be transient and self-resolving. Cell dispersal promotes resolution, but it remains unclear how transition from clustering to dispersal is regulated. Here we show, using quantitative live imaging in zebrafish, that differential ligand-induced trafficking of chemokine receptors such as Cxcr1 and Cxcr2 orchestrates the state of neutrophil congregation at sites of tissue damage. Through receptor mutagenesis and biosensors, we show that Cxcr1 promotes clustering at wound sites, but is promptly desensitized and internalized, which prevents excess congregation. By contrast, Cxcr2 promotes bidirectional motility and is sustained at the plasma membrane. Persistent plasma membrane residence of Cxcr2 prolongs downstream signaling and is required for sustained exploratory motion conducive to dispersal. Thus, differential trafficking of two chemokine receptors allows coordination of antagonistic cell behaviors, promoting a self-resolving migratory response.

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 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., Competing Interests: The authors have declared that no competing interests exist.

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., (Copyright © 2019 Elsevier B.V. All rights reserved.)

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., Competing Interests: RY, TH, FC, HS, QS, LL, CW, BC, JL, EA, AK, IS, JR, CT, IB, MA, EB, GG, SW No competing interests declared, (© 2019, Young et al.)

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., Competing Interests: ML, DD, NK, SB, HB, MS, TG, AR, EB, AP, KH, AP, JH, TP No competing interests declared, (© 2019, Lush et al.)

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., (Copyright © 2018 Minchin et al.)

Large-scale phenotyping efforts have demonstrated that approximately 25-30% of mouse gene knockouts cause intrauterine lethality. Analysis of these mutants has largely focused on the embryo and not the placenta, despite the crucial role of this extraembryonic organ for developmental progression. Here we screened 103 embryonic lethal and sub-viable mouse knockout lines from the Deciphering the Mechanisms of Developmental Disorders program for placental phenotypes. We found that 68% of knockout lines that are lethal at or after mid-gestation exhibited placental dysmorphologies. Early lethality (embryonic days 9.5-14.5) is almost always associated with severe placental malformations. Placental defects correlate strongly with abnormal brain, heart and vascular development. Analysis of mutant trophoblast stem cells and conditional knockouts suggests that a considerable number of factors that cause embryonic lethality when ablated have primary gene function in trophoblast cells. Our data highlight the hugely under-appreciated importance of placental defects in contributing to abnormal embryo development and suggest key molecular nodes that govern placenta formation.

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.

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.

Large-scale forward genetic screens have been instrumental for identifying genes that regulate development, homeostasis, and regeneration, as well as the mechanisms of disease. The zebrafish, Danio rerio , is an established genetic and developmental model used in genetic screens to uncover genes necessary for early development. However, the regulation of postembryonic development has received less attention as these screens are more labor intensive and require extensive resources. The lack of systematic interrogation of late development leaves large aspects of the genetic regulation of adult form and physiology unresolved. To understand the genetic control of postembryonic development, we performed a dominant screen for phenotypes affecting the adult zebrafish. In our screen, we identified 72 adult viable mutants showing changes in the shape of the skeleton as well as defects in pigmentation. For efficient mapping of these mutants and mutation identification, we devised a new mapping strategy based on identification of mutant-specific haplotypes. Using this method in combination with a candidate gene approach, we were able to identify linked mutations for 22 out of 25 mutants analyzed. Broadly, our mutational analysis suggests that there are key genes and pathways associated with late development. Many of these pathways are shared with humans and are affected in various disease conditions, suggesting constraint in the genetic pathways that can lead to change in adult form. Taken together, these results show that dominant screens are a feasible and productive means to identify mutations that can further our understanding of gene function during postembryonic development and in disease., (Copyright © 2017 by the Genetics Society of America.)

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., Competing Interests: Competing interests: No competing interests were disclosed.

KDM2A is a histone demethylase associated with transcriptional silencing, however very little is known about its in vivo role in development and disease. Here we demonstrate that loss of the orthologue kdm2aa in zebrafish causes widespread transcriptional disruption and leads to spontaneous melanomas at a high frequency. Fish homozygous for two independent premature stop codon alleles show reduced growth and survival, a strong male sex bias, and homozygous females exhibit a progressive oogenesis defect. kdm2aa mutant fish also develop melanomas from early adulthood onwards which are independent from mutations in braf and other common oncogenes and tumour suppressors as revealed by deep whole exome sequencing. In addition to effects on translation and DNA replication gene expression, high-replicate RNA-seq in morphologically normal individuals demonstrates a stable regulatory response of epigenetic modifiers and the specific de-repression of a group of zinc finger genes residing in constitutive heterochromatin. Together our data reveal a complex role for Kdm2aa in regulating normal mRNA levels and carcinogenesis. These findings establish kdm2aa mutants as the first single gene knockout model of melanoma biology.

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.

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.

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.

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.

The zebrafish mutation project (ZMP) aims to generate a loss of function allele for every protein-coding gene, but importantly to also characterise the phenotypes of these alleles during the first five days of development. Such a large-scale screen requires a systematic approach both to identifying phenotypes, and also to linking those phenotypes to specific mutations. This phenotyping pipeline simultaneously assesses the consequences of multiple alleles in a two-step process. First, mutations that do not produce a visible phenotype during the first five days of development are identified, while a second round of phenotyping focuses on detailed analysis of those alleles that are suspected to cause a phenotype. Allele-specific PCR single nucleotide polymorphism (SNP) assays are used to genotype F2 parents and individual F3 fry for mutations known to be present in the F1 founder. With this method specific phenotypes can be linked to induced mutations. In addition a method is described for cryopreserving sperm samples of mutagenised males and their subsequent use for in vitro fertilisation to generate F2 families for phenotyping. Ultimately this approach will lead to the functional annotation of the zebrafish genome, which will deepen our understanding of gene function in development and disease., (Copyright © 2013 The Authors. Published by Elsevier Inc. All rights reserved.)

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.