The mutational hazard hypothesis of organelle genome evolution: 10 years on

Why is there such a large variation in size and noncoding DNA content among organelle genomes? One explanation is that this genomic variation results from differences in the rates of organelle mutation and random genetic drift, as opposed to being the direct product of natural selection. Along these lines, the mutational hazard hypothesis (MHH) holds that ‘excess’ DNA is a mutational liability (because it increases the potential for harmful mutations) and, thus, has a greater tendency to accumulate in an organelle system with a low mutation rate as opposed to one with a high rate of mutation. Various studies have explored this hypothesis and, more generally, the relationship between organelle genome architecture and the mode and efficiency of organelle DNA repair. Although some of these investigations are in agreement with the MHH, others have contradicted it; nevertheless, they support a central role of mutation, DNA maintenance pathways and random genetic drift in fashioning organelle chromosomes. Arguably, one of the most important contributions of the MHH is that it has sparked crucial, widespread discussions about the importance of nonadaptive processes in genome evolution.     

SPONTANEOUS MUTATION ACCUMULATION IN MULTIPLE STRAINS OF THE GREEN ALGA,CHLAMYDOMONAS REINHARDTII

Estimates of mutational parameters, such as the average fitness effect of a new mutation and the rate at which new genetic variation for fitness is created by mutation, are important for the understanding of many biological processes. However, the causes of interspecific variation in mutational parameters and the extent to which they vary within species remain largely unknown. We maintained multiple strains of the unicellular eukaryote Chlamydomonas reinhardtii, for approximately 1000 generations under relaxed selection by transferring a single cell every ～10 generations. Mean fitness of the lines tended to decline with generations of mutation accumulation whereas mutational variance increased. We did not find any evidence for differences among strains in any of the mutational parameters estimated. The overall change in mean fitness per cell division and rate of input of mutational variance per cell division were more similar to values observed in multicellular organisms than to those in other single‐celled microbes. However, after taking into account differences in genome size among species, estimates from multicellular organisms and microbes, including our new estimates from C. reinhardtii, become substantially more similar. Thus, we suggest that variation in genome size is an important determinant of interspecific variation in mutational parameters.     

Behavioral degradation under mutation accumulation in Caenorhabditis elegans

Spontaneous mutations play a fundamental role in the maintenance of genetic variation in natural populations, the nature of inbreeding depression, the evolution of sexual reproduction, and the conservation of endangered species. Using long-term mutation-accumulation lines of the nematode Caenorhabditis elegans, we estimate the rate and magnitude of mutational effects for a suite of behaviors characterizing individual chemosensory responses to a repellant stimulus. In accordance with evidence that the vast majority of mutations are deleterious, we find that behavioral responses degrade over time as a result of spontaneous mutation accumulation. The rate of mutation for behavioral traits is roughly of the same order or slightly smaller than those previously estimated for reproductive traits and the average size of the mutational effects is also comparable. These results have important implications for the maintenance of genetic variation for behavior in natural populations as well as for expectations for behavioral change within endangered species and captive populations.     

Assessing Substitution Variation Across Sites in Grass Chloroplast DNA

We assess the similarity of base substitution processes, described by empirically derived 4 × 4 matrices, using chi-square homogeneity tests. Such significance analyses allow us to assess variation in sequence evolution across sites and we apply them to matrices derived from noncoding sites in different contexts in grass chloroplast DNA. We show that there is statistically significant variation in rates and patterns of mutation among noncoding sites in different contexts and then demonstrate a similar and significant influence of context on substitutions at fourfold degenerate sites of coding regions from grass chloroplast DNA. These results show that context has the same general effect on substitution bias in coding and noncoding DNA: the A+T content of flanking bases is correlated with rate of substitution, transition bias, and GC → AT pressure, while the number of flanking pyrimidines on a single strand is correlated with a mutational bias, or skew, toward pyrimidines. Despite the similarity in general trends, however, when we compare coding and noncoding matrices we find that there is a statistically significant difference between them even when we control for context. Most noticeably, fourfold degenerate sites in coding sequences are undergoing substitution at a higher rate and there are also significant differences in the relationship between pyrimidines skew and the number of flanking pyrimidines. Possible reasons for the differences between coding and noncoding sites are discussed. Furthermore, our analysis illustrates a simple statistical way for comparing substitution processes across sites allowing us to better study variation in evolutionary processes across a genome. [Reviewing Editor: Dr. Martin Kreitman]     

Evolution of mutational robustness in an RNA virus

Mutational (genetic) robustness is phenotypic constancy in the face of mutational changes to the genome. Robustness is critical to the understanding of evolution because phenotypically expressed genetic variation is the fuel of natural selection. Nonetheless, the evidence for adaptive evolution of mutational robustness in biological populations is controversial. Robustness should be selectively favored when mutation rates are high, a common feature of RNA viruses. However, selection for robustness may be relaxed under virus co-infection because complementation between virus genotypes can buffer mutational effects. We therefore hypothesized that selection for genetic robustness in viruses will be weakened with increasing frequency of co-infection. To test this idea, we used populations of RNA phage φ6 that were experimentally evolved at low and high levels of co-infection and subjected lineages of these viruses to mutation accumulation through population bottlenecking. The data demonstrate that viruses evolved under high co-infection show relatively greater mean magnitude and variance in the fitness changes generated by addition of random mutations, confirming our hypothesis that they experience weakened selection for robustness. Our study further suggests that co-infection of host cells may be advantageous to RNA viruses only in the short term. In addition, we observed higher mutation frequencies in the more robust viruses, indicating that evolution of robustness might foster less-accurate genome replication in RNA viruses.     

Epistasis in polygenic traits and the evolution of genetic architecture under stabilizing selection

We consider the effects of epistasis in a polygenic trait in the balance of mutation and stabilizing selection. The main issues are the genetic variation maintained in equilibrium and the evolution of the mutational effect distribution. The model assumes symmetric mutation and a continuum of alleles at all loci. Epistasis is modeled proportional to pairwise products of the single-locus effects. A general analytical formalism is developed. Assuming linkage equilibrium, we derive results for the equilibrium mutation load and the genetic and mutational variance in the house of cards and the Gaussian approximation. The additive genetic variation maintained in mutation-selection balance is reduced by any pattern of the epistatic interactions. The mutational variance, in contrast, is often increased. Large differences in mutational effects among loci emerge, and a negative correlation among (standard mean) locus mutation effects and mutation rates is predicted. Contrary to the common view since Waddington, we find that stabilizing selection in general does not lead to canalization of the trait. We propose that canalization as a target of selection instead occurs at the genic level. Here, primarily genes with a high mutation rate are buffered, often at the cost of decanalization of other genes. An intuitive interpretation of this view is given in the discussion.     

Frequent recovery of triplet mutations in UVB-exposed skin epidermis of Xpc-knockout mice

Mutations of the Xpc gene cause a deficiency in global genome repair, a subpathway of nucleotide excision repair (NER), in mammalian cells. We used transgenic mice harboring the lambda-phage-based lacZ mutational reporter gene to study the effect of an Xpc null mutation (Xpc-/-) on damage induction, repair and mutagenesis in mouse skin epidermis after UVB irradiation. UVB induced equal amounts of cyclobutane pyrimidine dimers (CPDs) and pyrimidine(6-4)pyrimidone photoproducts (64PPs) in mouse skin epidermis of Xpc-/- and wild-type mice. CPDs were not significantly removed in either of the mouse genotypes by 12h after irradiation, whereas removal of 64PPs was observed in the wild-type. Irradiation with 300 and 400J/m2 UVB increased the lacZ mutant frequency in the Xpc-/- epidermis to at least twice as high as in the wild-type. Ninety-nine lacZ mutants isolated from the UVB-exposed epidermis of Xpc(-/-)mice were analyzed and compared with mutant sequences from irradiated wild-type mice. The spectra of the mutations in the two genotypes were both highly UV-specific and similar in the dominance of C-->T transitions at dipyrimidine sites; however, Xpc-/- mice had a higher frequency of two-base tandem substitutions, including CC-->TT mutations, three-base tandem substitutions and double base substitutions that were separated by one unchanged base in a three-base sequence (alternating mutations). These tandem/alternating mutations included a remarkably large number of triplet mutations, a recently reported, novel type of UV-specific mutation, characterized by multiple base substitutions or frameshifts within a three-nucleotide sequence containing a dipyrimidine. We concluded that the triplet mutation is a UV-specific mutation that preferably occurs in NER deficient genetic backgrounds.     

Epistasis can lead to fragmented neutral spaces and contingency in evolution

In evolution, the effects of a single deleterious mutation can sometimes be compensated for by a second mutation which recovers the original phenotype. Such epistatic interactions have implications for the structure of genome space--namely, that networks of genomes encoding the same phenotype may not be connected by single mutational moves. We use the folding of RNA sequences into secondary structures as a model genotype-phenotype map and explore the neutral spaces corresponding to networks of genotypes with the same phenotype. In most of these networks, we find that it is not possible to connect all genotypes to one another by single point mutations. Instead, a network for a phenotypic structure with n bonds typically fragments into at least 2(n) neutral components, often of similar size. While components of the same network generate the same phenotype, they show important variations in their properties, most strikingly in their evolvability and mutational robustness. This heterogeneity implies contingency in the evolutionary process.     

Mutational History of a Human Cell Lineage from Somatic to Induced Pluripotent Stem Cells

The accuracy of replicating the genetic code is fundamental. DNA repair mechanisms protect the fidelity of the genome ensuring a low error rate between generations. This sustains the similarity of individuals whilst providing a repertoire of variants for evolution. The mutation rate in the human genome has recently been measured to be 50–70 de novo single nucleotide variants (SNVs) between generations. During development mutations accumulate in somatic cells so that an organism is a mosaic. However, variation within a tissue and between tissues has not been analysed. By reprogramming somatic cells into induced pluripotent stem cells (iPSCs), their genomes and the associated mutational history are captured. By sequencing the genomes of polyclonal and monoclonal somatic cells and derived iPSCs we have determined the mutation rates and show how the patterns change from a somatic lineage in vivo through to iPSCs. Somatic cells have a mutation rate of 14 SNVs per cell per generation while iPSCs exhibited a ten-fold lower rate. Analyses of mutational signatures suggested that deamination of methylated cytosine may be the major mutagenic source in vivo, whilst oxidative DNA damage becomes dominant in vitro. Our results provide insights for better understanding of mutational processes and lineage relationships between human somatic cells. Furthermore it provides a foundation for interpretation of elevated mutation rates and patterns in cancer.     

Phylogenomic analysis of the uracil-DNA glycosylase superfamily

The spontaneous deamination of cytosine produces uracil mispaired with guanine in DNA, which will produce a mutation, unless repaired. In all domains of life, uracil-DNA glycosylases (UDGs) are responsible for the elimination of uracil from DNA. Thus, UDGs contribute to the integrity of the genetic information and their loss results in mutator phenotypes. We are interested in understanding the role of UDG genes in the evolutionary variation of the rate and the spectrum of spontaneous mutations. To this end, we determined the presence or absence of the five main UDG families in more than 1,000 completely sequenced genomes and analyzed their patterns of gene loss and gain in eubacterial lineages. We observe nonindependent patterns of gene loss and gain between UDG families in Eubacteria, suggesting extensive functional overlap in an evolutionary timescale. Given that UDGs prevent transitions at G:C sites, we expected the loss of UDG genes to bias the mutational spectrum toward a lower equilibrium G + C content. To test this hypothesis, we used phylogenetically independent contrasts to compare the G + C content at intergenic and 4-fold redundant sites between lineages where UDG genes have been lost and their sister clades. None of the main UDG families present in Eubacteria was associated with a higher G + C content at intergenic or 4-fold redundant sites. We discuss the reasons of this negative result and report several features of the evolution of the UDG superfamily with implications for their functional study. uracil-DNA glycosylase, mutation rate evolution, mutational bias, GC content, DNA repair, mutator gene.     

DNA indels in coding regions reveal selective constraints on protein evolution in the human lineage

Background  

Insertions and deletions of DNA segments (indels) are together with substitutions the major mutational processes that generate genetic variation. Here we focus on recent DNA insertions and deletions in protein coding regions of the human genome to investigate selective constraints on indels in protein evolution.     

Elevated temperature increases genome-wide selection on de novo mutations

Adaptation in new environments depends on the amount of genetic variation available for evolution, and the efficacy by which natural selection discriminates among this variation. However, whether some ecological factors reveal more genetic variation, or impose stronger selection pressures than others, is typically not known. Here, we apply the enzyme kinetic theory to show that rising global temperatures are predicted to intensify natural selection throughout the genome by increasing the effects of DNA sequence variation on protein stability. We test this prediction by (i) estimating temperature-dependent fitness effects of induced mutations in seed beetles adapted to ancestral or elevated temperature, and (ii) calculate 100 paired selection estimates on mutations in benign versus stressful environments from unicellular and multicellular organisms. Environmental stress per se did not increase mean selection on de novo mutation, suggesting that the cost of adaptation does not generally increase in new ecological settings to which the organism is maladapted. However, elevated temperature increased the mean strength of selection on genome-wide polymorphism, signified by increases in both mutation load and mutational variance in fitness. These results have important implications for genetic diversity gradients and the rate and repeatability of evolution under climate change.     

Heterogeneous Nature and Distribution of Interruptions in Dinucleotides May Indicate the Existence of Biased Substitutions Underlying Microsatellite Evolution

Some aspects of microsatellite evolution, such as the role of base substitutions, are far from being fully understood. To examine the significance of base substitutions underlying the evolution of microsatellites we explored the nature and the distribution of interruptions in dinucleotide repeats from the human genome. The frequencies that we inferred in the repetitive sequences were statistically different from the frequencies observed in other noncoding sequences. Additionally, we detected that the interruptions tended to be towards the ends of the microsatellites and 5'-3' asymmetry. In all the estimates nucleotides forming the same repetitive motif seem to be affected by different base substitution rates in AC and AG. This tendency itself could generate patterning and similarity in flanking sequences and reconcile these phenomena with the high mutation rate found in flanking sequences without invoking convergent evolution. Nevertheless, our data suggest that there is a regional bias in the substitution pattern of microsatellites. The accumulation of random substitutions alone cannot explain the heterogeneity and the asymmetry of interruptions found in this study or the relative frequency of different compound microsatellites in the human genome. Therefore, we cannot rule out the possibility of a mutational bias leading to convergent or parallel evolution in flanking sequences.     

Human Microsatellites: Mutation and Evolution

Microsatellites (MSs) are short tandem DNA repeats with the repetitive motif of two to six nucleotides, forming tracts up to hundreds of nucleotides long. Notwithstanding the active use of MSs in genetic studies of various biological problems, the reasons for their wide occurrence in the genome, their possible functions, and mutational behavior are still unclear. The mutation rate in MS repeats is on average several orders of magnitude higher than in the remaining DNA, which allows for direct estimation of evolutionary transformation rate in nucleotide sequences of the genome. Mutation process in MSs is highly heterogeneous, with distinct differences between species; furthermore, within a species it differs among loci with different repeat size, among alleles of one locus, and among individuals of different sex and age. Most MS mutations are caused by DNA slippage during replication but the probability of this event depends on the locus. In this review, a number of models of MS evolution are discussed, which account for the relationship between mutation rate and allele size, different mutation direction in alleles of different size, and the appearance of point mutations within repeat tracts restricting allele size. The MS evolution is considered mainly in the context of selective neutrality, although there is evidence showing functional significance of some variants of tandem repeats and thus their possible selective value.     

Rates of genomic divergence in humans,chimpanzees and their lice

The rate of DNA mutation and divergence is highly variable across the tree of life. However, the reasons underlying this variation are not well understood. Comparing the rates of genetic changes between hosts and parasite lineages that diverged at the same time is one way to begin to understand differences in genetic mutation and substitution rates. Such studies have indicated that the rate of genetic divergence in parasites is often faster than that of their hosts when comparing single genes. However, the variation in this relative rate of molecular evolution across different genes in the genome is unknown. We compared the rate of DNA sequence divergence between humans, chimpanzees and their ectoparasitic lice for 1534 protein-coding genes across their genomes. The rate of DNA substitution in these orthologous genes was on average 14 times faster for lice than for humans and chimpanzees. In addition, these rates were positively correlated across genes. Because this correlation only occurred for substitutions that changed the amino acid, this pattern is probably produced by similar functional constraints across the same genes in humans, chimpanzees and their ectoparasites.     

Competition between Transposable Elements and Mutator Genes in Bacteria

Although both genotypes with elevated mutation rate (mutators) and mobilization of insertion sequence (IS) elements have substantial impact on genome diversification, their potential interactions are unknown. Moreover, the evolutionary forces driving gradual accumulation of these elements are unclear: Do these elements spread in an initially transposon-free bacterial genome as they enable rapid adaptive evolution? To address these issues, we inserted an active IS1 element into a reduced Escherichia coli genome devoid of all other mobile DNA. Evolutionary laboratory experiments revealed that IS elements increase mutational supply and occasionally generate variants with especially large phenotypic effects. However, their impact on adaptive evolution is small compared with mismatch repair mutator alleles, and hence, the latter impede the spread of IS-carrying strains. Given their ubiquity in natural populations, such mutator alleles could limit early phase of IS element evolution in a new bacterial host. More generally, our work demonstrates the existence of an evolutionary conflict between mutation-promoting mechanisms.     

Integrating genomics,bioinformatics, and classical genetics to study the effects of recombination on genome evolution

This study presents compelling evidence that recombination significantly increases the silent GC content of a genome in a selectively neutral manner, resulting in a highly significant positive correlation between recombination and "GC3s" in the yeast Saccharomyces cerevisiae. Neither selection nor mutation can explain this relationship. A highly significant GC-biased mismatch repair system is documented for the first time in any member of the Kingdom Fungi. Much of the variation in the GC3s within yeast appears to result from GC-biased gene conversion. Evidence suggests that GC-biased mismatch repair exists in numerous organisms spanning six kingdoms. This transkingdom GC mismatch repair bias may have evolved in response to a ubiquitous AT mutational bias. A significant positive correlation between recombination and GC content is found in many of these same organisms, suggesting that the processes influencing the evolution of the yeast genome may be a general phenomenon. Nonrecombining regions of the genome and nonrecombining genomes would not be subject to this type of molecular drive. It is suggested that the low GC content characteristic of many nonrecombining genomes may be the result of three processes (1) a prevailing AT mutational bias, (2) random fixation of the most common types of mutation, and (3) the absence of the GC-biased gene conversion which, in recombining organisms, permits the reversal of the most common types of mutation. A model is proposed to explain the observation that introns, intergenic regions, and pseudogenes typically have lower GC content than the silent sites of corresponding open reading frames. This model is based on the observation that the greater the heterology between two sequences, the less likely it is that recombination will occur between them. According to this "Constraint" hypothesis, the formation and propagation of heteroduplex DNA is expected to occur, on average, more frequently within conserved coding and regulatory regions of the genome. In organisms possessing GC-biased mismatch repair, this would enhance the GC content of these regions through biased gene conversion. These findings have a number of important implications for the way we view genome evolution and suggest a new model for the evolution of sex.     

Comparison of the Ames test and a newly developed assay for detection of mutagenic pollution of marine environments

Mismatch repair (MMR) genes, such as Msh2, are classified as "mutator" genes, responsible for the microsatellite instability identified in many tumors. In the current study, the mutation frequency and mutational spectrum in thymic lymphoma arising in Msh2 deficient mice are investigated. Thymic lymphoma developed in Msh2-/- background displayed an eight to nine-fold increase in mutation frequency compared to the normal thymi in Msh2 deficient animals. Sequencing demonstrated significantly different mutational spectra between normal thymus tissue and thymic lymphomas in Msh2-/- mice (P=0.02). The tumor mutational spectrum is characterized by an increase in base substitutions occurring at A:T sites, and multiple mutations, as well as a minor increase in -1 frameshifts. We analyzed mutations in different parts of the tumors, and different regional hotspots could be identified. Several hotspot mutations that are a rare event in normal tissues were identified in the tumor tissues. We conclude that thymic lymphomas arising in Msh2 deficient genetic background are hypermutable and the altered mutational spectrum might be an indication of infidelity of DNA replication during tumorigenesis.     

Hydroxymethylated Cytosines Are Associated with Elevated C to G Transversion Rates

It has long been known that methylated cytosines deaminate at higher rates than unmodified cytosines and constitute mutational hotspots in mammalian genomes. The repertoire of naturally occurring cytosine modifications, however, extends beyond 5-methylcytosine to include its oxidation derivatives, notably 5-hydroxymethylcytosine. The effects of these modifications on sequence evolution are unknown. Here, we combine base-resolution maps of methyl- and hydroxymethylcytosine in human and mouse with population genomic, divergence and somatic mutation data to show that hydroxymethylated and methylated cytosines show distinct patterns of variation and evolution. Surprisingly, hydroxymethylated sites are consistently associated with elevated C to G transversion rates at the level of segregating polymorphisms, fixed substitutions, and somatic mutations in tumors. Controlling for multiple potential confounders, we find derived C to G SNPs to be 1.43-fold (1.22-fold) more common at hydroxymethylated sites compared to methylated sites in human (mouse). Increased C to G rates are evident across diverse functional and sequence contexts and, in cancer genomes, correlate with the expression of Tet enzymes and specific components of the mismatch repair pathway (MSH2, MSH6, and MBD4). Based on these and other observations we suggest that hydroxymethylation is associated with a distinct mutational burden and that the mismatch repair pathway is implicated in causing elevated transversion rates at hydroxymethylated cytosines.