Targeted recovery of mutations in Drosophila

Reverse genetic techniques will be necessary to take full advantage of the genomic sequence data for Drosophila and other experimental organisms. To develop a method for the targeted recovery of mutations, we combined an EMS chemical mutagenesis regimen with mutation detection by denaturing high performance liquid chromatography (DHPLC). We recovered mutant strains at the high rate of approximately 4.8 mutations/kb for every 1000 mutagenized chromosomes from a screen for new mutations in the Drosophila awd gene. Furthermore, we observed that the EMS mutational spectrum in Drosophila germ cells shows a strong preference for 5'-PuG-3' sites, and for G/C within a stretch of three or more G/C base pairs. Our method should prove useful for targeted mutagenesis screens in Drosophila and other genetically tractable organisms and for more precise studies of mutagenesis and DNA repair mechanisms.     

Genome-wide analysis of mutagenesis bias and context sensitivity of N-methyl-N'-nitro-N-nitrosoguanidine (NTG)

We have analyzed the mutation spectrum of N-methyl-N'-nitro-N-nitrosoguanidine (NTG) from a set of 4099 mutations identified from whole-genome sequencing of 32 E. coli strains mutagenized with NTG. These data permit precise measurement of NTG's bias for G/C to A/T transitions (96.6% of all mutations) and also show that NTG mutagenesis is strongly sensitive to context, favoring guanine residues preceded by purines by five-fold over those preceded by pyrimidines. These data give confident estimates for the GC bias and transition/transversion ratios of NTG mutagenesis, which could not be estimated confidently from previous, much smaller datasets.     

A comparison of directed evolution approaches using the beta-glucuronidase model system

Protein engineers can alter the properties of enzymes by directing their evolution in vitro. Many methods to generate molecular diversity and to identify improved clones have been developed, but experimental evolution remains as much an art as a science. We previously used DNA shuffling (sexual recombination) and a histochemical screen to direct the evolution of Escherichia coli beta-glucuronidase (GUS) variants with improved beta-galactosidase (BGAL) activity. Here, we employ the same model evolutionary system to test the efficiencies of several other techniques: recursive random mutagenesis (asexual), combinatorial cassette mutagenesis (high-frequency recombination) and a versatile high-throughput microplate screen. GUS variants with altered specificity evolved in each trial, but different combinations of mutagenesis and screening techniques effected the fixation of different beneficial mutations. The new microplate screen identified a broader set of mutations than the previously employed X-gal colony screen. Recursive random mutagenesis produced essentially asexual populations, within which beneficial mutations drove each other into extinction (clonal interference); DNA shuffling and combinatorial cassette mutagenesis led instead to the accumulation of beneficial mutations within a single allele. These results explain why recombinational approaches generally increase the efficiency of laboratory evolution.     

From mouse to man: generating megabase chromosome rearrangements

Experimental approaches for deciphering the function of human genes rely heavily on our ability to generate mutations in model organisms such as the mouse. However, because recessive mutations are masked by the wild-type allele in the diploid context, conventional mutagenesis and screening is often laborious and costly. Chromosome engineering combines the power of gene targeting in embryonic stem (ES) cells with Cre--loxP technology to create mice that are functionally haploid in discrete portions of the genome. Chromosome deletions, duplications and inversions can be tagged with visible markers, facilitating strain maintenance. These approaches allow for more refined mutagenesis screens that will greatly accelerate functional mouse genomics and generate mammalian models for developmental processes and cancer.     

A general method for introducing a series of mutations into cloned DNA using the polymerase chain reaction

A simple and fast method for introducing a series of mutations in cloned DNA has been developed. The polymerase chain reaction (PCR) has been used for site-directed mutagenesis. Because mutations can be introduced only within the primer sequences used for PCR, a suitable restriction site in the vicinity of the mutated nucleotide is required to permit recloning. Several methods have been devised to overcome this limitation. Our present method is a modification of the overlap extension method [Ho et al., Gene 77 (1989) 51-57], and has some advantages over this and other published methods. In our method, three common primers and a series of primers specific for various mutations are chemically synthesized. Once the proper oligodeoxyribonucleotides are selected as common primers, each mutation requires only one additional primer. Therefore, this method is very useful for introducing many mutations in various sites of the target DNA. We describe our protocol for the site-directed mutagenesis and an example of the introduction of several mutations in the hen egg-white lysozyme-encoding gene.     

Mutational Analysis Defines a C-Terminal Tail Domain of Rap1 Essential for Telomeric Silencing in Saccharomyces Cerevisiae

Alleles specifically defective in telomeric silencing were generated by in vitro mutagenesis of the yeast RAP1 gene. The most severe phenotypes occur with three mutations in the C-terminal 28 amino acids. Two of the alleles are nonsense mutations resulting in truncated repressor/activator protein 1 (RAP1) species lacking the C-terminal 25-28 amino acids; the third allele is a missense mutation within this region. These alleles define a novel 28-amino acid region, termed the C-terminal tail domain, that is essential for telomeric and HML silencing. Using site-directed mutagenesis, an 8-amino acid region (amino acids 818-825) that is essential for telomeric silencing has been localized within this domain. Further characterization of these alleles has indicated that the C-terminal tail domain also plays a role in telomere size control. The function of the C-terminal tail in telomere maintenance is not mediated through the RAP1 interacting factor RIF1: rap1 alleles defective in both the C-terminal tail and RIF1 interaction domains have additive effects on telomere length. Overproduction of SIR3, a dose-dependent enhancer of telomeric silencing, suppresses the telomeric silencing, but not length, phenotypes of a subset of C-terminal tail alleles. In contrast, an allele that truncates the terminal 28 amino acids of RAP1 is refractory to SIR3 overproduction. These results indicate that the C-terminal tail domain is required for SIR3-dependent enhancement of telomeric silencing. These data also suggest a distinct set of C-terminal requirements for telomere size control and telomeric silencing.     

Three different frameshift mutations of the tyrosinase gene in type IA oculocutaneous albinism.

Mutations in the gene for the pigment-producing enzyme tyrosinase are responsible for type IA (tyrosinase-negative) oculocutaneous albinism (OCA). Most reported mutations have been single base substitutions. We now report three different frameshift mutations in three unrelated individuals with type IA OCA. The first individual has a single base deletion within a series of five guanidines, resulting in a premature stop codon in exon I on one allele and a missense mutation at codon 382 in exon III on the homologous allele. The second individual is a genetic compound of two separate frameshift mutations, including both the same exon I single base deletion found in the first individual and a deletion of a thymidine-guanidine pair, within the sequence GTGTG, forming a termination codon (TAG) in exon I on the homologous allele. The third individual has a single base insertion in exon I on one allele and a missense mutation at codon 373 in exon III on the homologous allele. The two missense mutations occur within the copper Bbinding region and may interfere with either copper binding to the enzyme or oxygen binding to the copper. These five different mutations disrupt tyrosinase function and are associated with a total lack of melanin biosynthesis.     

Allelic losses in mutations at the aprt locus of human lymphoblastoid cells.

We analyzed the nature of mutations at the autosomal locus coding for adenine phosphoribosyltransferase (aprt) in human cells to elucidate the process(es) governing mutagenesis at autosomal loci. A human lymphoblastoid cell line, WR10, was found to be heterozygous for mutated allele at the aprt locus, and was used for mutation analyses. By the use of a restriction fragment length polymorphism associated with the aprt locus in WR10 cells, the molecular characteristics of mutations arising spontaneously or induced by gamma-rays were investigated. Eighty-five percent (22/26) of the spontaneous mutant clones and 93% (64/69) of the gamma-ray-induced mutant clones resulted from loss of one of the two aprt alleles. Determination of the dosage of aprt genes in those mutants with allelic losses revealed that approximately half of them retained two copies of the mutated allele. These data suggest that the mutational events leading to APRT deficiency are analogous to those reported for tumor suppressor genes in malignancies.     

Metabolic adaptations of Leishmania donovani in relation to differentiation,drug resistance,and drug pressure

Antimonial (sodium stibogluconate, SSG) resistance and differentiation have been shown to be closely linked in Leishmania donovani, with SSG‐resistant strains showing an increased capacity to generate infectious (metacyclic) forms. This is the first untargeted LC‐MS metabolomics study which integrated both phenomena in one experimental design and provided insights into metabolic differences between three clinical L. donovani strains with a similar genetic background but different SSG‐susceptibilities. We performed this analysis at different stages during promastigote growth and in the absence or presence of drug pressure. When comparing SSG‐resistant and SSG‐sensitive strains, a number of metabolic changes appeared to be constitutively present in all growth stages, pointing towards a clear link with SSG‐resistance, whereas most metabolic changes were only detected in the stationary stage. These changes reflect the close intertwinement between SSG‐resistance and an increased metacyclogenesis in resistant parasites. The metabolic changes suggest that SSG‐resistant parasites have (i) an increased capacity for protection against oxidative stress; (ii) a higher fluidity of the plasma membrane; and (iii) a metabolic survival kit to better endure infection. These changes were even more pronounced in a resistant strain kept under Sb^III drug pressure.     

Single-step overlap-primer-walk polymerase chain reaction for multiple mutagenesis without overlap extension

We present a simple, single-step, single-tube, and rapid method for introducing a series of mutations into cloned DNA. Polymerase chain reaction (PCR)-based mutagenesis methods have become very prevalent due to their simplicity and efficiency for introducing mutations. Our method, overlap-primer-walk PCR, has several advantages over other published methods. It uses two common oligodeoxyribonucleotides and a series of overlapping primers specific for various mutations. Once common flanking primers are selected, two to three mutations require only one additional primer. Therefore, this method is very useful for introduction of multiple mutations in various sites of the target DNA. We illustrate the usefulness of the method by introducing several mutations into the human TNF-α encoding gene.     

Redox imbalance and mutagenesis in spleens of mice harboring a hypomorphic allele of Gpdx(a) encoding glucose 6-phosphate dehydrogenase

Mice harboring the activity-attenuated Gpdx(a-m2Neu) allele and also harboring a chromosomally integrated lacZ reporter gene to study mutagenesis (pUR288) were used to demonstrate that moderate glucose 6-phosphate dehydrogenase (G6PD) deficiency causes elevated mutagenesis and endogenous oxidative stress in the spleen. G6PD-deficient spleens with a residual enzyme activity of 22% exhibited a dramatic shift in the mutational pattern of lacZ (4.6-fold increase in the prevalence of recombination mutations of lacZ) together with a 1.8-fold increase in mutant frequencies in lacZ. A concomitant 3-fold reduction in catalase activity (dependent upon NADPH) indicated that the in vivo supply of G6PD-generated NADPH was insufficient. An additional 3-fold increase in oxidized glutathione suggested that redox control was disturbed in G6PD-deficient spleens. These findings indicate that G6PD is required for limiting oxidative mutagenesis in the mouse spleen. Gpdx(a-m2Neu) is the first hypomorphic allele of a mouse housekeeping gene associated with elevated somatic mutagenesis in vivo.     

Unlinked Noncomplementation: Isolation of New Conditional-Lethal Mutations in Each of the Tubulin Genes of Saccharomyces Cerevisiae

Mutations in genes of Saccharomyces cerevisiae that code for proteins that interact with beta-tubulin were sought by screening for unlinked mutations that fail to complement mutations in the single beta-tubulin-encoding gene (TUB2). Among the first three noncomplementing mutations examined, two are linked to TUB2 while one is unlinked. The unlinked mutation was shown to be a conditional-lethal allele of the major alpha-tubulin-encoding gene (TUB1) and represents the first such mutation in that gene. The tub1-1 mutation itself causes a cold-sensitive cell-cycle arrest, and confers supersensitivity to the antimicrotubule drug benomyl. These phenotypes occur in the presence of a wild-type copy of the minor alpha-tubulin-encoding gene, TUB3; the combination of tub1-1 and a tub3 null mutation is inviable in haploids. Through further application of this method, new mutations in TUB2 and TUB3 were isolated as unlinked noncomplementers of tub1-1. The noncomplementation between tub1 and tub2 mutations is gene specific and allele specific, suggesting that the phenotype is due to an interaction at the protein level. We conclude that isolation of unlinked noncomplementing mutations is likely to be a generally useful method for isolating mutations in interacting gene products.     

Multiple Site Directed Mutagenesis Strategy Based on Total RNA and RT-PCR Method

Site-directed PCR-based mutagenesis methods are widely used to generate mutations. All published methods work on DNA clones carrying the target sequence. However, DNA clones are not always available. We have previously published a RT-PCR-based site-directed mutagenesis method starting from total RNA to overcome this problem. In this article, we report an improvement of our previous method to facilitate introduction of multiple mutations into a target sequence. We demonstrate the efficacy and feasibility of this strategy by mutation of the human β-actin gene. BamHI restriction endonuclease cleavage sites were generated within the gene to assist screening. Using three mutagenic primers in a single RT-PCR reaction, seven different clones were produced carrying three single and four multiple mutations. An investigation of the effect of the cycle number and elongation time of the PCR reactions revealed that both have an influence on the ratio of clones carrying single and multiple mutations. An optimized protocol was established for efficient multiple site-directed mutagenesis.     

Isolation and characterization of noncleavable (Ind-) mutants of the LexA repressor of Escherichia coli K-12.

The LexA repressor of Escherichia coli represses a set of genes that are expressed in the response to DNA damage. After inducing treatments, the repressor is inactivated in vivo by a specific cleavage reaction which requires an activated form of RecA protein. In vitro, specific cleavage requires activated RecA at neutral pH and proceeds spontaneously at alkaline pH. We have isolated and characterized a set of lexA mutants that are deficient in in vivo RecA-mediated cleavage but retain significant repressor function. Forty-six independent mutants, generated by hydroxylamine and formic acid mutagenesis, were isolated by a screen involving the use of operon fusions. DNA sequence analysis identified 20 different mutations. In a recA mutant, all but four of the mutant proteins functioned as repressor as well as wild-type LexA. In a strain carrying a constitutively active recA allele, recA730, all the mutant proteins repressed a sulA::lacZ fusion more efficiently than the wild-type repressor, presumably because they were cleaved poorly or not at all by the activated RecA protein. These 20 mutations resulted in amino acid substitutions in 12 positions, most of which are conserved between LexA and four other cleavable proteins. All the mutations were located in the hinge region or C-terminal domain of the protein, portions of LexA previously implicated in the specific cleavage reactions. Furthermore, these mutations were clustered in three regions, around the cleavage site (Ala-84-Gly-85) and in blocks of conserved amino acids around two residues, Ser-119 and Lys-156, which are believed essential for the cleavage reactions. These three regions of the protein thus appear to play important roles in the cleavage reaction.     

Genetic analysis reveals a role for the C terminus of the Saccharomyces cerevisiae GTPase Snu114 during spliceosome activation

Snu114 is the only GTPase required for mRNA splicing. As a homolog of elongation factor G, it contains three domains (III-V) predicted to undergo a large rearrangement following GTP hydrolysis. To assess the functional importance of the domains of Snu114, we used random mutagenesis to create conditionally lethal alleles. We identified three main classes: (1) mutations that are predicted to affect GTP binding and hydrolysis, (2) mutations that are clustered in 10- to 20-amino-acid stretches in each of domains III-V, and (3) mutations that result in deletion of up to 70 amino acids from the C terminus. Representative mutations from each of these classes blocked the first step of splicing in vivo and in vitro. The growth defects caused by most alleles were synthetically exacerbated by mutations in PRP8, a U5 snRNP protein that physically interacts with Snu114, as well as in genes involved in snRNP biogenesis, including SAD1 and BRR1. The allele snu114-60, which truncates the C terminus, was synthetically lethal with factors required for activation of the spliceosome, including the DExD/H-box ATPases BRR2 and PRP28. We propose that GTP hydrolysis results in a rearrangement between Prp8 and the C terminus of Snu114 that leads to release of U1 and U4, thus activating the spliceosome for catalysis.     

The polyanion-binding domain of cytoplasmic Lys-tRNA synthetase from Saccharomyces cerevisiae is not essential for cell viability.

Cytoplasmic Lys-tRNA synthetase (LysRS) from Saccharomyces cerevisiae is a dimeric enzyme made up of identical subunits of 68 kDa. By limited proteolysis, this enzyme can be converted to a truncated dimer without loss of activity. Whereas the native enzyme strongly interacts with polyanionic carriers, the modified form displays reduced binding properties. KRS1 is the structural gene for yeast cytoplasmic LysRS. It encodes a polypeptide with an amino-terminal extension composed of about 60-70 amino acid residues, compared to its prokaryotic counterpart. This segment, containing 13 lysine residues, is removed upon proteolytic treatment of the native enzyme. The aim of the present study was to probe in vivo the significance of this amino-terminal extension. We have constructed derivatives of the KRS1 gene, encoding enzymes lacking 58 or 69 amino-terminal residues and, by site-directed mutagenesis, we have changed four or eight lysine residues from the amino-terminal segment of LysRS into glutamic acids. Engineered proteins were expressed in vivo after replacement of the wild-type KRS1 allele. The mutant enzymes displayed reduced specific activities (2-100-fold). A series of carboxy-terminal deletions, encompassing 3, 10 or 15 amino acids, were introduced into the LysRS mutants with modified amino-terminal extensions. The removal of three residues led to a 2-7-fold increase in the specific activity of the mutant enzymes. This partial compensatory effect suggests that interactions between the two extreme regions of yeast LysRS are required for a proper conformation of the native enzyme. All KRS1 derivatives were able to sustain growth of yeast cells, although the mutant cell lines displaying a low LysRS activity grew more slowly. The expression, as single-copy genes, of mutant enzymes with a complete deletion of the amino-terminal extension or with four Lys----Glu mutations, that displayed specific activities close to that of the wild-type LysRS, had no discernable effect on cell growth. We conclude that the polycationic extensions of eukaryotic aminoacyl-tRNA synthetases are dispensable, in vivo, for aminoacylation activities. The results are discussed in relation to the triggering role in in situ compartmentalization of protein synthesis that has been ascribed to the polypeptide-chain extensions that characterize most, if not all, eukaryotic aminoacyl-tRNA synthetases.     

Evaluation of structural and evolutionary contributions to deleterious mutation prediction

Methods for automated prediction of deleterious protein mutations have utilized both structural and evolutionary information but the relative contribution of these two factors remains unclear. To address this, we have used a variety of structural and evolutionary features to create simple deleterious mutation models that have been tested on both experimental mutagenesis and human allele data. We find that the most accurate predictions are obtained using a solvent-accessibility term, the C(beta) density, and a score derived from homologous sequences, SIFT. A classification tree using these two features has a cross-validated prediction error of 20.5% on an experimental mutagenesis test set when the prior probability for deleterious and neutral cases is equal, whereas this prediction error is 28.8% and 22.2% using either the C(beta) density or SIFT alone. The improvement imparted by structure increases when fewer homologs are available: when restricted to three homologs the prediction error improves from 26.9% using SIFT alone to 22.4% using SIFT and the C(beta) density, or 24.8% using SIFT and a noisy C(beta) density term approximating the inaccuracy of ab initio structures modeled by the Rosetta method. We conclude that methods for deleterious mutation prediction should include structural information when fewer than five to ten homologs are available, and that ab initio predicted structures may soon be useful in such cases when high-resolution structures are unavailable.