Mutational spectrum in the recent human genome inferred by single nucleotide polymorphisms

So far, there is no genome-wide estimation of the mutational spectrum in humans. In this study, we systematically examined the directionality of the point mutations and maintenance of GC content in the human genome using approximately 1.8 million high-quality human single nucleotide polymorphisms and their ancestral sequences in chimpanzees. The frequency of C-->T (G-->A) changes was the highest among all mutation types and the frequency of each type of transition was approximately fourfold that of each type of transversion. In intergenic regions, when the GC content increased, the frequency of changes from G or C increased. In exons, the frequency of G:C-->A:T was the highest among the genomic categories and contributed mainly by the frequent mutations at the CpG sites. In contrast, mutations at the CpG sites, or CpG-->TpG/CpA mutations, occurred less frequently in the CpG islands relative to intergenic regions with similar GC content. Our results suggest that the GC content is overall not in equilibrium in the human genome, with a trend toward shifting the human genome to be AT rich and shifting the GC content of a region to approach the genome average. Our results, which differ from previous estimates based on limited loci or on the rodent lineage, provide the first representative and reliable mutational spectrum in the recent human genome and categorized genomic regions.     

High prevalence of SLC6A8 deficiency in X-linked mental retardation

A novel X-linked mental retardation (XLMR) syndrome was recently identified, resulting from creatine deficiency in the brain caused by mutations in the creatine transporter gene, SLC6A8. We have studied the prevalence of SLC6A8 mutations in a panel of 290 patients with nonsyndromic XLMR archived by the European XLMR Consortium. The full-length open reading frame and splice sites of the SLC6A8 gene were investigated by DNA sequence analysis. Six pathogenic mutations, of which five were novel, were identified in a total of 288 patients with XLMR, showing a prevalence of at least 2.1% (6/288). The novel pathogenic mutations are a nonsense mutation (p.Y317X) and four missense mutations. Three missense mutations (p.G87R, p.P390L, and p.P554L) were concluded to be pathogenic on the basis of conservation, segregation, chemical properties of the residues involved, as well as the absence of these and any other missense mutation in 276 controls. For the p.C337W mutation, additional material was available to biochemically prove (i.e., by increased urinary creatine : creatinine ratio) pathogenicity. In addition, we found nine novel polymorphisms (IVS1+26G-->A, IVS7+37G-->A, IVS7+87A-->G, IVS7-35G-->A, IVS12-3C-->T, IVS2+88G-->C, IVS9-36G-->A, IVS12-82G-->C, and p.Y498) that were present in the XLMR panel and/or in the control panel. Two missense variants (p.V629I and p.M560V) that were not highly conserved and were not associated with increased creatine : creatinine ratio, one translational silent variant (p.L472), and 10 intervening sequence variants or untranslated region variants (IVS6+9C-->T, IVS7-151_152delGA, IVS7-99C-->A, IVS8-35G-->A, IVS8+28C-->T, IVS10-18C-->T, IVS11+21G-->A, IVS12+15C-->T, *207G-->C, IVS12+32C-->A) were found only in the XLMR panel but should be considered as unclassified variants or as a polymorphism (p.M560V). Our data indicate that the frequency of SLC6A8 mutations in the XLMR population is close to that of CGG expansions in FMR1, the gene responsible for fragile-X syndrome.     

Determination of the gene structure of human oligophrenin-1 and identification of three novel polymorphisms by screening of DNA from 164 patients with non-specific X-linked mental retardation

We have recently shown that mutations in oligophrenin-1 (OPHN1) are responsible for non-specific X-linked mental retardation (MRX). The structure of the gene encoding the OPHN1 protein was determined by isolation of genomic DNA clones from the human cosmid library. Genomic fragments containing exons were sequenced, and the sequences of the exons and flanking introns were defined. Knowledge of the genomic structure of the OPHN1 gene, which spans at least 500 kb and consists of 25 exons, will facilitate the search for additional mutations in OPHN1. OPHN1 was screened for mutations in 164 subjects with non-specific mental retardation. Three nucleotide substitutions were identified, one of which was a silent mutation in the codon threonine 301 at position 903 (G-->C). The other substitutions were located in exon 2, a G-->A substitution at position 133 (A45T), and in exon 10, a C-->T substitution at position 902 (T301M), but these are common polymorphisms rather than disease-causing mutations.     

p53 gene mutations in neoplastic transformation of C3H 10T1/2 and severe combined immunodeficiency fibroblasts.

The relevance of p53 mutations to the neoplastic malignant transformation of rodent fibroblasts by genotoxic physical and chemical agents is not clear. In the present study, we investigated p53 mutations (in exons 5-8) in non-transformed and neoplastically transformed C3H 10T1/2 and severe combined immunodeficiency (SCID) cells. No p53 mutations were detected in 15 neoplastically transformed (two spontaneous, one 3-methylcholanthrene-induced, seven gamma-ray-induced and five 'hot particle'-induced) and two non-transformed 10T1/2 cells. Wild-type p53 gene was also detected in all non-transformed (immortalized) SCID cell lines analyzed (four lines) whereas all three neoplastically transformed (two spontaneous, one gamma-ray-induced) cell lines displayed missense mutations in the p53 gene. These mutations were all transitions: A > G in codon 123, G > A in codon 152, and C > T in codon 238. We conclude that mutation in the p53 gene appears to be an infrequent event in 10T1/2 cells regardless of the transforming agent, but a frequent event in the neoplastic transformation of immortalized SCID cells. Non-transformed SCID cells are deficient in repair of DNA double-strand breaks, and neoplastically transformed cells are assumed to be deficient as well.     

Methylglyoxal induces G:C to C:G and G:C to T:A transversions in the supF gene on a shuttle vector plasmid replicated in mammalian cells

We previously reported that the majority of base-pair substitutions induced by an endogenous mutagen, methylglyoxal, were G:C-->T:A transversions and G:C-->A:T transitions in wild-type and nucleotide excision repair (NER)-deficient (uvrA or uvrC) Escherichia coli strains. To investigate the mutation spectrum of methylglyoxal in mammalian cells and to compare the spectrum with those detected in other experimental systems, we analyzed mutations in a bacterial suppressor tRNA (supF) gene in the shuttle vector plasmid pMY189. We treated pMY189 with methylglyoxal and immediately transfected it into simian COS-7 cells. The cytotoxicity and the mutation frequency (MF) increased according to the dose of methylglyoxal. In the mutants induced by methylglyoxal, multi-base deletions were predominant (50%), followed by base-pair substitutions (35%), in which 89% of the substitutions occurred at G:C sites. Among them, G:C-->C:G and G:C-->T:A transversions were predominant. The overall distribution of methylglyoxal-induced mutations detected in the supF gene was different from that for the spontaneous mutations. These results suggest that methylglyoxal may take part in causing G:C-->C:G and G:C-->T:A transversions in vivo.     

Conserved quartets near 5' intron junctions in primate nuclear pre-mRNA

Analysis of a 1000 nucleotide span around 664 primate 5' exon/intron junctions revealed frequent recurrences of G-rich runs downstream of the 5' splice sites. In particular, AGGG, GGGA, GGGG, GGGT and TGGG are frequent at this site. Some C-rich quarters are frequent upstream of the 5' splice site. Similar behaviour of these G- and C-rich quartets is indicated for the 587 rodent introns and for a combined eukaryotic file containing 1688 introns. (A)GGG(A) is also frequent in the introns 60 nucleotides upstream of the 3' splice site, and (A)CCC(A) is frequently found in the exons downstream of the 3' site. The same consistent behaviour of the 3' splice sites is obtained as for the 5' sites, for the primates, rodents and combined eukaryotic file. These results suggest that in addition to the well-conserved 5' and 3' splice sequences, exon as well as intron sequences may play a role in nuclear pre-mRNA splicing.     

Deficient nucleotide excision repair increases base-pair substitutions but decreases TGGC frameshifts induced by methylglyoxal in Escherichia coli.

To investigate the mutation spectrum of a well-known mutagen, methylglyoxal, and the influence of nucleotide excision repair (NER) on methylglyoxal-induced mutations, we treated wild-type and NER-deficient (uvrA or uvrC) Escherichia coli strains with methylglyoxal, and analyzed mutations in the chromosomal lacI gene. In the three strains, the cell death and the mutation frequency increased according to the dose of methylglyoxal added to the culture medium. The frequencies of methylglyoxal-induced base-pair substitutions were higher in the NER-deficient strains than in the wild-type strain, in the presence and absence of mucAB gene. Paradoxically, the frequency of methylglyoxal-induced TGGC frameshifts was higher in the wild-type strain than in the NER-deficient strains. When the methylglyoxal-induced mutation spectra in the presence and absence of mucAB gene are compared, the ratios of base-pair substitutions to frameshifts were increased by the effects of mucAB gene. In the three strains, more than 75% of the base-pair substitutions occurred at G:C sites, independent of the mucAB gene. When the mucAB gene was present, G:C-->T:A transversions were predominant, followed by G:C-->A:T transitions. When the mucAB gene was absent, the predominant mutations differed in the three strains: in the wild-type and uvrC strains, G:C-->A:T transitions were predominant, followed by G:C-->T:A transversions, while in the uvrA strains, G:C-->T:A transversions were predominant, followed by G:C-->A:T transitions. These results suggest that NER may be involved in both the repair and the fixation of methylglyoxal-induced mutations.     

Cytogenetic damage and ras p21 oncoprotein levels from patients with chronic obstructive pulmonary disease (COPD), untreated lung cancer and healthy controls

Human population has been continually exposed to benzene which is present in our environment as an essential component of petroleum. p-Benzoquinone (p-BQ) is one of the benzene metabolites and is thought to be an ultimate toxic or carcinogenic substance. For molecular analysis of carcinogen-induced mutations in mouse cells, we constructed a new shuttle vector plasmid pNY200 that has supF gene as a target of the mutations and replicates in mouse and in Escherichia coli cells. In p-BQ-treated pNY200 propagated in mouse cells, base substitutions were induced predominantly at G:C sites, and the major mutation was G:C-->A:T transition. Many tandem base substitutions were also induced at CC:GG sequences. By a postlabeling analysis and a polymerase stop assay, we confirmed that p-BQ adducts formed in DNA and mutation sites roughly correspond to the sites where the adducts were formed. Comparing data of pNY200 in mouse cells with those of the similar shuttle vector plasmid pMY189 in human cells should be important for extrapolation of data from mouse to human, because carcinogenicity of chemicals is tested in mice.     

Specificity of spontaneous mutations induced in mutA mutator cells

Escherichia coli cells expressing the mutA allele of a glyV (glycine tRNA) gene express a strong mutator phenotype. The mutA allele differs from the wild type glyV gene by a base substitution in the anticodon such that the resulting tRNA misreads certain aspartate codons as glycine, resulting in random, low-level Asp-->Gly substitutions in proteins. Subsequent work showed that many types of mistranslation can lead to a very similar phenotype, named TSM for translational stress-induced mutagenesis. Here, we have determined the specificity of forward mutations occurring in the lacI gene in mutA cells as well as in wild type cells. Our results show that in comparison to wild type cells, base substitutions are elevated 23-fold in mutA cells, as against a eight-fold increase in insertions and a five-fold increase in deletions. Among base substitutions, transitions are elevated 13-fold, with both G:C-->A:T and A:T-->G:C mutations showing roughly similar increases. Transversions are elevated 35-fold, with G:C-->T:A, G:C-->C:G and A:T-->C:G elevated 28-, 13- and 27-fold, respectively. A:T-->T:A mutations increase a striking 348-fold over parental cells, with most occurring at two hotspot sequences that share the G:C-rich sequence 5'-CCGCGTGG. The increase in transversion mutations is similar to that observed in cells defective for dnaQ, the gene encoding the proofreading function of DNA polymerase III. In particular, the relative proportions and sites of occurrence of A:T-->T:A transversions are similar in mutA and mutD5 (an allele of dnaQ) cells. Interestingly, transversions are also the predominant base substitutions induced in dnaE173 cells in which a missense mutation in the alpha subunit of polymerase III abolishes proofreading without affecting the 3'-->5' exonuclease activity of the epsilon subunit.     

Mutation spectra of chemical mutagens determined by Lac+ reversion assay with Escherichia coli WP3101P-WP3106P tester strains

We previously reported the development of mutation-specific Escherichia coli B tester strains WP3101 to WP3106 from strain WP2uvrA. In this study we constructed their pKM101-containing derivatives WP3101P to WP3106P, and further isolated their rfa derivatives WP4101-WP4106 and WP4101P-WP4106P. The six kinds of F' plasmids (lacI-, lacZ-, proAB+), each of which carries a different lacZ allele, contained in the above strains were originally derived from E. coli K-12 strains CC101-CC106. All the tester strains show Lac- and Trp- phenotype. Assays for transitions and transversions are based upon Lac+ reversion of a specific mutation located within the lacZ gene on an F' plasmid. The trpE65(ochre) allele in the same strains enables them to be used for Trp+ reversion assays as well. In the present paper, we evaluated the sensitivity, specificity, and usefulness of the newly developed tester strains. Strains WP3101P-WP3106P were highly sensitive to determine mutational profile of heterocyclic amines with S9 mix-mediated metabolic activation and most of the oxidative mutagens and free radical generators tested. Every type of base-pair substitutions induced by 2-amino-3,4-dimethylimidazo[4,5-f]quinoline (MeIQ) or 5-diazouracil were detected in strains WP3101P-WP3106P, while A:T-->C:G and G:C-->A:T mutations induced by MeIQ, and A:T-->C:G, G:C-->A:T, and G:C-->C:G by 5-diazouracil were not detected in pKM101-free tester strains. In pKM101-carrying strains, cumene hydroperoxide induced all types of base substitutions, while formaldehyde preferentially induced G:C-->T:A transversions. Phenazine methosulfate induced predominantly G:C-->A:T transitions and G:C-->T:A transversions, while H2O2 induced predominantly G:C-->T:A and A:T-->T:A transversions. Introduction of the rfa mutation considerably enhanced sensitivity to bulky mutagens such as polycyclic aromatic compounds. All six possible base substitutions induced by 9, 10-dimethyl-1,2-benzanthracene (DMBA) were detected in tester strains WP4101P-WP4106P. In conclusion, our tester strains WP3101P-WP3106P and WP4101P-WP4106P permitted rapid and simple detection of specific mutations induced by variety of mutagens.     

Comparison of the mutational spectra of the lacZ transgene in four organs of the MutaMouse treated with benzo[a]pyrene: target organ specificity

We recently demonstrated that not all organs with a high rate of induction of mutation in the lacZ transgene develop tumors in the lambdalacZ transgenic mice (MutaMouse) used for a long-term carcinogenicity study with benzo[a]pyrene (BP). To better understand the role of chemical-induced in vivo mutations in carcinogenesis, we compared the mutational spectra of the lacZ transgene in four organs of the MutaMouse obtained 2 weeks after five daily consecutive oral treatments with 125 mg/kg/day BP. lacZ transgenes were analyzed in two target organs (forestomach and spleen) and two non-target organs (colon and glandular stomach) for BP-induced carcinogenesis in MutaMouse, and all of these organs were highly mutated in the lacZ transgene. The sequence data showed similar mutational spectra of the lacZ transgene between the two target organs; the predominant mutations were G:C-->T:A transversions (55% and 50% for forestomach and spleen, respectively), followed by deletions (20% and 21% for forestomach and spleen, respectively) mainly at G:C site. The frequent G:C-->T:A transversions are consistent with reports of the mutational spectra produced in the p53 gene in tumors generated in rats and mice exposed to BP. In contrast, the mutational spectra of the lacZ transgene in the two non-target organs are different from those in the target organs, and are also suggested to differ from one another. These findings suggest an organ/tissue-specific mechanism of mutagenesis.     

An in vivo approach to identifying sequence context of 8-oxoguanine mutagenesis

Base substitution mutations are not distributed randomly in that most are located at a few specific hotspots sites. We have been studying 7,8-dihydro-8-oxoguanine mutagenesis in Escherichia coli in the supF gene carried in a plasmid. Among hotspots, guanine within the 5'-AGA-3' located in the anticodon site was susceptible to the induction of G:C-->T:A transversion. In this study, we constructed variants of the supF gene in which the hotspot 5'-AGA-3' was modified to 5'-AGT-3', 5'-AGG-3' and 5'-AGC-3' to determine the influence of 3' neighboring base on G:C-->T:A mutational activity. Using these variant supF genes propagated in a 7,8-dihydro-8-oxoguanine repair-deficient host, we found that guanine within 5'-AGA-3' and 5'-AGG-3' produce G:C-->T:A, but guanine within 5'-AGT-3' and 5'-AGC-3' reduce the formation of G:C-->T:A. These changes were thus due to the effect of sequence context on the efficiency of mutation formation at the sites of 7,8-dihydro-8-oxoguanine. We also observed a longer range base-pair effect on hotspot formation.     

Analysis of alkaptonuria (AKU) mutations and polymorphisms reveals that the CCC sequence motif is a mutational hot spot in the homogentisate 1,2 dioxygenase gene (HGO)

We recently showed that alkaptonuria (AKU) is caused by loss-of-function mutations in the homogentisate 1,2 dioxygenase gene (HGO). Herein we describe haplotype and mutational analyses of HGO in seven new AKU pedigrees. These analyses identified two novel single-nucleotide polymorphisms (INV4+31A-->G and INV11+18A-->G) and six novel AKU mutations (INV1-1G-->A, W60G, Y62C, A122D, P230T, and D291E), which further illustrates the remarkable allelic heterogeneity found in AKU. Reexamination of all 29 mutations and polymorphisms thus far described in HGO shows that these nucleotide changes are not randomly distributed; the CCC sequence motif and its inverted complement, GGG, are preferentially mutated. These analyses also demonstrated that the nucleotide substitutions in HGO do not involve CpG dinucleotides, which illustrates important differences between HGO and other genes for the occurrence of mutation at specific short-sequence motifs. Because the CCC sequence motifs comprise a significant proportion (34.5%) of all mutated bases that have been observed in HGO, we conclude that the CCC triplet is a mutational hot spot in HGO.     

Sequence specificity of Hprt lymphocyte mutation in rats fed the hepatocarcinogen 2-acetylaminofluorene

Rats fed the hepatocarcinogen 2-acetylaminofluorene (2-AAF) have a low, but significantly increased, frequency of lymphocyte Hprt mutants. In this study, mutants from 2-AAF-fed and control F344 rats were examined for mutations in the Hprt gene in order to determine if the 2-AAF treatment resulted in an agent-specific mutation profile. The most common mutation from 2-AAF-treated rats was G:C-->T:A transversion (32% of all mutations) followed by 1-basepair (bp) deletion (19%); there were very few (5%) G:C-->A:T transitions. Among mutations from control rats, G:C-->A:T transition was the most common (43%), and there were very few G:C-->T:A transversions (5%) and no 1-bp deletions. The profile of mutations from 2-AAF-fed rats was significantly different from control rats (P = 0.003) and was consistent with the types of mutations produced by 2-AAF in vitro. The results of this study indicate that even weak mutational responses in the lymphocyte Hprt assay are capable of producing mutation profiles that reflect the DNA damage inducing them.     

Conserved signals around the 5' splice sites in eukaryotic nuclear precursor mRNAs: G-runs are frequent in the introns and C in the exons near both 5' and 3' splice sites

It is known that the GT doublet is well conserved at the 5' exon/intron splice junction and is frequently embedded in the AGGT quartet. Although only the underlined G is invariable, splicing and ligation are accurately executed. In this work we search for additional conserved potential signals which may aid in 5' splice site recognition. Extensive searches which are not limited to a preconceived consensus sequence are carried out. We investigate the distributions of the 256 quartets in a 1000 nucleotide span around the 5' splice sites in approximately 1700 eukaryotic nuclear precursor mRNAs. Several potential signals are noted. Of particular interest are quartets containing runs of G, e.g., G4, G3T, G3C, G3A and AG3 in the intron immediately downstream and some C-containing quartets in the exon upstream of the 5' splice site. In an analogous calculation, (A)GGG(A) has also been found to be frequent in the intron, 60 nucleotides upstream and (A)CCC(A) in the exon downstream of the 3' splice site. These results are consistent with the recent indications that exon sequences may play a role in efficient splicing. Some models are proposed.     

Mismatch repair deficient human cells: spontaneous and MNNG-induced mutational spectra in the HPRT gene

We have determined both the spontaneous and N-methyl-N′-nitro-N-nitrosoguanidine (MNNG)-induced mutational spectra in the HPRT gene of human cells (MT1) defective in the mismatch repair gene hMSH6 (GTBP). Eight of nine exons and nine of sixteen intronic flanking sequences were scanned, encompassing >900 bp of the HPRT gene. Mutant hotspots were detected and separated by differences in their melting temperatures using constant denaturant capillary electrophoresis (CDCE) or denaturing gradient gel electrophoresis (DGGE).

A key finding of this work is that a high proportion of all HPRT inactivating mutations is represented by a small number of hotspots distributed over the exons and mRNA splice sites. Thirteen spontaneous hotspots and sixteen MNNG-induced hotspots accounted for 55% and 48% of all 6TG^R point mutations, respectively. MNNG-induced hotspots were predominantly G:C→A:T transitions. The spontaneous spectrum of cells deficient in hMSH6 contained transversions (A:T→T:A, G:C→T:A, A:T→C:G), transitions (A:T→G:C), a plus-one insertion, and a minus-one deletion. Curiously, G:C→A:T transitions, which dominate human germinal and somatic point mutations were absent from the spontaneous hMSH6 spectra.     

Mutations induced by glyoxal and methylglyoxal in mammalian cells.

To investigate the mutation spectra of glyoxal and methylglyoxal in mammalian cells, we analyzed mutations in a bacterial suppressor tRNA (supF) gene in the shuttle vector plasmid pMY189. The cytotoxicity and the mutation frequency increased according to the doses of glyoxal and methylglyoxal. The majority of glyoxal-induced mutations (65%) were base-pair substitutions, in which G:C-->C:G transversions were predominant. In the mutants induced by methylglyoxal, multi-base deletions were predominant (50%), followed by base-pair substitutions (35%), in which G:C-->C:G and G:C-->T:A transversions were predominant.