Real-time resolution of point mutations that cause phenovariance in mice

A viable hypomorphic Arnt2 mutation causes hyperphagic obesity, diabetes and hepatic steatosis

Aryl hydrocarbon receptor nuclear translocator 2 (ARNT2) is a member of the basic helix-loop-helix/PER-ARNT-SIM (bHLH/PAS) transcription factor family. ARNT2 heterodimerizes with several members of the family, including single-minded homolog-1 (SIM1) and neuronal PAS domain protein 4 (NPAS4), primarily in neurons of the central nervous system. We screened 64,424 third-generation germline mutant mice derived from N-ethyl-N-nitrosourea (ENU)-mutagenized great-grandsires for weight abnormalities. Among 17 elevated body weight phenotypes identified and mapped, one strongly correlated with an induced missense mutation in Arnt2 using a semidominant model of inheritance. Causation was confirmed by CRISPR/Cas9 gene targeting to recapitulate the original ENU allele, specifying Arg74Cys (R74C). The CRISPR/Cas9-targeted (Arnt2^R74C/R74C) mice demonstrated hyperphagia and increased adiposity as well as hepatic steatosis and abnormalities in glucose homeostasis. The mutant ARNT2 protein showed decreased transcriptional activity when coexpressed with SIM1. These findings establish a requirement for ARNT2-dependent genes in the maintenance of the homeostatic feeding response, necessary for prevention of obesity and obesity-related diseases.

Summary: A forward genetic approach identified genes involved in maintaining normal body weight. A mutation in Arnt2 caused obesity by regulating food intake.

Excessive endosomal TLR signaling causes inflammatory disease in mice with defective SMCR8-WDR41-C9ORF72 complex function

The SMCR8-WDR41-C9ORF72 complex is a regulator of autophagy and lysosomal function. Autoimmunity and inflammatory disease have been ascribed to loss-of-function mutations of Smcr8 or C9orf72 in mice. In humans, autoimmunity has been reported to precede amyotrophic lateral sclerosis caused by mutations of C9ORF72 However, the cellular and molecular mechanisms underlying autoimmunity and inflammation caused by C9ORF72 or SMCR8 deficiencies remain unknown. Here, we show that splenomegaly, lymphadenopathy, and activated circulating T cells observed in Smcr8 ^-/- mice were rescued by triple knockout of the endosomal Toll-like receptors (TLRs) TLR3, TLR7, and TLR9. Myeloid cells from Smcr8 ^-/- mice produced excessive inflammatory cytokines in response to endocytosed TLR3, TLR7, or TLR9 ligands administered in the growth medium and in response to TLR2 or TLR4 ligands internalized by phagocytosis. These defects likely stem from prolonged TLR signaling caused by accumulation of LysoTracker-positive vesicles and by delayed phagosome maturation, both of which were observed in Smcr8 ^-/- macrophages. Smcr8 ^-/- mice also showed elevated susceptibility to dextran sodium sulfate-induced colitis, which was not associated with increased TLR3, TLR7, or TLR9 signaling. Deficiency of WDR41 phenocopied loss of SMCR8. Our findings provide evidence that excessive endosomal TLR signaling resulting from prolonged ligand-receptor contact causes inflammatory disease in SMCR8-deficient mice.

Predicting variant deleteriousness in non-human species: applying the CADD approach in mouse

Background

Predicting the deleteriousness of observed genomic variants has taken a step forward with the introduction of the Combined Annotation Dependent Depletion (CADD) approach, which trains a classifier on the wealth of available human genomic information. This raises the question whether it can be done with less data for non-human species. Here, we investigate the prerequisites to construct a CADD-based model for a non-human species.

Results

Performance of the mouse model is competitive with that of the human CADD model and better than established methods like PhastCons conservation scores and SIFT. Like in the human case, performance varies for different genomic regions and is best for coding regions. We also show the benefits of generating a species-specific model over lifting variants to a different species or applying a generic model. With fewer genomic annotations, performance on the test set as well as on the three validation sets is still good.

Conclusions

It is feasible to construct species-specific CADD models even when annotations such as epigenetic markers are not available. The minimal requirement for these models is the availability of a set of genomes of closely related species that can be used to infer an ancestor genome and substitution rates for the data generation.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2337-5) contains supplementary material, which is available to authorized users.

The class I myosin MYO1D binds to lipid and protects against colitis

Myosin ID (MYO1D) is a member of the class I myosin family. We screened 48,649 third generation (G3) germline mutant mice derived from N-ethyl-N-nitrosourea-mutagenized grandsires for intestinal homeostasis abnormalities after oral administration of dextran sodium sulfate (DSS). We found and validated mutations in Myo1d as a cause of increased susceptibility to DSS-induced colitis. MYO1D is produced in the intestinal epithelium, and the colitis phenotype is dependent on the nonhematopoietic compartment of the mouse. Moreover, MYO1D appears to couple cytoskeletal elements to lipid in an ATP-dependent manner. These findings demonstrate that MYO1D is needed to maintain epithelial integrity and protect against DSS-induced colitis.

Summary: Using random germline mutagenesis and screening of mice, we determined that loss of MYO1D function in nonhematopoietic tissues renders mice susceptible to colitis induced by dextran sodium sulfate challenge.

Whole exome sequencing of ENU-induced thrombosis modifier mutations in the mouse

Although the Factor V Leiden (FVL) gene variant is the most prevalent genetic risk factor for venous thrombosis, only 10% of FVL carriers will experience such an event in their lifetime. To identify potential FVL modifier genes contributing to this incomplete penetrance, we took advantage of a perinatal synthetic lethal thrombosis phenotype in mice homozygous for FVL (F5L/L) and haploinsufficient for tissue factor pathway inhibitor (Tfpi+/-) to perform a sensitized dominant ENU mutagenesis screen. Linkage analysis conducted in the 3 largest pedigrees generated from the surviving F5L/L Tfpi+/- mice ('rescues') using ENU-induced coding variants as genetic markers was unsuccessful in identifying major suppressor loci. Whole exome sequencing was applied to DNA from 107 rescue mice to identify candidate genes enriched for ENU mutations. A total of 3,481 potentially deleterious candidate ENU variants were identified in 2,984 genes. After correcting for gene size and multiple testing, Arl6ip5 was identified as the most enriched gene, though not reaching genome-wide significance. Evaluation of CRISPR/Cas9 induced loss of function in the top 6 genes failed to demonstrate a clear rescue phenotype. However, a maternally inherited (not ENU-induced) de novo mutation (Plcb4R335Q) exhibited significant co-segregation with the rescue phenotype (p = 0.003) in the corresponding pedigree. Thrombosis suppression by heterozygous Plcb4 loss of function was confirmed through analysis of an independent, CRISPR/Cas9-induced Plcb4 mutation (p = 0.01).

Large-scale forward genetics screening identifies Trpa1 as a chemosensor for predator odor-evoked innate fear behaviors

Innate behaviors are genetically encoded, but their underlying molecular mechanisms remain largely unknown. Predator odor 2,4,5-trimethyl-3-thiazoline (TMT) and its potent analog 2-methyl-2-thiazoline (2MT) are believed to activate specific odorant receptors to elicit innate fear/defensive behaviors in naive mice. Here, we conduct a large-scale recessive genetics screen of ethylnitrosourea (ENU)-mutagenized mice. We find that loss of Trpa1, a pungency/irritancy receptor, diminishes TMT/2MT and snake skin-evoked innate fear/defensive responses. Accordingly, Trpa1 -/- mice fail to effectively activate known fear/stress brain centers upon 2MT exposure, despite their apparent ability to smell and learn to fear 2MT. Moreover, Trpa1 acts as a chemosensor for 2MT/TMT and Trpa1-expressing trigeminal ganglion neurons contribute critically to 2MT-evoked freezing. Our results indicate that Trpa1-mediated nociception plays a crucial role in predator odor-evoked innate fear/defensive behaviors. The work establishes the first forward genetics screen to uncover the molecular mechanism of innate fear, a basic emotion and evolutionarily conserved survival mechanism.

Genetics of Resistant Hypertension: the Missing Heritability and Opportunities

PURPOSE OF THE REVIEW

Blood pressure regulation in humans has long been known to be a genetically determined trait. The identification of causal genetic modulators for this trait has been unfulfilling at the least. Despite the recent advances of genome-wide genetic studies, loci associated with hypertension or blood pressure still explain a very low percentage of the overall variation of blood pressure in the general population. This has precluded the translation of discoveries in the genetics of human hypertension to clinical use. Here, we propose the combined use of resistant hypertension as a trait for mapping genetic determinants in humans and the integration of new large-scale technologies to approach in model systems the multidimensional nature of the problem.

RECENT FINDINGS

New large-scale efforts in the genetic and genomic arenas are paving the way for an increased and granular understanding of genetic determinants of hypertension. New technologies for whole genome sequence and large-scale forward genetic screens can help prioritize gene and gene-pathways for downstream characterization and large-scale population studies, and guided pharmacological design can be used to drive discoveries to the translational application through better risk stratification and new therapeutic approaches. Although significant challenges remain in the mapping and identification of genetic determinants of hypertension, new large-scale technological approaches have been proposed to surpass some of the shortcomings that have limited progress in the area for the last three decades. The incorporation of these technologies to hypertension research may significantly help in the understanding of inter-individual blood pressure variation and the deployment of new phenotyping and treatment approaches for the condition.

Probability of phenotypically detectable protein damage by ENU-induced mutations in the Mutagenetix database

Computational inference of mutation effects is necessary for genetic studies in which many mutations must be considered as etiologic candidates. Programs such as PolyPhen-2 predict the relative severity of damage caused by missense mutations, but not the actual probability that a mutation will reduce/eliminate protein function. Based on genotype and phenotype data for 116,330 ENU-induced mutations in the Mutagenetix database, we calculate that putative null mutations, and PolyPhen-2-classified “probably damaging”, “possibly damaging”, or “probably benign” mutations have, respectively, 61%, 17%, 9.8%, and 4.5% probabilities of causing phenotypically detectable damage in the homozygous state. We use these probabilities in the estimation of genome saturation and the probability that individual proteins have been adequately tested for function in specific genetic screens. We estimate the proportion of essential autosomal genes in Mus musculus (C57BL/6J) and show that viable mutations in essential genes are more likely to induce phenotype than mutations in non-essential genes.

Programs such as PolyPhen-2 predict the relative severity of damage by missense mutations. Here, Wang et al estimate probabilities that putative null or missense alleles would reduce protein function to cause detectable phenotype by analyzing data from ENU-induced mouse mutations.

Skin-specific regulation of SREBP processing and lipid biosynthesis by glycerol kinase 5

The recessive N-ethyl-N-nitrosourea-induced phenotype toku is characterized by delayed hair growth, progressive hair loss, and excessive accumulation of dermal cholesterol, triglycerides, and ceramides. The toku phenotype was attributed to a null allele of Gk5, encoding glycerol kinase 5 (GK5), a skin-specific kinase expressed predominantly in sebaceous glands. GK5 formed a complex with the sterol regulatory element-binding proteins (SREBPs) through their C-terminal regulatory domains, inhibiting SREBP processing and activation. In Gk5^toku/toku mice, transcriptionally active SREBPs accumulated in the skin, but not in the liver; they were localized to the nucleus and led to elevated lipid synthesis and subsequent hair growth defects. Similar defective hair growth was observed in kinase-inactive GK5 mutant mice. Hair growth defects of homozygous toku mice were partially rescued by treatment with the HMG-CoA reductase inhibitor simvastatin. GK5 exists as part of a skin-specific regulatory mechanism for cholesterol biosynthesis, independent of cholesterol regulation elsewhere in the body.

A mutagenesis-derived Lrp5 mouse mutant with abnormal retinal vasculature and low bone mineral density

PURPOSE

Familial exudative vitreoretinopathy (FEVR) is caused by mutations in the genes encoding low-density lipoprotein receptor-related protein (LRP5) or its interacting partners, namely frizzled class receptor 4 (FZD4) and norrin cystine knot growth factor (NDP). Mouse models for Lrp5, Fzd4, and Ndp have proven to be important for understanding the retinal pathophysiology underlying FEVR and systemic abnormalities related to defective Wnt signaling. Here, we report a new mouse mutant, tvrm111B, which was identified by electroretinogram (ERG) screening of mice generated in the Jackson Laboratory Translational Vision Research Models (TVRM) mutagenesis program.

METHODS

ERGs were used to examine outer retinal physiology. The retinal vasculature was examined by in vivo retinal imaging, as well as by histology and immunohistochemistry. The tvrm111B locus was identified by genetic mapping of mice generated in a cross to DBA/2J, and subsequent sequencing analysis. Gene expression was examined by real-time PCR of retinal RNA. Bone mineral density (BMD) was examined by peripheral dual-energy X-ray absorptiometry.

RESULTS

The tvrm111B allele is inherited as an autosomal recessive trait. Genetic mapping of the decreased ERG b-wave phenotype of tvrm111B mice localized the mutation to a region on chromosome 19 that included Lrp5. Sequencing of Lrp5 identified the insertion of a cytosine (c.4724_4725insC), which is predicted to cause a frameshift that disrupts the last three of five conserved PPPSPxS motifs in the cytoplasmic domain of LRP5, culminating in a premature termination. In addition to a reduced ERG b-wave, Lrp5tvrm111B homozygotes have low BMD and abnormal features of the retinal vasculature that have been reported previously in Lrp5 mutant mice, including persistent hyaloid vessels, leakage on fluorescein angiography, and an absence of the deep retinal capillary bed.

CONCLUSIONS

The phenotype of the Lrp5tvrm111B mutant includes abnormalities of the retinal vasculature and of BMD. This model may be a useful resource to further our understanding of the biological role of LRP5 and to evaluate experimental therapies for FEVR or other conditions associated with LRP5 dysfunction.

Creatine maintains intestinal homeostasis and protects against colitis

Creatine, a nitrogenous organic acid, replenishes cytoplasmic ATP at the expense of mitochondrial ATP via the phosphocreatine shuttle. Creatine levels are maintained by diet and endogenous synthesis from arginine and glycine. Glycine amidinotransferase (GATM) catalyzes the rate-limiting step of creatine biosynthesis: the transfer of an amidino group from arginine to glycine to form ornithine and guanidinoacetate. We screened 36,530 third-generation germline mutant mice derived from N-ethyl-N-nitrosourea-mutagenized grandsires for intestinal homeostasis abnormalities after oral administration of dextran sodium sulfate (DSS). Among 27 colitis susceptibility phenotypes identified and mapped, one was strongly correlated with a missense mutation in Gatm in a recessive model of inheritance, and causation was confirmed by CRISPR/Cas9 gene targeting. Supplementation of homozygous Gatm mutants with exogenous creatine ameliorated the colitis phenotype. CRISPR/Cas9-targeted (Gatm^c/c ) mice displayed a normal peripheral immune response and immune cell homeostasis. However, the intestinal epithelium of the Gatm^c/c mice displayed increased cell death and decreased proliferation during DSS treatment. In addition, Gatm^c/c colonocytes showed increased metabolic stress in response to DSS with higher levels of phospho-AMPK and lower levels of phosphorylation of mammalian target of rapamycin (phospho-mTOR). These findings establish an in vivo requirement for rapid replenishment of cytoplasmic ATP within colonic epithelial cells in the maintenance of the mucosal barrier after injury.

IgD class switching is initiated by microbiota and limited to mucosa-associated lymphoid tissue in mice

Class-switch recombination (CSR) alters the Ig isotype to diversify antibody effector functions. IgD CSR is a rare event, and its regulation is poorly understood. We report that deficiency of 53BP1, a DNA damage-response protein, caused age-dependent overproduction of secreted IgD resulting from increased IgD CSR exclusively within B cells of mucosa-associated lymphoid tissues. IgD overproduction was dependent on activation-induced cytidine deaminase, hematopoietic MyD88 expression, and an intact microbiome, against which circulating IgD, but not IgM, was reactive. IgD CSR occurred via both alternative nonhomologous end-joining and homologous recombination pathways. Microbiota-dependent IgD CSR also was detected in nasal-associated lymphoid tissue of WT mice. These results identify a pathway, present in WT mice and hyperactivated in 53BP1-deficient mice, by which microbiota signal via Toll-like receptors to elicit IgD CSR.

Innate immunity and the new forward genetics

As it is a hard-wired system for responses to microbes, innate immunity is particularly susceptible to classical genetic analysis. Mutations led the way to the discovery of many of the molecular elements of innate immune sensing and signaling pathways. In turn, the need for a faster way to find the molecular causes of mutation-induced phenotypes triggered a huge transformation in forward genetics. During the 1980s and 1990s, many heritable phenotypes were ascribed to mutations through positional cloning. In mice, this required three steps. First, a genetic mapping step was used to show that a given phenotype emanated from a circumscribed region of the genome. Second, a physical mapping step was undertaken, in which all of the region was cloned and its gene content determined. Finally, a concerted search for the mutation was performed. Such projects usually lasted for several years, but could produce breakthroughs in our understanding of biological processes. Publication of the annotated mouse genome sequence in 2002 made physical mapping unnecessary. More recently we devised a new technology for automated genetic mapping, which eliminated both genetic mapping and the search for mutations among candidate genes. The cause of phenotype can now be determined instantaneously. We have created more than 100,000 coding/splicing mutations. And by screening for defects of innate and adaptive immunity we have discovered many "new" proteins needed for innate immune function.

Insulin resistance and diabetes caused by genetic or diet-induced KBTBD2 deficiency in mice

We describe a metabolic disorder characterized by lipodystrophy, hepatic steatosis, insulin resistance, severe diabetes, and growth retardation observed in mice carrying N-ethyl-N-nitrosourea (ENU)-induced mutations. The disorder was ascribed to a mutation of kelch repeat and BTB (POZ) domain containing 2 (Kbtbd2) and was mimicked by a CRISPR/Cas9-targeted null allele of the same gene. Kbtbd2 encodes a BTB-Kelch family substrate recognition subunit of the Cullin-3-based E3 ubiquitin ligase. KBTBD2 targeted p85α, the regulatory subunit of the phosphoinositol-3-kinase (PI3K) heterodimer, causing p85α ubiquitination and proteasome-mediated degradation. In the absence of KBTBD2, p85α accumulated to 30-fold greater levels than in wild-type adipocytes, and excessive p110-free p85α blocked the binding of p85α-p110 heterodimers to IRS1, interrupting the insulin signal. Both transplantation of wild-type adipose tissue and homozygous germ line inactivation of the p85α-encoding gene Pik3r1 rescued diabetes and hepatic steatosis phenotypes of Kbtbd2^-/- mice. Kbtbd2 was down-regulated in diet-induced obese insulin-resistant mice in a leptin-dependent manner. KBTBD2 is an essential regulator of the insulin-signaling pathway, modulating insulin sensitivity by limiting p85α abundance.

Identification of mutations through dominant screening for obesity using C57BL/6 substrains

The discovery of leptin substantiated the usefulness of a forward genetic approach in elucidating the molecular network regulating energy metabolism. However, no successful dominant screening for obesity has been reported, which may be due to the influence of quantitative trait loci between the screening and counter strains and the low fertility of obese mice. Here, we performed a dominant screening for obesity using C57BL/6 substrains, C57BL/6J and C57BL/6N, with the routine use of in vitro fertilization. The screening of more than 5000 mutagenized mice established two obese pedigrees in which single nucleotide substitutions in Mc4r and Sim1 genes were identified through whole-exome sequencing. The mutation in the Mc4r gene produces a premature stop codon, and the mutant SIM1 protein lacks transcriptional activity, showing that the haploinsufficiency of SIM1 and MC4R results in obesity. We further examined the hypothalamic neuropeptide expressions in the mutant pedigrees and mice with diet-induced obesity, which showed that each obesity mouse model has distinct neuropeptide expression profiles. This forward genetic screening scheme is useful and applicable to any research field in which mouse models work.

Novel gene function revealed by mouse mutagenesis screens for models of age-related disease

Determining the genetic bases of age-related disease remains a major challenge requiring a spectrum of approaches from human and clinical genetics to the utilization of model organism studies. Here we report a large-scale genetic screen in mice employing a phenotype-driven discovery platform to identify mutations resulting in age-related disease, both late-onset and progressive. We have utilized N-ethyl-N-nitrosourea mutagenesis to generate pedigrees of mutagenized mice that were subject to recurrent screens for mutant phenotypes as the mice aged. In total, we identify 105 distinct mutant lines from 157 pedigrees analysed, out of which 27 are late-onset phenotypes across a range of physiological systems. Using whole-genome sequencing we uncover the underlying genes for 44 of these mutant phenotypes, including 12 late-onset phenotypes. These genes reveal a number of novel pathways involved with age-related disease. We illustrate our findings by the recovery and characterization of a novel mouse model of age-related hearing loss.

What Lies Ahead?

Viral pathogenesis is a field in rapid evolution, reflecting the dynamic development of systems biology and the continuing introduction of new or improved methodologies. Therefore, this final chapter is dedicated to “futurism,” a look at what lies ahead for this field. We have recruited a number of scientists to write short pieces where they are free to speculate on future developments in their respective areas of expertise.

NLRP3 activation and mitosis are mutually exclusive events coordinated by NEK7, a new inflammasome component

The NLRP3 inflammasome responds to microbes and danger signals by processing and activating proinflammatory cytokines, including interleukin 1β (IL-1β) and IL-18. We found here that activation of the NLRP3 inflammasome was restricted to interphase of the cell cycle by NEK7, a serine-threonine kinase previously linked to mitosis. Activation of the NLRP3 inflammasome required NEK7, which bound to the leucine-rich repeat domain of NLRP3 in a kinase-independent manner downstream of the induction of mitochondrial reactive oxygen species (ROS). This interaction was necessary for the formation of a complex containing NLRP3 and the adaptor ASC, oligomerization of ASC and activation of caspase-1. NEK7 promoted the NLRP3-dependent cellular inflammatory response to intraperitoneal challenge with monosodium urate and the development of experimental autoimmune encephalitis in mice. Our findings suggest that NEK7 serves as a cellular switch that enforces mutual exclusivity of the inflammasome response and cell division.

An ultra-dense library resource for rapid deconvolution of mutations that cause phenotypes in Escherichia coli

With the wide availability of whole-genome sequencing (WGS), genetic mapping has become the rate-limiting step, inhibiting unbiased forward genetics in even the most tractable model organisms. We introduce a rapid deconvolution resource and method for untagged causative mutations after mutagenesis, screens, and WGS in Escherichia coli. We created Deconvoluter—ordered libraries with selectable insertions every 50 kb in the E. coli genome. The Deconvoluter method uses these for replacement of untagged mutations in the genome using a phage-P1-based gene-replacement strategy. We validate the Deconvoluter resource by deconvolution of 17 of 17 phenotype-altering mutations from a screen of N-ethyl-N-nitrosourea-induced mutants. The Deconvoluter resource permits rapid unbiased screens and gene/function identification and will enable exploration of functions of essential genes and undiscovered genes/sites/alleles not represented in existing deletion collections. This resource for unbiased forward-genetic screens with mapping-by-sequencing (‘forward genomics’) demonstrates a strategy that could similarly enable rapid screens in many other microbes.

Variant mapping and mutation discovery in inbred mice using next-generation sequencing

Background

The development of powerful new methods for DNA sequencing enable the discovery of sequence variants, their utilization for the mapping of mutant loci, and the identification of causal variants in a single step. We have applied this approach for the analysis of ENU-mutagenized mice maintained on an inbred background.

Results

We ascertained ENU-induced variants in four different phenotypically mutant lines. These were then used as informative markers for positional cloning of the mutated genes. We tested both whole genome (WGS) and whole exome (WES) datasets.

Conclusion

Both approaches were successful as a means to localize a region of homozygosity, as well as identifying mutations of candidate genes, which could be individually assessed. As expected, the WGS strategy was more reliable, since many more ENU-induced variants were ascertained.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2173-1) contains supplementary material, which is available to authorized users.

Current strategies for mutation detection in phenotype-driven screens utilising next generation sequencing

Mutagenesis-based screens in mice are a powerful discovery platform to identify novel genes or gene functions associated with disease phenotypes. An N-ethyl-N-nitrosourea (ENU) mutagenesis screen induces single nucleotide variants randomly in the mouse genome. Subsequent phenotyping of mutant and wildtype mice enables the identification of mutated pathways resulting in phenotypes associated with a particular ENU lesion. This unbiased approach to gene discovery conducts the phenotyping with no prior knowledge of the functional mutations. Before the advent of affordable next generation sequencing (NGS), ENU variant identification was a limiting step in gene characterization, akin to ‘finding a needle in a haystack’. The emergence of a reliable reference genome alongside advances in NGS has propelled ENU mutation discovery from an arduous, time-consuming exercise to an effective and rapid form of mutation discovery. This has permitted large mouse facilities worldwide to use ENU for novel mutation discovery in a high-throughput manner, helping to accelerate basic science at the mechanistic level. Here, we describe three different strategies used to identify ENU variants from NGS data and some of the subsequent steps for mutation characterisation.