Evidence for mutation showers

Analyzing somatic mutations by single-cell whole-genome sequencing

Somatic mutations are the cause of cancer and have been implicated in other, noncancerous diseases and aging. While clonally expanded mutations can be studied by deep sequencing of bulk DNA, very few somatic mutations expand clonally, and most are unique to each cell. We describe a detailed protocol for single-cell whole-genome sequencing to discover and analyze somatic mutations in tissues and organs. The protocol comprises single-cell multiple displacement amplification (SCMDA), which ensures efficiency and high fidelity in amplification, and the SCcaller software tool to call single-nucleotide variations and small insertions and deletions from the sequencing data by filtering out amplification artifacts. With SCMDA and SCcaller at its core, this protocol describes a complete procedure for the comprehensive analysis of somatic mutations in a single cell, covering (1) single-cell or nucleus isolation, (2) single-cell or nucleus whole-genome amplification, (3) library preparation and sequencing, and (4) computational analyses, including alignment, variant calling, and mutation burden estimation. Methods are also provided for mutation annotation, hotspot discovery and signature analysis. The protocol takes 12-15 h from single-cell isolation to library preparation and 3-7 d of data processing. Compared with other single-cell amplification methods or single-molecular sequencing, it provides high genomic coverage, high accuracy in single-nucleotide variation and small insertions and deletion calling from the same single-cell genome, and fewer processing steps. SCMDA and SCcaller require basic experience in molecular biology and bioinformatics. The protocol can be utilized for studying mutagenesis and genome mosaicism in normal and diseased human and animal tissues under various conditions.

Identification of a signature of evolutionarily conserved stress-induced mutagenesis in cancer

The clustering of mutations observed in cancer cells is reminiscent of the stress-induced mutagenesis (SIM) response in bacteria. Bacteria deploy SIM when faced with DNA double-strand breaks in the presence of conditions that elicit an SOS response. SIM employs DinB, the evolutionary precursor to human trans-lesion synthesis (TLS) error-prone polymerases, and results in mutations concentrated around DNA double-strand breaks with an abundance that decays with distance. We performed a quantitative study on single nucleotide variant calls for whole-genome sequencing data from 1950 tumors, non-inherited mutations from 129 normal samples, and acquired mutations in 3 cell line models of stress-induced adaptive mutation. We introduce statistical methods to identify mutational clusters, quantify their shapes and tease out the potential mechanism that produced them. Our results show that mutations in both normal and cancer samples are indeed clustered and have shapes indicative of SIM. Clusters in normal samples occur more often in the same genomic location across samples than in cancer suggesting loss of regulation over the mutational process during carcinogenesis. Additionally, the signatures of TLS contribute the most to mutational cluster formation in both patient samples as well as experimental models of SIM. Furthermore, a measure of cluster shape heterogeneity was associated with cancer patient survival with a hazard ratio of 5.744 (Cox Proportional Hazard Regression, 95% CI: 1.824-18.09). Our results support the conclusion that the ancient and evolutionary-conserved adaptive mutation response found in bacteria is a source of genomic instability in cancer. Biological adaptation through SIM might explain the ability of tumors to evolve in the face of strong selective pressures such as treatment and suggests that the conventional 'hit it hard' approaches to therapy could prove themselves counterproductive.

Postzygotic Somatic Mutations in the Human Brain Expand the Threshold-Liability Model of Schizophrenia

The search for what causes schizophrenia has been onerous. This research has included extensive assessment of a variety of genetic and environmental factors using ever emerging high-resolution technologies and traditional understanding of the biology of the brain. These efforts have identified a large number of schizophrenia-associated genes, some of which are altered by mutational and epi-mutational mechanisms in a threshold liability model of schizophrenia development. The results, however, have limited predictability and the actual cause of the disease remains unknown. This current state asks for conceptualizing the problem differently in light of novel insights into the nature of mutations, the biology of the brain and the fine precision and resolution of emerging technologies. There is mounting evidence that mutations acquired during postzygotic development are more common than germline mutations. Also, the postzygotic somatic mutations including epimutations (PZMs), which often lead to somatic mosaicism, are relatively common in the mammalian brain in comparison to most other tissues and PZMs are more common in patients with neurodevelopmental mental disorders, including schizophrenia. Further, previously inaccessible, detection of PZMs is becoming feasible with the advent of novel technologies that include single-cell genomics and epigenomics and the use of exquisite experimental designs including use of monozygotic twins discordant for the disease. These developments allow us to propose a working hypothesis and expand the threshold liability model of schizophrenia that already encompasses familial genetic, epigenetic and environmental factors to include somatic de novo PZMs. Further, we offer a test for this expanded model using currently available genome sequences and methylome data on monozygotic twins discordant for schizophrenia (MZD) and their parents. The results of this analysis argue that PZMs play a significant role in the development of schizophrenia and explain extensive heterogeneity seen across patients. It also offers the potential to convincingly link PZMs to both nervous system health and disease, an area that has remained challenging to study and relatively under explored.

Single-cell whole-genome sequencing reveals the functional landscape of somatic mutations in B lymphocytes across the human lifespan

Accumulation of mutations in somatic cells has been implicated as a cause of aging since the 1950s. However, attempts to establish a causal relationship between somatic mutations and aging have been constrained by the lack of methods to directly identify mutational events in primary human tissues. Here we provide genome-wide mutation frequencies and spectra of human B lymphocytes from healthy individuals across the entire human lifespan using a highly accurate single-cell whole-genome sequencing method. We found that the number of somatic mutations increases from <500 per cell in newborns to >3,000 per cell in centenarians. We discovered mutational hotspot regions, some of which, as expected, were located at Ig genes associated with somatic hypermutation (SHM). B cell-specific mutation signatures associated with development, aging, or SHM were found. The SHM signature strongly correlated with the signature found in human B cell tumors, indicating that potential cancer-causing events are already present even in B cells of healthy individuals. We also identified multiple mutations in sequence features relevant to cellular function (i.e., transcribed genes and gene regulatory regions). Such mutations increased significantly during aging, but only at approximately one-half the rate of the genome average, indicating selection against mutations that impact B cell function. This full characterization of the landscape of somatic mutations in human B lymphocytes indicates that spontaneous somatic mutations accumulating with age can be deleterious and may contribute to both the increased risk for leukemia and the functional decline of B lymphocytes in the elderly.

Homologous Recombination and the Formation of Complex Genomic Rearrangements

The maintenance of genome integrity involves multiple independent DNA damage avoidance and repair mechanisms. However, the origin and pathways of the focal chromosomal reshuffling phenomena collectively referred to as chromothripsis remain mechanistically obscure. We discuss here the role, mechanisms, and regulation of homologous recombination (HR) in the formation of simple and complex chromosomal rearrangements. We emphasize features of the recently characterized multi-invasion (MI)-induced rearrangement (MIR) pathway which uniquely amplifies the initial DNA damage. HR intermediates and cellular contexts that endanger genomic stability are discussed as well as the emerging roles of various classes of nucleases in the formation of genome rearrangements. Long-read sequencing and improved mapping of repeats should enable better appreciation of the significance of recombination in generating genomic rearrangements.

Spatial statistical tools for genome-wide mutation cluster detection under a microarray probe sampling system

Mutation cluster analysis is critical for understanding certain mutational mechanisms relevant to genetic disease, diversity, and evolution. Yet, whole genome sequencing for detection of mutation clusters is prohibitive with high cost for most organisms and population surveys. Single nucleotide polymorphism (SNP) genotyping arrays, like the Mouse Diversity Genotyping Array, offer an alternative low-cost, screening for mutations at hundreds of thousands of loci across the genome using experimental designs that permit capture of de novo mutations in any tissue. Formal statistical tools for genome-wide detection of mutation clusters under a microarray probe sampling system are yet to be established. A challenge in the development of statistical methods is that microarray detection of mutation clusters is constrained to select SNP loci captured by probes on the array. This paper develops a Monte Carlo framework for cluster testing and assesses test statistics for capturing potential deviations from spatial randomness which are motivated by, and incorporate, the array design. While null distributions of the test statistics are established under spatial randomness via the homogeneous Poisson process, power performance of the test statistics is evaluated under postulated types of Neyman-Scott clustering processes through Monte Carlo simulation. A new statistic is developed and recommended as a screening tool for mutation cluster detection. The statistic is demonstrated to be excellent in terms of its robustness and power performance, and useful for cluster analysis in settings of missing data. The test statistic can also be generalized to any one dimensional system where every site is observed, such as DNA sequencing data. The paper illustrates how the informal graphical tools for detecting clusters may be misleading. The statistic is used for finding clusters of putative SNP differences in a mixture of different mouse genetic backgrounds and clusters of de novo SNP differences arising between tissues with development and carcinogenesis.

Two de novo novel mutations in one SHANK3 allele in a patient with autism and moderate intellectual disability

SHANK3 encodes for a scaffolding protein that links neurotransmitter receptors to the cytoskeleton and is enriched in postsynaptic densities of excitatory synapses. Deletions or mutations in one copy of the SHANK3 gene cause Phelan-McDermid syndrome, also called 22q13.3 deletion syndrome, a neurodevelopmental disorder with common features including global developmental delay, absent to severely impaired language, autistic behavior, and minor dysmorphic features. By whole exome sequencing, we identified two de novo novel variants including one frameshift pathogenic variant and one missense variant of unknown significance in a 14-year-old boy with delayed motor milestones, delayed language acquisition, autism, intellectual disability, ataxia, progressively worsening spasticity of the lower extremities, dysmorphic features, short stature, microcephaly, failure to thrive, chronic constipation, intrauterine growth restriction, and bilateral inguinal hernias. Both changes are within the CpG island in exon 21, separated by a 375 bp sequence. Next generation sequencing of PCR products revealed that the two variants are most frequently associated with each other. Sanger sequencing of the cloned PCR products further confirmed that both changes were on a single allele. The clinical presentation in this individual is consistent with other patients with a truncating mutation in exon 21, suggesting that the missense change contributes none or minimally to the phenotypes. This is the first report of two de novo mutations in one SHANK3 allele.

Contributions of Zea mays subspecies mexicana haplotypes to modern maize

Maize was domesticated from lowland teosinte (Zea mays ssp. parviglumis), but the contribution of highland teosinte (Zea mays ssp. mexicana, hereafter mexicana) to modern maize is not clear. Here, two genomes for Mo17 (a modern maize inbred) and mexicana are assembled using a meta-assembly strategy after sequencing of 10 lines derived from a maize-teosinte cross. Comparative analyses reveal a high level of diversity between Mo17, B73, and mexicana, including three Mb-size structural rearrangements. The maize spontaneous mutation rate is estimated to be 2.17 × 10^-8 ~3.87 × 10^-8 per site per generation with a nonrandom distribution across the genome. A higher deleterious mutation rate is observed in the pericentromeric regions, and might be caused by differences in recombination frequency. Over 10% of the maize genome shows evidence of introgression from the mexicana genome, suggesting that mexicana contributed to maize adaptation and improvement. Our data offer a rich resource for constructing the pan-genome of Zea mays and genetic improvement of modern maize varieties.

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.

AID/APOBEC deaminases and cancer

Mutations are the basis for evolution and the development of genetic diseases. Especially in cancer, somatic mutations in oncogenes and tumor suppressor genes alongside the occurrence of passenger mutations have been observed by recent deep-sequencing approaches. While mutations have long been considered random events induced by DNA-replication errors or by DNA damaging agents, genome sequencing led to the discovery of non-random mutation signatures in many human cancer. Common non-random mutations comprise DNA strand-biased mutation showers and mutations restricted to certain DNA motifs, which recently have become attributed to the activity of the AID/APOBEC family of DNA deaminases. Hence, APOBEC enzymes, which have evolved as key players in natural and adaptive immunity, have been proposed to contribute to cancer development and clonal evolution of cancer by inducing collateral genomic damage due to their DNA deaminating activity. This review focuses on how mutagenic events through AID/APOBEC deaminases may contribute to cancer development.

APOBEC3 genes: retroviral restriction factors to cancer drivers

The APOBEC3 cytosine deaminases play key roles in innate immunity through their ability to mutagenize viral DNA and restrict viral replication. Recent advances in cancer genomics, together with biochemical characterization of the APOBEC3 enzymes, have now implicated at least two family members in somatic mutagenesis during tumor development. We review the evidence linking these enzymes to carcinogenesis and highlight key questions, including the potential mechanisms that misdirect APOBEC3 activity to the host genome, the links to viral infection, and the association between a common APOBEC3 polymorphism and cancer risk.

Hypermutation in human cancer genomes: footprints and mechanisms

A role for somatic mutations in carcinogenesis is well accepted, but the degree to which mutation rates influence cancer initiation and development is under continuous debate. Recently accumulated genomic data have revealed that thousands of tumour samples are riddled by hypermutation, broadening support for the idea that many cancers acquire a mutator phenotype. This major expansion of cancer mutation data sets has provided unprecedented statistical power for the analysis of mutation spectra, which has confirmed several classical sources of mutation in cancer, highlighted new prominent mutation sources (such as apolipoprotein B mRNA editing enzyme catalytic polypeptide-like (APOBEC) enzymes) and empowered the search for cancer drivers. The confluence of cancer mutation genomics and mechanistic insight provides great promise for understanding the basic development of cancer through mutations.

Do mutator mutations fuel tumorigenesis?

The mutator phenotype hypothesis proposes that the mutation rate of normal cells is insufficient to account for the large number of mutations found in human cancers. Consequently, human tumors exhibit an elevated mutation rate that increases the likelihood of a tumor acquiring advantageous mutations. The hypothesis predicts that tumors are composed of cells harboring hundreds of thousands of mutations, as opposed to a small number of specific driver mutations, and that malignant cells within a tumor therefore constitute a highly heterogeneous population. As a result, drugs targeting specific mutated driver genes or even pathways of mutated driver genes will have only limited anticancer potential. In addition, because the tumor is composed of such a diverse cell population, tumor cells harboring drug-resistant mutations will exist prior to the administration of any chemotherapeutic agent. We present recent evidence in support of the mutator phenotype hypothesis, major arguments against this concept, and discuss the clinical consequences of tumor evolution fueled by an elevated mutation rate. We also consider the therapeutic possibility of altering the rate of mutation accumulation. Most significantly, we contend that there is a need to fundamentally reconsider current approaches to personalized cancer therapy. We propose that targeting cellular pathways that alter the rate of mutation accumulation in tumors will ultimately prove more effective than attempting to identify and target mutant driver genes or driver pathways.

Stress-induced cellular adaptive strategies: ancient evolutionarily conserved programs as new anticancer therapeutic targets

Despite the remarkable achievements of novel targeted anti-cancer drugs, most therapies only produce remission for a limited time, resistance to treatment, and relapse, often being the ultimate outcome. Drug resistance is due to highly efficient adaptive strategies utilized by cancer cells. Exogenous and endogenous stress stimuli are known to induce first-line responses, capable of re-establishing cellular homeostasis and determining cell fate decisions. Cancer cells may also mount second-line adaptive strategies, such as the mutator response. Hypermutable subpopulations of cells may expand under severe selective stress, thereby accelerating the emergence of adapted clones. As with first-line protective responses, these strategies appear highly conserved, and are found in yeasts and bacteria. We hypothesize that evolutionarily conserved programs rheostatically regulate mutability in fluctuating environments, and contribute to drug resistance in cancer cells. Elucidating the conserved genetic and molecular mechanisms may present novel opportunities to increase the effectiveness of cancer therapies.

Clustered and genome-wide transient mutagenesis in human cancers: Hypermutation without permanent mutators or loss of fitness

The gain of a selective advantage in cancer as well as the establishment of complex traits during evolution require multiple genetic alterations, but how these mutations accumulate over time is currently unclear. There is increasing evidence that a mutator phenotype perpetuates the development of many human cancers. While in some cases the increased mutation rate is the result of a genetic disruption of DNA repair and replication or environmental exposures, other evidence suggests that endogenous DNA damage induced by AID/APOBEC cytidine deaminases can result in transient localized hypermutation generating simultaneous, closely spaced (i.e. "clustered") multiple mutations. Here, we discuss mechanisms that lead to mutation cluster formation, the biological consequences of their formation in cancer and evidence suggesting that APOBEC mutagenesis can also occur genome-wide. This raises the possibility that dysregulation of these enzymes may enable rapid malignant transformation by increasing mutation rates without the loss of fitness associated with permanent mutators.

Endogenous Mutagenesis in Recombinant Sulfolobus Plasmids

Low rates of replication errors in chromosomal genes of Sulfolobus spp. demonstrate that these extreme thermoacidophiles can maintain genome integrity in environments with high temperature and low pH. In contrast to this genetic stability, we observed unusually frequent mutation of the β- d -glycosidase gene ( lacS ) of a shuttle plasmid (pJ lacS ) propagated in Sulfolobus acidocaldarius . The resulting Lac − mutants also grew faster than the Lac + parent, thereby amplifying the impact of the frequent lacS mutations on the population. We developed a mutant accumulation assay and corrections for the effects of copy number and differential growth for this system; the resulting measurements and calculations yielded a corrected rate of 5.1 × 10 −4 mutational events at the lacS gene per plasmid replication. Analysis of independent lacS mutants revealed three types of mutations: (i) G·C-to-A·T transitions, (ii) slipped-strand events, and (iii) deletions. These mutations were frequent in plasmid-borne lacS expressed at a high level but not in single-copy lacS in the chromosome or at lower levels of expression in a plasmid. Substitution mutations arose at only two of 12 potential priming sites of the DNA primase of the pRN1 replicon, but nearly all these mutations created nonsense (chain termination) codons. The spontaneous mutation rate of plasmid-borne lacS was 175-fold higher under high-expression than under low-expression conditions. The results suggest that important DNA repair or replication fidelity functions are impaired or overwhelmed in pJ lacS , with results analogous to those of the “transcription-associated mutagenesis” seen in bacteria and eukaryotes.

The rate of spontaneous mutations in human myeloid cells

The mutation rate (μ) is likely to be a key parameter in leukemogenesis, but historically, it has been difficult to measure in humans. The PIG-A gene has some advantages for the detection of spontaneous mutations because it is X-linked, and therefore only one mutation is required to disrupt its function. Furthermore, the PIG-A-null phenotype is readily detected by flow cytometry. Using PIG-A, we have now provided the first in vitro measurement of μ in myeloid cells, using cultures of CD34+ cells that are transduced with either the AML-ETO or the MLL-AF9 fusion genes and expanded with cytokines. For the AML-ETO cultures, the median μ value was ∼9.4×10(-7) (range ∼3.6-23×10(-7)) per cell division. In contrast, few spontaneous mutations were observed in the MLL-AF9 cultures. Knockdown of p53 or introduction of mutant NRAS or FLT3 alleles did not have much of an effect on μ. Based on these data, we provide a model to predict whether hypermutability must occur in the process of leukemogenesis.

Whole-genome sequencing in autism identifies hot spots for de novo germline mutation

De novo mutation plays an important role in autism spectrum disorders (ASDs). Notably, pathogenic copy number variants (CNVs) are characterized by high mutation rates. We hypothesize that hypermutability is a property of ASD genes and may also include nucleotide-substitution hot spots. We investigated global patterns of germline mutation by whole-genome sequencing of monozygotic twins concordant for ASD and their parents. Mutation rates varied widely throughout the genome (by 100-fold) and could be explained by intrinsic characteristics of DNA sequence and chromatin structure. Dense clusters of mutations within individual genomes were attributable to compound mutation or gene conversion. Hypermutability was a characteristic of genes involved in ASD and other diseases. In addition, genes impacted by mutations in this study were associated with ASD in independent exome-sequencing data sets. Our findings suggest that regional hypermutation is a significant factor shaping patterns of genetic variation and disease risk in humans.

Prevalence of multinucleotide replacements in evolution of primates and Drosophila

Evolution of sequences mostly involves independent changes at different sites. However, substitutions at neighboring sites may co-occur as multinucleotide replacement events (MNRs). Here, we compare noncoding sequences of several species of primates, and of three species of Drosophila fruit flies, in a phylogenetic analysis of the replacements that occurred between species at nearby nucleotide sites. Both in primates and in Drosophila, the frequency of single-nucleotide replacements is substantially elevated within 10 nucleotides from other replacements that occurred on the same lineage but not on another lineage. The data imply that dinucleotide replacements (DNRs) affecting sites at distances of up to 10 nucleotides from each other are responsible for 2.3% of single-nucleotide replacements in primate genomes and for 5.6% in Drosophila genomes. Among these DNRs, 26% and 69%, respectively, are in fact parts of replacements of three or more trinucleotide replacements (TNRs). The plurality of MNRs affect nearby nucleotides, so that at least six times as many DNRs affect two adjacent nucleotide sites than sites 10 nucleotides apart. Still, approximately 60% of DNRs, and approximately 90% of TNRs, span distances more than two (or three) nucleotides. MNRs make a major contribution to the observed clustering of substitutions: In the human–chimpanzee comparison, DNRs are responsible for 50% of cases when two nearby replacements are observed on the human lineage, and TNRs are responsible for 83% of cases when three replacements at three immediately adjacent sites are observed on the human lineage. The prevalence of MNRs matches that is observed in data on de novo mutations and is also observed in the regions with the lowest sequence conservation, suggesting that MNRs mainly have mutational origin; however, epistatic selection and/or gene conversion may also play a role.

Two mechanisms produce mutation hotspots at DNA breaks in Escherichia coli

Mutation hotspots and showers occur across phylogeny and profoundly influence genome evolution, yet the mechanisms that produce hotspots remain obscure. We report that DNA double-strand breaks (DSBs) provoke mutation hotspots via stress-induced mutation in Escherichia coli. With tet reporters placed 2 kb to 2 Mb (half the genome) away from an I-SceI site, RpoS/DinB-dependent mutations occur maximally within the first 2 kb and decrease logarithmically to ∼60 kb. A weak mutation tail extends to 1 Mb. Hotspotting occurs independently of I-site/tet-reporter-pair position in the genome, upstream and downstream in the replication path. RecD, which allows RecBCD DSB-exonuclease activity, is required for strong local but not long-distance hotspotting, indicating that double-strand resection and gap-filling synthesis underlie local hotspotting, and newly illuminating DSB resection in vivo. Hotspotting near DSBs opens the possibility that specific genomic regions could be targeted for mutagenesis, and could also promote concerted evolution (coincident mutations) within genes/gene clusters, an important issue in the evolution of protein functions.

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.

Mutational processes molding the genomes of 21 breast cancers

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed “kataegis,” was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

PaperClip

Clustered mutations in yeast and in human cancers can arise from damaged long single-strand DNA regions

Mutations are typically perceived as random, independent events. We describe here nonrandom clustered mutations in yeast and in human cancers. Genome sequencing of yeast grown under chronic alkylation damage identified mutation clusters that extend up to 200 kb. A predominance of "strand-coordinated" changes of either cytosines or guanines in the same strand, mutation patterns, and genetic controls indicated that simultaneous mutations were generated by base alkylation in abnormally long single-strand DNA (ssDNA) formed at double-strand breaks (DSBs) and replication forks. Significantly, we found mutation clusters with analogous features in sequenced human cancers. Strand-coordinated clusters of mutated cytosines or guanines often resided near chromosome rearrangement breakpoints and were highly enriched with a motif targeted by APOBEC family cytosine-deaminases, which strongly prefer ssDNA. These data indicate that hypermutation via multiple simultaneous changes in randomly formed ssDNA is a general phenomenon that may be an important mechanism producing rapid genetic variation.

Two closely spaced nonsense mutations in the DMD gene in a Malaysian family

In Duchenne muscular dystrophy (DMD), identification of one nonsense mutation in the DMD gene has been considered an endpoint of genetic diagnosis. Here, we identified two closely spaced nonsense mutations in the DMD gene. In a Malaysian DMD patient two nonsense mutations (p.234S>X and p.249Q>X, respectively) were identified within exon 8. The proband's mother carried both mutations on one allele. Multiple mutations may explain the occasional discrepancies between genotype and phenotype in dystrophinopathy.

Peto's Paradox: evolution's prescription for cancer prevention

The evolution of multicellularity required the suppression of cancer. If every cell has some chance of becoming cancerous, large, long-lived organisms should have an increased risk of developing cancer compared with small, short-lived organisms. The lack of correlation between body size and cancer risk is known as Peto's paradox. Animals with 1000 times more cells than humans do not exhibit an increased cancer risk, suggesting that natural mechanisms can suppress cancer 1000 times more effectively than is done in human cells. Because cancer has proven difficult to cure, attention has turned to cancer prevention. In this review, similar to pharmaceutical companies mining natural products, we seek to understand how evolution has suppressed cancer to develop ultimately improved cancer prevention in humans.

"Patchy-tachy" leads to false positives for recombination

Indirect tests have detected recombination in mitochondrial DNA (mtDNA) from many animal lineages, including mammals. However, it is possible that features of the molecular evolutionary process without recombination could be incorrectly inferred by indirect tests as being due to recombination. We have identified one such example, which we call "patchy-tachy" (PT), where different partitions of sequences evolve at different rates, that leads to an excess of false positives for recombination inferred by indirect tests. To explore this phenomena, we characterized the false positive rates of six widely used indirect tests for recombination using simulations of general models for mtDNA evolution with PT but without recombination. All tests produced 30-99% false positives for recombination, although the conditions that produced the maximal level of false positives differed between the tests. To evaluate the degree to which conditions that exacerbate false positives are found in published sequence data, we turned to 20 animal mtDNA data sets in which recombination is suggested by indirect tests. Using a model where different regions of the sequences were free to evolve at different rates in different lineages, we demonstrated that PT is prevalent in many data sets in which recombination was previously inferred using indirect tests. Taken together, our results argue that PT without recombination is a viable alternative explanation for detection of widespread recombination in animal mtDNA using indirect tests.

Damage-induced localized hypermutability

Genome instability continuously presents perils of cancer, genetic disease and death of a cell or an organism. At the same time, it provides for genome plasticity that is essential for development and evolution. We address here the genome instability confined to a small fraction of DNA adjacent to free DNA ends at uncapped telomeres and double-strand breaks. We found that budding yeast cells can tolerate nearly 20 kilobase regions of subtelomeric single-strand DNA that contain multiple UV-damaged nucleotides. During restoration to the double-strand state, multiple mutations are generated by error-prone translesion synthesis. Genome-wide sequencing demonstrated that multiple regions of damage-induced localized hypermutability can be tolerated, which leads to the simultaneous appearance of multiple mutation clusters in the genomes of UV- irradiated cells. High multiplicity and density of mutations suggest that this novel form of genome instability may play significant roles in generating new alleles for evolutionary selection as well as in the incidence of cancer and genetic disease.

Correlated evolution of nearby residues in Drosophilid proteins

Here we investigate the correlations between coding sequence substitutions as a function of their separation along the protein sequence. We consider both substitutions between the reference genomes of several Drosophilids as well as polymorphisms in a population sample of Zimbabwean Drosophila melanogaster. We find that amino acid substitutions are “clustered” along the protein sequence, that is, the frequency of additional substitutions is strongly enhanced within ≈10 residues of a first such substitution. No such clustering is observed for synonymous substitutions, supporting a “correlation length” associated with selection on proteins as the causative mechanism. Clustering is stronger between substitutions that arose in the same lineage than it is between substitutions that arose in different lineages. We consider several possible origins of clustering, concluding that epistasis (interactions between amino acids within a protein that affect function) and positional heterogeneity in the strength of purifying selection are primarily responsible. The role of epistasis is directly supported by the tendency of nearby substitutions that arose on the same lineage to preserve the total charge of the residues within the correlation length and by the preferential cosegregation of neighboring derived alleles in our population sample. We interpret the observed length scale of clustering as a statistical reflection of the functional locality (or modularity) of proteins: amino acids that are near each other on the protein backbone are more likely to contribute to, and collaborate toward, a common subfunction.

Genes are templates for proteins, yet evolutionary studies of genes and proteins often bear little resemblance. Analyses of gene evolution typically treat each codon independently, quantifying gene evolution by summing over the constituent codons. In contrast, studies of protein evolution generally incorporate protein structure and interactions between amino acids explicitly. We investigate correlations in the evolution of codons as a function of their distance from each other along the protein coding sequence. This approach is motivated by the expectation that codons near each other in sequence often encode amino acids belonging to the same functional unit. Consequently, these amino acids are more likely to interact and/or experience similar selective regimes, introducing correlation between the evolution of the underlying codons. We find codon evolution in Drosophilids to be correlated over a characteristic length scale of ≈10 codons. Specifically, the presence of a non-synonymous substitution substantially increases the probability of further such substitutions nearby, particularly within that lineage. Further analysis suggests both functional interactions between amino acids and correlation in the strength of selection contribute to this effect. These findings are relevant for understanding the relative importance of different modes of selection, and particularly the role of epistasis, in gene and protein evolution.

Break-induced replication is highly inaccurate

DNA replication initiated by one-ended homologous recombination at a double-strand break is highly inaccurate, as it greatly stimulates frameshift mutations over the entire path of the replication fork.

DNA must be synthesized for purposes of genome duplication and DNA repair. While the former is a highly accurate process, short-patch synthesis associated with repair of DNA damage is often error-prone. Break-induced replication (BIR) is a unique cellular process that mimics normal DNA replication in its processivity, rate, and capacity to duplicate hundreds of kilobases, but is initiated at double-strand breaks (DSBs) rather than at replication origins. Here we employed a series of frameshift reporters to measure mutagenesis associated with BIR in Saccharomyces cerevisiae. We demonstrate that BIR DNA synthesis is intrinsically inaccurate over the entire path of the replication fork, as the rate of frameshift mutagenesis during BIR is up to 2,800-fold higher than during normal replication. Importantly, this high rate of mutagenesis was observed not only close to the DSB where BIR is less stable, but also far from the DSB where the BIR replication fork is fast and stabilized. We established that polymerase proofreading and mismatch repair correct BIR errors. Also, dNTP levels were elevated during BIR, and this contributed to BIR-related mutagenesis. We propose that a high level of DNA polymerase errors that is not fully compensated by error-correction mechanisms is largely responsible for mutagenesis during BIR, with Pol δ generating many of the mutagenic errors. We further postulate that activation of BIR in eukaryotic cells may significantly contribute to accumulation of mutations that fuel cancer and evolution.

Accurate transmission of genetic information requires the precise replication of parental DNA. Mutations (which can be beneficial or deleterious) arise from errors that remain uncorrected. DNA replication occurs during S-phase of the cell cycle and is extremely accurate due to highly selective DNA polymerases coupled with effective error-correction mechanisms. In contrast, DNA synthesis associated with short-patch DNA repair is often error-prone. Break-induced replication (BIR) presents an interesting case of large-scale DNA duplication that occurs in the context of DNA repair. In this study we employed a yeast-based system to investigate the level of mutagenesis associated with BIR compared to mutagenesis during normal DNA replication. We report that frameshifts, which are the most deleterious kind of point mutation, are much more frequent during BIR than during normal DNA replication. Surprisingly, we observed that the majority of mutations associated with BIR were created by polymerases responsible for normal DNA replication, which are assumed to be highly precise. Overall, we propose that BIR is a novel source of mutagenesis that may contribute to disease genesis and evolution.

Effects of calorie restriction on the age-dependent accumulation of mutations in the small intestine of lacZ-transgenic mice

To understand the effect of calorie restriction on genome maintenance systems, the age-dependent accumulation of mutations in animals maintained on high and low calorie diets was examined using lacZ-transgenic mice. Mice were fed a diet of 95 kcal/w or 65 kcal/w from 2 to 17 months of age. The mutation frequencies in the lacZ gene in epithelial tissues from the small intestine were examined at 12 and 17 months. Mutation frequencies were found to be lower in mice fed with a low calorie diet than in mice fed with a high calorie diet at the two age points. The molecular nature of the mutations was examined with DNA sequencing. It showed a predominance of transversions from G:C to T:A, and this is a typical type of mutation induced by reactive oxygen species. The fraction of this type of mutation among the different types of mutations detected was not affected by calorie restriction. The percentage of the other types of mutation was not influenced either. These results suggest that calorie restriction reduces the age-dependent accumulation of mutations by stimulating or inducing various types of DNA protection and repair systems rather than protecting cells against any specific type of DNA alteration.

Massive genomic rearrangement acquired in a single catastrophic event during cancer development

Cancer is driven by somatically acquired point mutations and chromosomal rearrangements, conventionally thought to accumulate gradually over time. Using next-generation sequencing, we characterize a phenomenon, which we term chromothripsis, whereby tens to hundreds of genomic rearrangements occur in a one-off cellular crisis. Rearrangements involving one or a few chromosomes crisscross back and forth across involved regions, generating frequent oscillations between two copy number states. These genomic hallmarks are highly improbable if rearrangements accumulate over time and instead imply that nearly all occur during a single cellular catastrophe. The stamp of chromothripsis can be seen in at least 2%–3% of all cancers, across many subtypes, and is present in ∼25% of bone cancers. We find that one, or indeed more than one, cancer-causing lesion can emerge out of the genomic crisis. This phenomenon has important implications for the origins of genomic remodeling and temporal emergence of cancer.

PaperClip

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.

Genesis of ancestral haplotypes: RNA modifications and reverse transcription-mediated polymorphisms

Understanding the genesis of the block haplotype structure of the genome is a major challenge. With the completion of the sequencing of the Human Genome and the initiation of the HapMap project the concept that the chromosomes of the mammalian genome are a mosaic, or patchwork, of conserved extended block haplotype sequences is now accepted by the mainstream genomics research community. Ancestral Haplotypes (AHs) can be viewed as a recombined string of smaller Polymorphic Frozen Blocks (PFBs). How have such variant extended DNA sequence tracts emerged in evolution? Here the relevant literature on the problem is reviewed from various fields of molecular and cell biology particularly molecular immunology and comparative and functional genomics. Based on our synthesis we then advance a testable molecular and cellular model. A critical part of the analysis concerns the origin of the strand biased mutation signatures in the transcribed regions of the human and higher primate genome, A-to-G versus T-to-C (ratio ∼ 1.5 fold) and C-to-T versus G-to-A (≥ 1.5 fold). A comparison and evaluation of the current state of the fields of immunoglobulin Somatic Hypermutation (SHM) and Transcription-Coupled DNA Repair focused on how mutations in newly synthesized RNA might be copied back to DNA thus accounting for some of the genome-wide strand biases (e.g., the A-to-G vs T-to-C component of the strand biased spectrum). We hypothesize that the genesis of PFBs and extended AHs occurs during mutagenic episodes in evolution (e.g., retroviral infections) and that many of the critical DNA sequence diversifying events occur first at the RNA level, e.g., recombination between RNA strings resulting in tandem and dispersed RNA duplications (retroduplications), RNA mutations via adenosine-to-inosine pre-mRNA editing events as well as error prone RNA synthesis. These are then copied back into DNA by a cellular reverse transcription process (also likely to be error-prone) that we have called "reverse transcription-mediated long DNA conversion." Finally we suggest that all these activities and others can be envisaged as being brought physically under the umbrella of special sites in the nucleus involved in transcription known as "transcription factories."

Localized hypermutation and associated gene losses in legume chloroplast genomes

Point mutations result from errors made during DNA replication or repair, so they are usually expected to be homogeneous across all regions of a genome. However, we have found a region of chloroplast DNA in plants related to sweetpea (Lathyrus) whose local point mutation rate is at least 20 times higher than elsewhere in the same molecule. There are very few precedents for such heterogeneity in any genome, and we suspect that the hypermutable region may be subject to an unusual process such as repeated DNA breakage and repair. The region is 1.5 kb long and coincides with a gene, ycf4, whose rate of evolution has increased dramatically. The product of ycf4, a photosystem I assembly protein, is more divergent within the single genus Lathyrus than between cyanobacteria and other angiosperms. Moreover, ycf4 has been lost from the chloroplast genome in Lathyrus odoratus and separately in three other groups of legumes. Each of the four consecutive genes ycf4-psaI-accD-rps16 has been lost in at least one member of the legume "inverted repeat loss" clade, despite the rarity of chloroplast gene losses in angiosperms. We established that accD has relocated to the nucleus in Trifolium species, but were unable to find nuclear copies of ycf4 or psaI in Lathyrus. Our results suggest that, as well as accelerating sequence evolution, localized hypermutation has contributed to the phenomenon of gene loss or relocation to the nucleus.

Mutator phenotypes due to DNA replication infidelity

This article considers the fidelity of DNA replication performed by eukaryotic DNA polymerases involved in replicating the nuclear genome. DNA replication fidelity can vary widely depending on the DNA polymerase, the composition of the error, the flanking sequence, the presence of DNA damage and the ability to correct errors. As a consequence, defects in processes that determine DNA replication fidelity can confer strong mutator phenotypes whose specificity can help determine the molecular nature of the defect.

A single-strand specific lesion drives MMS-induced hyper-mutability at a double-strand break in yeast

Localized hyper-mutability (LHM) can be important in evolution, immunity, and genetic diseases. We previously reported that single-strand DNA (ssDNA) can be an important source of damage-induced LHM in yeast. Here, we establish that the generation of LHM by methyl methanesulfonate (MMS) during repair of a chromosomal double-strand break (DSB) can result in over 0.2 mutations/kb, which is approximately 20,000-fold higher than the MMS-induced mutation density without a DSB. The MMS-induced mutations associated with DSB repair were primarily due to substitutions via translesion DNA synthesis at damaged cytosines, even though there are nearly 10 times more MMS-induced lesions at other bases. Based on this mutation bias, the promutagenic lesion dominating LHM is likely 3-methylcytosine, which is single-strand specific. Thus, the dramatic increase in mutagenesis at a DSB is concluded to result primarily from the generation of non-repairable lesions in ssDNA associated with DSB repair along with efficient induction of highly mutagenic ssDNA-specific lesions. These findings with MMS-induced LHM have broad biological implications for unrepaired damage generated in ssDNA and possibly ssRNA.

Mutational clusters generated by non-processive polymerases: A case study using DNA polymerase betain vitro

Available DNA mutational spectra reveal that the number of mutants with multiple mutations ("multiples") is usually greater than expected from a random distribution of mutations among mutants. These overloads imply the occurrence of non-random clusters of mutations, probably generated during episodes of low-fidelity DNA synthesis. Excess multiples have been reported not only for viruses, bacteria, and eukaryotic cells but also for the DNA polymerases of phages T4 and RB69 in vitro. In the simplest case of a purified polymerase, non-random clusters may be generated by a subfraction of phenotypic variants able to introduce more errors per cycle of DNA synthesis than the normal enzyme. According to this hypothesis, excess multiples are not expected with non-processive polymerases even if they harbor rare mutator variants. DNA polymerase beta (Pol beta) is a mammalian DNA-repair polymerase with very low processivity. Although several Pol beta mutational spectra have been described, there is conflicting evidence on whether or not excess multiples occur, with spectra based on the HSV-tk system tending to show excess multiples. Excess multiples generated by Pol beta or any of its mutants might imply that the excesses of multiples observed in numerous other systems, especially those with processive polymerases, could be artifactual. Here, the distributions of mutations generated by native and recombinant rat Pol beta and by the Pol beta(Y265C) mutator were analyzed in the M13mp2 lacZalpha system. Our results present no evidence for a significant excess of multiples over the expected numbers with any of the Pol beta enzymes tested in this system. The reported excess of Pol beta-generated multiples in the HSV-tk system may reflect a reduced efficiency of detection of base substitutions that cause weak phenotypes, which in turn may artifactually increase the frequency of multiples.

Small variable segments constitute a major type of diversity of bacterial genomes at the species level

BACKGROUND

Analysis of large scale diversity in bacterial genomes has mainly focused on elements such as pathogenicity islands, or more generally, genomic islands. These comprise numerous genes and confer important phenotypes, which are present or absent depending on strains. We report that despite this widely accepted notion, most diversity at the species level is composed of much smaller DNA segments, 20 to 500 bp in size, which we call microdiversity.

RESULTS

We performed a systematic analysis of the variable segments detected by multiple whole genome alignments at the DNA level on three species for which the greatest number of genomes have been sequenced: Escherichia coli, Staphylococcus aureus, and Streptococcus pyogenes. Among the numerous sites of variability, 62 to 73% were loci of microdiversity, many of which were located within genes. They contribute to phenotypic variations, as 3 to 6% of all genes harbor microdiversity, and 1 to 9% of total genes are located downstream from a microdiversity locus. Microdiversity loci are particularly abundant in genes encoding membrane proteins. In-depth analysis of the E. coli alignments shows that most of the diversity does not correspond to known mobile or repeated elements, and it is likely that they were generated by illegitimate recombination. An intriguing class of microdiversity includes small blocks of highly diverged sequences, whose origin is discussed.

CONCLUSIONS

This analysis uncovers the importance of this small-sized genome diversity, which we expect to be present in a wide range of bacteria, and possibly also in many eukaryotic genomes.

The rate of establishment of complex adaptations

A central problem in evolutionary theory concerns the mechanisms by which adaptations requiring multiple mutations emerge in natural populations. We develop a series of expressions that clarify the scaling of the time to establishment of complex adaptations with population size, mutation rate, magnitude of the selective disadvantage of intermediate-state alleles, and the complexity of the adaptation. In general, even in the face of deleterious intermediate steps, the time to establishment is minimized in populations with very large size. Under a broad range of conditions, the time to establishment also scales by no more than the square of the mutation rate, regardless of the number of sites contributing to the adaptive change, demonstrating that the emergence of complex adaptations is only weakly constrained by the independent acquisition of mutations at the underlying sites. Mutator alleles with deleterious side effects have only moderate effects on the rate of adaptation in large populations but can cause a quantum decrease in the time to establishment of some adaptive alleles in small populations, although probably not at a high enough rate to offset the increased deleterious mutation load. Transient hypermutability, whereby a subset of gamete-producing cells mutate at an elevated rate in a nonheritable manner, may also elevate the rate of adaptation, although the effect is modest and appears to result from a simple increase in the rate of transitions between intermediate states rather than from the saltational production of doublet mutations. Taken together, these results illustrate the plausibility of the relatively rapid emergence of specific complex adaptations by conventional population genetic mechanisms and provide insight into the relative incidences of various paths of allelic adaptation in organisms with different population genetic features.

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.

Is evolution Darwinian or/and Lamarckian?

Background

The year 2009 is the 200^th anniversary of the publication of Jean-Bapteste Lamarck's Philosophie Zoologique and the 150^th anniversary of Charles Darwin's On the Origin of Species. Lamarck believed that evolution is driven primarily by non-randomly acquired, beneficial phenotypic changes, in particular, those directly affected by the use of organs, which Lamarck believed to be inheritable. In contrast, Darwin assigned a greater importance to random, undirected change that provided material for natural selection.

The concept

The classic Lamarckian scheme appears untenable owing to the non-existence of mechanisms for direct reverse engineering of adaptive phenotypic characters acquired by an individual during its life span into the genome. However, various evolutionary phenomena that came to fore in the last few years, seem to fit a more broadly interpreted (quasi)Lamarckian paradigm. The prokaryotic CRISPR-Cas system of defense against mobile elements seems to function via a bona fide Lamarckian mechanism, namely, by integrating small segments of viral or plasmid DNA into specific loci in the host prokaryote genome and then utilizing the respective transcripts to destroy the cognate mobile element DNA (or RNA). A similar principle seems to be employed in the piRNA branch of RNA interference which is involved in defense against transposable elements in the animal germ line. Horizontal gene transfer (HGT), a dominant evolutionary process, at least, in prokaryotes, appears to be a form of (quasi)Lamarckian inheritance. The rate of HGT and the nature of acquired genes depend on the environment of the recipient organism and, in some cases, the transferred genes confer a selective advantage for growth in that environment, meeting the Lamarckian criteria. Various forms of stress-induced mutagenesis are tightly regulated and comprise a universal adaptive response to environmental stress in cellular life forms. Stress-induced mutagenesis can be construed as a quasi-Lamarckian phenomenon because the induced genomic changes, although random, are triggered by environmental factors and are beneficial to the organism.

Conclusion

Both Darwinian and Lamarckian modalities of evolution appear to be important, and reflect different aspects of the interaction between populations and the environment.

Reviewers

this article was reviewed by Juergen Brosius, Valerian Dolja, and Martijn Huynen. For complete reports, see the Reviewers' reports section.

XPC is involved in genome maintenance through multiple pathways in different tissues

In an attempt to evaluate the role of the Xpc gene in maintaining genomic stability in vivo under normal conditions, the age-dependent accumulation of spontaneous mutations in different tissues was analyzed in Xpc-deficient lacZ-transgenic mice. Brain, testis, and small intestine revealed no effects from the Xpc-deficiency, whereas liver, spleen, heart, and lung showed an enhanced age-related accumulation of mutations in Xpc-deficient mice. In the spleen, the effect was not obvious at 2 and 12 months of age, but became apparent at 23 months. The magnitude of the observed effect at an advanced age was similar in the liver, spleen and heart, but was comparatively smaller in the lung. Haploinsufficiency was observed in liver and spleen but not in heart and lung. Analysis of DNA sequences in the mutants revealed that the frequency of G:C to T:A changes were elevated in the liver and heart of Xpc-deficient aged mice, supporting the possible involvement of XPC in base excision repair of oxidized guanine. The occurrence of two or more mutations within a single lacZ gene was termed a multiple mutation and was also elevated in old Xpc-deficient mice. Among the clones examined, two mutant clones showed as many as four mutations within a short stretch of DNA. This is the first demonstration to support suggestions for the existence of a role for XPC in the suppression of multiple mutations. These multiple mutations could conceivably be generated by error-prone trans-lesional DNA synthesis. Overall, these results indicate that there may be diverse roles or mechanisms through which XPC participates in genome maintenance in different tissues.

Genomic mutation rates: what high-throughput methods can tell us

High-throughput DNA analyses are increasingly being used to detect rare mutations in moderately sized genomes. These methods have yielded genome mutation rates that are markedly higher than those obtained using pre-genomic strategies. Recent work in a variety of organisms has shown that mutation rate is strongly affected by sequence context and genome position. These observations suggest that high-throughput DNA analyses will ultimately allow researchers to identify trans-acting factors and cis sequences that underlie mutation rate variation. Such work should provide insights on how mutation rate variability can impact genome organization and disease progression.

Analysis of cancer mutation signatures in blood by a novel ultra-sensitive assay: monitoring of therapy or recurrence in non-metastatic breast cancer

BACKGROUND

Tumor DNA has been shown to be present both in circulating tumor cells in blood and as fragments in the plasma of metastatic cancer patients. The identification of ultra-rare tumor-specific mutations in blood would be the ultimate marker to measure efficacy of cancer therapy and/or early recurrence. Herein we present a method for detecting microinsertions/deletions/indels (MIDIs) at ultra-high analytical selectivity. MIDIs comprise about 15% of mutations.

METHODS AND FINDINGS

We describe MIDI-Activated Pyrophosphorolysis (MAP), a method of ultra-high analytical selectivity for detecting MIDIs. The high analytical selectivity of MAP is putatively due to serial coupling of two rare events: heteroduplex slippage and mis-pyrophosphorolysis. MAP generally has an analytical selectivity of one mutant molecule per >1 billion wild type molecules and an analytical sensitivity of one mutant molecule per reaction. The analytical selectivity of MAP is about 100,000-fold better than that of our previously described method of Pyrophosphorolysis Activated Polymerization-Allele specific amplification (PAP-A) for detecting MIDIs. The utility of this method is illustrated in two ways. 1) We demonstrate that two EGFR deletions commonly found in lung cancers are not present in tissue from four normal human lungs (10(7) copies of gDNA each) or in blood samples from 10 healthy individuals (10(7) copies of gDNA each). This is inconsistent, at least at an analytical sensitivity of 10(-7), with the hypotheses of (a) hypermutation or (b) strong selection of these growth factor-mutated cells during normal lung development leads to accumulation of pre-neoplastic cells with these EGFR mutations, which sometimes can lead to lung cancer in late adulthood. Moreover, MAP was used for large scale, high throughput "gene pool" analysis. No germline or early embryonic somatic mosaic mutation was detected (at a frequency of >0.3%) for the 15/18 bp EGFR deletion mutations in 6,400 individuals, suggesting that early embryonic EGFR somatic mutation is very rare, inconsistent with hypermutation or strong selection of these deletions in the embryo. 2) The second illustration of MAP utility is in personalized monitoring of therapy and early recurrence in cancer. Tumor-specific p53 mutations identified at diagnosis in the plasma of six patients with stage II and III breast cancer were undetectable after therapy in four women, consistent with clinical remission, and continued to be detected after treatment in two others, reflecting tumor progression.

CONCLUSIONS

MAP has an analytical selectivity of one part per billion for detection of MIDIs and an analytical sensitivity of one molecule. MAP provides a general tool for monitoring ultra-rare mutations in tissues and blood. As an example, we show that the personalized cancer signature in six out of six patients with non-metastatic breast cancer can be detected and that levels over time are correlated with the clinical course of disease.

High mutability of the tumor suppressor genes RASSF1 and RBSP3 (CTDSPL) in cancer

Background

Many different genetic alterations are observed in cancer cells. Individual cancer genes display point mutations such as base changes, insertions and deletions that initiate and promote cancer growth and spread. Somatic hypermutation is a powerful mechanism for generation of different mutations. It was shown previously that somatic hypermutability of proto-oncogenes can induce development of lymphomas.

Methodology/Principal Findings

We found an exceptionally high incidence of single-base mutations in the tumor suppressor genes RASSF1 and RBSP3 (CTDSPL) both located in 3p21.3 regions, LUCA and AP20 respectively. These regions contain clusters of tumor suppressor genes involved in multiple cancer types such as lung, kidney, breast, cervical, head and neck, nasopharyngeal, prostate and other carcinomas. Altogether in 144 sequenced RASSF1A clones (exons 1–2), 129 mutations were detected (mutation frequency, MF = 0.23 per 100 bp) and in 98 clones of exons 3–5 we found 146 mutations (MF = 0.29). In 85 sequenced RBSP3 clones, 89 mutations were found (MF = 0.10). The mutations were not cytidine-specific, as would be expected from alterations generated by AID/APOBEC family enzymes, and appeared de novo during cell proliferation. They diminished the ability of corresponding transgenes to suppress cell and tumor growth implying a loss of function. These high levels of somatic mutations were found both in cancer biopsies and cancer cell lines.

Conclusions/Significance

This is the first report of high frequencies of somatic mutations in RASSF1 and RBSP3 in different cancers suggesting it may underlay the mutator phenotype of cancer. Somatic hypermutations in tumor suppressor genes involved in major human malignancies offer a novel insight in cancer development, progression and spread.

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.