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.

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.

Placentation defects are highly prevalent in embryonic lethal mouse mutants.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Multi-allelic phenotyping--a systematic approach for the simultaneous analysis of multiple induced mutations.

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.

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.

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.

Zebrafish null mutants of Sept6 and Sept15 are viable but more susceptible to Shigella infection.

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.

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.

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.

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.