Influences of bacterial DNA base-ratios and amino acid composition of the gene product on the consequences of frameshift mutations

Using as an example the E. coli lac I (repressor) gene product, it can be shown that DNA base-ratio is a major determinant of the detailed outcomes of both +1 (?2) and ?1 (+2) types of frameshift mutations. Potential reinitiation codons (AUG or GUG) and premature stop codons (UAA, UAG and UGA) occur in very different proportions depending on the type of frameshift and the DNA base-ratio.A comparison of the H. halobium bacteriorhodopsin precursor gene with other actual and theoretical genes reveals that the amino acid composition of the gene product is a second, important, determinant of the detailed outcome of frameshift mutations. Rules are formulated for the occurrence of particular codon, and hence amino acid, doublets in AT-GC-rich or intermediate base-ratio DNA as these affect frameshift-generated reinitiation and premature stop codons.     

The influence of DNA base ratios on intragenic mutational spectra in prokaryotes

Because of the degeneracy of the Genetic Code, a given amino acid sequence can be written in terms of either GC- or AT-rich DNA-specified mRNA codons. It is shown that the genetic consequences in terms of detectable single base pair substitution mutations are very different in such GC- or AT-rich sequences. Furthermore the occurrences of potential frameshift mutation hot-spots, in runs of repeated or reiterated base pairs, and sites of mutationally important base methylation, also differ markedly between AT- and GC-rich sequences. Thus the evolutionary pathways which can be followed by using single or multiple changes in the amino acid sequence of a given polypeptide will be very different in prokaryotes with GC-rich or AT-rich DNA.     

Resolving a fidelity paradox: why Escherichia coli DNA polymerase II makes more base substitution errors in AT- compared with GC-rich DNA.

The activity of DNA polymerase-associated proofreading 3'-exonucleases is generally enhanced in less stable DNA regions leading to a reduction in base substitution error frequencies in AT- versus GC-rich sequences. Unexpectedly, however, the opposite result was found for Escherichia coli DNA polymerase II (pol II). Nucleotide misincorporation frequencies for pol II were found to be 3-5-fold higher in AT- compared with GC-rich DNA, both in the presence and absence of polymerase processivity subunits, beta dimer and gamma complex. In contrast, E. coli pol III holoenzyme, behaving "as expected," exhibited 3-5-fold lower misincorporation frequencies in AT-rich DNA. A reduction in fidelity in AT-rich regions occurred for pol II despite having an associated 3'-exonuclease proofreading activity that preferentially degrades AT-rich compared with GC-rich DNA primer-template in the absence of DNA synthesis. Concomitant with a reduction in fidelity, pol II polymerization efficiencies were 2-6-fold higher in AT-rich DNA, depending on sequence context. Pol II paradoxical fidelity behavior can be accounted for by the enzyme's preference for forward polymerization in AT-rich sequences. The more efficient polymerization suppresses proofreading thereby causing a significant increase in base substitution error rates in AT-rich regions.     

Codon cassette mutagenesis: a general method to insert or replace individual codons by using universal mutagenic cassettes.

We describe codon cassette mutagenesis, a simple method of mutagenesis that uses universal mutagenic cassettes to deposit single codons at specific sites in double-stranded DNA. A target molecule is first constructed that contains a blunt, double-strand break at the site targeted for mutagenesis. A double-stranded mutagenic codon cassette is then inserted at the target site. Each mutagenic codon cassette contains a three base pair direct terminal repeat and two head-to-head recognition sequences for the restriction endonuclease Sapl, an enzyme that cleaves outside of its recognition sequence. The intermediate molecule containing the mutagenic cassette is then digested with Sapl, thereby removing most of the mutagenic cassette, leaving only a three base cohesive overhang that is ligated to generate the final insertion or substitution mutation. A general method for constructing blunt-end target molecules suitable for this approach is also described. Because the mutagenic cassette is excised during this procedure and alters the target only by introducing the desired mutation, the same cassette can be used to introduce a particular codon at all target sites. Each cassette can deposit two different codons, depending on the orientation in which it is inserted into the target molecule. Therefore, a series of eleven cassettes is sufficient to insert all possible amino acids at any constructed target site. Thus codon cassettes are 'off-the-shelf' reagents, and this methodology should be a particularly useful and inexpensive approach for subjecting multiple different positions in a protein sequence to saturation mutagenesis.     

Evolutionary implications of potential frameshift hot-spots, methylation sites, and codon usage differences between genes and between organisms in Escherichia coli and Salmonella typhimurium

The published nucleotide sequences of the E. coli and S. typhimurium trp A and trp B genes show a high degree of similarity between homologous genes of the two organisms, and an even greater degree of similarity between the amino acid sequences of the gene products. In spite of this, analysis of the nucleotide sequences reveals that there are marked differences between E. coli and S. typhimurium genes with respect to potential frameshift mutation hot-spots and dam and mec, mutationally important, methylation sites. Such existing differences may well lead to divergent evolution of these two, presently closely related, bacteria. Codon usage patterns in the trp A and trp B genes of E. coli and S. typhimurium, and the lac I gene of E. coli, have been re-analysed in terms of AT-rich, GC-rich, neutral, or unique codons and marked preferences found. In some cases particular amino acids are most often specified by AT-rich, in others by the GC-rich, alternative codons. In still other cases the codon preference depends on the gene studied. These patterns can be interpreted in terms of enteric bacterial evolution, via hybridizations, from ancestral bacteria with AT- or GC-rich DNA.     

DNA base changes and RNA levels in N-acetoxy-2-acetylaminofluorene-induced dihydrofolate reductase mutants of Chinese hamster ovary cells

Formerly, we isolated a series of dihydrofolate reductase-deficient Chinese hamster ovary cell mutants that were induced by N-acetoxy-2-acetylaminofluorene. Deletions and complex gene rearrangements were detected in 28% of these mutants; 72% contained putative point mutations. In the present study, we have localized the putative point mutations in the 25,000 base dhfr gene by RNase heteroduplex mapping. Assignment of a position for each mutation was successful in 16 of 19 mutants studied. We cloned DNA fragments containing the mapped mutations from nine mutants into a bacteriophage lambda vector. In the case of 11 other mutants, DNA was amplified by the polymerase chain reaction procedure. Sequence analysis of cloned and amplified DNA confirmed the presence of point mutations. Most mutants (90%) carried base substitutions; the rest contained frameshift mutations. Of the point mutations, 75% were G.C to T.A transversions in either the dhfr coding sequence or at splice sites; transition G.C to A.T mutations were found in two mutants (10%). In one of these transition mutants, the base substitution occurred at the fifth base of the third intron. Of the frameshift mutations, one was a deletion of G.C pair and the other was an insertion of an A.T pair. Of the mapped mutants, 38% exhibited greatly reduced (approximately 10-fold) steady-state levels of dhfr mRNA. All eight sequenced mutants displaying this phenotype contained premature chain termination codons. Normal levels of dhfr mRNA were observed in five missense mutants and in five mutants carrying nonsense codons in the translated portion of exon VI. Taken together with the results of other mutagens at this locus, we conclude that the low dhfr mRNA phenotype is correlated with the presence of nonsense codons in exons II to V but not in the last exon of the dhfr gene.     

Functional mutation of DNA polymerase beta found in human gastric cancer--inability of the base excision repair in vitro.

DNA polymerase beta (polbeta) is one of mammalian DNA polymerases and is known to be involved in a G:T/G:U mismatch repair. In order to investigate an involvement of this enzyme in a base excision repair, we searched a mutation of human polbeta in human gastric cancer and studied a function of the mutation. We observed cancer-specific missense mutations in 6 of 20 samples. All of these mutations were, however, heterozygous. We further analyzed the base excision repair activity of these mutants to know whether these mutants cause an error of mismatch repair. One of these mutants, which resulted in an amino acid substitution of Glu for Lys at codon 295, showed an inhibitory effect by in vitro base excision repair assay, suggesting that this mutation might play some role in carcinogenesis of the gastric mucosa.     

Characterization of mutations in the factor VIII gene by direct sequencing of amplified genomic DNA

In order to search for mutations resulting in hemophilia A that are not detectable by restriction analysis, three regions of the factor VIII gene were chosen for direct sequence analysis. Short segments of genomic DNA of 127 unrelated patients with hemophilia A were amplified by polymerase chain reaction. A total of 136,017 nucleotides were sequenced, and four mutations leading to the disease were found: a frameshift at codon 360 due to deletion of two nucleotides (GA), a nonsense codon 1705 due to a C----T transition, and two missense codons at positions 1699 and 1708. The first missense mutation (A----T) results in a Tyr----Phe substitution at a putative von Willebrand factor binding site. The second results in an Arg----Cys substitution at a thrombin cleavage site. In addition, we identified three rare sequence variants: a silent C----T transition at codon 34 which does not result in an amino acid change, a G----C change at codon 345 (Val----Leu), and an A----G change at the third nucleotide of intron 14. Direct sequence analysis of amplified DNA is a powerful but labor-intensive method of identifying mutations in large genes such as the human factor VIII gene.     

Molecular evolution of the histone 3 multigene family in the Drosophila melanogaster species subgroup

Molecular evolution of the histone multigene family was studied by cloning and sequencing regions of the histone 3 gene in the Drosophila melanogaster species subgroup. Analysis of the nucleotide substitution pattern showed that in the coding region synonymous changes occurred more frequently to A or T in contrast to the GC-rich base composition, while in the 3' region the nucleotide substitutions were most likely in equilibrium. These results suggested that the base composition at the third codon position of the H3 gene, i.e., codon usage, has been changing to A or T in the Drosophila melanogaster species subgroup.     

New Methods for Detecting Positive Selection at Single Amino Acid Sites

Inferring positive selection at single amino acid sites is of particular importance for studying evolutionary mechanisms of a protein. For this purpose, Suzuki and Gojobori (1999) developed a method (SG method) for comparing the rates of synonymous and nonsynonymous substitutions at each codon site in a protein-coding nucleotide sequence, using ancestral codons at interior nodes of the phylogenetic tree as inferred by the maximum parsimony method. In the SG method, however, selective neutrality of nucleotide substitutions cannot be tested at codon sites, where only termination codons are inferred at any interior node or the number of equally parsimonious inferences of ancestral codons at all interior nodes exceeds 10,000. Here I present a modified SG method which is free from these problems. Specifically, I use the distance-based Bayesian method for inferring the single most likely ancestral codon from 61 sense codons at each interior node. In the computer simulation and real data analysis, the modified SG method showed a higher overall efficiency of detecting positive selection than the original SG method, particularly at highly polymorphic codon sites. These results indicate that the modified SG method is useful for inferring positive selection at codon sites where neutrality cannot be tested by the original SG method. I also discuss that the p-distance is preferable to the number of synonymous substitutions for inferring the phylogenetic tree in the SG method, and present a maximum likelihood method for detecting positive selection at single amino acid sites, which produced reasonable results in the real data analysis.     

Molecular basis of beta (0)-thalassemia/HbE disease in Thailand

The molecular basis of beta(0)-thalassemia/HbE disease in 30 Thai patients was investigated using DNA amplification and dot-blot hybridization with a number of allele specific oligonucleotide probes. The mutations identified were 17 cases of 4 base-pair deletion at codons 41-42, 4 cases of amber mutation at codon 17, and one case each of an ochre mutation at codon 35, a single base substitution at position 5 of IVS-1, and a single base substitution at position 654 of IVS-2.     

A relationship between GC content and coding-sequence length

Since base composition of translational stop codons (TAG, TAA, and TGA) is biased toward a low G+C content, a differential density for these termination signals is expected in random DNA sequences of different base compositions. The expected length of reading frames (DNA segments of sense codons flanked by in-phase stop codons) in random sequences is thus a function of GC content. The analysis of DNA sequences from several genome databases stratified according to GC content reveals that the longest coding sequences—exons in vertebrates and genes in prokaryotes—are GC-rich, while the shortest ones are GC-poor. Exon lengthening in GC-rich vertebrate regions does not result, however, in longer vertebrate proteins, perhaps because of the lower number of exons in the genes located in these regions. The effects on coding-sequence lengths constitute a new evolutionary meaning for compositional variations in DNA GC content. Correspondence to: J. L. Oliver     

Reversion of trpA nonsense mutations by deletion of the chain-termination codons

This paper describes a novel mechanism for reversion of nonsense mutations in the trpA gene of Escherichia coli. This mechanism, deletion of the nonsense codon, was discovered in the course of selecting for missense revertants of trpA(UGA211) and for catalytically active tryptophan synthetase alpha chain revertants of trpA(UAA234) and trpA(UAG234). Each type of revertant trpA was cloned and its DNA sequence determined. trpA(UGA211) gave rise to two previously unidentified types of missense revertant. The first type was expected, namely trpA(CGA211), the result of a base substitution event. The other type, representing approximately 1% of the missense revertants, was unexpected on the basis of single base substitutions and an understanding of which amino acids are functional at alpha chain position 211. It was found to be the result of a 21 base-pair deletion of a region containing codon 211. The tryptophan-independent revertants of both position 234 nonsense mutants occurred at a frequency of approximately 2 per 10(9) viable cells. They were identical in that they both resulted from a 3 base-pair deletion, namely deletion of the chain-terminating codon at position 234. One of them, however, also displayed an A instead of the normal G in the third position of codon 235. The revertants were characterized according to growth in different media and tryptophan synthetase assays performed on crude extracts. These types of mutants should prove interesting and important for the elucidation of alpha chain structure-function relationships, for insight into the assembly and interaction of subunits in this model multienzyme complex, and for the study of mechanisms by which deletions can be generated.     

Decoding the genome: a modified view

Transfer RNA’s role in decoding the genome is critical to the accuracy and efficiency of protein synthesis. Though modified nucleosides were identified in RNA 50 years ago, only recently has their importance to tRNA’s ability to decode cognate and wobble codons become apparent. RNA modifications are ubiquitous. To date, some 100 different posttranslational modifications have been identified. Modifications of tRNA are the most extensively investigated; however, many other RNAs have modified nucleosides. The modifications that occur at the first, or wobble position, of tRNA’s anticodon and those 3′-adjacent to the anticodon are of particular interest. The tRNAs most affected by individual and combinations of modifications respond to codons in mixed codon boxes where distinction of the third codon base is important for discriminating between the correct cognate or wobble codons and the incorrect near-cognate codons (e.g. AAA/G for lysine versus AAU/C asparagine). In contrast, other modifications expand wobble codon recognition, such as U·U base pairing, for tRNAs that respond to multiple codons of a 4-fold degenerate codon box (e.g. GUU/A/C/G for valine). Whether restricting codon recognition, expanding wobble, enabling translocation, or maintaining the messenger RNA, reading frame modifications appear to reduce anticodon loop dynamics to that accepted by the ribosome. Therefore, we suggest that anticodon stem and loop domain nucleoside modifications allow a limited number of tRNAs to accurately and efficiently decode the 61 amino acid codons by selectively restricting some anticodon–codon interactions and expanding others.     

Texture analysis of fluorescence lifetime images of AT- and GC-rich regions in nuclei.

We used intensity and fluorescence lifetime microscopy (FLIM) of 3T3 nuclei to investigate the existence of AT-rich and GC-rich regions of the nuclear DNA. Hoechst 33258 (Ho) and 7-aminoactinomycin D (7-AAD) were used as fluorescence probes specific for AT and GC base pairs, respectively. YOYO-1 (Yo) was used as a dye that displays distinct fluorescence lifetimes when bound to AT or GC base pairs. We combined fluorescence imaging of Ho and 7-AAD with time-resolved measurements of Yo and took advantage of an additional information content of the time-resolved fluorescence. Because a single nucleus could not be stained and measured with all three dyes, we used texture analysis to compare the spatial distribution of AT-rich and GC-rich DNA in 100 nuclei in different phases of the cell cycle. The fluorescence intensity-based analysis of Ho- or 7-AAD-stained images indicates increased number and larger size of the DNA condensation centers in the G2/M-phases compared to G0/1-phases. The lifetime-based study of Yo-stained images suggests spatial separation of the AT- or GC-rich DNA regions in the G2/M-phase. Texture analysis of fluorescence intensity and lifetime images was used to quantitatively study the spatial change of condensation and separation of AT- and GC-rich DNA during the cell cycle.     

Amino Acid Usage Is Asymmetrically Biased in AT- and GC-Rich Microbial Genomes

Introduction

Genomic base composition ranges from less than 25% AT to more than 85% AT in prokaryotes. Since only a small fraction of prokaryotic genomes is not protein coding even a minor change in genomic base composition will induce profound protein changes. We examined how amino acid and codon frequencies were distributed in over 2000 microbial genomes and how these distributions were affected by base compositional changes. In addition, we wanted to know how genome-wide amino acid usage was biased in the different genomes and how changes to base composition and mutations affected this bias. To carry this out, we used a Generalized Additive Mixed-effects Model (GAMM) to explore non-linear associations and strong data dependences in closely related microbes; principal component analysis (PCA) was used to examine genomic amino acid- and codon frequencies, while the concept of relative entropy was used to analyze genomic mutation rates.

Results

We found that genomic amino acid frequencies carried a stronger phylogenetic signal than codon frequencies, but that this signal was weak compared to that of genomic %AT. Further, in contrast to codon usage bias (CUB), amino acid usage bias (AAUB) was differently distributed in AT- and GC-rich genomes in the sense that AT-rich genomes did not prefer specific amino acids over others to the same extent as GC-rich genomes. AAUB was also associated with relative entropy; genomes with low AAUB contained more random mutations as a consequence of relaxed purifying selection than genomes with higher AAUB.

Conclusion

Genomic base composition has a substantial effect on both amino acid- and codon frequencies in bacterial genomes. While phylogeny influenced amino acid usage more in GC-rich genomes, AT-content was driving amino acid usage in AT-rich genomes. We found the GAMM model to be an excellent tool to analyze the genomic data used in this study.     

Nonrandom frequency patterns of synonymous substitutions in homologous mammalian genes

All 69 homologous coding sequences that are currently available in four mammalian orders were aligned and the synonymous positions of quartet and duet (fourfold and twofold degenerate) codons were divided into three classes (that will be called conserved, intermediate, and variable) according to whether they show no change, one change, or more than one change, respectively. We observed (1) that the frequencies of conserved, intermediate, and variable positions of quartet and duet codons are different in different genes; (2) that the frequencies of the three classes are significantly different from expectations based on a random substitution process in the majority of genes (especially for GC-rich genes) for quartet codons and in a minority of genes for doublet codons; and (3) that the frequencies of the three classes of positions of quartet codons are correlated with those of duet codons, the conserved positions of quartet and duet codons being, in addition, correlated with the degree of amino acid conservation. Our main conclusions are that synonymous substitution frequencies: (1) are gene-specific; (2) are not simply the result of a stochastic process in which nucleotide substitutions accumulate at random, over time; and (3) are correlated in quartet and duet codons.