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

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

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

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

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, dicer1sa9205, which has a single point mutation induced by N-ethyl-N-nitrosourea mutagenesis. Heterozygous dicer1sa9205 developed normally, being phenotypically indistinguishable from wild-type siblings. Homozygous dicer1sa9205 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 dicer1sa9205 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

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)

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

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

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