Predicting mutation frequencies in stem-loop structures of derepressed genes: implications for evolution

L’impact des mutations neutres sur l’évolvabilité et l’évolution des génomes

Les mutations bénéfiques à forts effets sont rares et les mutations délétères sont éliminées par la sélection naturelle. La majorité des mutations qui s’accumulent dans les génomes ont donc des effets sélectifs très faibles, voire nuls ; elles sont alors appelées mutations neutres. Au cours des deux dernières décennies, il a été montré que les mutations, même en l’absence d’effet sur la valeur sélective des organismes, affectent leur évolvabilité, en donnant accès à de nouveaux phénotypes par le biais de mutations apparaissant ultérieurement, et qui n’auraient pas été disponibles autrement. En plus de cet effet, de nombreuses mutations neutres – indépendamment de leurs effets sélectifs – peuvent affecter la mutabilité de séquences d’ADN voisines, et moduler l’efficacité de la recombinaison homologue. De telles mutations ne modifient pas le spectre des phénotypes accessibles, mais plutôt la vitesse à laquelle de nouveaux phénotypes seront produits, un processus qui a des conséquences à long terme mais aussi potentiellement à court terme, en lien avec l’émergence de cancers.

A mutational hotspot that determines highly repeatable evolution can be built and broken by silent genetic changes

Mutational hotspots can determine evolutionary outcomes and make evolution repeatable. Hotspots are products of multiple evolutionary forces including mutation rate heterogeneity, but this variable is often hard to identify. In this work, we reveal that a near-deterministic genetic hotspot can be built and broken by a handful of silent mutations. We observe this when studying homologous immotile variants of the bacteria Pseudomonas fluorescens, AR2 and Pf0-2x. AR2 resurrects motility through highly repeatable de novo mutation of the same nucleotide in >95% lines in minimal media (ntrB A289C). Pf0-2x, however, evolves via a number of mutations meaning the two strains diverge significantly during adaptation. We determine that this evolutionary disparity is owed to just 6 synonymous variations within the ntrB locus, which we demonstrate by swapping the sites and observing that we are able to both break (>95% to 0%) and build (0% to 80%) a deterministic mutational hotspot. Our work reveals a key role for silent genetic variation in determining adaptive outcomes.

The Impact of Neutral Mutations on Genome Evolvability

Beneficial mutations are rare and deleterious mutations are purged by natural selection. As a result, the vast majority of mutations that accumulate in genomes belong to the class of neutral mutations. Over the last two decades, neutral mutations, despite their null effect on fitness, have been shown to affect evolvability by providing access to new phenotypes through subsequent mutations that would not have been available otherwise. Here we propose that in addition, many mutations - independent of their selective effects - can affect the mutability of neighboring DNA sequences and modulate the efficacy of homologous recombination. Such mutations do not change the spectrum of accessible phenotypes, but rather the rate at which new phenotypes will be produced. Therefore, neutral mutations that accumulate in genomes have an important long-term impact on the evolutionary fate of genomes.

Spontaneous Mutations in the Nitrate Reductase Gene napC Drive the Emergence of Eco-friendly Low-N2O-Emitting Alfalfa Rhizobia in Regions with Different Climates

We have recently shown that commercial alfalfa inoculants (e.g., Sinorhizobium meliloti B399), which are closely related to the denitrifier model strain Sinorhizobium meliloti 1021, have conserved nitrate, nitrite, and nitric oxide reductases associated with the production of the greenhouse gas nitrous oxide (N₂O) from nitrate but lost the N₂O reductase related to the degradation of N₂O to gas nitrogen. Here, we screened a library of nitrogen-fixing alfalfa symbionts originating from different ecoregions and containing N₂O reductase genes and identified novel rhizobia (Sinorhizobium meliloti INTA1-6) exhibiting exceptionally low N₂O emissions. To understand the genetic basis of this novel eco-friendly phenotype, we sequenced and analyzed the genomes of these strains, focusing on their denitrification genes, and found mutations only in the nitrate reductase structural gene napC. The evolutionary analysis supported that, in these natural strains, the denitrification genes were inherited by vertical transfer and that their defective nitrate reductase napC alleles emerged by independent spontaneous mutations. In silico analyses showed that mutations in this gene occurred in ssDNA loop structures with high negative free energy (-ΔG) and that the resulting mutated stem-loop structures exhibited increased stability, suggesting the occurrence of transcription-associated mutation events. In vivo assays supported that at least one of these ssDNA sites is a mutational hot spot under denitrification conditions. Similar benefits from nitrogen fixation were observed when plants were inoculated with the commercial inoculant B399 and strains INTA4-6, suggesting that the low-N₂O-emitting rhizobia can be an ecological alternative to the current inoculants without resigning economic profitability.

Therapeutic effect of green tea extract on alcohol induced hepatic mitochondrial DNA damage in albino wistar rats

The present study principally sought to investigate the effect of green tea extract (GTE) supplementation on hepatic mitochondrial DNA (mtDNA) damage in alcohol receiving rats. MtDNA was isolated from hepatic tissues of albino wistar rats after alcohol treatment with and without GTE supplementation. Entire displacement loop (D-loop) of mtDNA was screened by PCR-Sanger's sequencing method. In addition, mtDNA deletions and antioxidant activity were measured in hepatic tissue of all rats. Results showed increased frequency of D-loop mutations in alcoholic rats (ALC). DNA mfold analysis predicted higher free energy for 15507C and 16116C alleles compared to their corresponding wild alleles which represents less stable secondary structures with negative impact on overall mtDNA function. Interestingly, D-loop mutations observed in ALC rats were successfully restored on GTE supplementation. MtDNA deletions were observed in ALC rats, but intact native mtDNA was found in ALC + GTE group suggesting alcohol induced oxidative damage of mtDNA and ameliorative effect of GTE. Furthermore, markedly decreased activities of glutathione peroxidise, superoxide dismutase, catalase and glutathione content were identified in ALC rats; however, GTE supplementation significantly (P < 0.05) restored these levels close to normal. In conclusion, green tea could be used as an effective nutraceutical against alcohol induced mitochondrial DNA damage.

Cell-SELEX-based selection and characterization of a G-quadruplex DNA aptamer against mouse dendritic cells

Targeting of dendritic cells (DCs) by aptamers increases antigen capture and presentation to the immune system. Our aim was to produce aptamers against DC molecules using the cell-SELEX procedure. For this purpose, 18 rounds of cell-SELEX were performed on mouse macrophage J774A.1 and CT26 as target and control cells, respectively. The selected aptamers were truncated and their binding to mouse macrophages, and immature and mature DCs analyzed. Two macrophage-specific aptamers, Seq6 and Seq7, were identified. A truncated form of Seq7, Seq7-4, 33 nucleotides in length and containing the G-quadruplex, bound macrophages and immature DCs with KD values in the nanomolar range. We anticipate that Seq7-4 has potential as a therapeutic tool in targeting of mouse macrophages and immature DCs to efficiently improve different immunotherapy approaches.

Does transcription-associated DNA damage limit lifespan?

Small mammals undergo an aging process similar to that of larger mammals, but aging occurs at a dramatically faster rate. This phenomenon is often assumed to be the result of damage caused by reactive oxygen species generated in mitochondria. An alternative explanation for the phenomenon is suggested here. The rate of RNA synthesis is dramatically elevated in small mammals and correlates quantitatively with the rate of aging among different mammalian species. The rate of RNA synthesis is reduced by caloric restriction and inhibition of TOR pathway signaling, two perturbations that increase lifespan in multiple metazoan species. From bacteria to man, the transcription of a gene has been found to increase the rate at which it is damaged, and a number of lines of evidence suggest that DNA damage is sufficient to induce multiple symptoms associated with normal aging. Thus, the correlations frequently found between the rate of RNA synthesis and the rate of aging could potentially reflect an important role for transcription-associated DNA damage in the aging process.

Genetic instability in budding and fission yeast-sources and mechanisms

Cells are constantly confronted with endogenous and exogenous factors that affect their genomes. Eons of evolution have allowed the cellular mechanisms responsible for preserving the genome to adjust for achieving contradictory objectives: to maintain the genome unchanged and to acquire mutations that allow adaptation to environmental changes. One evolutionary mechanism that has been refined for survival is genetic variation. In this review, we describe the mechanisms responsible for two biological processes: genome maintenance and mutation tolerance involved in generations of genetic variations in mitotic cells of both Saccharomyces cerevisiae and Schizosaccharomyces pombe. These processes encompass mechanisms that ensure the fidelity of replication, DNA lesion sensing and DNA damage response pathways, as well as mechanisms that ensure precision in chromosome segregation during cell division. We discuss various factors that may influence genome stability, such as cellular ploidy, the phase of the cell cycle, transcriptional activity of a particular region of DNA, the proficiency of DNA quality control systems, the metabolic stage of the cell and its respiratory potential, and finally potential exposure to endogenous or environmental stress.

The stability of budding and fission yeast genomes is influenced by two contradictory factors: (1) the need to be fully functional, which is ensured through the replication fidelity pathways of nuclear and mitochondrial genomes through sensing and repairing DNA damage, through precise chromosome segregation during cell division; and (2) the need to acquire changes for adaptation to environmental challenges.

Transcription-associated mutagenesis

Transcription requires unwinding complementary DNA strands, generating torsional stress, and sensitizing the exposed single strands to chemical reactions and endogenous damaging agents. In addition, transcription can occur concomitantly with the other major DNA metabolic processes (replication, repair, and recombination), creating opportunities for either cooperation or conflict. Genetic modifications associated with transcription are a global issue in the small genomes of microorganisms in which noncoding sequences are rare. Transcription likewise becomes significant when one considers that most of the human genome is transcriptionally active. In this review, we focus specifically on the mutagenic consequences of transcription. Mechanisms of transcription-associated mutagenesis in microorganisms are discussed, as is the role of transcription in somatic instability of the vertebrate immune system.

Frequency-dependent selection can lead to evolution of high mutation rates

Theoretical and experimental studies have shown that high mutation rates can be advantageous, especially in novel or fluctuating environments. Here we examine how frequency-dependent competition may lead to fluctuations in trait frequencies that exert upward selective pressure on mutation rates. We use a mathematical model to show that cyclical trait dynamics generated by "rock-paper-scissors" competition can cause the mutation rate in a population to converge to a high evolutionarily stable mutation rate, reflecting a trade-off between generating novelty and reproducing past success. Introducing recombination lowers the evolutionarily stable mutation rate but allows stable coexistence between mutation rates above and below the evolutionarily stable rate. Even considering strong mutational load and ignoring the costs of faithful replication, evolution favors positive mutation rates if the selective advantage of prevailing in competition exceeds the ratio of recombining to nonrecombining offspring. We discuss a number of genomic mechanisms that may meet our theoretical requirements for the adaptive evolution of mutation. Overall, our results suggest that local mutation rates may be higher on genes influencing cyclical competition and that global mutation rates in asexual species may be higher in populations subject to strong cyclical competition.

Conservation of mRNA secondary structures may filter out mutations in Escherichia coli evolution

Recent reports indicate that mutations in viral genomes tend to preserve RNA secondary structure, and those mutations that disrupt secondary structural elements may reduce gene expression levels, thereby serving as a functional knockout. In this article, we explore the conservation of secondary structures of mRNA coding regions, a previously unknown factor in bacterial evolution, by comparing the structural consequences of mutations in essential and nonessential Escherichia coli genes accumulated over 40 000 generations in the course of the ‘long-term evolution experiment’. We monitored the extent to which mutations influence minimum free energy (MFE) values, assuming that a substantial change in MFE is indicative of structural perturbation. Our principal finding is that purifying selection tends to eliminate those mutations in essential genes that lead to greater changes of MFE values and, therefore, may be more disruptive for the corresponding mRNA secondary structures. This effect implies that synonymous mutations disrupting mRNA secondary structures may directly affect the fitness of the organism. These results demonstrate that the need to maintain intact mRNA structures imposes additional evolutionary constraints on bacterial genomes, which go beyond preservation of structure and function of the encoded proteins.

Kinetic models reveal the in vivo mechanisms of mutagenesis in microbes and man

This review summarizes the evidence indicating that mutagenic mechanisms in vivo are essentially the same in all living cells. Unique metabolic reactions to a particular environmental stress apparently target specific genes for increased rates of transcription and mutation, resulting in higher mutation rates for those genes most likely to solve the problem. Kinetic models which have demonstrated predictive value are described and are shown to simulate mutagenesis in vivo in Escherichia coli, the p53 tumor suppressor gene, and somatic hypermutation. In all three models, direct correlations are seen between mutation frequencies and transcription rates. G and C nucleosides in single-stranded DNA (ssDNA) are intrinsically mutable, and G and C silent mutations in p53 and in VH framework regions provide compelling evidence for intrinsic mechanisms of mutability, since mutation outcomes are neutral and are not selected. During transcription, the availability of unpaired bases in the ssDNA of secondary structures is rate-limiting for, and determines the frequency of mutations in vivo. In vitro analyses also verify the conclusion that intrinsically mutable bases are in fact located in ssDNA loops of predicted stem-loop structures (SLSs).

Gene Conversion-Like Events in the Diversification of Human Rearranged IGHV3-23*01 Gene Sequences

Gene conversion (GCV), a mechanism mediated by activation-induced cytidine deaminase (AID) is well established as a mechanism of immunoglobulin diversification in a few species. However, definitive evidence of GCV-like events in human immunoglobulin genes is scarce. The lack of evidence of GCV in human rearranged immunoglobulin gene sequences is puzzling given the presence of highly similar germline donors and the presence of all the enzymatic machinery required for GCV. In this study, we undertook a computational analysis of rearranged IGHV3-23(*)01 gene sequences from common variable immunodeficiency (CVID) patients, AID-deficient patients, and healthy individuals to survey "GCV-like" activities. We analyzed rearranged IGHV3-23(*)01 gene sequences obtained from total PBMC RNA and single-cell polymerase chain reaction of individual B cell lysates. Our search identified strong evidence of GCV-like activity. We observed that GCV-like tracts are flanked by AID hotspot motifs. Structural modeling of IGHV3-23(*)01 gene sequence revealed that hypermutable bases flanking GCV-like tracts are in the single stranded DNA (ssDNA) of stable stem-loop structures (SLSs). ssDNA is inherently fragile and also an optimal target for AID. We speculate that GCV could have been initiated by the targeting of hypermutable bases in ssDNA state in stable SLSs, plausibly by AID. We have observed that the frequency of GCV-like events is significantly higher in rearranged IGHV3-23-(*)01 sequences from healthy individuals compared to that of CVID patients. We did not observe GCV-like events in rearranged IGHV3-23-(*)01 sequences from AID-deficient patients. GCV, unlike somatic hypermutation (SHM), can result in multiple base substitutions that can alter many amino acids. The extensive changes in antibody affinity by GCV-like events would be instrumental in protecting humans against pathogens that diversify their genome by antigenic shift.

Genetic variability of human respiratory syncytial virus A strains circulating in Ontario: a novel genotype with a 72 nucleotide G gene duplication

Human respiratory syncytial virus (HRSV) is the main cause of acute lower respiratory infections in children under 2 years of age and causes repeated infections throughout life. We investigated the genetic variability of RSV-A circulating in Ontario during 2010–2011 winter season by sequencing and phylogenetic analysis of the G glycoprotein gene.

Among the 201 consecutive RSV isolates studied, RSV-A (55.7%) was more commonly observed than RSV-B (42.3%). 59.8% and 90.1% of RSV-A infections were among children ≤12 months and ≤5 years old, respectively. On phylogenetic analysis of the second hypervariable region of the 112 RSV-A strains, 110 (98.2%) clustered within or adjacent to the NA1 genotype; two isolates were GA5 genotype. Eleven (10%) NA1-related isolates clustered together phylogenetically as a novel RSV-A genotype, named ON1, containing a 72 nucleotide duplication in the C-terminal region of the attachment (G) glycoprotein. The predicted polypeptide is lengthened by 24 amino acids and includes a23 amino acid duplication. Using RNA secondary structural software, a possible mechanism of duplication occurrence was derived. The 23 amino acid ON1 G gene duplication results in a repeat of 7 potential O-glycosylation sites including three O-linked sugar acceptors at residues 270, 275, and 283. Using Phylogenetic Analysis by Maximum Likelihood analysis, a total of 19 positively selected sites were observed among Ontario NA1 isolates; six were found to be codons which reverted to the previous state observed in the prototype RSV-A2 strain. The tendency of codon regression in the G-ectodomain may infer a decreased avidity of antibody to the current circulating strains. Further work is needed to document and further understand the emergence, virulence, pathogenicity and transmissibility of this novel RSV-A genotype with a72 nucleotide G gene duplication.

Stabilised DNA secondary structures with increasing transcription localise hypermutable bases for somatic hypermutation in IGHV3-23

Somatic hypermutation (SHM) mediated by activation-induced cytidine deaminase (AID) is a transcription-coupled mechanism most responsible for generating high affinity antibodies. An issue remaining enigmatic in SHM is how AID is preferentially targeted during transcription to hypermutable bases in its substrates (WRC motifs) on both DNA strands. AID targets only single stranded DNA. By modelling the dynamical behaviour of IGHV3-23 DNA, a commonly used human variable gene segment, we observed that hypermutable bases on the non-transcribed strand are paired whereas those on transcribed strand are mostly unpaired. Hypermutable bases (both paired and unpaired) are made accessible to AID in stabilised secondary structures formed with increasing transcription levels. This observation provides a rationale for the hypermutable bases on both the strands of DNA being targeted to a similar extent despite having differences in unpairedness. We propose that increasing transcription and RNAP II stalling resulting in the formation and stabilisation of stem-loop structures with AID hotspots in negatively supercoiled region can localise the hypermutable bases of both strands of DNA, to AID-mediated SHM.

Transcription as a source of genome instability

Alterations in genome sequence and structure contribute to somatic disease, affect the fitness of subsequent generations and drive evolutionary processes. The crucial roles of highly accurate replication and efficient repair in maintaining overall genome integrity are well-known, but the more localized stability costs that are associated with transcribing DNA into RNA molecules are less appreciated. Here we review the diverse ways in which the essential process of transcription alters the underlying DNA template and thereby modifies the genetic landscape.

Evolution of coordinated mutagenesis and somatic hypermutation in VH5

The VH5 human antibody gene was analyzed using a computer program (mfg) which simulates transcription, to better understand transcription-driven mutagenesis events that occur during "phase 1" of somatic hypermutation. Results show that the great majority of mutations in the non-transcribed strand occur within loops of two predicted high-stability stem-loop structures, termed SLSs 14.9 and 13.9. In fact, 89% of the 2505 mutations reported are within the encoded complementarity-determining region (CDR) and occur in loops of these high-stability structures. In vitro studies were also done and verified the existence of SLS 14.9. Following the formation of SLSs 14.9 and 13.9, a sustained period of transcriptional activity occurs within a window size of 60-70 nucleotides. During this period, the stability of these two SLSs does not change, and may provide the substrate for base exchanges and mutagenesis. The data suggest that many mutable bases are exposed simultaneously at pause sites, allowing for coordinated mutagenesis.

Self-catalyzed site-specific depurination of G residues mediated by cruciform extrusion in closed circular DNA plasmids

A major variety of "spontaneous" genomic damage is endogenous generation of apurinic sites. Depurination rates vary widely across genomes, occurring with higher frequency at "depurination hot spots." Recently, we discovered a site-specific self-catalyzed depurinating activity in short (14-18 nucleotides) DNA stem-loop-forming sequences with a 5'-G(T/A)GG-3' loop and T·A or G·C as the first base pair at the base of the loop; the 5'-G residue of the loop self-depurinates at least 10(5)-fold faster than random "spontaneous" depurination at pH 5. Formation of the catalytic intermediate for self-depurination in double-stranded DNA requires a stem-loop to extrude as part of a cruciform. In this study, evidence is presented for self-catalyzed depurination mediated by cruciform formation in plasmid DNA in vitro. Cruciform extrusion was confirmed, and its extent was quantitated by digestion of the plasmid with single strand-specific mung bean endonuclease, followed by restriction digestion and sequencing of resulting mung bean-generated fragments. Appearance of the apurinic site in the self-depurinating stem-loop was confirmed by digestion of plasmid DNA with apurinic endonuclease IV, followed by primer extension and/or PCR amplification to detect the endonuclease-generated strand break and identify its location. Self-catalyzed depurination was contingent on the plasmid being supercoiled and was not observed in linearized plasmids, consistent with the presence of the extruded cruciform in the supercoiled plasmid and not in the linear one. These results indicate that self-catalyzed depurination is not unique to single-stranded DNA; rather, it can occur in stem-loop structures extruding from double-stranded DNA and therefore could, in principle, occur in vivo.

The roles of transcription and genotoxins underlying p53 mutagenesis in vivo

Transcription drives supercoiling which forms and stabilizes single-stranded (ss) DNA secondary structures with loops exposing G and C bases that are intrinsically mutable and vulnerable to non-enzymatic hydrolytic reactions. Since many studies in prokaryotes have shown direct correlations between the frequencies of transcription and mutation, we conducted in silico analyses using the computer program, mfg, which simulates transcription and predicts the location of known mutable bases in loops of high-stability secondary structures. Mfg analyses of the p53 tumor suppressor gene predicted the location of mutable bases and mutation frequencies correlated with the extent to which these mutable bases were exposed in secondary structures. In vitro analyses have now confirmed that the 12 most mutable bases in p53 are in fact located in predicted ssDNA loops of these structures. Data show that genotoxins have two independent effects on mutagenesis and the incidence of cancer: Firstly, they activate p53 transcription, which increases the number of exposed mutable bases and also increases mutation frequency. Secondly, genotoxins increase the frequency of G-to-T transversions resulting in a decrease in G-to-A and C mutations. This precise compensatory shift in the 'fate' of G mutations has no impact on mutation frequency. Moreover, it is consistent with our proposed mechanism of mutagenesis in which the frequency of G exposure in ssDNA via transcription is rate limiting for mutation frequency in vivo.

Functional and metabolic effects of adaptive glycerol kinase (GLPK) mutants in Escherichia coli

Herein we measure the effect of four adaptive non-synonymous mutations to the glycerol kinase (glpK) gene on catalytic function and regulation, to identify changes that correlate to increased fitness in glycerol media. The mutations significantly reduce affinity for the allosteric inhibitor fructose-1,6-bisphosphate (FBP) and formation of the tetramer, which are structurally related, in a manner that correlates inversely with imparted fitness during growth on glycerol, which strongly suggests that these enzymatic parameters drive growth improvement. Counterintuitively, the glpK mutations also increase glycerol-induced auto-catabolite repression that reduces glpK transcription in a manner that correlates to fitness. This suggests that increased specific GlpK activity is attenuated by negative feedback on glpK expression via catabolite repression, possibly to prevent methylglyoxal toxicity. We additionally report that glpK mutations were fixed in 47 of 50 independent glycerol-adapted lineages. By far the most frequently mutated locus (nucleotide 218) was mutated in 20 lineages, strongly suggesting this position has an elevated mutation rate. This study demonstrates that fitness correlations can be used to interrogate adaptive processes at the protein level and to identify the regulatory constraints underlying selection and improved growth.

Transcription-associated mutagenesis increases protein sequence diversity more effectively than does random mutagenesis in Escherichia coli

Background

During transcription, the nontranscribed DNA strand becomes single-stranded DNA (ssDNA), which can form secondary structures. Unpaired bases in the ssDNA are less protected from mutagens and hence experience more mutations than do paired bases. These mutations are called transcription-associated mutations. Transcription-associated mutagenesis is increased under stress and depends on the DNA sequence. Therefore, selection might significantly influence protein-coding sequences in terms of the transcription-associated mutability per transcription event under stress to improve the survival of Escherichia coli.

Methodology/Principal Findings

The mutability index (MI) was developed by Wright et al. to estimate the relative transcription-associated mutability of bases per transcription event. Using the most stable fold of each ssDNA that have an average length n, MI was defined as (the number of folds in which the base is unpaired)/n×(highest –ΔG of all n folds in which the base is unpaired), where ΔG is the free energy. The MI values show a significant correlation with mutation data under stress but not with spontaneous mutations in E. coli. Protein sequence diversity is preferred under stress but not under favorable conditions. Therefore, we evaluated the selection pressure on MI in terms of the protein sequence diversity for all the protein-coding sequences in E. coli. The distributions of the MI values were lower at bases that could be substituted with each of the other three bases without affecting the amino acid sequence than at bases that could not be so substituted. Start codons had lower distributions of MI values than did nonstart codons.

Conclusions/Significance

Our results suggest that the majority of protein-coding sequences have evolved to promote protein sequence diversity and to reduce gene knockout under stress. Consequently, transcription-associated mutagenesis increases protein sequence diversity more effectively than does random mutagenesis under stress. Nonrandom transcription-associated mutagenesis under stress should improve the survival of E. coli.

Evidence for variable selective pressures at a large secondary structure of the human mitochondrial DNA control region

A combined effect of functional constraints and random mutational events is responsible for the sequence evolution of the human mitochondrial DNA (mtDNA) control region. Most studies targeting this noncoding segment usually focus on its primary sequence information disregarding other informative levels such as secondary or tertiary DNA conformations. In this work, we combined the most recent developments in DNA folding calculations with a phylogenetic comparative approach in order to investigate the formation of intrastrand secondary structures in the human mtDNA control region. Our most striking results are those regarding a new cloverleaf-like secondary structure predicted for a 93-bp stretch of the control region 5'-peripheral domain. Randomized sequences indicated that this structure has a more negative folding energy than the average of random sequences with the same nucleotide composition. In addition, a sliding window scan across the complete mitochondrial genome revealed that it stands out as having one of the highest folding potential. Moreover, we detected several lines of evidence of both negative and positive selection on this structure with high levels of conservation at the structure-relevant stem regions and the occurrence of compensatory base changes in the primate lineage. In the light of previous data, we discuss the possible involvement of this structure in mtDNA replication and/or transcription. We conclude that maintenance of this structure is responsible for the observed heterogeneity in the rate of substitution among sites in part of the human hypervariable region I and that it is a hot spot for the 3' end of human mtDNA deletions.

Analyzing the evolution of RNA secondary structures in vertebrate introns using Kimura's model of compensatory fitness interactions

Previous studies have shown that splicing efficiency, and thus maturation of pre-mRNA, depends on the correct folding of the RNA molecule into a secondary or higher order structure. When disrupted by a mutation, aberrant folding may result in a lower splicing efficiency. However, the structure can be restored by a second, compensatory mutation. Here, we present a logistic regression approach to analyze the evolutionary dynamics of RNA secondary structures. We apply our approach to a set of computationally predicted RNA secondary structures in vertebrate introns. Our results are consistent with the hypothesis of a negative influence of the physical distance between pairing nucleotides on the occurrence of covariations, as predicted by Kimura's model of compensatory evolution. We also confirm the hypothesis that longer local secondary structure elements (helices) can accommodate a larger number of covariations, wobbles, and mismatches. Furthermore, we find that wobbles and mismatches are more frequent in the middle of a helix, whereas covariations occur preferentially at the helix ends. The GC content is a major determinant of this pattern.

Evolutionary origins of invasive populations

What factors shape the evolution of invasive populations? Recent theoretical and empirical studies suggest that an evolutionary history of disturbance might be an important factor. This perspective presents hypotheses regarding the impact of disturbance on the evolution of invasive populations, based on a synthesis of the existing literature. Disturbance might select for life-history traits that are favorable for colonizing novel habitats, such as rapid population growth and persistence. Theoretical results suggest that disturbance in the form of fluctuating environments might select for organismal flexibility, or alternatively, the evolution of evolvability. Rapidly fluctuating environments might favor organismal flexibility, such as broad tolerance or plasticity. Alternatively, longer fluctuations or environmental stress might lead to the evolution of evolvability by acting on features of the mutation matrix. Once genetic variance is generated via mutations, temporally fluctuating selection across generations might promote the accumulation and maintenance of genetic variation. Deeper insights into how disturbance in native habitats affects evolutionary and physiological responses of populations would give us greater capacity to predict the populations that are most likely to tolerate or adapt to novel environments during habitat invasions. Moreover, we would gain fundamental insights into the evolutionary origins of invasive populations.

I. VH gene transcription creates stabilized secondary structures for coordinated mutagenesis during somatic hypermutation

During the adaptive immune response, antigen challenge triggers a million-fold increase in mutation rates in the variable-region antibody genes. The frequency of mutation is causally and directly linked to transcription, which provides ssDNA and drives supercoiling that stabilizes secondary structures containing unpaired, intrinsically mutable bases. Simulation analysis of transcription in VH5 reveals a dominant 65nt secondary structure in the non-transcribed strand containing six sites of mutable ssDNA that have also been identified independently in human B cell lines and in primary mouse B cells. This dominant structure inter-converts briefly with less stable structures and is formed repeatedly during transcription, due to periodic pauses and backtracking. In effect, this creates a stable yet dynamic "mutability platform" consisting of ever-changing patterns of unpaired bases that are simultaneously exposed and therefore able to coordinate mutagenesis. Such a complex of secondary structures may be the source of ssDNA for enzyme-based diversification, which ultimately results in high affinity antibodies.

II. Correlations between secondary structure stability and mutation frequency during somatic hypermutation

The role of secondary structures and base mutability at different levels of transcription and supercoiling is analyzed in variable region antibody genes VH5, VH94 and VH186.2. The data are consistent with a model of somatic hypermutation in which increasing levels of transcription and secondary structure stability correlate with the initial formation of successive mutable sites. Encoded differences exist in stem length and the number of GC pairs at low versus high levels of transcription in CDRs. These circumstances simplify the complexities of coordinating mutagenesis by confining this process to each mutable site successively, as they form in response to increasing levels of transcription during affinity maturation.

Ectopic expression of AID in a non-B cell line triggers A:T and G:C point mutations in non-replicating episomal vectors

Somatic hypermutation (SHM) of immunoglobulin genes is currently viewed as a two step process initiated by the deamination of deoxycytidine (C) to deoxyuridine (U), catalysed by the activation induced deaminase (AID). Phase 1 mutations arise from DNA replication across the uracil residue or the abasic site, generated by the uracil-DNA glycosylase, yielding transitions or transversions at G:C pairs. Phase 2 mutations result from the recognition of the U∶G mismatch by the Msh2/Msh6 complex (MutS Homologue), followed by the excision of the mismatched nucleotide and the repair, by the low fidelity DNA polymerase η, of the gap generated by the exonuclease I. These mutations are mainly focused at A∶T pairs. Whereas in activated B cells both G:C and A∶T pairs are equally targeted, ectopic expression of AID was shown to trigger only G:C mutations on a stably integrated reporter gene. Here we show that when using non-replicative episomal vectors containing a GFP gene, inactivated by the introduction of stop codons at various positions, a high level of EGFP positive cells was obtained after transient expression in Jurkat cells constitutively expressing AID. We show that mutations at G:C and A∶T pairs are produced. EGFP positive cells are obtained in the absence of vector replication demonstrating that the mutations are dependent only on the mismatch repair (MMR) pathway. This implies that the generation of phase 1 mutations is not a prerequisite for the expression of phase 2 mutations.

Secondary structures as predictors of mutation potential in the lacZ gene of Escherichia coli

Four independent nonsense mutations were engineered into the Escherichia coli chromosomal lacZ gene, and reversion rates back to LacZ(+) phenotypes were determined. The mutation potential of bases within putative DNA secondary structures formed during transcription was predicted by a sliding-window analysis that simulates successive folding of the ssDNA creating these structures. The relative base mutabilities predicted by the mfg computer program correlated with experimentally determined reversion rates in three of the four mutants analysed. The nucleotide changes in revertants at one nonsense codon site consisted of a triple mutation, presumed to occur by a templated repair mechanism. Additionally, the effect of supercoiling on mutation was investigated and, in general, reversion rates increased with higher levels of negative supercoiling. Evidence indicates that predicted secondary structures are in fact formed in vivo and that directed mutation in response to starvation stress is dependent upon the exposure of particular bases, the stability of the structures in which these bases are unpaired and the level of DNA supercoiling within the cell.

The effect of promoter strength, supercoiling and secondary structure on mutation rates in Escherichia coli

Four mutations resulting in opal stop codons were individually engineered into a plasmid-borne chloramphenicol-resistance (cat) gene driven by the lac promoter. These four mutations were located at different sites in secondary structures. The mutations were analysed with the computer program mfg, which predicted their relative reversion frequencies. Reversion frequencies determined experimentally correlated with the mutability of the bases as predicted by mfg. To examine the effect of increased transcription on reversion frequencies, the lac promoter was replaced with the stronger tac promoter, which resulted in 12- to 30-fold increases in reversion rates. The effect of increased and decreased supercoiling was also investigated. The cat mutants had higher reversion rates in a topA mutant strain with increased negative supercoiling compared with wild-type levels, and the cat reversion rates were lower in a topA gyrB mutant strain with decreased negative supercoiling, as predicted.

Mechanisms of genotoxin-induced transcription and hypermutation in p53

It is widely assumed that genotoxin-induced damage (e.g., G-to-T transversions) to the tumor suppressor gene, p53, is a direct cause of cancer. However, genotoxins also induce the stress response, which upregulates p53 transcription and the formation of secondary structures from ssDNA. Since unpaired bases are thermodynamically unstable and intrinsically mutable, increased transcription could be the cause of hypermutation, and thus cancer. Support for this hypothesis has been obtained by analyzing 6662 mutations in all types of cancer compared to lung and colon cancers, using the p53 mutation database. The data suggest that genotoxins have two independent effects: first, they induce p53 transcription, which increases the number of mutable bases that determine the incidence of cancer. Second, genotoxins may alter the fate, or ultimate mutation of a mutable base, for example, by causing more of the available mutable Gs to mutate to T, leaving fewer to mutate to A. Such effects on the fate of mutable bases have no impact on the incidence of cancer, as both types of mutations lead to cancer.

Secondary structure formations of conotoxin genes: A possible role in mediating variability

Small venomous peptides called conotoxins produced by the predatory marine snail (genus Conus) present an interesting case for mutational studies. They have a high degree of amino acid variability among them yet they possess highly conserved structural elements that are defined by cysteine residues forming disulfide bridges along the length of the mature peptide. It has been observed that codons specifying these cysteines are also highly conserved. It is unknown how such codon conservation is maintained within the mature conotoxin gene since this entire region undergoes an accelerated rate of mutation. There is evidence suggesting that nucleic acids wield some influence in mechanisms that dictate the region and frequency where mutations occur in DNA. Nucleic acids exert this effect primarily through secondary structures that bring about local peaks and troughs in the energy relief of these transient formations. Secondary structure predictions of several conotoxin genes were analyzed to see if there was any correspondence between the highly variable regions of the conotoxin. Regions of the DNA encompassing the conserved Cys codons (and several other conserved amino acid codons) have been found to correspond to predicted secondary structures of higher stabilities. In stark contrast the regions of the conotoxin that have a higher degree of variation correlate to regions of lower stability. This striking co-relation allows for a simple model of inaccessibility of a mutator to these highly conserved regions of the conotoxin gene allowing them a relative degree of resistance towards change.

Selection acts on DNA secondary structures to decrease transcriptional mutagenesis

Single-stranded DNA is more subject to mutation than double stranded. During transcription, DNA is transiently single stranded and therefore subject to higher mutagenesis. However, if local intra-strand secondary structures are formed, some bases will be paired and therefore less sensitive to mutation than unpaired bases. Using complete genome sequences of Escherichia coli, we show that local intra-strand secondary structures can, as a consequence, be used to define an index of transcription-driven mutability. At gene level, we show that natural selection has favoured a reduced transcription-driven mutagenesis via the higher than expected frequency of occurrence of intra-strand secondary structures. Such selection is stronger in highly expressed genes and suggests a sequence-dependent way to control mutation rates and a novel form of selection affecting the evolution of synonymous mutations.

Genome sequence evolution results from the interplay between mutagenesis and natural selection. Mutations occur as the result of biochemical or physical alteration of DNA and/or from the errors made by polymerases while replicating DNA. As many mutations tend to be detrimental to the organism's fitness, natural selection favours a decrease in mutation rate. Hence, many mechanisms have evolved to control mutation rate. The mechanisms described to date have relied on (i) the existence of enzymes repairing the damaged DNA or correcting mismatched bases, which are mechanisms having an effect on whole genome mutation rate, and (ii) the avoidance in the sequence of repetition that could be misread by the polymerases, which is a sequence-dependent local control of mutation rate. In the present paper, the authors suggest that another sequence-dependent control of mutation exists and shapes the overall evolution of the genome. Using a comparative analysis of Escherichia coli genomes, they show that local secondary structures that are formed during the transcription of genes into RNA can modulate the base-to-base mutation rate. Moreover, the authors show that natural selection seems to have favoured the occurrence of such structures to minimise mutability, especially in the most expressed genes. This paper proposes a new way in which gene sequences can be constrained by natural selection.

Transcription promotes guanine to thymine mutations in the non-transcribed strand of an Escherichia coli gene

Transcription of DNA opens the chromatin, causes topological changes in DNA and transiently exposes the two strands to different biochemical environments. Consequently, it has long been argued that transcription may promote damage to DNA and there are data in Escherichia coli and yeast supporting a correlation between high transcription and mutations. We examined the transcription-dependence of the reversion of a nonsense codon (TGA) in E. coli and found that there was a strong dependence of mutations on transcription in strains defective in the repair of 8-oxoguanine in DNA. Under conditions of high transcription there was a three to five-fold increase in mutations that changed TGA in the non-transcribed strand to a sense codon. Furthermore, in both mutY and mutM mutY backgrounds the mutations were overwhelmingly G:C to T:A. In contrast, when the TGA was in the transcribed strand in relation with the inducible promoter, high transcription decreased the rate of reversion. Similar results were obtained in a strain defective in the transcription-repair coupling factor, Mfd, suggesting that transcription dependent increase in base substitutions does not require transcription-dependent DNA repair. However, Mfd does modulate the magnitude of the mutagenic effect of transcription. These data are consistent with a model in which the non-transcribed strand is more susceptible to oxidative damage during transcription than the transcribed strand. These results suggest that the magnitudes of individual base substitutions and their relative numbers in other studies of mutational spectra may also be affected by transcription.

Hypothesis: biological role for J-C intronic matrix attachment regions in the molecular mechanism of antigen-driven somatic hypermutation

A major function of J-C intronic matrix attachment regions (MAR) during immune diversification via somatic hypermutation (SHM) at immunoglobulin loci may be to manipulate the topology of DNA within the upstream target domain. The suggestion that SHM induction requires MAR-induced torsional strain, in conjunction with DNA remodelling at the J-C intron, completes the definition of a cogent paradigm within which all extant molecular data on the issue may be interpreted. Moreover, the suggestion that a mutagenic mechanism relieves MAR-generated superhelicity could provide an indication as to the evolutionary basis of SHM.

Mechanisms by which transcription can regulate somatic hypermutation

Mechanisms for somatic hypermutation (SHM) have proven elusive. An actively transcribed substrate was analyzed to elucidate the role of stem-loop structures (SLSs) in SHM. Analysis with a new computer algorithm indicates that the location and mutability of a base are regulated by: (a) the extent to which it is unpaired, (b) the degree to which it is exposed by stabilization of SLSs containing and flanking it, and (c) the level of transcription that drives supercoiling, which creates and stabilizes SLSs containing unpaired bases vulnerable to mutation. New mechanisms are described by which transcription can differentially stabilize SLSs harboring targeted bases and establish specific base exposure patterns. Assuming that transcription levels correlate with the magnitude of superhelicity induced and the lengths of ssDNA forming SLSs, this analysis accounts for the location of all mutable bases during SHM.

Increased transcription rates correlate with increased reversion rates in leuB and argH Escherichia coli auxotrophs

Escherichia coli auxotrophs of leuB and argH were examined to determine if higher rates of transcription in derepressed genes were correlated with increased reversion rates. Rates of leuB and argH mRNA synthesis were determined using half-lives and concentrations, during exponential growth and at several time points during 30 min of amino acid starvation. Changes in mRNA concentration were primarily due to increased mRNA synthesis and not to increased stability. Four strains of E. coli amino acid auxotrophs, isogenic except for relA and argR, were examined. In both the leuB and argH genes, rates of transcription and mutation were compared. In general, strains able to activate transcription with guanosine tetraphosphate (ppGpp) had higher rates of mRNA synthesis and mutation than those lacking ppGpp (relA2 mutants). argR knockout strains were constructed in relA(+) and relA mutant strains, and rates of both argH reversion and mRNA synthesis were significantly higher in the argR knockouts than in the regulated strains. A statistically significant linear correlation between increased rates of transcription and mutation was found for data from both genes. In general, changes in mRNA half-lives were less than threefold, whereas changes in rates of mRNA synthesis were often two orders of magnitude. The results suggest that specific starvation conditions target the biosynthetic genes for derepression and increased rates of transcription and mutation.

Stress-directed adaptive mutations and evolution

Comparative biochemistry demonstrates that the metabolites, complex biochemical networks, enzymes and regulatory mechanisms essential to all living cells are conserved in amazing detail throughout evolution. Thus, in order to evolve, an organism must overcome new adverse conditions without creating different but equally dangerous alterations in its ongoing successful metabolic relationship with its environment. Evidence suggests that stable long-term acquisitive evolution results from minor increases in mutation rates of genes related to a particular stress, with minimal disturbance to the balanced and resilient metabolism critical for responding to an unpredictable environment. Microorganisms have evolved specific biochemical feedback mechanisms that direct mutations to genes derepressed by starvation or other stressors in their environment. Transcription of the activated genes creates localized supercoiling and DNA secondary structures with unpaired bases vulnerable to mutation. The resulting mutants provide appropriate variants for selection by the stress involved, thus accelerating evolution with minimal random damage to the genome. This model has successfully predicted mutation frequencies in genes of E. coli and humans. Stressed cells observed in the laboratory over hundreds of generations accumulate mutations that also arise by this mechanism. When this occurs in repair-deficient mutator strains with high rates of random mutation, the specific stress-directed mutations are also enhanced.