Evidence that proximal multiple mutations in Big Blue transgenic mice are dependent events

A Framework for Investigating Rules of Life by Establishing Zones of Influence

The incredible complexity of biological processes across temporal and spatial scales hampers defining common underlying mechanisms driving the patterns of life. However, recent advances in sequencing, big data analysis, machine learning, and molecular dynamics simulation have renewed the hope and urgency of finding potential hidden rules of life. There currently exists no framework to develop such synoptic investigations. Some efforts aim to identify unifying rules of life across hierarchical levels of time, space, and biological organization, but not all phenomena occur across all the levels of these hierarchies. Instead of identifying the same parameters and rules across levels, we posit that each level of a temporal and spatial scale and each level of biological organization has unique parameters and rules that may or may not predict outcomes in neighboring levels. We define this neighborhood, or the set of levels, across which a rule functions as the zone of influence. Here, we introduce the zone of influence framework and explain using three examples: (Smocovitis, 1992) randomness in biology, where we use a Poisson process to describe processes from protein dynamics to DNA mutations to gene expressions, (Leroi, 2014) island biogeography, and (Gropp, 2016) animal coloration. The zone of influence framework may enable researchers to identify which levels are worth investigating for a particular phenomenon and reframe the narrative of searching for a unifying rule of life to the investigation of how, when, and where various rules of life operate.

Evolutionary switches between two serine codon sets are driven by selection

Serine is the only amino acid that is encoded by two disjoint codon sets so that a tandem substitution of two nucleotides is required to switch between the two sets. Previously published evidence suggests that, for the most evolutionarily conserved serines, the codon set switch occurs by simultaneous substitution of two nucleotides. Here we report a genome-wide reconstruction of the evolution of serine codons in triplets of closely related species from diverse prokaryotes and eukaryotes. The results indicate that the great majority of codon set switches proceed by two consecutive nucleotide substitutions, via a threonine or cysteine intermediate, and are driven by selection. These findings imply a strong pressure of purifying selection in protein evolution, which in the case of serine codon set switches occurs via an initial deleterious substitution quickly followed by a second, compensatory substitution. The result is frequent reversal of amino acid replacements and, at short evolutionary distances, pervasive homoplasy.

Clusters of Multiple Mutations: Incidence and Molecular Mechanisms

It has been long understood that mutation distribution is not completely random across genomic space and in time. Indeed, recent surprising discoveries identified multiple simultaneous mutations occurring in tiny regions within chromosomes while the rest of the genome remains relatively mutation-free. Mechanistic elucidation of these phenomena, called mutation showers, mutation clusters, or kataegis, in parallel with findings of abundant clustered mutagenesis in cancer genomes, is ongoing. So far, the combination of factors most important for clustered mutagenesis is the induction of DNA lesions within unusually long and persistent single-strand DNA intermediates. In addition to being a fascinating phenomenon, clustered mutagenesis also became an indispensable tool for identifying a previously unrecognized major source of mutation in cancer, APOBEC cytidine deaminases. Future research on clustered mutagenesis may shed light onto important mechanistic details of genome maintenance, with potentially profound implications for human health.

Concurrent nucleotide substitution mutations in the human genome are characterized by a significantly decreased transition/transversion ratio

There is accumulating evidence that the number of multiple-nucleotide substitutions (MNS) occurring in closely spaced sites in eukaryotic genomes is significantly higher than would be predicted from the random accumulation of independently generated single-nucleotide substitutions (SNS). Although this excess can in principle be accounted for by the concept of transient hypermutability, a general mutational signature of concurrent MNS mutations has not so far been evident. Employing a dataset (N = 449) of "concurrent" double MNS mutations causing human inherited disease, we have identified just such a mutational signature: concurrently generated double MNS mutations exhibit a >twofold lower transition/transversion ratio (termed RTs/Tv ) than independently generated de novo SNS mutations (<0.80 vs. 2.10; P = 2.69 × 10(-14) ). We replicated this novel finding through a similar analysis employing two double MNS variant datasets with differing abundances of concurrent events (150,521 variants with both substitutions on the same haplotypic lineage vs. 94,875 variants whose component substitutions were on different haplotypic lineages) plus 5,430,874 SNS variants, all being derived from the whole-genome sequencing of seven Chinese individuals. Evaluation of the newly observed mutational signature in diverse contexts provides solid support for the postulated role of translesion synthesis DNA polymerases in transient hypermutability.

Double SMCHD1 variants in FSHD2: the synergistic effect of two SMCHD1 variants on D4Z4 hypomethylation and disease penetrance in FSHD2

Facioscapulohumeral muscular dystrophy (FSHD) predominantly affects the muscles in the face, trunk and upper extremities and is marked by large clinical variability in disease onset and progression. FSHD is associated with partial chromatin relaxation of the D4Z4 repeat array on chromosome 4 and the somatic expression of the D4Z4 encoded DUX4 gene. The most common form, FSHD1, is caused by a contraction of the D4Z4 repeat array on chromosome 4 to a size of 1-10 units. FSHD2, the less common form of FSHD, is most often caused by heterozygous variants in the chromatin modifier SMCHD1, which is involved in the maintenance of D4Z4 methylation. We identified three families in which the proband carries two potentially damaging SMCHD1 variants. We investigated whether these variants were located in cis or in trans and determined their functional consequences by detailed clinical information and D4Z4 methylation studies. In the first family, both variants in trans were shown to act synergistically on D4Z4 hypomethylation and disease penetrance, in the second family both SMCHD1 function-affecting variants were located in cis while in the third family one of the two variants did not affect function. This study demonstrates that having two SMCHD1 missense variants that affect function is compatible with life in males and females, which is remarkable considering its role in X inactivation in mice. The study also highlights the variability in SMCHD1 variants underlying FSHD2 and the predictive value of D4Z4 methylation analysis in determining the functional consequences of SMCHD1 variants of unknown significance.

Oncogenic potential is related to activating effect of cancer single and double somatic mutations in receptor tyrosine kinases

Aberrant activation of receptor tyrosine kinases (RTKs) is a common feature of many cancer cells. It was previously suggested that the mechanisms of kinase activation in cancer might be linked to transitions between active and inactive states. Here, we estimate the effects of single and double cancer mutations on the stability of active and inactive states of the kinase domains from different RTKs. We show that singleton cancer mutations destabilize active and inactive states; however, inactive states are destabilized more than the active ones, leading to kinase activation. We show that there exists a relationship between the estimate of oncogenic potential of cancer mutation and kinase activation. Namely, more frequent mutations have a higher activating effect, which might allow us to predict the activating effect of the mutations from the mutation spectra. Independent evolutionary analysis of mutation spectra complements this observation and finds the same frequency threshold defining mutation hotspots. We analyze double mutations and report a positive epistasis and additional advantage of doublets with respect to cancer cell fitness. The activation mechanisms of double mutations differ from those of single mutations and double mutation spectrum is found to be dissimilar to the mutation spectrum of singletons.

Transient hypermutability, chromothripsis and replication-based mechanisms in the generation of concurrent clustered mutations

Clustered mutations may be broadly defined as the presence of two or more mutations within a spatially localized genomic region on a single chromosome. Known instances vary in terms of both the number and type of the component mutations, ranging from two closely spaced point mutations to tens or even hundreds of genomic rearrangements. Although clustered mutations can represent the observable net result of independent lesions sequentially acquired over multiple cell cycles, they can also be generated in a simultaneous or quasi-simultaneous manner within a single cell cycle. This review focuses on those mechanisms known to underlie the latter type. Both gene conversion and transient hypermutability are capable of generating closely spaced multiple mutations. However, a recently described phenomenon in human cancer cells, known as 'chromothripsis', has provided convincing evidence that tens to hundreds of genomic rearrangements can sometimes be generated simultaneously via a single catastrophic event. The distinctive genomic features observed in the derivative chromosomes, together with the highly characteristic junction sequences, point to non-homologous end joining (NHEJ) as being the likely underlying mutational mechanism. By contrast, replication-based mechanisms such as microhomology-mediated break-induced replication (MMBIR) which involves serial replication slippage or serial template switching probably account for those complex genomic rearrangements that comprise multiple duplications and/or triplications.

Gene conversion in human genetic disease

Gene conversion is a specific type of homologous recombination that involves the unidirectional transfer of genetic material from a ‘donor’ sequence to a highly homologous ‘acceptor’. We have recently reviewed the molecular mechanisms underlying gene conversion, explored the key part that this process has played in fashioning extant human genes, and performed a meta-analysis of gene-conversion events known to have caused human genetic disease. Here we shall briefly summarize some of the latest developments in the study of pathogenic gene conversion events, including (i) the emerging idea of minimal efficient sequence homology (MESH) for homologous recombination, (ii) the local DNA sequence features that appear to predispose to gene conversion, (iii) a mechanistic comparison of gene conversion and transient hypermutability, and (iv) recently reported examples of pathogenic gene conversion events.

HO* radicals induce an unexpected high proportion of tandem base lesions refractory to repair by DNA glycosylases

Reaction of HO(*) radicals with double-stranded calf thymus DNA produces high levels of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodGuo) and, to a minor extent, 8-oxo-7,8-dihydro-2'-deoxyadenosine (8-oxodAdo). Formation of the hydroxylated purine lesions is explained by addition of HO(*) to the C8 position of the purine moiety. It has been reported that tandem lesions containing a formylamine residue neighboring 8-oxodGuo could be produced through addition of a transiently generated pyrimidine peroxyl radical onto the C8 of an adjacent purine base. Formation of such tandem lesions accounted for approximately 10% of the total 8-oxodGuo. In the present work we show that addition of HO(*) onto the C8 of purine accounts for only approximately 5% of the generated 8-oxodGuo. About 50% of the 8-hydroxylated purine lesions, including 8-oxodGuo and 8-oxodAdo, are involved in tandem damage and are produced by peroxyl addition onto the C8 of a vicinal purine base. In addition, the remaining 45% of the 8-oxodGuo are produced by an electron transfer reaction, providing an explanation for the higher yield of formation of 8-oxodGuo compared to 8-oxodAdo. Interestingly, we show that >40% of the 8-oxodGuo involved in tandem lesions is refractory to excision by DNA glycosylases. Altogether our results demonstrate that, subsequently to a single oxidation event, peroxidation reactions significantly increase the yield of formation of hydroxylated purine modifications, generating a high proportion of tandem lesions partly refractory to base excision repair.

Closely spaced multiple mutations as potential signatures of transient hypermutability in human genes

Data from diverse organisms suggests that transient hypermutability is a general mutational mechanism with the potential to generate multiple synchronous mutations, a phenomenon probably best exemplified by closely spaced multiple mutations (CSMMs). Here we have attempted to extend the concept of transient hypermutability from somatic cells to the germline, using human inherited disease-causing multiple mutations as a model system. Employing stringent criteria for data inclusion, we have retrospectively identified numerous potential examples of pathogenic CSMMs that exhibit marked similarities to the CSMMs reported in other systems. These examples include (1) eight multiple mutations, each comprising three or more components within a sequence tract of <100 bp; (2) three possible instances of "mutation showers"; and (3) numerous highly informative "homocoordinate" mutations. Using the proportion of CpG substitution as a crude indicator of the relative likelihood of transient hypermutability, we present evidence to suggest that CSMMs comprising at least one pair of mutations separated by < or =100 bp may constitute signatures of transient hypermutability in human genes. Although this analysis extends the generality of the concept of transient hypermutability and provides new insights into what may be considered a novel mechanism of mutagenesis underlying human inherited disease, it has raised serious concerns regarding current practices in mutation screening.

The rate of mutation of a single gene

The rate of mutation refers to the probability that a unit length of DNA (generally a base pair) mutates with time. Fluctuation analysis or mutant accumulation assays applied to phenotypic changes measure mutation rates of cells. However, only a few phenotypic changes indicative of mutations are known thus limiting the analysis to those rare genes. Direct sequencing overcomes the limitations imposed by phenotypic analysis but is limited by the extensive number of clones or cells that have to be analyzed in fluctuation or mutant accumulation assays. We propose a strategy to determine the rate of mutation of a gene by limited direct sequencing of a few single cells of a defined lineage. To accomplish this, we determined the average number of mutations per position in each DNA length sequenced from the proportion of the non-mutated positions, according to the Poisson process and/or the Taylor series. Measuring the rate of mutation by direct sequencing of genes does not require ascertaining a phenotype and can be applied to any area of the genome in a cell. The approach avoids fluctuation errors.

A novel germline mutation in Big Blue mice

The Big Blue lacI mutation detection assay is well validated and has permitted detailed analysis of spontaneous mutations in individual tissues over the lifespan of the mouse. In a recent assay of spontaneous mutations, a novel lacI mutation (C354T) recurred in six of seven mutants with a second mutation. The frequency of spontaneous doublets (mutants with two nontandem mutations) was elevated 2.7-fold over that previously reported (Hill KA et al., [2004b]: Mutat Res 554:223-240) for normal tissues (6.3 x 10(-7) herein vs. 2.36 x 10(-7)). The average spacing between mutations in the doublets (237 bp) was greater than previously reported for spontaneous doublets. The frequency of C354T as a "hitchhiker" mutation in doublets was consistent with a germline mutation in one of 38 mutation targets in the Big Blue mouse genome. C354T is a missense mutation at a CpG dinucleotide producing a conservative amino acid change (Ala109Val) and a very light blue mutant phenotype. Mutant phenotypes of doublets with C354T were typical of the second mutation. C354T was observed in mutants from five tissues of five Big Blue mice. A bidirectional-PCR amplification of specific alleles (Bi-PASA) assay detected C354T in genomic DNA from multiple tissues of five Big Blue mice. These observations are consistent with a novel lacI C354T germline mutation in Big Blue mice that introduces a significant artifact in the analysis of spontaneous mutations. This finding reiterates the importance of identifying all mutations and examining new mutations in the context of our increasingly detailed knowledge of features of spontaneous murine mutations.

Epidemiology of doublet/multiplet mutations in lung cancers: evidence that a subset arises by chronocoordinate events

Background

Evidence strongly suggests that spontaneous doublet mutations in normal mouse tissues generally arise from chronocoordinate events. These chronocoordinate mutations sometimes reflect “mutation showers”, which are multiple chronocoordinate mutations spanning many kilobases. However, little is known about mutagenesis of doublet and multiplet mutations (domuplets) in human cancer. Lung cancer accounts for about 25% of all cancer deaths. Herein, we analyze the epidemiology of domuplets in the EGFR and TP53 genes in lung cancer. The EGFR gene is an oncogene in which doublets are generally driver plus driver mutations, while the TP53 gene is a tumor suppressor gene with a more typical situation in which doublets derive from a driver and passenger mutation.

Methodology/Principal Findings

EGFR mutations identified by sequencing were collected from 66 published papers and our updated EGFR mutation database (www.egfr.org). TP53 mutations were collected from IARC version 12 (www-p53.iarc.fr). For EGFR and TP53 doublets, no clearly significant differences in race, ethnicity, gender and smoking status were observed. Doublets in the EGFR and TP53 genes in human lung cancer are elevated about eight- and three-fold, respectively, relative to spontaneous doublets in mouse (6% and 2.3% versus 0.7%).

Conclusions/Significance

Although no one characteristic is definitive, the aggregate properties of doublet and multiplet mutations in lung cancer are consistent with a subset derived from chronocoordinate events in the EGFR gene: i) the eight frameshift doublets (present in 0.5% of all patients with EGFR mutations) are clustered and produce a net in-frame change; ii) about 32% of doublets are very closely spaced (≤30 nt); and iii) multiplets contain two or more closely spaced mutations. TP53 mutations in lung cancer are very closely spaced (≤30 nt) in 33% of doublets, and multiplets generally contain two or more very closely spaced mutations. Work in model systems is necessary to confirm the significance of chronocoordinate events in lung and other cancers.

EGFR somatic doublets in lung cancer are frequent and generally arise from a pair of driver mutations uncommonly seen as singlet mutations: one-third of doublets occur at five pairs of amino acids

Doublet mutations in cancer are not well studied. We find that allelic somatic doublet mutations are present at high frequency in the epidermal growth factor receptor (EGFR) tyrosine kinase (TK) domain in lung cancers. When doublets from the literature are added, a total of 96 doublets are available for analysis. The frequency of doublets overall is 6%, which is sevenfold greater than that observed in normal tissue in mouse. All characterized doublets are allelic, and silent mutations occur rarely. About half of all doublets contain one or two of 12 distinct missense mutations at five amino acids: E709, G719, S768, T790 and L861. The mutations at these five amino acids are seldom reported as singlets. Moreover, when the common L858 target is included, more than one-third of EGFR doublets are one of five specific missense pairs: G719/E709, G719/S768, G719/L861, L858/E709 and L858/T790. Structure suggests function: The data imply that most EGFR doublets are NOT consistent with a 'driver and passenger' mutation mechanism. EGFR doublets are highly skewed relative to singlets, consistent with functional selection of two individually suboptimal mutations that, in combination, have enhanced oncogenic potential.

Too many mutants with multiple mutations

It has recently become clear that the classical notion of the random nature of mutation does not hold for the distribution of mutations among genes: most collections of mutants contain more isolates with two or more mutations than predicted by the mutant frequency on the assumption of a random distribution of mutations. Excesses of multiples are seen in a wide range of organisms, including riboviruses, DNA viruses, prokaryotes, yeasts, and higher eukaryotic cell lines and tissues. In addition, such excesses are produced by DNA polymerases in vitro. These "multiples" appear to be generated by transient, localized hypermutation rather than by heritable mutator mutations. The components of multiples are sometimes scattered at random and sometimes display an excess of smaller distances between mutations. As yet, almost nothing is known about the mechanisms that generate multiples, but such mutations have the capacity to accelerate those evolutionary pathways that require multiple mutations where the individual mutations are neutral or deleterious. Examples that impinge on human health may include carcinogenesis and the adaptation of microbial pathogens as they move between individual hosts.

Database of somatic mutations in EGFR with analyses revealing indel hotspots but no smoking-associated signature

We created an Epidermal Growth Factor Receptor (EGFR) Mutation Database (http://www.cityofhope.org/cmdl/egfr_db) that curates a convenient compilation of somatic EGFR mutations in non-small-cell lung cancer (NSCLC) and associated epidemiological and methodological data, including response to the tyrosine kinase inhibitors Gefitinib and Erlotinib. Herein, we analyze 809 mutations collected from 26 publications. Four super hotspots account for 70% of reported mutations while two-thirds of 131 unique mutations have been reported only once and account for only 11% of reported mutations. Consistent with strong biological selection for gain of function, the reported mutations are virtually all missense substitutions or in-frame microdeletions, microinsertions, or microindels (colocalized insertion and deletion with a net gain or loss of 1-50 nucleotides). Microdeletions and microindels are common in a region of exon 19. Microindels, which account for 8% of mutations, have smaller inserted sequences (95% are 1 to 5 bp) and are elevated 16-fold relative to mouse somatic microindels and to human germline microindels. Microdeletions/microindels are significantly more frequent in responders to Gefitinib or Erlotinib (P = 0.003). In addition, EGFR mutations in smokers do not carry signatures of mutagens in cigarette smoke. Otherwise, the mutation pattern does not differ significantly with respect to gender, age, or tumor histology. The EGFR Mutation Database is a central resource of EGFR sequence variant data for clinicians, geneticists, and other researchers. Authors are encouraged to submit new publications with EGFR sequence variants to be included in the database or to provide direct submissions via The WayStation submission and publication process (http://www.centralmutations.org).

Evidence for mutation showers

Mutants in the Big Blue transgenic mouse system show spontaneous clustered multiple mutations with unexpectedly high frequency, consistent with chronocoordinate events. We tested the prediction that the multiple mutations seen within the lacI mutation target sometimes occur in the context of chronocoordinate multiple mutations spanning multiple kilobases (mutation showers). Additional sequencing of mutants was performed in regions immediately flanking the lacI region (total of 10.7 kb). Nineteen additional mutations were found outside the lacI region ("ectomutations") from 10 mutants containing two or more lacI mutations, whereas only one ectomutation was found in 130 mutants with a single mutation (P < 0.0001). The mutation showers had an average of approximately one mutation per 3 kb. Four mutants showed closely spaced double mutations in the new sequence, and analysis of the spacing between these mutations revealed significant clustering (P = 0.0098). To determine the extent of the mutation showers, regions (8.5 kb total) remote from the lacI region (approximately 16-17 kb away) were sequenced. Only two additional ectomutations were found in these remote regions, consistent with mutation showers that generally do not extend more than approximately 30 kb. We conclude that mutation showers exist and that they constitute at least 0.2% and possibly 1% or more of mutational events observed in this system. The existence of mutation showers has implications for oncogenesis and evolution, raising the possibilities of "cancer in an instant" and "introns as sponges to reduce the deleterious impact of mutation showers."

Somatic microindels: analysis in mouse soma and comparison with the human germline

Microindels, defined as mutations that result in a colocalized microinsertion and microdeletion with a net gain or loss of between 1 and 50 nucleotides, may be an important contributor to cancer. We report the first comprehensive analysis of somatic microindels. Our large database of mutations in the lacI transgene of Big Blue((R)) mice contains 0.5% microindels, 2.8% pure microinsertions, and 11.5% pure microdeletions. There appears to be no age, gender, or tissue-type specificity in the frequency of microindels. Of the independent somatic mutations that result in a net in-frame insertion or deletion, microindels are responsible for 13% of protein expansions and 6% of protein contractions. These in-frame microindels may play a crucial role in oncogenesis and evolution via "protein tinkering" (i.e., modest expansion or contraction of proteins). Four characteristics suggest that microindels are caused by unique mechanisms, not just simple combinations of the same mechanisms that cause pure microinsertions and pure microdeletions. First, microinsertions and microdeletions commonly occur at hotspots, but none of the 30 microindels are recurrent. Second, the sizes of the deletions and insertions in microindels are larger and more varied than in pure microdeletions and pure microinsertions. Third, microinsertions overwhelmingly repeat the adjacent base (97%) while the insertions in microindels do so only infrequently (17%). Fourth, analysis of the sequence contexts of microindels is consistent with unique mechanisms including recruitment of translesion DNA synthesis polymerases. The mouse somatic microindels have characteristics similar to those of human germline microindels, consistent with similar causative mechanisms in mouse and human, and in soma and germline.

Preferential occurrence of 1-2 microindels

Microindels are unique, infrequent mutations that result in inserted and deleted sequences of different sizes (between one and 50 nucleotides) at the same nucleotide position. Little is known about the mutational mechanisms that are responsible for these mutations. From our database of 6,016 independent somatic mutational events in the lacI gene in Big Blue mice, we assembled the 30 microindels (0.5%) for analysis. Microindels with one nucleotide inserted and two nucleotides deleted (1-2 microindels) accounted for seven (23%) of the microindels observed, with the remaining microindels distributed among 21 other combinations of insertion and deletion sizes. A preferential occurrence of 1-2 microindels (20%) was also observed in human germline transmitted mutations in the Human Gene Mutation Database (HGMD). An examination of the sequence flanking the mouse 1-2 microindels did not reveal obvious site specificity or associated secondary structure. A detailed examination of 1-2 microindels did not reveal the features typical of pure microinsertion and microdeletion events, but rather suggested a unique mutational mechanism. The 1 bp insertion in 1-2 microinsertions, and pure 1 bp insertions show distinct features. The mechanism for 1-2 microindels is not obviously a simple combination of pure microinsertion and microdeletion events. The dramatic enhancement of 1-2 microindels requires explanation. We speculate that certain error-prone polymerases may be responsible for the preferential occurrence of 1-2 microindels in both somatic tissues and germ cells. It is estimated that a human adult carries roughly 400 billion somatic 1-2 microindels with the potential to predispose to cancer.

Clusters of mutations from transient hypermutability

Collections of mutants usually contain more mutants bearing multiple mutations than expected from the mutant frequency and a random distribution of mutations. This excess is seen in a variety of organisms and also after DNA synthesis in vitro. The excess is unlikely to originate in mutator mutants but rather from transient hypermutability resulting from a perturbation of one of the many transactions that maintain genetic fidelity. The multiple mutations are sometimes clustered and sometimes randomly distributed. We model some spectra as populations comprising a majority with a low mutation frequency and a minority with a high mutation frequency. In the case of mutants produced in vitro by a bacteriophage RB69 mutator DNA polymerase, mutants with two mutations are in approximately 10-fold excess and mutants with three mutations are in even greater excess. However, phenotypically undetectable mutations seen only as hitchhikers with detectable mutations are approximately 5-fold more frequent than mutants bearing detectable mutations, indicating that they arose in a subpopulation with a higher mutation frequency. Excess multiple mutations may contribute critically to carcinogenesis and to adaptive mutation, including the adaptations of pathogens as they move from host to host. In the case of the rapidly mutating riboviruses, the viral population appears to be composed of a majority with a mutation frequency substantially lower than the average and a minority with a huge mutational load.

Tissue-specific time courses of spontaneous mutation frequency and deviations in mutation pattern are observed in middle to late adulthood in Big Blue mice

To better define the time course of spontaneous mutation frequency in middle to late adulthood of the mouse, measurements were made at 10, 14, 17, 23, 25, and 30 months of age in samples of adipose tissue, liver, cerebellum (90% neurons), and the male germline (95% germ cells). A total of 46 million plaque-forming units (pfus) were screened at the six time points and 1,450 circular blue plaques were harvested and sequenced. These data improve resolution and confirm the previously observed occurrence of at least two tissue-specific profiles of spontaneous mutation frequency (elevation with age in adipose tissue and liver, and constancy with age in neurons and male germ cells), a low mutation frequency in the male germline, and a mutation pattern unchanged with age within a tissue. These findings appear to extend to very old age (30 months). Additional findings include interanimal variation in spontaneous mutation frequency is larger in adipose tissues and liver compared with neurons and male germ cells, and subtle but significant differences in the mutation pattern among tissues, consistent with a minor effect of tissue-specific metabolism. The presumptive unaltered balance of DNA damage and repair with age in the male germline has evolutionary consequences. It is of particular interest given the controversy over whether or not increasing germline mutation frequency with paternal age underlies the reports associating older males with a higher incidence of some types of genetic disease. These most detailed measurements available to date regarding the time course of spontaneous mutation frequency and pattern in individual tissues help to constrain hypotheses regarding the role of mutational mechanisms in DNA repair and aging.

Spontaneous multiple mutations show both proximal spacing consistent with chronocoordinate events and alterations with p53-deficiency

Analysis of spontaneous multiple mutations in normal and tumor cells may constrain hypotheses about the mechanisms responsible for multiple mutations and provide insight into the mutator phenotype. In a previous study, spontaneous doublets in Big Blue mice were dramatically more frequent than expected by chance and exhibited a mutation pattern similar to that observed for single mutations [Mutat. Res. 452 (2000) 219]. The spacing between mutations in doublets was generally closer than expected by chance and the distribution of mutation spacing fit an exponential, albeit with substantial scatter. We now analyze 2658 additional mutants and confirm that doublets are enhanced dramatically relative to chance expectation. The spacing, frequency and pattern of spontaneous doublets and multiplets (domuplets) are examined as a function of age, tissue type, p53-deficiency and neoplasia in the new and combined data. The new and combined data confirm that the distribution of the spacing between mutations in doublets is non-random with the mutations more closely spaced than expected by chance (P < 0.0005; combined data), consistent with temporally coordinate (chronocoordinate) events. An exponential provides an excellent fit to the distribution (R2 = 0.98) and estimates that half of doublets have mutations separated by 120 nucleotides or less (the "half-life of mutation spacing"). We make several novel observations: (i) singlets and doublets show similar overall increases in frequency with age (ii) doublet frequency may be lower in the male germline, consistent with the generally reduced mutation frequency in the male germline (iii) doublet frequencies are elevated in somatic tissues of p53-deficient mice (Li-Fraumini cancer syndrome model; P = 0.005) and (iv) doublets and singlets in tumors from p53-deficient mice have a different mutation pattern (P = 0.007). The observations are consistent with chronocoordinate occurrence of spontaneous doublets and multiplets due to a transient error-prone condition and do not suggest a major role for the recently discovered Y family of error-prone polymerases. The enhancement of doublets in p53-deficient mice may contribute to cancer risk.

In vivo mutation in gene A of splenic lymphocytes from phiX174 transgenic mice

Single-burst analysis was applied to a forward assay for gene A mutation in splenic lymphocytes of phiX174 transgenic mice for the purpose of optimizing analytical parameters for identifying in vivo mutations. The effect of varying the cutoff value for an in vivo burst on induced mutant frequency, fold increase, and the significance of the difference between control and N-ethyl-N-nitrosourea (ENU)-treated mice was calculated by two different methods. The plating density was reduced to an average of less than 10 background mutant plaques per aliquot in order to separate in vitro bursts. The spectrum of mutations contributing < 60 plaques per aliquot from control animals was not significantly different from the control spectra from E. coli or transgenic phiX174 cells in culture. The mutant spectra from ENU-treated animals was highly different between mutant bursts of > 80 plaques per aliquot compared to mutations contributing < 60 plaques per aliquot (P < 0.000001), the former fitting the spectrum expected for ENU-induced mutations. The latter spectrum was also different from control animals and E. coli (P < 0.000001), suggesting the difference was caused by ex vivo mutation. With the mutations found in this study, the total number of reported target sites for gene A is now 33. The results support the interpretation that, in contrast to results for the lacI transgene, 100% of mutants isolated in gene A from control animals and cells were fixed in E. coli. We attribute the difference between the two genes to hot-spot sites for mutation in gene A and to a testable hypothesis that the mosaic plaque assay for the lacI transgene underestimates the frequency of ex vivo mutants.

Spontaneous tandem-base mutations (TBM) show dramatic tissue, age, pattern and spectrum specificity

To supplement a previous analysis of spontaneous tandem-base mutations (TBM) in the lacI gene of Big Blue((R)) mice, 2658 additional mutants were sequenced from 13 tissues and 44 spontaneous TBM were identified (tripling the sample size). Previous findings were confirmed and generalized and several new observations were made. TBM differ from single and other double mutations in that TBM frequency varies dramatically with tissue type. In certain tissues, most notably male germ cells, no TBM are observed despite screening as many as 26 million plaque forming units. TBM are most frequent in kidney and liver (3.45 and 2x10(-6), respectively), accounting for 7.6 and 4.8% of all mutational events in kidney and liver, respectively. There is a trend for elevated TBM frequency in thymic lymphomas in p53-deficient mice. TBM are more frequent in old age in both liver and kidney. TBM differ from single mutations and other double mutations because they display a marked difference in pattern and dramatic tissue specificity for target sequence. Five of the 78 possible TBM outcomes comprise 79% of those observed, and mutations at GG/CC predominate. TBM in mice were compared with TBM found in human mutation databases. TBM are also rare in the human germline (one in 5133 germline mutations reported in five human mutation databases). In general, the types of somatic TBM are similar in mice and humans except for an excess of TG/CA to CA/TG TBM in humans (TBM related to ultraviolet light-induced skin cancer were excluded). TBM may be the result of unknown mechanisms that may have some similarities in mice and humans.

Exposure to chronic, low-dose rate gamma-radiation at Chornobyl does not induce point mutations in Big Blue mice

Empirical genetic effects resulting from low-dose rate irradiation and chronic, cumulative exposure are poorly characterized. Expected effects are based on epidemiological studies and downward, linear extrapolations from nonthreshold models derived from acute, high-dose exposures. These extrapolations and their associated risk coefficients have no experimental support, and because of their inherent uncertainty they are the subject of considerable debate. The expectation of deleterious genetic effects resulting from low-dose rate irradiation and chronic exposure is in need of empirical assessment because this type of exposure is typical of those encountered in occupational, residential, and environmental settings. Recent acute low-dose (<10 cGy) studies using cytogenetic and point mutation endpoints indicate that observed effects range from those lower than spontaneous to an increase in the frequency of point mutations. Using the Big Blue assay, we examined the ability of chronic, continuous gamma-irradiation (2.3 x 10(-3) cGy/min) in the Chornobyl environment to induce point mutations. This system has demonstrated a significant point mutation sensitivity (4.5-fold increase) to acute, high-dose (1-3 Gy) gamma-radiation. Mutant frequencies and the mutation spectra were examined in exposed and reference samples of Big Blue mice following 90 days exposure (cumulative absorbed dose = 3 Gy) to the Chornobyl environment. No significant increase in the mutant frequency or bias in the mutational spectrum was observed in exposed individuals. This finding suggests that low-dose rate gamma-irradiation at Chornobyl does not induce point mutations and that cumulative, chronically absorbed doses do not induce the same genetic effects as acute doses of the same magnitude.

The rate and character of spontaneous mutation in an RNA virus

Estimates of spontaneous mutation rates for RNA viruses are few and uncertain, most notably due to their dependence on tiny mutation reporter sequences that may not well represent the whole genome. We report here an estimate of the spontaneous mutation rate of tobacco mosaic virus using an 804-base cognate mutational target, the viral MP gene that encodes the movement protein (MP). Selection against newly arising mutants was countered by providing MP function from a transgene. The estimated genomic mutation rate was on the lower side of the range previously estimated for lytic animal riboviruses. We also present the first unbiased riboviral mutational spectrum. The proportion of base substitutions is the same as that in a retrovirus but is lower than that in most DNA-based organisms. Although the MP mutant frequency was 0.02-0.05, 35% of the sequenced mutants contained two or more mutations. Therefore, the mutation process in populations of TMV and perhaps of riboviruses generally differs profoundly from that in populations of DNA-based microbes and may be strongly influenced by a subpopulation of mutator polymerases.

Multiple mutations are common at mouse Aprt in genotoxin-exposed mismatch repair deficient cells

Mismatch repair deficiency is known to contribute to elevated rates of mutations, particularly at mono- and dinucleotide repeat sequences. However, such repeats are often missing from the coding regions of endogenous genes. To determine the types of mutations that can occur within an endogenous gene lacking highly susceptible repeat sequences, we examined mutagenic events at the 2.3 kb mouse Aprt gene in kidney cell lines derived from mice deficient for the PMS2 and MLH1 mismatch repair proteins. The Aprt mutation rate was increased 33-fold and 3.6-20-fold for Mlh1 and Pms2 null cell lines, respectively, when compared with a wild-type kidney cell line. For the Pms2 null cells this increase resulted from both intragenic events, which were predominantly base-pairs substitutions, and loss of heterozygosity events. Almost all mutations in the Mlh1 null cells were due to base-pair substitutions. A:T-->G:C transitions (54% of small events) were predominant in the Pms2 null cells whereas G:C-->A:T transitions (36%) were the most common base-pair change in the Mlh1 null cells. Interestingly, 4-9% of the spontaneous mutant alleles in the mismatch repair deficient cells exhibited two well-separated base-pair substitution events. The percentage of mutant alleles with two and occasionally three base-pair substitutions increased when the Pms2 and Mlh1 null cells were treated with ultraviolet radiation (15-21%) and when the Mlh1 null cells were treated with hydrogen peroxide (35%). In most cases the distance separating the multiple base-pair substitutions on a given allele was in excess of 100 base-pairs, suggesting that the two mutational events were not linked directly to a single DNA lesion. The significance of these results is discussed with regards to the roles for the PMS2 and MLH1 proteins in preventing spontaneous and genotoxin-related mutations.

The unexpected landscape of in vivo somatic mutation in a human epithelial cell lineage

Few data exist on somatic mutation in the epithelial cell lineages that play a central role in human biology and disease. To delineate the "landscape" of somatic mutation in a human epithelial cell lineage, we determined the frequency and molecular nature of somatic mutations occurring in vivo in the X-linked HPRT gene of kidney tubular epithelial cells. Kidney epithelial mutants were frequent (range 0.5 to 4.2 x 10(-4)) and contained a high proportion of unreported HPRT base substitutions, -1-bp deletions and multiple mutations. This spectrum of somatic mutation differed from HPRT mutations identified in human peripheral blood T lymphocytes and from germ-line HPRT mutations identified in Lesch-Nyhan syndrome or hyperuricemia patients. Our results indicate that DNA damage and mutagenesis may have unusual or mechanistically interesting features in kidney tubular epithelium, and that somatic mutation may play a more important role in human kidney disease than has been previously appreciated.

Characterization of mutant spectra generated by a forward mutational assay for gene A of Phi X174 from ENU-treated transgenic mouse embryonic cell line PX-2

The sensitivity of in vivo transgenic mutation assays benefits from the sequencing of mutations, although the large number of possible mutations hinders high throughput sequencing. A forward mutational assay exists for Phi X174 that requires an altered, functional Phi X174 protein and therefore should have fewer targets (sense, base-pair substitutions) than forward assays that inactivate a protein. We investigated this assay to determine the number of targets and their suitability for detecting a known mutagen, N-ethyl-N-nitrosourea (ENU). We identified 25 target sites and 33 different mutations in Phi X174 gene A after sequencing over 350 spontaneous and ENU-induced mutants, mostly from mouse embryonic cell line PX-2 isolated from mice transgenic for Phi X174 am3, cs70 (line 54). All six types of base-pair substitution were represented among both the spontaneous and ENU-treated mutant spectra. The mutant spectra from cells treated with 200 and 400 microg/ml ENU were both highly different from the spontaneous spectrum (P < 0.000001) but not from each other. The dose trend was significant (P < 0.0001) for a linear regression of mutant frequencies (R(2) = 0.79), with a ninefold increase in mutant frequency at the 400 microg/ml dose. The spontaneous mutant frequency was 1.9 x 10(-5) and the spontaneous spectrum occurred at 11 target base pairs with 15 different mutations. Thirteen mutations at 12 targets were identified only from ENU-treated cells. Seven mutations had highly significant increases with ENU treatment (P < 0.0001) and 15 showed significant increases. The results suggest that the Phi X174 forward assay might be developed into a sensitive, inexpensive in vivo mutagenicity assay.

Spontaneous microdeletions and microinsertions in a transgenic mouse mutation detection system: analysis of age, tissue, and sequence specificity

A total of 3497 independent spontaneous mutations were examined using the Big Blue transgenic mouse mutation detection system. Base substitutions predominate, although 16% of somatic and germline mutations are microdeletions, microinsertions, or deletions combined with insertions. The pattern of microdeletions and microinsertions is similar in both the lacI transgene and the human p53 gene. Single-base deletions (D1) and insertions (I1) are evenly distributed in the lacI transgene, whereas microdeletions from 2 to 50 bp are clustered at two regions (bp 129-228 and 529-628). The pattern of microdeletions and microinsertions is similar between young (< or =3 months) and old (25 months) mice. Brain tissue has a paucity of deletions combined with insertions when compared with that of thymus and nine other tissues (P = 0.01). A 16-bp deletion at lacI base position 272 is a tissue-specific hotspot preferentially occurring in brain. Approximately 68 and 93% of D1 and I1, respectively, occur at mononucleotide repeats. The frequencies of D1 and I1 in mononucleotide repeats increase in an exponential manner with the length of the repeat. The lacI transgene shows similarity to the human p53 gene in the pattern of microdeletions and microinsertions and the size distribution of microdeletions.