The effect of the Primary Structure of DNA on Induction of Mutations by Alkylating Agents

Analysis and comparison of mutation spectra is one of the major tasks of molecular biology, since mutation spectra often reveal important properties of various mutagens and proteins involved in the repair/replication systems. Mutability is known to vary significantly along the nucleotide sequence. Mutations are abundant at certain positions (mutation hotspots). In this work, we applied regression analysis based on the basic logic patterns to understand the role of the nucleotide sequence context in mutation induction. The spectra of mutations induced by various alkylating agents were studied. The nucleotide bases at positions –2, –1, +1 and +2 were shown to have the most significant effect in G : C A : T replacements.     

Computational analysis of mutation spectra

Mutation frequencies vary along a nucleotide sequence, and nucleotide positions with an exceptionally high mutation frequency are called hotspots. Mutation hotspots in DNA often reflect intrinsic properties of the mutation process, such as the specificity with which mutagens interact with nucleic acids and the sequence-specificity of DNA repair/replication enzymes. They might also reflect structural and functional features of target protein or RNA sequences in which they occur. The determinants of mutation frequency and specificity are complex and there are many analytical methods for their study. This paper discusses computational approaches to analysing mutation spectra (distribution of mutations along the target genes) that include many detectable (mutable) positions. The following methods are reviewed: mutation hotspot prediction; pairwise and multiple comparisons of mutation spectra; derivation of a consensus sequence; and analysis of correlation between nucleotide sequence features and mutation spectra. Spectra of spontaneous and induced mutations are used for illustration of the complexities and pitfalls of such analyses. In general, the DNA sequence context of mutation hotspots is a fingerprint of interactions between DNA and DNA repair/replication/modification enzymes, and the analysis of hotspot context provides evidence of such interactions.     

Theoretical analysis of mutation hotspots and their DNA sequence context specificity

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of "context-free" characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.     

Regression trees for analysis of mutational spectra in nucleotide sequences.

MOTIVATION: The study and comparison of mutational spectra is an important problem in molecular biology, because these spectra often reveal important features of the action of various mutagens and the functioning of repair/replication enzymes. As is known, mutability varies significantly along nucleotide sequences: mutations often concentrate at certain positions in a sequence, otherwise termed 'hotspots'. RESULTS: Herein, we propose a regression analysis method based on the use of regression trees in order to analyse the influence of nucleotide context on the occurrence of such hotspots. The REGRT program developed has been tested on simulated and real mutational spectra. For the G:C-->T:A mutational spectra induced by Sn1 alkylating agents (nine spectra), the prediction accuracy was 0. 99. AVAILABILITY: The REGRT program is available upon request from V.Berikov.     

Mutation and polymorphism analysis of the human homogentisate 1, 2-dioxygenase gene in alkaptonuria patients.

Alkaptonuria (AKU), a rare hereditary disorder of phenylalanine and tyrosine catabolism, was the first disease to be interpreted as an inborn error of metabolism. AKU patients are deficient for homogentisate 1,2 dioxygenase (HGO); this deficiency causes homogentisic aciduria, ochronosis, and arthritis. We cloned the human HGO gene and characterized two loss-of-function mutations, P230S and V300G, in the HGO gene in AKU patients. Here we report haplotype and mutational analysis of the HGO gene in 29 novel AKU chromosomes. We identified 12 novel mutations: 8 (E42A, W97G, D153G, S189I, I216T, R225H, F227S, and M368V) missense mutations that result in amino acid substitutions at positions conserved in HGO in different species, 1 (F10fs) frameshift mutation, 2 intronic mutations (IVS9-56G-->A, IVS9-17G-->A), and 1 splice-site mutation (IVS5+1G-->T). We also report characterization of five polymorphic sites in HGO and describe the haplotypic associations of alleles at these sites in normal and AKU chromosomes. One of these sites, HGO-3, is a variable dinucleotide repeat; IVS2+35T/A, IVS5+25T/C, and IVS6+46C/A are intronic sites at which single nucleotide substitutions (dimorphisms) have been detected; and c407T/A is a relatively frequent nucleotide substitution in the coding sequence, exon 4, resulting in an amino acid change (H80Q). These data provide insight into the origin and evolution of the various AKU alleles.     

Recognition sites of eukaryotic DNA topoisomerase I: DNA nucleotide sequencing analysis of topo I cleavage sites on SV40 DNA.

Eukaryotic DNA topoisomerase I introduces transient single-stranded breaks on double-stranded DNA and spontaneously breaks down single-stranded DNA. The cleavage sites on both single and double-stranded SV40 DNA have been determined by DNA sequencing. Consistent with other reports, the eukaryotic enzymes, in contrast to prokaryotic type I topoisomerases, links to the 3'-end of the cleaved DNA and generates a free 5'-hydroxyl end on the other half of the broken DNA strand. Both human and calf enzymes cleave SV40 DNA at the identical and specific sites. From 827 nucleotides sequenced, 68 cleavage sites were mapped. The majority of the cleavage sites were present on both double and single-stranded DNA at exactly the same nucleotide positions, suggesting that the DNA sequence is essential for enzyme recognition. By analyzing all the cleavage sequences, certain nucleotides are found to be less favored at the cleavage sites. There is a high probability to exclude G from positions -4, -2, -1 and +1, T from position -3, and A from position -1. These five positions (-4 to +1 oriented in the 5' to 3' direction) around the cleavage sites must interact intimately with topo I and thus are essential for enzyme recognition. One topo I cleavage site which shows atypical cleavage sequence maps in the middle of a palindromic sequence near the origin of SV40 DNA replication. It occurs only on single-stranded SV40 DNA, suggesting that the DNA hairpin can alter the cleavage specificity. The strongest cleavage site maps near the origin of SV40 DNA replication at nucleotide 31-32 and has a pentanucleotide sequence of 5'-TGACT-3'.     

Promoter Sequence for Stringent Control of Bacterial Ribonucleic Acid Synthesis

In bacterial promoters subject to stringent control the heptanucleotide sequence spanning nucleotide positions -5 to +2 (defining the -10 region [D. Pribnow, J. Mol. Biol. 99:419-443, 1975] as positions -12 to -6) is highly conserved. This conserved region, for which the best consensus sequence is [Formula: see text] is absent from promoters which are not under stringent control.     

Site-directed cross-linking of mRNA analogues to 16S ribosomal RNA; a complete scan of cross-links from all positions between '+1' and '+16' on the mRNA, downstream from the decoding site.

mRNA analogues containing 4-thiouridine residues at selected sites were used to extend our analysis of photo-induced cross-links between mRNA and 16S RNA to cover the entire downstream range between positions +1 and +16 on the mRNA (position +1 is the 5'-base of the P-site codon). No tRNA-dependent cross-links were observed from positions +1, +2, +3 or +5. Position +4 on the mRNA was cross-linked in a tRNA-dependent manner to 16S RNA at a site between nucleotides ca 1402-1415 (most probably to the modified residue 1402), and this was absolutely specific for the +4 position. Similarly, the previously observed cross-link to nucleotide 1052 was absolutely specific for the +6 position. The previously observed cross-links from +7 to nucleotide 1395 and from +11 to 532 were however seen to a lesser extent with certain types of mRNA sequence from neighbouring positions (+6 to +10, and +10 to +13, respectively); no tRNA-dependent cross-links to other sites on 16S RNA were found from these positions, and no cross-linking was seen from positions +14 to +16. In each case the effect of a second cognate tRNA (at the ribosomal A-site) on the level of cross-linking was studied, and the specificity of each cross-link was confirmed by translocation experiments with elongation factor G, using appropriate mRNA analogues.     

Congenital end-plate acetylcholinesterase deficiency caused by a nonsense mutation and an A-->G splice-donor-site mutation at position +3 of the collagenlike-tail-subunit gene (COLQ): how does G at position +3 result in aberrant splicing?

Congenital end-plate acetylcholinesterase (AChE) deficiency (CEAD), the cause of a disabling myasthenic syndrome, arises from defects in the COLQ gene, which encodes the AChE triple-helical collagenlike-tail subunit that anchors catalytic subunits of AChE to the synaptic basal lamina. Here we describe a patient with CEAD with a nonsense mutation (R315X) and a splice-donor-site mutation at position +3 of intron 16 (IVS16+3A-->G) of COLQ. Because both A and G are consensus nucleotides at the +3 position of splice-donor sites, we constructed a minigene that spans exons 15-17 and harbors IVS16+3A-->G for expression in COS cells. We found that the mutation causes skipping of exon 16. The mutant splice-donor site of intron 16 harbors five discordant nucleotides (at -3, -2, +3, +4, and +6) that do not base-pair with U1 small-nuclear RNA (snRNA), the molecule responsible for splice-donor-site recognition. Versions of the minigene harboring, at either +4 or +6, nucleotides complementary to U1 snRNA restore normal splicing. Analysis of 1,801 native splice-donor sites reveals that presence of a G nucleotide at +3 is associated with preferential usage, at positions +4 to +6, of nucleotides concordant to U1 snRNA. Analysis of 11 disease-associated IVS+3A-->G mutations indicates that, on average, two of three nucleotides at positions +4 to +6 fail to base-pair, and that the nucleotide at +4 never base-pairs, with U1 snRNA. We conclude that, with G at +3, normal splicing generally depends on the concordance that residues at +4 to +6 have with U1 snRNA, but other cis-acting elements may also be important in assuring the fidelity of splicing.     

Frameshift mutagenesis in Escherichia coli by reversible DNA intercalators: sequence specificity

The mutagenic potency of the simple reversible intercalators isopropyl-OPC (iPr-OPC) and 9-aminoacridine (9-AA) is assessed in E. coli using reversion assays based on plasmids derived from pBR322 carrying various frameshift mutations within the tetracycline resistance gene in repetitive sequences: +/- 2 frameshift mutations within alternating GC sequences; +/- 1 frameshift mutation at runs of guanines. The results obtained show that iPr-OPC and 9-AA have a sequence specificity for mutagenesis: they revert +1 and -1 frameshift mutations within runs of monotonous G:C base pairs. The precise determination of the size of a small restriction fragment which contains the mutation allowed us to demonstrate that reversion occurred by -1 deletions for the +1 frameshift mutations and by +1 additions for the -1 frameshift mutations. The possible relations of this specific reversion with the base sequence specificity of the mutagenesis are briefly discussed.     

Phylogeny of Greya (Lepidoptera: Prodoxidae), based on nucleotide sequence variation in mitochondrial cytochrome oxidase I and II: congruence with morphological data

The phylogeny of Greya Busck (Lepidoptera: Prodoxidae) was inferred from nucleotide sequence variation across a 765-bp region in the cytochrome oxidase I and II genes of the mitochondrial genome. Most parsimonious relationships of 25 haplotypes from 16 Greya species and two outgroup genera (Tetragma and Prodoxus) showed substantial congruence with the species relationships indicated by morphological variation. Differences between mitochondrial and morphological trees were found primarily in the positions of two species, G. variabilis and G. pectinifera, and in the branching order of the three major species groups in the genus. Conflicts between the data sets were examined by comparing levels of homoplasy in characters supporting alternative hypotheses. The phylogeny of Greya species suggests that host-plant association at the family level and larval feeding mode are conservative characters. Transition/transversion ratios estimated by reconstruction of nucleotide substitutions on the phylogeny had a range of 2.0-9.3, when different subsets of the phylogeny were used. The decline of this ratio with the increase in maximum sequence divergence among taxa indicates that transitions are masked by transversions along deeper internodes or long branches of the phylogeny. Among transitions, substitutions of A-->G and T-->C outnumbered their reciprocal substitutions by 2-6 times, presumably because of the approximately 4:1 (77%) A+T-bias in nucleotide base composition. Of all transversions, 73%-80% were A<-->T substitutions, 85% of which occurred at third positions of codons; these estimates did not decrease with an increase in maximum sequence divergence of taxa included in the analysis. The high frequency of A<-->T substitutions is either a reflection or an explanation of the 92% A+T bias at third codon positions.      

Several regions of a tRNA precursor determine the Escherichia coli RNase P cleavage site.

The RNase P cleavage reaction was studied as a function of the number of base-pairs in the acceptor-stem and/or T-stem of a natural tRNA precursor, the tRNA(Tyr)Su3 precursor. Our data suggest that the location of the Escherichia coli RNase P cleavage site does not depend merely on the lengths of the acceptor-stem and T-stem as previously suggested. Surprisingly, we find that precursors with only four base-pairs in the acceptor-stem are cleaved by M1 RNA and by holoenzyme. Furthermore, we show that both disruption of base-pairing, and alteration of the nucleotide sequence (without disruption of base-pairing) proximal to the cleavage site result in aberrant cleavage. Thus, the identity of the nucleotides near the cleavage site is important for recognition of the cleavage site rather than base-pairing. The important nucleotides are those at positions -2, -1, +1, +72, +73 and +74. We propose that the nucleotide at position +1 functions as a guiding nucleotide. These results raise the possibility that Mg2+ binding near the cleavage site is dependent on the identity of the nucleotides at these positions. In addition, we show that disruption of base-pairing in the acceptor-stem affects both Michaelis-Menten constants, Km and kcat.     

Correlation between sequence conservation of the 5' untranslated region and codon usage bias in Mus musculus genes

The codon adaptation index (CAI) values of all protein-coding sequences of the full-length cDNA libraries of Mus musculus were computed based on the RIKEN mouse full-length cDNA library. We have also computed the extent of consensus in flanking sequences of the initiator ATG codon based on the 'relative entropy' values of respective nucleotide positions (from -20 to +12 bp relative to the initiator ATG codon) for each group of genes classified by CAI values. With regard to the two nucleotides positions (-3 and +4) known to be highly conserved in Kozak's consensus sequence, a clear correlation between CAI values and relative entropy values was observed at position -3 but this was not significant at position +4, although a significant correlation was found at position -1 of the consensus sequence. Further, although no correlation was observed at any additional positions, relative entropy values were very high at positions -4, -6, and -8 in genes with high CAI values. These findings suggest that the extent of conservation in the flanking sequence of the initiator ATG codon including Kozak's consensus sequence was an important factor in modulation of the translation efficiency as well as synonymous codon usage bias particularly in highly expressed genes.     

Mutations That Affect the Efficiency of Translation of mRNA for the cII Gene of Coliphage Lambda

Starting with the lambda pRE-strain lambda ctr1 cy3008, which forms clear plaques, we have isolated two mutant strains, lambda dya2 ctr1 cy3008 and lambda dya3 ctr1 cy3008, that form plaques with very slightly turbid centers. The dya2 and dya3 mutations lie in the region of overlap between the PRE promoter and the ribosome recognition region of the cII gene, and have nucleotide alterations at positions -1 and +5 of pRE, and alterations in cII mRNA at -16 and -21 nucleotides before the initial AUG codon of the gene. Both mutations destabilize a stem structure that may be formed by cII mRNA, and dya2 also changes the sequence on cII mRNA that is complementary to the 3'-end of 16 S rRNA from 5'-UAAGGA-3' to 5'-UGAGGA-3'. --The dya2 and dya3 mutations, along with the ctr1 mutation, which destabilizes either of two alternate stem structures which may be formed by cII mRNA (these being more stable stem structures than the one affected by dya2 and dya3), were tested for their ability to reverse two cII-mutations that are characterized by inefficient translation of cII mRNA. These are cII3088, an A----G mutation four bases before the initial AUG codon, and cII3059, a GUU----GAU (Val2----Asp) second codon mutation. It was found that ctr1 completely reverses the translation defects of these two mutations, while dya2 partially reverses these translation defects. The dya3 mutation has no effect on translation efficiency under any condition tested.(ABSTRACT TRUNCATED AT 250 WORDS)