Expression of miR-34a is a sensitive biomarker for exposure to genotoxic agents in human lymphoblastoid TK6 cells.

miR-34a has been identified as a tumor suppressor microRNA (miRNA) involved in the P53 network. Its expression levels correlate to carcinogenesis, which are generally lower in tumor tissue and higher in response to DNA damage. In this study, the response of miR-34a from exposure to genotoxic agents in human lymphoblastoid TK6 cells was evaluated to assess whether the expression of this miRNA could be used as an early indicator for genotoxic damage in mammalian cells. TK6 cells were treated with seven genotoxic agents with different mode-of-actions (cisplatin, N-ethyl-N-nitrosourea, etoposide, mitomycin C, methyl methanesulphonate, taxol, and X-ray radiation) and a non-genetic toxin (usnic acid) at different concentrations for four hours (except for X-rays) and the expression levels of miR-34a were measured 24 h after the beginning of the treatments. The expression levels of miR-34a were significantly increased by these genetic toxins in a dose-dependent manner, while no significant change in miRNA expression was found in the usnic acid-treated cells. These results suggest that miR-34a can respond to genotoxic insults sensitively; thus, miR-34a expression has the potential to be used to evaluate genotoxicity of agents.

Mutational Analysis of N-Ethyl-N-Nitrosourea (ENU) in the Fission Yeast Schizosaccharomyces pombe.

Forward genetics in model organisms has boosted our knowledge of the genetic bases of development, aging, and human diseases. In this experimental pipeline, it is crucial to start by inducing a large number of random mutations in the genome of the model organism to search for phenotypes of interest. Many chemical mutagens are used to this end because most of them display particular reactivity properties and act differently over DNA. Here we report the use of N-ethyl-N-nitrosourea (ENU) as a mutagen in the fission yeast Schizosaccharomyces pombe As opposed to many other alkylating agents, ENU only induces an S N 1-type reaction with a low s constant (s = 0.26), attacking preferentially O2 and O4 in thymine and O6 deoxyguanosine, leading to base substitutions rather than indels, which are extremely rare in its resulting mutagenic repertoire. Using ENU, we gathered a collection of 13 temperature-sensitive mutants and 80 auxotrophic mutants including two deleterious alleles of the human ortholog ATIC. Defective alleles of this gene cause AICA-ribosiduria, a severe genetic disease. In this screen, we also identified 13 aminoglycoside-resistance inactivating mutations in APH genes. Mutations reported here may be of interest for metabolism related diseases and antibiotic resistance research fields.

IL-15 and GM-CSF stimulated macrophages enhances phagocytic activity in ENU induced leukemic mice.

Murine splenic macrophage plays a decisive role in host immunity through phagocytosis against pathogens. It was reported that, macrophages also involves in phagocytosis of some tumour cells upon its activation initiated by certain cytokines produced by other immune cell or by indigenously treated. In this study, we have investigated the killing of leukemic blast cells by macrophages upon stimulated with IL-15 and GM-CSF alone or in combination in ENU challenged leukemic murine model. Along with, the release of TNF-α, IL-12 and IFN-γ by macrophages were assayed by ELISA. NO production by macrophages was also investigated. The molecular expressions like GM-CSF and TLRs were investigated for better understand of macrophage-leukemic cell interaction. Result shows that in disease condition macrophages have poor phagocytic activities which may be due to less release of TNF-α, IL-12 and IFN-γ by macrophages. This impaired phagocytic activity in leukemic mice was increase upon stimulation with IL-15 and GM-CSF.

EGFR-Mutated Lung Cancers Resistant to Osimertinib through EGFR C797S Respond to First-Generation Reversible EGFR Inhibitors but Eventually Acquire EGFR T790M/C797S in Preclinical Models and Clinical Samples.

INTRODUCTION: Osimertinib is approved for advanced EGFR-mutated NSCLC, and identification of on-target mechanisms of resistance (i.e., EGFR C797S) to this third-generation EGFR inhibitor are evolving. Whether durable control of subsequently osimertinib-resistant NSCLC with the EGFR-sensitizing mutation (SM)/C797S is possible with first-generation EGFR inhibitors (such as gefitinib or erlotinib) remains underreported, as does the resultant acquired resistance profile. METHODS: We used N-ethyl-N-nitrosourea mutagenesis to determine the profile of EGFR SM/C797S preclinical models exposed to reversible EGFR inhibitors. In addition, we retrospectively probed a case of EGFR SM lung adenocarcinoma treated with first-line osimertinib, followed by second-line erlotinib in the setting of EGFR SM/C797S. RESULTS: Use of N-ethyl-N-nitrosourea mutagenesis against the background of EGFR L858R/C797S in conjunction with administration of gefitinib revealed preferential outgrowth of cells with EGFR L858R/T790M/C797S. A patient with EGFR delE746_T751insV NSCLC was treated with osimertinib with sustained response for 10 months before acquiring EGFR C797S. The patient was subsequently treated with erlotinib, with response for a period of 4 months, but disease progression ensued. Liquid biopsy disclosed EGFR delE746_T751insV with T790M and C797S present in cis. CONCLUSION: EGFR SM NSCLC can acquire resistance to osimertinib through development of the EGFR C797S mutation. In this clinical scenario, the tumor may respond transiently to reversible first-generation EGFR inhibitors (gefitinib or erlotinib), but evolving mechanisms of on-target resistance-in clinical specimens and preclinical systems-indicate that EGFR C797S along with EGFR T790M can evolve. This report adds to the growing understanding of tumor evolution or adaptability to sequential EGFR inhibition and augments support for exploring combination therapies to delay or prevent on-target resistance.

An essential role for α4A-tubulin in platelet biogenesis.

During platelet biogenesis, microtubules (MTs) are arranged into submembranous structures (the marginal band) that encircle the cell in a single plane. This unique MT array has no equivalent in any other mammalian cell, and the mechanisms responsible for this particular mode of assembly are not fully understood. One possibility is that platelet MTs are composed of a particular set of tubulin isotypes that carry specific posttranslational modifications. Although ß1-tubulin is known to be essential, no equivalent roles of α-tubulin isotypes in platelet formation or function have so far been reported. Here, we identify α4A-tubulin as a predominant α-tubulin isotype in platelets. Similar to ß1-tubulin, α4A-tubulin expression is up-regulated during the late stages of megakaryocyte differentiation. Missense mutations in the α4A-tubulin gene cause macrothrombocytopenia in mice and humans. Defects in α4A-tubulin lead to changes in tubulin tyrosination status of the platelet tubulin pool. Ultrastructural defects include reduced numbers and misarranged MT coils in the platelet marginal band. We further observed defects in megakaryocyte maturation and proplatelet formation in Tuba4a-mutant mice. We have, thus, discovered an α-tubulin isotype with specific and essential roles in platelet biogenesis.

Oxidant stress induction and signalling in xenografted (human breast cancer-tissues) plus estradiol treated or N-ethyl-N-nitrosourea treated female rats via altered estrogen sulfotransferase (rSULT1E1) expressions and SOD1/catalase regulations.

N-ethyl-N-nitrosourea (ENU) is highly used in rodent models of tumerogenesis/carcinogenesis. Xenografting human-cancer tissues/cells with estradiol (E2) treatment is also used to generate rodent-models of gynaecological cancers. The altered metabolic-redox environment leading to establishment of pre-tumorigenesis condition and their mechanism are less studied. Here, female Wister rats were treated with these drugs at their pre-tumerogenic dosage (one group ENU single intra-peritoneal dose of 90 mg/kg b.w. and another group were implanted with human breast tumor (stage-IIIB) and fed with 2.5 mg of 17ß-estradiol once in a week for 4 months). After 4 months, animals were sacrificed; their serum and liver tissues were tested. A brief comparison was made with a rat model (regarded as positive control) of toxicity induced by mutagenic environmental pollutant arsenic (0.6 ppm daily/4 weeks). The increase in serum alkaline phosphatase and glutamate-pyruvate transaminase suggests the possible organ toxicity is favoured by the increase in hepatic/systemic free radicals and oxidative stress in all drug application models. But the increase in the serum E2 level as noted in the ELISA data with impairment in the hepatic estrogen sulfotransferase (SULT1E1) protein expression (immuno-blot data) were noticed with interfered hepatic free-thiols only in ENU and xenograft-E2 group compared to arsenic group. It is also evident in the in vitro result from E2/GSH/NAC added hepatic slices with altered antioxidant regulations. Moreover, impairment in hepatic SOD1, catalase and glutathiole peroxidase activities (PAGEzymographic data), especially in the ENU-treated group makes them more vulnerable to the oxidative threat in creating pre-tumerogenic microenvironment. This is evident in the result of their higher DNA-damage and histological abnormalities. The Bioinformatics study revealed an important role of rSULT1E1 in the regulations of E2 metabolism. This study is important for the exploration of the pre-tumerogenic condition by ENU and E2 by impairing SULT1E1 expression and E2 regulations via oxidant-stress signalling. The finding may help to find new therapeutic-targets to treat gynaecological-cancers more effectively.

Genome-wide mutation detection by interclonal genetic variation.

Genetic toxicology assays estimate mutation frequencies by phenotypically screening for the activation or inactivation of endogenous or exogenous reporter genes. These reporters can only detect mutations in narrow areas of the genome and their use is often restricted to certain in vitro and in vivo models. Here, we show that Interclonal Genetic Variation (ICGV) can directly identify mutations genome-wide by comparing sequencing data of single-cell clones derived from the same source or organism. Upon ethyl methanesulfonate (EMS) exposure, ICGV detected greater levels of mutation in a dose- and time-dependent manner in E. coli. In addition, ICGV was also able to identify a ∼20-fold increase in somatic mutations in T-cell clones derived from an N-ethyl-N-nitrosourea (ENU)-treated rat vs. a vehicle-treated rat. These results demonstrate that the genetic differences of single-cell clones can be used for genome-wide mutation detection.

Establishing a novel Pig-a gene mutation assay in L5178YTk+/- mouse lymphoma cells.

The X-linked Pig-a gene encodes an enzyme required for the biosynthesis of glycosyl phosphatidylinositol (GPI) anchors. Pig-a mutant cells fail to synthesize GPI and to express GPI-anchored protein markers (e.g., CD90) on their surface. Marker deficiency serves as a phenotypic indicator of Pig-a mutation in various in vivo assays. Here, we describe an in vitro Pig-a mutation assay in L5178YTk+/- mouse lymphoma cells, in which mutant-phenotype cells are measured by flow cytometry using a fluorescent anti-CD90 antibody. Increased frequencies of CD90-deficient mutants were detected in cells treated with benzo[a]pyrene (B[a]P), N-ethyl-N-nitrosourea (ENU), ethyl methanesulphonate, and 7,12-dimethylbenz[a]anthracene, with near maximum mutant frequencies measured eight days after treatment. The CD90 deficiency in mutant cells quantified by flow cytometry was shown to be due to loss of GPI anchors in a limiting-dilution cloning assay using proaerolysin selection. Individual CD90-deficient cells from cultures treated with ENU, B[a]P, and vehicle were sorted and clonally expanded for molecular analysis of their Pig-a gene. Pig-a mutations with agent-specific signatures were found in nearly all clones that developed from sorted CD90-deficient cells. These results indicate that a Pig-a mutation assay can be successfully conducted in L5178YTk+/- cells. The assay may be useful for mutagenicity screening of environmental agents as well as for testing hypotheses in vitro before committing to in vivo Pig-a assays. Environ. Mol. Mutagen. 59:4-17, 2018. Published 2017. This article is a US Government work and is in the public domain in the USA.

Cinacalcet corrects hypercalcemia in mice with an inactivating Gα11 mutation.

Loss-of-function mutations of GNA11, which encodes G-protein subunit α11 (Gα11), a signaling partner for the calcium-sensing receptor (CaSR), result in familial hypocalciuric hypercalcemia type 2 (FHH2). FHH2 is characterized by hypercalcemia, inappropriately normal or raised parathyroid hormone (PTH) concentrations, and normal or low urinary calcium excretion. A mouse model for FHH2 that would facilitate investigations of the in vivo role of Gα11 and the evaluation of calcimimetic drugs, which are CaSR allosteric activators, is not available. We therefore screened DNA from > 10,000 mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU) for GNA11 mutations and identified a Gα11 variant, Asp195Gly (D195G), which downregulated CaSR-mediated intracellular calcium signaling in vitro, consistent with it being a loss-of-function mutation. Treatment with the calcimimetic cinacalcet rectified these signaling responses. In vivo studies showed mutant heterozygous (Gna11+/195G) and homozygous (Gna11195G/195G) mice to be hypercalcemic with normal or increased plasma PTH concentrations and normal urinary calcium excretion. Cinacalcet (30mg/kg orally) significantly reduced plasma albumin-adjusted calcium and PTH concentrations in Gna11+/195G and Gna11195G/195G mice. Thus, our studies have established a mouse model with a germline loss-of-function Gα11 mutation that is representative for FHH2 in humans and demonstrated that cinacalcet can correct the associated abnormalities of plasma calcium and PTH.

An ENU-induced splice site mutation of mouse Col1a1 causing recessive osteogenesis imperfecta and revealing a novel splicing rescue.

GU-AG consensus sequences are used for intron recognition in the majority of cases of pre-mRNA splicing in eukaryotes. Mutations at splice junctions often cause exon skipping, short deletions, or insertions in the mature mRNA, underlying one common molecular mechanism of genetic diseases. Using N-ethyl-N-nitrosourea, a novel recessive mutation named seal was produced, associated with fragile bones and susceptibility to fractures (spine and limbs). A single nucleotide transversion (T → A) at the second position of intron 36 of the Col1a1 gene, encoding the type I collagen, α1 chain, was responsible for the phenotype. Col1a1 seal mRNA expression occurred at greatly reduced levels compared to the wild-type transcript, resulting in reduced and aberrant collagen fibers in tibiae of seal homozygous mice. Unexpectedly, splicing of Col1a1 seal mRNA followed the normal pattern despite the presence of the donor splice site mutation, likely due to the action of a putative intronic splicing enhancer present in intron 25, which appeared to function redundantly with the splice donor site of intron 36. Seal mice represent a model of human osteogenesis imperfecta, and reveal a previously unknown mechanism for splicing "rescue."

TKI-addicted ROS1-rearranged cells are destined to survival or death by the intensity of ROS1 kinase activity.

ROS1 rearrangement is observed in 1-2% of non-small cell lung cancers (NSCLC). The ROS1 tyrosine kinase inhibitor (TKI) crizotinib has induced marked tumour shrinkage in ROS1-rearranged cancers. However, emergence of acquired resistance to TKI is inevitable within a few years. Previous findings indicate that cabozantinib overcomes secondary mutation-mediated crizotinib-resistance in ROS1-fusion-positive cells. Here we attempted to establish cabozantinib-resistant cells by N-ethyl-N-nitrosourea mutagenesis screening using CD74-ROS1-expressing Ba/F3 cells. Two resistant cell lines with CD74-ROS1 F2004V or F2075C mutations, which are homologous to ALK F1174 or F1245 mutations, survived in the presence of a low dose of ROS1-TKI. Removal of ROS1-TKI from these TKI-addicted cells induced excessive activation of ROS1 tyrosine kinase followed by apoptosis. We succeeded in recapturing the TKI-addicted phenotype using doxycycline-inducible CD74-ROS1 mutant over-expression in Ba/F3 cells, suggesting that excessive ROS1 oncogenic signaling itself induced apoptosis instead of cell growth. Phosphoproteomic analysis and high-throughput inhibitor screening revealed that excessive ROS1 signaling in the TKI-addicted cells phosphorylated or activated apoptosis-related molecules such as FAF1 or p38. Collectively, our findings partly clarify molecular mechanisms of excessive ROS1 oncogenic signaling that mediates paradoxical induction of apoptosis.

The DNA methyltransferase DNMT3C protects male germ cells from transposon activity.

DNA methylation is prevalent in mammalian genomes and plays a central role in the epigenetic control of development. The mammalian DNA methylation machinery is thought to be composed of three DNA methyltransferase enzymes (DNMT1, DNMT3A, and DNMT3B) and one cofactor (DNMT3L). Here, we describe the discovery of Dnmt3C, a de novo DNA methyltransferase gene that evolved via a duplication of Dnmt3B in rodent genomes and was previously annotated as a pseudogene. We show that DNMT3C is the enzyme responsible for methylating the promoters of evolutionarily young retrotransposons in the male germ line and that this specialized activity is required for mouse fertility. DNMT3C reveals the plasticity of the mammalian DNA methylation system and expands the scope of the mechanisms involved in the epigenetic control of retrotransposons.

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.

Confirmation of Pig-a mutation in flow cytometry-identified CD48-deficient T-lymphocytes from F344 rats.

The Pig-a assay is used for monitoring somatic cell mutation in laboratory animals and humans. The assay detects haematopoietic cells deficient in glycosylphosphatidylinositol (GPI)-anchored protein surface markers using flow cytometry. However, given that synthesis of the protein markers (and the expression of their genes) is independent of the expression of the X-linked Pig-a gene and the function of its enzyme product, the deficiency of markers at the surface of the cells may be caused by a number of events (e.g. by mutation or epigenetic silencing in the marker gene itself or in any of about two dozen autosomal genes involved in the synthesis of GPI). Here we provide direct evidence that the deficiency of the GPI-anchored surface marker CD48 in rat T-cells is accompanied by mutation in the endogenous X-linked Pig-a gene. We treated male F344 rats with N-ethyl-N-nitrosourea (ENU), and established colonies from flow cytometry-identified and sorted CD48-deficient spleen T-lymphocytes. Molecular analysis confirmed that the expanded sorted cells have mutations in the Pig-a gene. The spectrum of Pig-a mutation in our model was consistent with the spectrum of ENU-induced mutation determined in other in vivo models, mostly base-pair substitutions at A:T with the mutated T on the non-transcribed strand of Pig-a genomic DNA. We also used next generation sequencing to derive a similar mutational spectrum from a pool of 64 clones developed from flow-sorted CD48-deficient lymphocytes. Our findings confirm that Pig-a assays detect what they are designed to detect-gene mutation in the Pig-a gene.

Single-cell transcriptogenomics reveals transcriptional exclusion of ENU-mutated alleles.

Recently, great progress has been made in single cell genomics and transcriptomics. Here, we present an integrative method, termed single-cell transcriptogenomics (SCTG), in which whole exome sequencing and RNA-seq is performed concurrently on single cells. This methodology enables one to track germline and somatic variants directly from the genome to the transcriptome in individual cells. Mouse embryonic fibroblasts were treated with the powerful mutagen ethylnitrosourea (ENU) and subjected to SCTG. Interestingly, while germline variants were found to be transcribed in an allelically balanced fashion, a significantly different pattern of allelic exclusion was observed for ENU-mutant variants. These results suggest that the adverse effects of induced mutations, in contrast to germline variants, may be mitigated by allelically biased transcription. They also illustrate how SCTG can be instrumental in the direct assessment of phenotypic consequences of genomic variants.

ENU mutagenesis identifies mice modeling Warburg Micro Syndrome with sensory axon degeneration caused by a deletion in Rab18.

Mutations in the gene of RAB18, a member of Ras superfamily of small G-proteins, cause Warburg Micro Syndrome (WARBM) which is characterized by defective neurodevelopmental and ophthalmological phenotypes. Despite loss of Rab18 had been reported to induce disruption of the endoplasmic reticulum structure and neuronal cytoskeleton organization, parts of the pathogenic mechanism caused by RAB18 mutation remain unclear. From the N-ethyl-N-nitrosourea (ENU)-induced mutagenesis library, we identified a mouse line whose Rab18 was knocked out. This Rab18(-/-) mouse exhibited stomping gait, smaller testis and eyes, mimicking several features of WARBM. Rab18(-/-) mice were obviously less sensitive to pain and touch than WT mice. Histological examinations on Rab18(-/-) mice revealed progressive axonal degeneration in the optic nerves, dorsal column of the spinal cord and sensory roots of the spinal nerves while the motor roots were spared. All the behavioral and pathological changes that resulted from abnormalities in the sensory axons were prevented by introducing an extra copy of Rab18 transgene in Rab18(-/-) mice. Our results reveal that sensory axonal degeneration is the primary cause of stomping gait and progressive weakness of the hind limbs in Rab18(-/-) mice, and optic nerve degeneration should be the major pathology of progressive optic atrophy in children with WARBM. Our results indicate that the sensory nervous system is more vulnerable to Rab18 deficiency and WARBM is not only a neurodevelopmental but also neurodegenerative disease.

Screen for alterations of iron related parameters in N-ethyl-N-nitrosourea-treated mice identified mutant lines with increased plasma ferritin levels.

Iron is essential for numerous cellular processes. For diagnostic purposes iron-related parameters in patients are assessed by clinical chemical blood analysis including the analysis of ferritin, transferrin and iron levels. Here, we retrospectively evaluated the use of these parameters in the phenotype-driven Munich N-ethyl-N-nitrosourea mouse mutagenesis project for the generation of novel animal models for human diseases. The clinical chemical blood analysis was carried out on more than 10,700 G1 and G3 offspring of chemically mutagenized inbred C3H mice to detect dominant and recessive mutations leading to deviations in the plasma levels of iron-related plasma parameters. We identified animals consistently exhibiting altered plasma ferritin or transferrin values. Transmission of the phenotypic deviations to the subsequent generations led to the successful establishment of three mutant lines with increased plasma ferritin levels. For two of these lines the causative mutations were identified in the Fth1gene and the Ireb2 gene, respectively. Thus, novel mouse models for the functional analysis of iron homeostasis were established by a phenotype-driven screen for mutant mice.

Role of TP53 mutations in the origin and evolution of therapy-related acute myeloid leukaemia.

Therapy-related acute myeloid leukaemia (t-AML) and therapy-related myelodysplastic syndrome (t-MDS) are well-recognized complications of cytotoxic chemotherapy and/or radiotherapy. There are several features that distinguish t-AML from de novo AML, including a higher incidence of TP53 mutations, abnormalities of chromosomes 5 or 7, complex cytogenetics and a reduced response to chemotherapy. However, it is not clear how prior exposure to cytotoxic therapy influences leukaemogenesis. In particular, the mechanism by which TP53 mutations are selectively enriched in t-AML/t-MDS is unknown. Here, by sequencing the genomes of 22 patients with t-AML, we show that the total number of somatic single-nucleotide variants and the percentage of chemotherapy-related transversions are similar in t-AML and de novo AML, indicating that previous chemotherapy does not induce genome-wide DNA damage. We identified four cases of t-AML/t-MDS in which the exact TP53 mutation found at diagnosis was also present at low frequencies (0.003-0.7%) in mobilized blood leukocytes or bone marrow 3-6 years before the development of t-AML/t-MDS, including two cases in which the relevant TP53 mutation was detected before any chemotherapy. Moreover, functional TP53 mutations were identified in small populations of peripheral blood cells of healthy chemotherapy-naive elderly individuals. Finally, in mouse bone marrow chimaeras containing both wild-type and Tp53(+/-) haematopoietic stem/progenitor cells (HSPCs), the Tp53(+/-) HSPCs preferentially expanded after exposure to chemotherapy. These data suggest that cytotoxic therapy does not directly induce TP53 mutations. Rather, they support a model in which rare HSPCs carrying age-related TP53 mutations are resistant to chemotherapy and expand preferentially after treatment. The early acquisition of TP53 mutations in the founding HSPC clone probably contributes to the frequent cytogenetic abnormalities and poor responses to chemotherapy that are typical of patients with t-AML/t-MDS.

Simulation and estimation of gene number in a biological pathway using almost complete saturation mutagenesis screening of haploid mouse cells.

BACKGROUND: Genome-wide saturation mutagenesis and subsequent phenotype-driven screening has been central to a comprehensive understanding of complex biological processes in classical model organisms such as flies, nematodes, and plants. The degree of "saturation" (i.e., the fraction of possible target genes identified) has been shown to be a critical parameter in determining all relevant genes involved in a biological function, without prior knowledge of their products. In mammalian model systems, however, the relatively large scale and labor intensity of experiments have hampered the achievement of actual saturation mutagenesis, especially for recessive traits that require biallelic mutations to manifest detectable phenotypes. RESULTS: By exploiting the recently established haploid mouse embryonic stem cells (ESCs), we present an implementation of almost complete saturation mutagenesis in a mammalian system. The haploid ESCs were mutagenized with the chemical mutagen N-ethyl-N-nitrosourea (ENU) and processed for the screening of mutants defective in various steps of the glycosylphosphatidylinositol-anchor biosynthetic pathway. The resulting 114 independent mutant clones were characterized by a functional complementation assay, and were shown to be defective in any of 20 genes among all 22 known genes essential for this well-characterized pathway. Ten mutants were further validated by whole-exome sequencing. The predominant generation of single-nucleotide substitutions by ENU resulted in a gene mutation rate proportional to the length of the coding sequence, which facilitated the experimental design of saturation mutagenesis screening with the aid of computational simulation. CONCLUSIONS: Our study enables mammalian saturation mutagenesis to become a realistic proposition. Computational simulation, combined with a pilot mutagenesis experiment, could serve as a tool for the estimation of the number of genes essential for biological processes such as drug target pathways when a positive selection of mutants is available.