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.

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.

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.

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

This work provides evidence that, during transcription, the mutability (propensity to mutate) of a base in a DNA secondary structure depends both on the stability of the structure and on the extent to which the base is unpaired. Zuker's DNA folding computer program reveals the most probable stem-loop structures (SLSs) and negative energies of folding (-DeltaG) for any given nucleotide sequence. We developed an interfacing program that calculates (i) the percentage of folds in which each base is unpaired during transcription; and (ii) the mutability index (MI) for each base, expressed as an absolute value and defined as -follows: MI = (% total folds in which the base is unpaired) x (highest -DeltaG of all folds in which it is unpaired). Thus, MIs predict the relative mutation or reversion frequencies of unpaired bases in SLSs. MIs for 16 mutable bases in auxotrophs, selected during starvation in derepressed genes, are compared with 70 background mutations in lacI and ebgR that were not derepressed during mutant selection. All the results are consistent with the location of known mutable bases in SLSs. Specific conclusions are: (i) Of 16 mutable bases in transcribing genes, 87% have higher MIs than the average base of the sequence analysed, compared with 50% for the 70 background mutations. (ii) In 15 of the mutable bases of transcribing genes, the correlation between MIs and relative mutation frequencies determined experimentally is good. There is no correlation for 35 mutable bases in the lacI gene. (iii) In derepressed auxotrophs, 100% of the codons containing the mutable bases are within one codon's length of a stem, compared with 53% for the background mutable bases in lacI. (iv) The data suggest that environmental stressors may cause as well as select mutations in derepressed genes. The implications of these results for evolution are discussed.

Hypermutable bases in the p53 cancer gene are at vulnerable positions in DNA secondary structures.

A DNA folding analysis indicates that the most hypermutable bases in exons 5, 7, and 8 of the p53 tumor suppressor gene are located immediately next to stems in stable DNA stem-loop structures. On the basis of the highest negative energy (-DeltaG) value of the structures containing each mutable bases and on the extent to which each base is unpaired during transcription, their relative mutabilities are calculated using a new computer algorithm. These predicted mutation frequencies correlate well with those observed in 14,000 human cancers (R(2) = 0.76), whereas there is no such correlation (R(2) = 0.0005) for nearby control bases. The correlation of hypermutable base frequencies with -DeltaG values is poor (R(2) = 0.19), indicating that the extent to which a base is unpaired during transcription is a significant contribution to predicting mutation frequencies.