Mutations In Essential Genes Research Articles

Fitness benefits of a synonymous substitution in an ancient EF-Tu gene depend on the genetic background.

Synonymous mutations are changes to DNA sequence, which occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation elongation factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels as well as global polysome abundance on RNA transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.IMPORTANCEThis study explores the degree to which synonymous mutations in essential genes can influence adaptation in bacteria. An experimental system whereby an Escherichia coli strain harboring an engineered translation protein elongation factor-Tu (EF-Tu) was subjected to laboratory evolution. We find that a synonymous mutation acquired on the gene encoding for EF-Tu is conditionally beneficial for bacterial fitness. Our findings provide insight into the importance of the genetic background when a synonymous substitution is favored by natural selection and how such changes have the potential to impact evolution when critical cellular processes are involved.

Genome-scale analysis of essential gene knockout mutants to identify an antibiotic target process.

We describe a genome-scale approach to identify the essential biological process targeted by a new antibiotic. The procedure is based on the identification of essential genes whose inactivation sensitizes a Gram-negative bacterium (Acinetobacter baylyi) to a drug and employs recently developed transposon mutant screening and single-mutant validation procedures. The approach, based on measuring the rates of loss of newly generated knockout mutants in the presence of antibiotic, provides an alternative to traditional procedures for studying essential functions using conditional expression or activity alleles. As a proof of principle study, we evaluated whether mutations enhancing sensitivity to the β-lactam antibiotic meropenem corresponded to the known essential target process of the antibiotic (septal peptidoglycan synthesis). We found that indeed mutations inactivating most genes needed for peptidoglycan synthesis and cell division strongly sensitized cells to meropenem. Additional classes of sensitizing mutations in essential genes were also identified, including those that inactivated capsule synthesis, DNA replication, or envelope stress response regulation. The essential capsule synthesis mutants appeared to enhance meropenem sensitivity by depleting a precursor needed for both capsule and peptidoglycan synthesis. The replication mutants may sensitize cells by impairing division. Nonessential gene mutations sensitizing cells to meropenem were also identified in the screen and largely corresponded to functions subordinately associated with the essential target process, such as in peptidoglycan recycling. Overall, these results help validate a new approach to identify the essential process targeted by an antibiotic and define the larger functional network determining sensitivity to it.

Strategies to identify and edit improvements in synthetic genome segments episomally.

Genome engineering projects often utilize bacterial artificial chromosomes (BACs) to carry multi-kilobase DNA segments at low copy number. However, all stages of whole-genome engineering have the potential to impose mutations on the synthetic genome that can reduce or eliminate the fitness of the final strain. Here, we describe improvements to a multiplex automated genome engineering (MAGE) protocol to improve recombineering frequency and multiplexability. This protocol was applied to recoding an Escherichia coli strain to replace seven codons with synonymous alternatives genome wide. Ten 44 402-47 179 bp de novo synthesized DNA segments contained in a BAC from the recoded strain were unable to complement deletion of the corresponding 33-61 wild-type genes using a single antibiotic resistance marker. Next-generation sequencing (NGS) was used to identify 1-7 non-recoding mutations in essential genes per segment, and MAGE in turn proved a useful strategy to repair these mutations on therecoded segment contained in the BAC when both the recoded and wild-type copies of the mutated genes had to exist by necessity during the repair process. Finally, two web-based tools were used to predict the impact of a subset of non-recoding missense mutations on strain fitness using protein structure and function calls.

Permissiveness and competition within and between Neurospora crassa syncytia.

A multinucleate syncytium is a common growth form in filamentous fungi. Comprehensive functions of the syncytial state remain unknown, but it likely allows for a wide range of adaptations to enable filamentous fungi to coordinate growth, reproduction, responses to the environment, and to distribute nuclear and cytoplasmic elements across a colony. Indeed, the underlying mechanistic details of how syncytia regulate cellular and molecular processes spatiotemporally across a colony are largely unexplored. Here, we implemented a strategy to analyze the relative fitness of different nuclear populations in syncytia of Neurospora crassa, including nuclei with loss-of-function mutations in essential genes, based on production of multinucleate asexual spores using flow cytometry of pairings between strains with differentially fluorescently tagged nuclear histones. The distribution of homokaryotic and heterokaryotic asexual spores in pairings was assessed between different auxotrophic and morphological mutants, as well as with strains that were defective in somatic cell fusion or were heterokaryon incompatible. Mutant nuclei were compartmentalized into both homokaryotic and heterokaryotic asexual spores, a type of bet hedging for maintenance and evolution of mutational events, despite disadvantages to the syncytium. However, in pairings between strains that were blocked in somatic cell fusion or were heterokaryon incompatible, we observed a "winner-takes-all" phenotype, where asexual spores originating from paired strains were predominantly one genotype. These data indicate that syncytial fungal cells are permissive and tolerate a wide array of nuclear functionality, but that cells/colonies that are unable to cooperate via syncytia formation actively compete for resources.

Evolutionary Analysis of TCGA Data Using Over- and Under- Mutated Genes Identify Key Molecular Pathways and Cellular Functions in Lung Cancer Subtypes.

We identify critical conserved and mutated genes through a theoretical model linking a gene’s fitness contribution to its observed mutational frequency in a clinical cohort. “Passenger” gene mutations do not alter fitness and have mutational frequencies determined by gene size and the mutation rate. Driver mutations, which increase fitness (and proliferation), are observed more frequently than expected. Non-synonymous mutations in essential genes reduce fitness and are eliminated by natural selection resulting in lower prevalence than expected. We apply this “evolutionary triage” principle to TCGA data from EGFR-mutant, KRAS-mutant, and NEK (non-EGFR/KRAS) lung adenocarcinomas. We find frequent overlap of evolutionarily selected non-synonymous gene mutations among the subtypes suggesting enrichment for adaptations to common local tissue selection forces. Overlap of conserved genes in the LUAD subtypes is rare suggesting negative evolutionary selection is strongly dependent on initiating mutational events during carcinogenesis. Highly expressed genes are more likely to be conserved and significant changes in expression (>20% increased/decreased) are common in genes with evolutionarily selected mutations but not in conserved genes. EGFR-mut cancers have fewer average mutations (89) than KRAS-mut (228) and NEK (313). Subtype-specific variation in conserved and mutated genes identify critical molecular components in cell signaling, extracellular matrix remodeling, and membrane transporters. These findings demonstrate subtype-specific patterns of co-adaptations between the defining driver mutation and somatically conserved genes as well as novel insights into epigenetic versus genetic contributions to cancer evolution.

Target Identification of Small Molecules Using Large-Scale CRISPR-Cas Mutagenesis Scanning of Essential Genes.

Target deconvolution of new bioactive agents identified from phenotypic screens remains a challenging task. The discovery of mutations that confer resistance to such agents is regarded as the gold standard proof of target identification. Here, we describe a method that exploits the error-prone repair of CRISPR-induced DNA double-strand breaks to enhance mutagenesis and increase the incidence of drug resistancemutations in essential genes. As each DNA double-strand break is introduced at a targeted genomic site predefined by the presence of a protospacer adjacent motif (PAM) and a particular CRISPR single guide RNA (sgRNA), the genetic location of drug resistance mutations can be easily uncovered through targeted sequencing of CRISPR sgRNAs. Moreover, the method allows for the identification of not only the drug target gene, but also the drug-binding domain within the target gene.

Identification of Essential Genes in Caenorhabditis elegans with Lethal Mutations Maintained by Genetic Balancers.

Genetic balancer systems, which allow effective capture and maintenance of lethal mutations stably, play an important role in identifying essential genes. Whole-genome sequencing (WGS) followed by bioinformatics analysis, combined with genetic mapping data analysis, allows for an efficient and economical means of identifying genomic mutations in essential genes. Using this approach, we successfully identified 104 essential genes on ChrI, ChrIII, and ChrV in C. elegans. In this report, we described a protocol that sequences the genome of prebalanced Caenorhabditis elegans (C. elegans) strains to carry lethal mutations and identifies candidate causal mutations and candidate essential genes using a robust bioinformatics procedure.

A CRISPR-Based Method for Constructing Conditional Mutations of Essential Genes in Cyanobacteria.

Cyanobacteria, a group of diverse bacteria capable of oxygenic photosynthesis, are excellent models for investigating many important cellular processes, such as photosynthesis, nitrogen fixation, and prokaryotic cell differentiation. They also have great potential to become the next-generation cell factories for sustainable biosynthesis of valuable products. However, genetic manipulation in cyanobacteria is not as convenient as in other model bacteria. Particularly, handling essential genes in cyanobacteria has been difficult due to the lack of appropriate tools, limiting our understanding of many important cellular functions encoded by them. We recently develop a CRISPR-based method for constructing the conditional mutants of cyanobacterial essential genes by engineering the ribosome binding site to a theophylline-responsive riboswitch. Here, we provide the details of this method. The principle of this method could be used to construct conditional mutants in a wide range of bacterial species.

Exploring a Local Genetic Interaction Network Using Evolutionary Replay Experiments

Understanding how genes interact is a central challenge in biology. Experimental evolution provides a useful, but underutilized, tool for identifying genetic interactions, particularly those that involve non-loss-of-function mutations or mutations in essential genes. We previously identified a strong positive genetic interaction between specific mutations in KEL1 (P344T) and HSL7 (A695fs) that arose in an experimentally evolved Saccharomyces cerevisiae population. Because this genetic interaction is not phenocopied by gene deletion, it was previously unknown. Using “evolutionary replay” experiments, we identified additional mutations that have positive genetic interactions with the kel1-P344T mutation. We replayed the evolution of this population 672 times from six timepoints. We identified 30 populations where the kel1-P344T mutation reached high frequency. We performed whole-genome sequencing on these populations to identify genes in which mutations arose specifically in the kel1-P344T background. We reconstructed mutations in the ancestral and kel1-P344T backgrounds to validate positive genetic interactions. We identify several genetic interactors with KEL1, we validate these interactions by reconstruction experiments, and we show these interactions are not recapitulated by loss-of-function mutations. Our results demonstrate the power of experimental evolution to identify genetic interactions that are positive, allele specific, and not readily detected by other methods, shedding light on an underexplored region of the yeast genetic interaction network.

Protein Kinase C Isozymes Associated With Relapse Free Survival in Non-Small Cell Lung Cancer Patients.

IntroductionProtein expression is deregulated in cancer, and the proteomic changes observed in lung cancer may be a consequence of mutations in essential genes. The purpose of this study was to identify protein expression associated with prognosis in lung cancers stratified by smoking status, molecular subtypes, and EGFR-, TP53-, and KRAS-mutations.MethodsWe performed profiling of 295 cancer-relevant phosphorylated and non-phosphorylated proteins, using reverse phase protein arrays. Biopsies from 80 patients with operable lung adenocarcinomas were analyzed for protein expression and association with relapse free survival (RFS) were studied.ResultsSpearman’s rank correlation analysis identified 46 proteins with significant association to RFS (p<0.05). High expression of protein kinase C (PKC)-α and the phosporylated state of PKC-α, PKC-β, and PKC-δ, showed the strongest positive correlation to RFS, especially in the wild type samples. This was confirmed in gene expression data from 172 samples. Based on protein expression, unsupervised hierarchical clustering separated the samples into four subclusters enriched with the molecular subtypes terminal respiratory unit (TRU), proximal proliferative (PP), and proximal inflammatory (PI) (p=0.0001). Subcluster 2 contained a smaller cluster (2a) enriched with samples of the subtype PP, low expression of the PKC isozymes, and associated with poor RFS (p=0.003) compared to the other samples. Low expression of the PKC isozymes in the subtype PP and a reduced relapse free survival was confirmed with The Cancer Genome Atlas (TCGA) lung adenocarcinoma (LUAD) samples.ConclusionThis study identified different proteins associated with RFS depending on molecular subtype, smoking- and mutational-status, with PKC-α, PKC-β, and PKC-δ showing the strongest correlation.

High prevalence of autosomal recessive congenital ichthyosis in a Mexican population caused by a new mutation in the TGM1 gene: epidemiological evidence of a founder effect.

Autosomal recessive congenital ichthyoses (ARCI) are inherited disorders produced by mutations in essential genes for the skin function. A low prevalence of this disease has been resportedworldwide; however, in a recent study, we identified a large cluster of ARCI families who resided in the High Mountains Region from the Veracruz State, Mexico. Thus, we aimed to identify the causative mutation of ARCI and describe the high prevalence of this disease in this region. We selected seven familiar trios and performed whole-exome sequencing to identify the mutation associated with ARCI. To validate the identified mutation, we performed Sanger sequencing in 62 patients, 30 unaffected relatives, and 100 healthy volunteers. Finally, we performed molecular modeling to investigate the possible functional consequences produced by the mutation. We identified a novel homozygous mutation (c.1054C>G [p.Pro352Ala]) in the exon 7 of the TGM1 gene in all the patients. We calculated a prevalence rate of ARCI of 74:100,000 (1:1,348) in the studied communities. Molecular modeling revealed that the mutation leads to a nonconservative amino acid substitution, which is very probably damaging to the protein structure/function. We report a novel mutation in the TGM1 gene in 62 Mexican patients. The unusually high frequency of this mutation suggests a founder effect; however, further haplotype analysis is necessary to corroborate this hypothesis. In this respect, to our knowledge, the prevalence of ARCI found in the studied communities is the highest observed worldwide.

Efficient targeted mutation of genomic essential genes in yeast Saccharomyces cerevisiae.

Targeted gene mutation by allelic replacement is important for functional genomic analysis and metabolic engineering. However, it is challenging in mutating the essential genes with the traditional method by using a selection marker, since the first step of essential gene knockout will result in a lethal phenotype. Here, we developed a two-end selection marker (Two-ESM) method for site-directed mutation of essential genes in Saccharomyces cerevisiae with the aid of the CRISPR/Cas9 system. With this method, single and double mutations of the essential gene ERG20 (encoding farnesyl diphosphate synthase) in S. cerevisiae were successfully constructed with high efficiencies of 100%. In addition, the Two-ESM method significantly improved the mutation efficiency and simplified the genetic manipulation procedure compared with traditional methods. The genome integration and mutation efficiencies were further improved by dynamic regulation of mutant gene expression and optimization of the integration modules. This Two-ESM method will facilitate the construction of genomic mutations of essential genes for functional genomic analysis and metabolic flux regulation in yeasts. KEY POINTS: • A Two-ESM strategy achieves mutations of essential genes with high efficiency of 100%. • The optimized three-module method improves the integration efficiency by more than three times. • This method will facilitate the functional genomic analysis and metabolic flux regulation.

Loss of UGP2 in brain leads to a severe epileptic encephalopathy, emphasizing that bi-allelic isoform-specific start-loss mutations of essential genes can cause genetic diseases

Developmental and/or epileptic encephalopathies (DEEs) are a group of devastating genetic disorders, resulting in early-onset, therapy-resistant seizures and developmental delay. Here we report on 22 individuals from 15 families presenting with a severe form of intractable epilepsy, severe developmental delay, progressive microcephaly, visual disturbance and similar minor dysmorphisms. Whole exome sequencing identified a recurrent, homozygous variant (chr2:64083454A > G) in the essential UDP-glucose pyrophosphorylase (UGP2) gene in all probands. This rare variant results in a tolerable Met12Val missense change of the longer UGP2 protein isoform but causes a disruption of the start codon of the shorter isoform, which is predominant in brain. We show that the absence of the shorter isoform leads to a reduction of functional UGP2 enzyme in neural stem cells, leading to altered glycogen metabolism, upregulated unfolded protein response and premature neuronal differentiation, as modeled during pluripotent stem cell differentiation in vitro. In contrast, the complete lack of all UGP2 isoforms leads to differentiation defects in multiple lineages in human cells. Reduced expression of Ugp2a/Ugp2b in vivo in zebrafish mimics visual disturbance and mutant animals show a behavioral phenotype. Our study identifies a recurrent start codon mutation in UGP2 as a cause of a novel autosomal recessive DEE syndrome. Importantly, it also shows that isoform-specific start-loss mutations causing expression loss of a tissue-relevant isoform of an essential protein can cause a genetic disease, even when an organism-wide protein absence is incompatible with life. We provide additional examples where a similar disease mechanism applies.

CRISPR-TSKO: A Technique for Efficient Mutagenesis in Specific Cell Types, Tissues, or Organs in Arabidopsis.

Detailed functional analyses of many fundamentally important plant genes via conventional loss-of-function approaches are impeded by the severe pleiotropic phenotypes resulting from these losses. In particular, mutations in genes that are required for basic cellular functions and/or reproduction often interfere with the generation of homozygous mutant plants, precluding further functional studies. To overcome this limitation, we devised a clustered regularly interspaced short palindromic repeats (CRISPR)-based tissue-specific knockout system, CRISPR-TSKO, enabling the generation of somatic mutations in particular plant cell types, tissues, and organs. In Arabidopsis (Arabidopsis thaliana), CRISPR-TSKO mutations in essential genes caused well-defined, localized phenotypes in the root cap, stomatal lineage, or entire lateral roots. The modular cloning system developed in this study allows for the efficient selection, identification, and functional analysis of mutant lines directly in the first transgenic generation. The efficacy of CRISPR-TSKO opens avenues for discovering and analyzing gene functions in the spatial and temporal contexts of plant life while avoiding the pleiotropic effects of system-wide losses of gene function.

Escherichia coli vectors having stringently repressible replication origins allow a streamlining of Crispr/Cas9 gene editing

Readily curable plasmids facilitate the construction of plasmid-free bacterial strains after the plasmid encoded genes are no longer needed. The most popular of these plasmids features a temperature-sensitive (Ts) pSC101 origin of replication which can readily revert during usage and cannot be used to construct Ts mutations in essential genes. Plasmid pAM34 which contains an IPTG-dependent origin of replication largely overcomes this issue but is limited by carrying the most commonly utilized antibiotic selection and replication origin. This study describes the construction of an expanded series of plasmid vectors having replication origins of p15a, RSF1030 or RSF1031 that like pAM34 have IPTG-dependent replication. Surprisingly, these plasmids can be cured in fewer generations than pAM34. Derivatives of pAM34 with alternative antibiotic selection markers were also constructed. The utility of these vectors is demonstrated in the construction of a CRISPR-Cas9 system consisting of an IPTG-dependent Cas9 plasmid and a curable guide RNA plasmid having a streptomycin counterselection marker. This system was successfully demonstrated by construction of point mutations, deletions and insertions in the E. coli genome with a very high efficiency and in a shorter timescale than extant methods. The plasmids themselves were readily cured either together or singly from the resultant strains with minimal effort.

KrCRISPR: an easy and efficient strategy for generating conditional knockout of essential genes in cells

sBackgroundCRISPR/Cas9 system is a powerful tool for knocking out genes in cells. However, genes essential for cell survival cannot be directly knocked out. Traditionally, generation of conditional knockout cells requires multiple steps.ResultsIn this study, we developed an easy and efficient strategy to generate conditional knockout cells by using double episomal vectors – one which expresses gRNA and Cas9 nuclease, and the other expresses an inducible rescue gene. Using this system which we named “krCRISPR” (knockout-rescue CRISPR), we showed that essential genes, HDAC3 and DNMT1, can be efficiently knocked out. When cells reach a desired confluency, the exogenous rescue genes can be silenced by the addition of doxycycline. Furthermore, the krCRISPR system enabled us to study the effects of the essential gene mutations on cells. We showed that the P507L mutation in DNMT1 led to downregulation of global DNA methylation in cells, indicating that it is a disease-causing mutation.ConclusionsThe krCRISPR system offers an easy and efficient platform that facilitates the study of essential genes’ function.

Loss of function mutations in essential genes cause embryonic lethality in pigs.

Lethal recessive alleles cause pre- or postnatal death in homozygous affected individuals, reducing fertility. Especially in small size domestic and wild populations, those alleles might be exposed by inbreeding, caused by matings between related parents that inherited the same recessive lethal allele from a common ancestor. In this study we report five relatively common (up to 13.4% carrier frequency) recessive lethal haplotypes in two commercial pig populations. The lethal haplotypes have a large effect on carrier-by-carrier matings, decreasing litter sizes by 15.1 to 21.6%. The causal mutations are of different type including two splice-site variants (affecting POLR1B and TADA2A genes), one frameshift (URB1), and one missense (PNKP) variant, resulting in a complete loss-of-function of these essential genes. The recessive lethal alleles affect up to 2.9% of the litters within a single population and are responsible for the death of 0.52% of the total population of embryos. Moreover, we provide compelling evidence that the identified embryonic lethal alleles contribute to the observed heterosis effect for fertility (i.e. larger litters in crossbred offspring). Together, this work marks specific recessive lethal variation describing its functional consequences at the molecular, phenotypic, and population level, providing a unique model to better understand fertility and heterosis in livestock.

Plexin-Semaphorin Signaling Modifies Neuromuscular Defects in a Drosophila Model of Peripheral Neuropathy.

Dominant mutations in GARS, encoding the ubiquitous enzyme glycyl-tRNA synthetase (GlyRS), cause peripheral nerve degeneration and Charcot-Marie-Tooth disease type 2D (CMT2D). This genetic disorder exemplifies a recurring paradigm in neurodegeneration, in which mutations in essential genes cause selective degeneration of the nervous system. Recent evidence suggests that the mechanism underlying CMT2D involves extracellular neomorphic binding of mutant GlyRS to neuronally-expressed proteins. Consistent with this, our previous studies indicate a non-cell autonomous mechanism, whereby mutant GlyRS is secreted and interacts with the neuromuscular junction (NMJ). In this Drosophila model for CMT2D, we have previously shown that mutant gars expression decreases viability and larval motor function, and causes a concurrent build-up of mutant GlyRS at the larval neuromuscular presynapse. Here, we report additional phenotypes that closely mimic the axonal branching defects of Drosophila plexin transmembrane receptor mutants, implying interference of plexin signaling in gars mutants. Individual dosage reduction of two Drosophila Plexins, plexin A (plexA) and B (plexB) enhances and represses the viability and larval motor defects caused by mutant GlyRS, respectively. However, we find plexB levels, but not plexA levels, modify mutant GlyRS association with the presynaptic membrane. Furthermore, increasing availability of the plexB ligand, Semaphorin-2a (Sema2a), alleviates the pathology and the build-up of mutant GlyRS, suggesting competition for plexB binding may be occurring between these two ligands. This toxic gain-of-function and subversion of neurodevelopmental processes indicate that signaling pathways governing axonal guidance could be integral to neuropathology and may underlie the non-cell autonomous CMT2D mechanism.

Target identification of small molecules using large-scale CRISPR-Cas mutagenesis scanning of essential genes

Unraveling the mechanism of action and molecular target of small molecules remains a major challenge in drug discovery. While many cancer drugs target genetic vulnerabilities, loss-of-function screens fail to identify essential genes in drug mechanism of action. Here, we report CRISPRres, a CRISPR-Cas-based genetic screening approach to rapidly derive and identify drug resistance mutations in essential genes. It exploits the local genetic variation created by CRISPR-Cas-induced non-homologous end-joining (NHEJ) repair to generate a wide variety of functional in-frame mutations. Using large sgRNA tiling libraries and known drug–target pairs, we validate it as a target identification approach. We apply CRISPRres to the anticancer agent KPT-9274 and identify nicotinamide phosphoribosyltransferase (NAMPT) as its main target. These results present a powerful and simple genetic approach to create many protein variants that, in combination with positive selection, can be applied to reveal the cellular target of small-molecule inhibitors.

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.