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 to study 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 hours post-infection). Consistent with its clearance from the host, the signature of late Shigella infection (24 hours 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

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

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.This project has received funding from the European Union’s Horizon 2020 432 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 641458. The 433 work carried out at the University of Edinburgh was partly funded by EEP, MRC HGU Programme 434 (MC_UU_00007/9), European Research Council (ZF-MEL-CHEMBIO-648489), and L'Oreal-Melanoma 435 Research Alliance (401181)

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.This project has received funding from the European Union’s Horizon 2020 432 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 641458. The 433 work carried out at the University of Edinburgh was partly funded by EEP, MRC HGU Programme 434 (MC_UU_00007/9), European Research Council (ZF-MEL-CHEMBIO-648489), and L'Oreal-Melanoma 435 Research Alliance (401181)

Loss of the chromatin modifier Kdm2aa causes BrafV600E-independent spontaneous melanoma in zebrafish.

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

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

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 unexposed and exposed to MnCl2. Differentially expressed genes mapped to the central nervous system and eye, and pathway analysis suggested that calcium 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 calcium 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 calcium dyshomeostasis and cellular stress

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

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

Large-scale mapping of mutations affecting zebrafish development

BACKGROUND: Large-scale mutagenesis screens in the zebrafish employing the mutagen ENU have isolated several hundred mutant loci that represent putative developmental control genes. In order to realize the potential of such screens, systematic genetic mapping of the mutations is necessary. Here we report on a large-scale effort to map the mutations generated in mutagenesis screening at the Max Planck Institute for Developmental Biology by genome scanning with microsatellite markers. RESULTS: We have selected a set of microsatellite markers and developed methods and scoring criteria suitable for efficient, high-throughput genome scanning. We have used these methods to successfully obtain a rough map position for 319 mutant loci from the Tübingen I mutagenesis screen and subsequent screening of the mutant collection. For 277 of these the corresponding gene is not yet identified. Mapping was successful for 80 % of the tested loci. By comparing 21 mutation and gene positions of cloned mutations we have validated the correctness of our linkage group assignments and estimated the standard error of our map positions to be approximately 6 cM. CONCLUSION: By obtaining rough map positions for over 300 zebrafish loci with developmental phenotypes, we have generated a dataset that will be useful not only for cloning of the affected genes, but also to suggest allelism of mutations with similar phenotypes that will be identified in future screens. Furthermore this work validates the usefulness of our methodology for rapid, systematic and inexpensive microsatellite mapping of zebrafish mutations