Mutations induced by bacteriophage T7 RNA polymerase and their effects on the composition of the T7 genome

Hypermutation of specific genomic loci of Pseudomonas putida for continuous evolution of target genes

The ability of T7 RNA polymerase (RNAPT7 ) fusions to cytosine deaminases (CdA) for entering C➔T changes in any DNA segment downstream of a T7 promoter was exploited for hyperdiversification of defined genomic portions of Pseudomonas putida KT2440. To this end, test strains were constructed in which the chromosomally encoded pyrF gene (the prokaryotic homologue of yeast URA3) was flanked by T7 transcription initiation and termination signals and also carried plasmids expressing constitutively either high-activity (lamprey's) or low-activity (rat's) CdA-RNAPT7 fusions. The DNA segment-specific mutagenic action of these fusions was then tested in strains lacking or not uracil-DNA glycosylase (UDG), that is ∆ung/ung+ variants. The resulting diversification was measured by counting single nucleotide changes in clones resistant to 5-fluoroorotic acid (5FOA), which otherwise is transformed by wild-type PyrF into a toxic compound. Although the absence of UDG dramatically increased mutagenic rates with both CdA-RNAPT7 fusions, the most active variant - pmCDA1 - caused extensive appearance of 5FOA-resistant colonies in the wild-type strain not limited to C➔T but including also a range of other changes. Furthermore, the presence/absence of UDG activity swapped cytosine deamination preference between DNA strands. These qualities provided the basis of a robust system for continuous evolution of preset genomic portions of P. putida and beyond.

Heat adaptation of phage T7 under an extended genetic code

While bacteriophages have previously been used as a model system to understand thermal adaptation, most adapted genomes observed to date contain very few modifications and cover a limited temperature range. Here, we set out to investigate genome adaptation to thermal stress by adapting six populations of T7 bacteriophage virions to increasingly stringent heat challenges. Further, we provided three of the phage populations’ access to a new genetic code in which Amber codons could be read as selenocysteine, potentially allowing the formation of more stable selenide-containing bonds. Phage virions responded to the thermal challenges with a greater than 10°C increase in heat tolerance and fixed highly reproducible patterns of non-synonymous substitutions and genome deletions. Most fixed mutations mapped to either the tail complex or to the three internal virion proteins that form a pore across the E. coli cell membrane during DNA injection. However, few global changes in Amber codon usage were observed, with only one natural Amber codon being lost. These results reinforce a model in which adaptation to thermal stress proceeds via the cumulative fixation of a small set of highly adaptive substitutions and that adaptation to new genetic codes proceeds only slowly, even with the possibility of potential phenotypic advantages.

RNA polymerase between lesion bypass and DNA repair

DNA damage leads to heritable changes in the genome via DNA replication. However, as the DNA helix is the site of numerous other transactions, notably transcription, DNA damage can have diverse repercussions on cellular physiology. In particular, DNA lesions have distinct effects on the passage of transcribing RNA polymerases, from easy bypass to almost complete block of transcription elongation. The fate of the RNA polymerase positioned at a lesion is largely determined by whether the lesion is structurally subtle and can be accommodated and eventually bypassed, or bulky, structurally distorting and requiring remodeling/complete dissociation of the transcription elongation complex, excision, and repair. Here we review cellular responses to DNA damage that involve RNA polymerases with a focus on bacterial transcription-coupled nucleotide excision repair and lesion bypass via transcriptional mutagenesis. Emphasis is placed on the explosion of new structural information on RNA polymerases and relevant DNA repair factors and the mechanistic models derived from it.

Topoisomerase 1 provokes the formation of short deletions in repeated sequences upon high transcription in Saccharomyces cerevisiae

High transcription is associated with genetic instability, notably increased spontaneous mutation rates, which is a phenomenon termed Transcription-Associated-Mutagenesis (TAM). In this study, we investigated TAM using the chromosomal CAN1 gene under the transcriptional control of two strong and inducible promoters (pGAL1 and pTET) in Saccharomyces cerevisiae. Both pTET- and pGAL1-driven high transcription at the CAN1 gene result in enhanced spontaneous mutation rates. Comparison of both promoters reveals differences in the type of mutagenesis, except for short (-2 and -3 nt) deletions, which depend only on the level of transcription. This mutation type, characteristic of TAM, is sequence dependent, occurring prefentially at di- and trinucleotides repeats, notably at two mutational hotspots encompassing the same 5'-ACATAT-3' sequence. To explore the mechanisms underlying the formation of short deletions in the course of TAM, we have determined Can(R) mutation spectra in yeast mutants affected in DNA metabolism. We identified topoisomerase 1-deficient strains (top1Δ) that specifically abolish the formation of short deletions under high transcription. The rate of the formation of (-2/-3nt) deletions is also reduced in the absence of RAD1 and MUS81 genes, involved in the repair of Top1p-DNA covalent complex. Furthermore ChIP analysis reveals an enrichment of trapped Top1p in the CAN1 ORF under high transcription. We propose a model, in which the repair of trapped Top1p-DNA complexes provokes the formation of short deletion in S. cerevisiae. This study reveals unavoidable conflicts between Top1p and the transcriptional machinery and their potential impact on genome stability.

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.

Transcription-induced mutational strand bias and its effect on substitution rates in human genes

If substitution rates are not the same on the two complementary DNA strands, a substitution is considered strand asymmetric. Such substitutional strand asymmetries are determined here for the three most frequent types of substitution on the human genome (C --> T, A --> G, and G --> T). Substitution rate differences between both strands are estimated for 4,590 human genes by aligning all repeats occurring within the introns with their ancestral consensus sequences. For 1,630 of these genes, both coding strand and noncoding strand rates could be compared with rates in gene-flanking regions. All three rates considered are found to be on average higher on the coding strand and lower on the transcribed strand in comparison to their values in the gene-flanking regions. This finding points to the simultaneous action of rate-increasing effects on the coding strand--such as increased adenine and cytosine deamination--and transcription-coupled repair as a rate-reducing effect on the transcribed strand. The common behavior of the three rates leads to strong correlations of the rate asymmetries: Whenever one rate is strand biased, the other two rates are likely to show the same bias. Furthermore, we determine all three rate asymmetries as a function of time: the A --> G and G --> T rate asymmetries are both found to be constant in time, whereas the C --> T rate asymmetry shows a pronounced time dependence, an observation that explains the difference between our results and those of an earlier work by Green et al. (2003. Transcription-associated mutational asymmetry in mammalian evolution. Nat Genet. 33:514-517.). Finally, we show that in addition to transcription also the replication process biases the substitution rates in genes.

Transcriptional pausing and stalling causes multiple clustered mutations by human activation-induced deaminase

Transcription of the rearranged immunoglobulin gene and expression of the enzyme activation-induced deaminase (AID) are essential for somatic hypermutations of this gene during antibody maturation. While AID acts as a single-strand DNA-cytosine deaminase creating U . G mispairs that lead to mutations, the role played by transcription in this process is less clear. We have used in vitro transcription of the kan gene by the T7 RNA polymerase (RNAP) in the presence of AID and a genetic reversion assay for kanamycin-resistance to investigate the causes of multiple clustered mutations (MCMs) during somatic hypermutations. We find that, depending on transcription conditions, AID can cause single-base substitutions or MCMs. When wild-type RNAP is used for transcription at physiologically relevant concentrations of ribonucleoside triphosphates (NTPs), few MCMs are found. In contrast, slowing the rate of elongation by reducing the NTP concentration or using a mutant RNAP increases several-fold the percent of revertants containing MCMs. Arresting the elongation complexes by a quick removal of NTPs leads to formation of RNA-DNA hybrids (R-loops). Treatment of these structures with AID results in a high percentage of Kan(R) revertants with MCMs. Furthermore, selecting for transcription elongation complexes stalled near the codon that suffers mutations during acquisition of kanamycin-resistance results in an overwhelming majority of revertants with MCMs. These results show that if RNAP II pauses or stalls during transcription of immunoglobulin gene, AID is likely to promote MCMs. As changes in physiological conditions such as occurrence of certain DNA primary or secondary structures or DNA adducts are known to cause transcriptional pausing and stalling in mammalian cells, this process may cause MCMs during somatic hypermutation.

Exploring tree-building methods and distinct molecular data to recover a known asymmetric phage phylogeny

An experimental phylogeny was constructed using bacteriophage T7 and a propagation protocol, in the presence of the mutagen N-methyl-N'-nitro-N'-nitrosoguanidine, based on Hillis et al. [Hillis, D.M., Bull, J.J., White, M.E., Badgett, M.R., Molineux, I.J., 1992. Experimental phylogenetics, generation of a known phylogeny. Science 255, 589-592]. The topology presented in this study has a considerable variation in branch lengths and is less symmetric than the one presented by Hillis et al. [Hillis, D.M., Bull, J.J., White, M.E., Badgett, M.R., Molineux, I.J., 1992. Experimental phylogenetics, generation of a known phylogeny. Science 255, 589-592]. These features are known to present additional difficulties to phylogenetic inference methods. The performance of several phylogenetic methods (conventional and less conventional) was tested using restriction site and nucleotide data. Only methods that encompassed a molecular clock or those based on sequence signatures recovered the true phylogeny. Nevertheless a likelihood ratio test rejected the hypothesis of the existence of a molecular clock when the whole sequence data set was considered. This fact or the particular substitution pattern (mainly G-->A and C-->T) may be related to the unexpected performance of distance methods based on sequence signatures. To test if the results could have been predicted by simulation studies we estimated the evolution parameters from the real phylogeny and used them to simulate evolution along the same tree (parametric bootstrap). We found that simulation could predict most but not all of the problems encountered by phylogenetic inference methods in the real phylogeny. Short interior branches may be more prone to error than predicted by theoretical studies.

Bifractality of human DNA strand-asymmetry profiles results from transcription

We use the wavelet transform modulus maxima method to investigate the multifractal properties of strand-asymmetry DNA walk profiles in the human genome. This study reveals the bifractal nature of these profiles, which involve two competing scale-invariant (up to repeat-masked distances less, or similar 40 kbp) components characterized by Hölder exponents h{1}=0.78 and h{2}=1, respectively. The former corresponds to the long-range-correlated homogeneous fluctuations previously observed in DNA walks generated with structural codings. The latter is associated with the presence of jumps in the original strand-asymmetry noisy signal S. We show that a majority of upward (downward) jumps co-locate with gene transcription start (end) sites. Here 7228 human gene transcription start sites from the refGene database are found within 2 kbp from an upward jump of amplitude DeltaS > or = 0.1 which suggests that about 36% of annotated human genes present significant transcription-induced strand asymmetry and very likely high expression rate.

Sequence-dependent enhancement of hydrolytic deamination of cytosines in DNA by the restriction enzyme PspGI

Hydrolytic deamination of cytosines in DNA creates uracil and, if unrepaired, these lesions result in C to T mutations. We have suggested previously that a possible way in which cells may prevent or reduce this chemical reaction is through the binding of proteins to DNA. We use a genetic reversion assay to show that a restriction enzyme, PspGI, protects cytosines within its cognate site, 5'-CCWGG (W is A or T), against deamination under conditions where no DNA cleavage can occur. It decreases the rate of cytosine deamination to uracil by 7-fold. However, the same protein dramatically increases the rate of deaminations within the site 5'-CCSGG (S is C or G) by approximately 15-fold. Furthermore, a similar increase in cytosine deaminations is also seen with a catalytically inactive mutant of the enzyme showing that endonucleolytic ability of the protein is dispensable for its mutagenic action. The sequences of the mutants generated in the presence of PspGI show that only one of the cytosines in CCSGG is predominantly converted to thymine. Our results are consistent with PspGI 'sensitizing' the cytosine in the central base pair in CCSGG for deamination. Remarkably, PspGI sensitizes this base for damage despite its inability to form stable complexes at CCSGG sites. These results can be explained if the enzyme has a transient interaction with this sequence during which it flips the central cytosine out of the helix. This prediction was validated by modeling the structure of PspGI-DNA complex based on the structure of the related enzyme Ecl18kI which is known to cause base-flipping.

Mutation-biased adaptation in a protein NK model

Evolutionary trends responsible for systematic differences in genome and proteome composition have been attributed to GC:AT mutation bias in the context of neutral evolution or to selection acting on genome composition. A possibility that has been ignored, presumably because it is part of neither the Modern Synthesis nor the Neutral Theory, is that mutation may impose a directional bias on adaptation. This possibility is explored here with simulations of the effect of a GC:AT bias on amino acid composition during adaptive walks on an abstract protein fitness landscape called an "NK" model. The results indicate that adaptation does not preclude mutation-biased evolution. In the complete absence of neutral evolution, a modest GC:AT bias of realistic magnitude can displace the trajectory of adaptation in a mutationally favored direction, to such a degree that amino acid composition is biased substantially and persistently. Thus, mutational explanations for evolved patterns need not presuppose neutral evolution.

Features of Arabidopsis genes and genome discovered using full-length cDNAs

Arabidopsis is currently the reference genome for higher plants. A new, more detailed statistical analysis of Arabidopsis gene structure is presented including intron and exon lengths, intergenic distances, features of promoters, and variant 5'-ends of mRNAs transcribed from the same transcription unit. We also provide a statistical characterization of Arabidopsis transcripts in terms of their size, UTR lengths, 3'-end cleavage sites, splicing variants, and coding potential. These analyses were facilitated by scrutiny of our collection of sequenced full-length cDNAs and much larger collection of 5'-ESTs, together with another set of full-length cDNAs from Salk/Stanford/Plant Gene Expression Center/RIKEN. Examples of alternative splicing are observed for transcripts from 7% of the genes and many of these genes display multiple spliced isoforms. Most splicing variants lie in non-coding regions of the transcripts. Non-canonical splice sites constitute less than 1% of all splice sites. Genes with fewer than four introns display reduced average mRNA levels. Putative alternative transcription start sites were observed in 30% of highly expressed genes and in more than 50% of the genes with low expression. Transcription start sites correlate remarkably well with a CG skew peak in the DNA sequences. The intergenic distances vary considerably, those where genes are transcribed towards one another being significantly shorter. New transcripts, missing in the current TIGR genome annotation and ESTs that are non-coding, including those antisense to known genes, are derived and cataloged in the Supplementary Material. They identify 148 new loci in the Arabidopsis genome. The conclusions drawn provide a better understanding of the Arabidopsis genome and how the gene transcripts are processed. The results also allow better predictions to be made for, as yet, poorly defined genes and provide a reference for comparisons with other plant genomes whose complete sequences are currently being determined. Some comparisons with rice are included in this paper.

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.

Characterization of the Type III restriction endonuclease PstII from Providencia stuartii

A new Type III restriction endonuclease designated PstII has been purified from Providencia stuartii. PstII recognizes the hexanucleotide sequence 5'-CTGATG(N)(25-26/27-28)-3'. Endonuclease activity requires a substrate with two copies of the recognition site in head-to-head repeat and is dependent on a low level of ATP hydrolysis ( approximately 40 ATP/site/min). Cleavage occurs at just one of the two sites and results in a staggered cut 25-26 nt downstream of the top strand sequence to generate a two base 5'-protruding end. Methylation of the site occurs on one strand only at the first adenine of 5'-CATCAG-3'. Therefore, PstII has characteristic Type III restriction enzyme activity as exemplified by EcoPI or EcoP15I. Moreover, sequence asymmetry of the PstII recognition site in the T7 genome acts as an historical imprint of Type III restriction activity in vivo. In contrast to other Type I and III enzymes, PstII has a more relaxed nucleotide specificity and can cut DNA with GTP and CTP (but not UTP). We also demonstrate that PstII and EcoP15I cannot interact and cleave a DNA substrate suggesting that Type III enzymes must make specific protein-protein contacts to activate endonuclease activity.

Replication-associated strand asymmetries in mammalian genomes: toward detection of replication origins

In the course of evolution, mutations do not affect both strands of genomic DNA equally. This imbalance mainly results from asymmetric DNA mutation and repair processes associated with replication and transcription. In prokaryotes, prevalence of G over C and T over A is frequently observed in the leading strand. The sign of the resulting TA and GC skews changes abruptly when crossing replication-origin and termination sites, producing characteristic step-like transitions. In mammals, transcription-coupled skews have been detected, but so far, no bias has been associated with replication. Here, analysis of intergenic and transcribed regions flanking experimentally identified human replication origins and the corresponding mouse and dog homologous regions demonstrates the existence of compositional strand asymmetries associated with replication. Multiscale analysis of human genome skew profiles reveals numerous transitions that allow us to identify a set of 1,000 putative replication initiation zones. Around these putative origins, the skew profile displays a characteristic jagged pattern also observed in mouse and dog genomes. We therefore propose that in mammalian cells, replication termination sites are randomly distributed between adjacent origins. Taken together, these analyses constitute a step toward genome-wide studies of replication mechanisms.

Differential selection and mutation between dsDNA and ssDNA phages shape the evolution of their genomic AT percentage

Background

Bacterial genomes differ dramatically in AT%. We have developed a model to show that the genomic AT% in rapidly replicating bacterial species can be used as an index of the availability of nucleotides A and T for DNA replication in cellular medium. This index is then used to (1) study the evolution and adaptation of the bacteriophage genomic AT% in response to the differential nucleotide availability of the host and (2) test the prediction that double-stranded DNA (dsDNA) phage should exhibit better adaptation than single-stranded DNA (ssDNA) phage because the rate of spontaneous deamination, which leads to C→T or C→U mutations depending on whether C is methylated or not, is about 100-fold greater in ssDNA than in dsDNA.

Results

We retrieved 79 dsDNA phage and 27 ssDNA phage genomes together with their host genomic sequences. The dsDNA phages have their genomic AT% better adapted to the host genomic AT% than ssDNA phage. The poorer adaptation of the ssDNA phage can be partially accounted for by the C→T(U) mutations mediated by the spontaneous deamination. For ssDNA phage, the genomic A% is more strongly correlated with their host genomic AT% than the genomic T%.

Conclusion

A significant fraction of variation in the genomic AT% in the dsDNA phage, and that in the genomic A% and T% of the ssDNA phage, can be explained by the difference in selection and mutation between them.

Agrobacterium T-DNA integration in Arabidopsis is correlated with DNA sequence compositions that occur frequently in gene promoter regions

Mobile insertion elements such as transposons and T-DNA generate useful genetic variation and are important tools for functional genomics studies in plants and animals. The spectrum of mutations obtained in different systems can be highly influenced by target site preferences inherent in the mechanism of DNA integration. We investigated the target site preferences of Agrobacterium T-DNA insertions in the chromosomes of the model plant Arabidopsis thaliana. The relative frequencies of insertions in genic and intergenic regions of the genome were calculated and DNA composition features associated with the insertion site flanking sequences were identified. Insertion frequencies across the genome indicate that T-strand integration is suppressed near centromeres and rDNA loci, progressively increases towards telomeres, and is highly correlated with gene density. At the gene level, T-DNA integration events show a statistically significant preference for insertion in the 5' and 3' flanking regions of protein coding sequences as well as the promoter region of RNA polymerase I transcribed rRNA gene repeats. The increased insertion frequencies in 5' upstream regions compared to coding sequences are positively correlated with gene expression activity and DNA sequence composition. Analysis of the relationship between DNA sequence composition and gene activity further demonstrates that DNA sequences with high CG-skew ratios are consistently correlated with T-DNA insertion site preference and high gene expression. The results demonstrate genomic and gene-specific preferences for T-strand integration and suggest that DNA sequences with a pronounced transition in CG- and AT-skew ratios are preferred targets for T-DNA integration.

The Myotonic Dystrophy Type 1 Triplet Repeat Sequence Induces Gross Deletions and Inversions

The capacity of (CTG.CAG)n and (GAA.TTC)n repeat tracts in plasmids to induce mutations in DNA flanking regions was evaluated in Escherichia coli. Long repeats of these sequences are involved in the etiology of myotonic dystrophy type 1 and Friedreich's ataxia, respectively. Long (CTG.CAG)n (where n = 98 and 175) caused the deletion of most, or all, of the repeats and the flanking GFP gene. Deletions of 0.6-1.8 kbp were found as well as inversions. Shorter repeat tracts (where n = 0 or 17) were essentially inert, as observed for the (GAA.TTC)176-containing plasmid. The orientation of the triplet repeat sequence (TRS) relative to the unidirectional origin of replication had a pronounced effect, signaling the participation of replication and/or repair systems. Also, when the TRS was transcribed, the level of deletions was greatly elevated. Under certain conditions, 30-50% of the products contained gross deletions. DNA sequence analyses of the breakpoint junctions in 47 deletions revealed the presence of 1-8-bp direct or inverted homologies in all cases. Also, the presence of non-B folded conformations (i.e. slipped structures, cruciforms, or triplexes) at or near the breakpoints was predicted in all cases. This genetic behavior, which was previously unrecognized for a TRS, may provide the basis for a new type of instability of the myotonic dystrophy protein kinase (DMPK) gene in patients with a full mutation.

Translocation of Escherichia coli RNA polymerase against a protein roadblock in vivo highlights a passive sliding mechanism for transcript elongation

Current models for transcription elongation infer that RNA polymerase (RNAP) moves along the template by a passive sliding mechanism that takes advantage of random lateral oscillations in which single basepair sliding movements interconvert the elongation complex between pre- and post-translocated states. Such passive translocational equilibrium was tested in vivo by a systematic change in the templated NTP that is to be incorporated by RNAP, which is temporarily roadblocked by the lac repressor. Our results show that, under these conditions that hinder the forward movement of the polymerase, the elongation complex is able to extend its RNA chain one nucleotide further when the incoming NTP is a kinetically favoured substrate (i.e. low K(m)). The addition of an extra nucleotide destabilizes the repressor-operator roadblock leading to an increase in transcriptional readthrough. Similar results are obtained when the incoming NTPs are less kinetically favoured substrates (i.e. high K(m)s) by specifically increasing their intracellular concentrations. Altogether, these in vivo data are consistent with a passive sliding model in which RNAP forward translocation is favoured by NTP binding. They also suggest that fluctuations in the intracellular NTP pools may play a key role in gene regulation at the transcript elongation level.

Identification of a distinctive mutation spectrum associated with high levels of transcription in yeast

High levels of transcription are associated with increased mutation rates in Saccharomyces cerevisiae, a phenomenon termed transcription-associated mutation (TAM). To obtain insight into the mechanism of TAM, we obtained LYS2 forward mutation spectra under low- versus high-transcription conditions in which LYS2 was expressed from either the low-level pLYS2 promoter or the strong pGAL1-10 promoter, respectively. Because of the large size of the LYS2 locus, forward mutations first were mapped to specific LYS2 subregions, and then those mutations that occurred within a defined 736-bp target region were sequenced. In the low-transcription strain base substitutions comprised the majority (64%) of mutations, whereas short insertion-deletion mutations predominated (56%) in the high-transcription strain. Most notably, deletions of 2 nucleotides (nt) comprised 21% of the mutations in the high-transcription strain, and these events occurred predominantly at 5'-(G/C)AAA-3' sites. No -2 events were present in the low-transcription spectrum, thus identifying 2-nt deletions as a unique mutational signature for TAM.

Mutational patterns correlate with genome organization in SARS and other coronaviruses

Focused efforts by several international laboratories have resulted in the sequencing of the genome of the causative agent of severe acute respiratory syndrome (SARS), novel coronavirus SARS-CoV, in record time. Using cumulative skew diagrams, I found that mutational patterns in the SARS-CoV genome were strikingly different from other coronaviruses in terms of mutation rates, although they were in general agreement with the model of the coronavirus lifecycle. These findings might be relevant for the development of sequence-based diagnostics and the design of agents to treat SARS.

Recombinogenic Effects of DNA-Damaging Agents Are Synergistically Increased by Transcription inSaccharomyces cerevisiae: New Insights Into Transcription-Associated Recombination

Homologous recombination of a particular DNA sequence is strongly stimulated by transcription, a phenomenon observed from bacteria to mammals, which we refer to as transcription-associated recombination (TAR). TAR might be an accidental feature of DNA chemistry with important consequences for genetic stability. However, it is also essential for developmentally regulated processes such as class switching of immunoglobulin genes. Consequently, it is likely that TAR embraces more than one mechanism. In this study we tested the possibility that transcription induces recombination by making DNA more susceptible to recombinogenic DNA damage. Using different plasmid-chromosome and direct-repeat recombination constructs in which transcription is driven from either the PGAL1- or the Ptet-regulated promoters, we have

shown that either 4-nitroquinoline-N-oxide (4-NQO) or methyl methanesulfonate (MMS) produces a synergistic increase of recombination when combined with transcription. 4-NQO and MMS stimulated recombination of a transcriptionally active DNA sequence up to 12,800- and 130-fold above the spontaneous levels observed in the absence of transcription, whereas 4-NQO and MMS alone increased recombination 193- and 4.5-fold, respectively. Our results provide evidence that TAR is due, at least in part, to the ability of transcription to enhance the accessibility of DNA to exogenous chemicals and internal metabolites responsible for recombinogenic lesions. We discuss possible parallelisms between the mechanisms of induction of recombination and mutation by transcription.

Mismatch repair in methylated DNA. Structure and activity of the mismatch-specific thymine glycosylase domain of methyl-CpG-binding protein MBD4

MBD4 is a member of the methyl-CpG-binding protein family. It contains two DNA binding domains, an amino-proximal methyl-CpG binding domain (MBD) and a C-terminal mismatch-specific glycosylase domain. Limited in vitro proteolysis of mouse MBD4 yields two stable fragments: a 139-residue fragment including the MBD, and the other 155-residue fragment including the glycosylase domain. Here we show that the latter fragment is active as a glycosylase on a DNA duplex containing a G:T mismatch within a CpG sequence context. The crystal structure confirmed the C-terminal domain is a member of the helix-hairpin-helix DNA glycosylase superfamily. The MBD4 active site is situated in a cleft that likely orients and binds DNA. Modeling studies suggest the mismatched target nucleotide will be flipped out into the active site where candidate residues for catalysis and substrate specificity are present.

Transcription-dependent increase in multiple classes of base substitution mutations in Escherichia coli

We showed previously that transcription in Escherichia coli promotes C. G-to-T. A transitions due to increased deamination of cytosines to uracils in the nontranscribed but not the transcribed strand (A. Beletskii and A. S. Bhagwat, Proc. Natl. Acad. Sci. USA 93:13919-13924, 1996). To study mutations other than that of C to T, we developed a new genetic assay that selects only base substitution mutations and additionally excludes C. G to T. A transitions. This novel genetic reversion system is based on mutations in a termination codon and involves positive selection for resistance to bleomycin or kanamycin. Using this genetic system, we show here that transcription from a strong promoter increases the level of non-C-to-T as well as C-to-T mutations. We find that high-level transcription increases the level of non-C-to-T mutations in DNA repair-proficient cells in three different sequence contexts in two genes and that the rate of mutation is higher by a factor of 2 to 4 under these conditions. These increases are not caused by a growth advantage for the revertants and are restricted to genes that are induced for transcription. In particular, high levels of transcription do not create a general mutator phenotype in E. coli. Sequence analysis of the revertants revealed that the frequency of several different base substitutions increased upon transcription of the bleomycin resistance gene and that G. C-to-T. A transversions dominated the spectrum in cells transcribing the gene. These results suggest that high levels of transcription promote many different spontaneous base substitutions in E. coli.

Lack of dependance of transcription-induced cytosine deaminations on protein synthesis

Transcription-induced mutations (TIM) is a phenomenon in Escherichia coli in which transcription promotes C to T and other mutations in a strand-specific manner. Because the processes of transcription and translation are coupled in prokaryotes and some models regarding creating a hypermutagenic state in E. coli require new protein synthesis, we tested the possibility that TIM was dependent on efficient synthesis of proteins. We used puromycin to reversibly inhibit protein synthesis and found that it had little effect on mRNA synthesis, plasmid copy-number or TIM. Our results show that TIM is not dependent on efficient translation of mRNA and this helps eliminate certain models concerning the mechanism underlying TIM.

Transcription and double-strand breaks induce similar mitotic recombination events in Saccharomyces cerevisiae

We have made a comparative analysis of double-strand-break (DSB)-induced recombination and spontaneous recombination under low- and high-transcription conditions in yeast. We constructed two different recombination substrates, one for the analysis of intermolecular gene conversions and the other for intramolecular gene conversions and inversions. Such substrates were based on the same leu2-HOr allele fused to the tet promoter and containing a 21-bp HO site. Gene conversions and inversions were differently affected by rad1, rad51, rad52, and rad59 single and double mutations, consistent with the actual view that such events occur by different recombination mechanisms. However, the effect of each mutation on each type of recombination event was the same, whether associated with transcription or induced by the HO-mediated DSB. Both the highly transcribed DNA and the HO-cut sequence acted as recipients of the gene conversion events. These results are consistent with the hypothesis that transcription promotes initiation of recombination along the DNA sequence being transcribed. The similarity between transcription-associated and DSB-induced recombination suggests that transcription promotes DNA breaks.

The connection between transcription and genomic instability

Transcription is a central aspect of DNA metabolism that takes place on the same substrate as replication, repair and recombination. Not surprisingly, therefore, there is a physical and functional connection between these processes. In recent years, transcription has proven to be a relevant player in the maintenance of genome integrity and in the induction of genetic instability and diversity. The aim of this review is to provide an integrative view on how transcription can control different aspects of genomic integrity, by exploring different mechanisms that might be responsible for transcription-associated mutation (TAM) and transcription-associated recombination (TAR).

Asymmetric directional mutation pressures in bacteria

BACKGROUND

When there are no strand-specific biases in mutation and selection rates (that is, in the substitution rates) between the two strands of DNA, the average nucleotide composition is theoretically expected to be A = T and G = C within each strand. Deviations from these equalities are therefore evidence for an asymmetry in selection and/or mutation between the two strands. By focusing on weakly selected regions that could be oriented with respect to replication in 43 out of 51 completely sequenced bacterial chromosomes, we have been able to detect asymmetric directional mutation pressures.

RESULTS

Most of the 43 chromosomes were found to be relatively enriched in G over C and T over A, and slightly depleted in G+C, in their weakly selected positions (intergenic regions and third codon positions) in the leading strand compared with the lagging strand. Deviations from A = T and G = C were highly correlated between third codon positions and intergenic regions, with a lower degree of deviation in intergenic regions, and were not correlated with overall genomic G+C content.

CONCLUSIONS

During the course of bacterial chromosome evolution, the effects of asymmetric directional mutation pressures are commonly observed in weakly selected positions. The degree of deviation from equality is highly variable among species, and within species is higher in third codon positions than in intergenic regions. The orientation of these effects is almost universal and is compatible in most cases with the hypothesis of an excess of cytosine deamination in the single-stranded state during DNA replication. However, the variation in G+C content between species is influenced by factors other than asymmetric mutation pressure.

Transcription-induced cytosine-to-thymine mutations are not dependent on sequence context of the target cytosine

We showed previously that transcription of a plasmid-borne kan allele increases C-to-T mutations in the nontranscribed strand. Using two new plasmid-borne kan alleles, one cmp allele, and a chromosomal kan allele, we found in this study that transcription-induced mutations are not limited to specific genes, alleles, or locations and are likely to be a general property of transcript elongation in Escherichia coli.

A relationship between gene expression and protein interactions on the proteome scale: analysis of the bacteriophage T7 and the yeast Saccharomyces cerevisiae

The relationship between the similarity of expression patterns for a pair of genes and interaction of the proteins they encode is demonstrated both for the simple genome of the bacteriophage T7 and the considerably more complex genome of the yeast Saccharomyces cerevisiae. Statistical analysis of large-scale gene expression and protein interaction data shows that protein pairs encoded by co-expressed genes interact with each other more frequently than with random proteins. Furthermore, the mean similarity of expression profiles is significantly higher for respective interacting protein pairs than for random ones. Such coupled analysis of gene expression and protein interaction data may allow evaluation of the results of large-scale gene expression and protein interaction screens as demonstrated for several publicly available datasets. The role of this link between expression and interaction in the evolution from monomeric to oligomeric protein structures is also discussed.

Bias in the introduction of variation as an orienting factor in evolution

According to New Synthesis doctrine, the direction of evolution is determined by selection and not by "internal causes" that act by way of propensities of variation. This doctrine rests on the theoretical claim that because mutation rates are small in comparison to selection coefficients, mutation is powerless to overcome opposing selection. Using a simple population-genetic model, this claim is shown to depend on assuming the prior availability of variation, so that mutation may act only as a "pressure" on the frequencies of existing alleles, and not as the evolutionary process that introduces novelty. As shown here, mutational bias in the introduction of novelty can strongly influence the course of evolution, even when mutation rates are small in comparison to selection coefficients. Recognizing this mode of causation provides a distinct mechanistic basis for an "internalist" approach to determining the contribution of mutational and developmental factors to evolutionary phenomena such as homoplasy, parallelism, and directionality.