Molecular mechanisms of substitution mutagenesis. An experimental test of the Watson-Crick and topal-fresco models of base mispairings

High frequency of transition to transversion ratio in the stem region of RNA secondary structure of untranslated region of SARS-CoV-2

Introduction

The propensity of nucleotide bases to form pairs, causes folding and the formation of secondary structure in the RNA. Therefore, purine (R): pyrimidine (Y) base-pairing is vital to maintain uniform lateral dimension in RNA secondary structure. Transversions or base substitutions between R and Y bases, are more detrimental to the stability of RNA secondary structure, than transitions derived from substitutions between A and G or C and T. The study of transversion and transition base substitutions is important to understand evolutionary mechanisms of RNA secondary structure in the 5'  and 3'  untranslated (UTR) regions of SARS-CoV-2. In this work, we carried out comparative analysis of transition and transversion base substitutions in the stem and loop regions of RNA secondary structure of SARS-CoV-2.

Methods

We have considered the experimentally determined and well documented stem and loop regions of 5' and 3' UTR regions of SARS-CoV-2 for base substitution analysis. The secondary structure comprising of stem and loop regions were visualized using the RNAfold web server. The GISAID repository was used to extract base sequence alignment of the UTR regions. Python scripts were developed for comparative analysis of transversion and transition frequencies in the stem and the loop regions.

Results

The results of base substitution analysis revealed a higher transition (ti) to transversion (tv) ratio (ti/tv) in the stem region of UTR of RNA secondary structure of SARS-CoV-2 reported during the early stage of the pandemic. The higher ti/tv ratio in the stem region suggested the influence of secondary structure in selecting the pattern of base substitutions. This differential pattern of ti/tv values between stem and loop regions was not observed among the Delta and Omicron variants that dominated the later stage of the pandemic. It is noteworthy that the ti/tv values in the stem and loop regions were similar among the later dominant Delta and Omicron variant strains which is to be investigated to understand the rapid evolution and global adaptation of SARS-CoV-2.

Conclusion

Our findings implicate the lower frequency of transversions than the transitions in the stem regions of UTRs of SARS-CoV-2. The RNA secondary structures are associated with replication, translation, and packaging, further investigations are needed to understand these base substitutions across different variants of SARS-CoV-2.

Stem Region of tRNA Genes Favors Transition Substitution Towards Keto Bases in Bacteria

Transversion and transition mutations have variable effects on the stability of RNA secondary structure considering that the former destabilizes the double helix geometry to a greater extent by introducing purine:purine (R:R) or pyrimidine:pyrimidine (Y:Y) base pairs. Therefore, transversion frequency is likely to be lower than that of transition in the secondary structure regions of RNA genes. Here, we performed an analysis of transition and transversion frequencies in tRNA genes defined well with secondary structure and compared with the intergenic regions in five bacterial species namely Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Staphylococcus aureus and Streptococcus pneumoniae using a large genome sequence data set. In general, the transversion frequency was observed to be lower than that of transition in both tRNA genes and intergenic regions. The transition to transversion ratio was observed to be greater in tRNA genes than that in the intergenic regions in all the five bacteria that we studied. Interestingly, the intraspecies base substitution analysis in tRNA genes revealed that non-compensatory substitutions were more frequent than compensatory substitutions in the stem region. Further, transition to transversion ratio in the loop region was observed to be significantly lesser than that among the non-compensatory substitutions in the stem region. This indicated that the transversion is more deleterious than transition in the stem regions. In addition, substitutions from amino bases (A/C) to keto bases (G/T) were also observed to be more than the reverse substitutions in the stem region. Substitution from amino bases to keto bases are likely to facilitate the stable G:U pairing unlike the reverse substitution that facilitates the unstable A:C pairing in the stem region of tRNA. This work provides additional support that the secondary structure of tRNA molecule is what drives the different substitutions in its gene sequence.

High-resolution mapping of DNA polymerase fidelity using nucleotide imbalances and next-generation sequencing

DNA polymerase fidelity is affected by both intrinsic properties and environmental conditions. Current strategies for measuring DNA polymerase error rate in vitro are constrained by low error subtype sensitivity, poor scalability, and lack of flexibility in types of sequence contexts that can be tested. We have developed the Magnification via Nucleotide Imbalance Fidelity (MagNIFi) assay, a scalable next-generation sequencing assay that uses a biased deoxynucleotide pool to quantitatively shift error rates into a range where errors are frequent and hence measurement is robust, while still allowing for accurate mapping to error rates under typical conditions. This assay is compatible with a wide range of fidelity-modulating conditions, and enables high-throughput analysis of sequence context effects on base substitution and single nucleotide deletion fidelity using a built-in template library. We validate this assay by comparing to previously established fidelity metrics, and use it to investigate neighboring sequence-mediated effects on fidelity for several DNA polymerases. Through these demonstrations, we establish the MagNIFi assay for robust, high-throughput analysis of DNA polymerase fidelity.

A novel constructed SPT15 mutagenesis library of Saccharomyces cerevisiae by using gTME technique for enhanced ethanol production

During the last few years, the global transcription machinery engineering (gTME) technique has gained more attention as an effective approach for the construction of novel mutants. Genetic strategies (molecular biology methods) were utilized to get mutational for both genes (SPT15 and TAF23) basically existed in the Saccharomyces cerevisiae genome via screening the gTME approach in order to obtain a new mutant S. cerevisiae diploid strain. The vector pYX212 was utilized to transform these genes into the control diploid strain S. cerevisiae through the process of mating between haploids control strains S. cerevisiae (MAT-a [CICC 1374]) and (MAT-α [CICC 31144]), by using the oligonucleotide primers SPT15-EcoRI-FW/SPT15-SalI-RV and TAF23-SalI-FW/TAF23-NheI-RV, respectively. The resultant mutants were examined for a series of stability tests. This study showed how strong the effect of using strategic gTME with the importance of the modification we conducted in Error Prone PCR protocol by increasing MnCl2 concentration instead of MgCl2. More than ninety mutants we obtained in this study were distinguished by a high level production of bio-ethanol as compared to the diploid control strain.

Can DNA-binding proteins of replisome tautomerize nucleotide bases? Ab initio model study

Ab initio quantum-chemical study of specific point contacts of replisome proteins with DNA modeled by acetic acid with canonical and mutagenic tautomers of DNA bases methylated at the glycosidic nitrogen atoms was performed in vacuo and continuum with a low dielectric constant (ϵ ∼ 4) corresponding to a hydrophobic interface of protein-nucleic acid interaction. All tautomerized complexes were found to be dynamically unstable, because the electronic energies of their back-reaction barriers do not exceed zero-point vibrational energies associated with the vibrational modes whose harmonic vibrational frequencies become imaginary in the transition states of the tautomerization reaction. Additionally, based on the physicochemical arguments, it was demonstrated that the effects of biomolecular environment cannot ensure dynamic stabilization. This result allows suggesting that hypothetically generated by DNA-binding proteins of replisome rare tautomers will have no impact on the total spontaneous mutation due to the low reverse barrier allowing a quick return to the canonical form.

Discovery of mutations in Saccharomyces cerevisiae by pooled linkage analysis and whole-genome sequencing

Many novel and important mutations arise in model organisms and human patients that can be difficult or impossible to identify using standard genetic approaches, especially for complex traits. Working with a previously uncharacterized dominant Saccharomyces cerevisiae mutant with impaired vacuole inheritance, we developed a pooled linkage strategy based on next-generation DNA sequencing to specifically identify functional mutations from among a large excess of polymorphisms, incidental mutations, and sequencing errors. The VAC6-1 mutation was verified to correspond to PHO81-R701S, the highest priority candidate reported by VAMP, the new software platform developed for these studies. Sequence data further revealed the large extent of strain background polymorphisms and structural alterations present in the host strain, which occurred by several mechanisms including a novel Ty insertion. The results provide a snapshot of the ongoing genomic changes that ultimately result in strain divergence and evolution, as well as a general model for the discovery of functional mutations in many organisms.

Evidence that localized variation in primate sequence divergence arises from an influence of nucleosome placement on DNA repair

Understanding the origins of localized substitution rate heterogeneity has important implications for identifying functional genomic sequences. Outside of gene regions, the origins of rate heterogeneity remain unclear. Experimental studies establish that chromatin compaction affects rates of both DNA lesion formation and repair. A functional association between chromatin status and 5-methyl-cytosine also exists. These suggest that both the total rate and the type of substitution will be affected by chromatin status. Regular positioning of nucleosomes, the building block of chromatin, further predicts that substitution rate and type should vary spatially in an oscillating manner. We addressed chromatin's influence on substitution rate and type in primates. Matched numbers of sites were sampled from Dnase I hypersensitive (DHS) and closed chromatin control flank (Flank). Likelihood ratio tests revealed significant excesses of total and of transition substitutions in Flank compared with matched DHS for both intergenic and intronic samples. An additional excess of CpG transitions was evident for the intergenic, but not intronic, regions. Fluctuation in substitution rate along approximately 1,800 primate promoters was measured using phylogenetic footprinting. Significant positive correlations were evident between the substitution rate and a nucleosome score from resting human T-cells, with up to approximately 50% of the variance in substitution rate accounted for. Using signal processing techniques, a dominant oscillation at approximately 200 bp was evident in both the substitution rate and the nucleosome score. Our results support a role for differential DNA repair rates between open and closed chromatin in the spatial distribution of rate heterogeneity.

Effect of ribavirin on the mutation rate and spectrum of hepatitis C virus in vivo

Their extremely error-prone replication makes RNA viruses targets for lethal mutagenesis. In the case of hepatitis C virus (HCV), the standard treatment includes ribavirin, a base analog with an in vitro mutagenic effect, but the in vivo mode of action of ribavirin remains poorly understood. Here, we test the mutagenic effects of ribavirin plus interferon treatment in vivo using a new method to estimate mutation rates based on the analysis of nonsense mutations. We apply this methodology to a large HCV sequence database containing over 15,000 reverse transcription-PCR molecular clone sequences from 74 patients infected with HCV. We obtained an estimate of the spontaneous mutation rate of ca. 10(-4) substitutions per site or lower, a value within the typically accepted range for RNA viruses. A roughly threefold increase in mutation rate and a significant shift in mutation spectrum were observed in samples from patients undergoing 6 months of interferon plus ribavirin treatment. This result is consistent with the known in vitro mutagenic effect of ribavirin and suggests that the antiviral effect of ribavirin plus interferon treatment is at least partly exerted through lethal mutagenesis.

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

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

Hypoxanthine incorporation is nonmutagenic in Escherichia coli

Endonuclease V, encoded by the nfi gene, initiates removal of the base analogs hypoxanthine and xanthine from DNA, acting to prevent mutagenesis from purine base deamination within the DNA. On the other hand, the RdgB nucleotide hydrolase in Escherichia coli is proposed to prevent hypoxanthine and xanthine incorporation into DNA by intercepting the noncanonical DNA precursors dITP and dXTP. Because many base analogs are mutagenic when incorporated into DNA, it is intuitive to think of RdgB as acting to prevent similar mutagenesis from deaminated purines in the DNA precursor pools. To test this idea, we used a set of Claire Cupples' strains to detect changes in spontaneous mutagenesis spectra, as well as in nitrous acid-induced mutagenesis spectra, in wild-type cells and in rdgB single, nfi single, and rdgB nfi double mutants. We found neither a significant increase in spontaneous mutagenesis in rdgB and nfi single mutants or the double mutant nor any changes in nitrous acid-induced mutagenesis for rdgB mutant strains. We conclude that incorporation of deaminated purines into DNA is nonmutagenic.

Role in tumorigenesis of silent mutations in the TP53 gene

Over 10,000 mutations in the TP53 suppressor gene have been recorded in the International Agency for Research on Cancer (IARC) tumor data base. About 4% of these mutations are silent. It is a question whether these mutations play a role in tumor development. In order to approach this question, we asked whether the reported silent mutations are randomly distributed throughout the TP53 gene. The p53 data base was searched exon by exon. From the frequency of codons with no silent mutations, the average number of silent mutations per codon for each exon was calculated using the Poisson distribution. The results indicate the distribution to be non-random. About one-third of all silent mutations occur in "hot-spots" and after subtraction of these hot-spots, the remaining silent mutations are randomly distributed. In addition, the percentage of silent mutations among the total in the silent mutation hot-spots is close to that expected for random mutation. We conclude that most of the silent mutations recorded in tumors play no role in tumor development and that the percentage of silent mutation is an indication of the amount of random mutation during tumorigenesis. Silent mutations occur to a significantly different extent in different tumor types. Tumors of the esophagus and colon have a low frequency of silent mutations, tumors of the prostate have a high frequency.

Expression of Trigonopsis variabilis D-amino acid oxidase gene in Escherichia coli and characterization of its inactive mutants

The D-amino acid oxidase cDNA gene (daao) of Trigonopsis variabilis was prepared by reverse transcriptase-polymerase chain reaction (PCR) and cloned into Escherichia coli expression vector, pTrc99A, under the control of tac promoter. Expression of daao gene significantly affected the growth and morphology of E. coli. The highest D-amino acid oxidase (DAAO) activity was 705 U (mg of protein)(-)(1), which was about 12-fold higher than that of D-alanine-induced T. variabilis. The DAAO protein exhibited activity on native-PAGE and had a M(r)value of 39.3 kDa. We also constructed an expression plasmid, pKm-DAAO, in which kanamycin instead of ampicillin was used as the selective marker. High-performance liquid chromatography (HPLC) analysis demonstrated that cephalosporin C could be converted to 7-glutarylcephalosporanic acid by cell-free extract of E. coli harboring pKm-DAAO. Four inactive DAAO mutants were obtained by error-prone PCR. Sequence analysis of these four DAAO mutants indicated the occurrence of mutations at Val-167, Pro-291, Pro-309, and Ala-343 residues. The His(6)-tagged DAAOs were expressed in E. coli and purified by nickel ion affinity chromatography. The results showed that all DAAO mutants lost their enzymatic activities and characteristic adsorption spectra for flavoenzyme. Based on the crystal structure of a homologous protein, pig DAAO, it is suggested that these four residues may play essential structural roles in DAAO conformation, thereby influencing DAAO's catalytic activity.

Heavy metal mutagenicity: insights from bioinorganic model chemistry

The mutagenicity of metal species may be the result of a direct interaction with the target molecule DNA. Possible scenarios leading to nucleobase mispairing are discussed, and selected examples are presented. They include changes in nucleobase selectivity as a consequence of alterations in acid-base properties of nucleobase atoms and groups involved in complementary H bond formation, guanine deprotonation, and stabilization of rare nucleobase tautomers by metal ions. Oxidative nucleobase damage brought about by metal species will not be considered.

Strand asymmetries in DNA evolution

The complementary strands of DNA differ with respect to replication and transcription. Both of these processes are asymmetric and can bias the occurrence of mutations between the strands: during replication, the discontinuous lagging strand undergoes certain errors at higher rates, and transcription overexposes the nontranscribed strand to DNA damage while targeting repair enzymes to the transcribed strand. While biases introduced during replication apparently have little impact on sequence evolution, the effects of transcription are observed in the asymmetric patterns of substitution in bacterial genes and might be influencing genome-wide patterns of base composition.

T4 phage evolution data in terms of a time-dependent Topal-Fresco mechanism

A physical interpretation of the Topal-Fresco [Nature 263, 285 (1976)] model for spontaneous base substitutions suggests that hydrogen-bonded DNA protons satisfy the criteria for a classical noninteracting isolated system. Accessible states for duplex G-C protons include the keto-amino state and the six complementary enol-imine isomers. Hydrogen-bonded enol and imine protons occupy symmetric double-minima created by the two sets of indistinguishable electron lone pairs and a single proton belonging to each enol-imine end group. These protons will consequently participate in coupled quantum mechanical flip-flop, tunneling back and forth between symmetric energy wells. This results in a quantum mixing of proton energy states where the lowest energy state will be a linear combination of available G-C isomers. The resulting conclusion is that metastable keto-amino G-C protons will populate accessible enol-imine stationary states at rates governed by quantum laws of statistical equilibrium, consistent with achieving the lowest energy condition for duplex G-C protons. Enol-imine G-C stationary states are bound more tightly, of the order of 3 to 12 kcal/mol, which requires a modified mode of Topal-Fresco replication that will inhibit reequilibration of enol and imine G and C template isomers and, thus, promote the formation of complementary mispairs. The model is demonstrated on time-dependent base substitutions expressed by T4 phage DNA systems where data are consistent with model explanations, including the prediction that time-dependent evolutionary transversion sites will exhibit both G-C-to-T-A and G-C-to-C-G transversions at replication, due to proton flip-flop alteration of G template genetic specificity. The observation that A-T sites are resistant to time-dependent evolutionary base substitutions, expressed exclusively at G-C sites, allows codons to be classified as either evolutionary sensitive (16 codons) or evolutionary resistant (8 codons). These criteria provide possible explanations for expansion properties of the CGG fragile X sequences. Enol-imine G-C stationary states appear to have been misdiagnosed as deamination of cytosine and oxidation of guanine to 8-hydroxy-guanine.

Base mismatches and mutagenesis: how important is tautomerism?

Tautomerism of nucleotides is an attractive model to explain transitional mutations; the structure of the mismatched base pairs in their enol or imino forms do not distort the DNA double helix. However recent structural data suggest that the mismatched nucleotides are in their normal keto and amino forms. 'Wobble' and other base pairings (for transversions) have C1'-C1' distances close to those found in the classical Watson-Crick base pairs. Kinetic studies show substrates are in their major tautomeric forms in their transition states. This suggests tautomerism may not be important for substitution mutations.

Biochemical basis of DNA replication fidelity

DNA polymerase is the critical enzyme maintaining genetic integrity during DNA replication. Individual steps in the replication process that contribute to DNA synthesis fidelity include nucleotide insertion, exonucleolytic proofreading, and binding to and elongation of matched and mismatched primer termini. Each process has been investigated using polyacrylamide gel electrophoresis (PAGE) to resolve 32P-labeled primer molecules extended by polymerase. We describe how integrated gel band intensities can be used to obtain site-specific velocities for addition of correct and incorrect nucleotides, extending mismatched compared to correctly matched primer termini and measuring polymerase dissociation rates and equilibrium DNA binding constants. The analysis is based on steady-state "single completed hit conditions", where polymerases encounter many DNA molecules but where each DNA encounters an enzyme at most once. Specific topics addressed include nucleotide misinsertion, mismatch extension, exonucleolytic proofreading, single nucleotide discrimination using PCR, promiscuous mismatch extension by HIV-1 and AMV reverse transcriptases, sequence context effects on fidelity and polymerase dissociation, structural and kinetic properties of mispairs relating to fidelity, error avoidance mechanisms, kinetics of copying template lesions, the "A-rule" for insertion at abasic template lesions, an interesting exception to the "A-rule", thermodynamic and kinetic determinants of base pair discrimination by polymerases.

Randomization of genes by PCR mutagenesis

A modified polymerase chain reaction (PCR) was developed to introduce random point mutations into cloned genes. The modifications were made to decrease the fidelity of Taq polymerase during DNA synthesis without significantly decreasing the level of amplification achieved in the PCR. The resulting PCR products can be cloned to produce random mutant libraries or transcribed directly if a T7 promoter is incorporated within the appropriate PCR primer. We used this method to mutagenize the gene that encodes the Tetrahymena ribozyme with a mutation rate of 0.66% +/- 0.13% (95% C.I.) per position per PCR, as determined by sequence analysis. There are no strong preferneces with respect to the type of base substituion. The number of mutations per DNA sequence follows a Poisson distribution and the mutations are randomly distributed throughout the amplified sequence.

HeLa cell single-stranded DNA-binding protein increases the accuracy of DNA synthesis by DNA polymerase alpha in vitro

To determine whether cellular replication factors can influence the fidelity of DNA replication, the effect of HeLa cell single-stranded DNA-binding protein (SSB) on the accuracy of DNA replication by HeLa cell DNA polymerase alpha has been examined. An in vitro gap-filling assay, in which the single-stranded gap contains the supF target gene, was used to measure mutagenesis. Addition of SSB to the in vitro DNA synthesis reaction increased the accuracy of DNA polymerase alpha by 2- to 8-fold. Analysis of the products of DNA synthesis indicated that SSB reduces pausing by the polymerase at specific sites in the single-stranded supF template. Sequence analysis of the types of errors resulting from synthesis in the absence or presence of SSB reveals that, while the errors are primarily base substitutions under both conditions, SSB reduces the number of errors found at 3 hotspots in the supF gene. Thus, a cellular replication factor (SSB) can influence the fidelity of a mammalian DNA polymerase in vitro, suggesting that the high accuracy of cellular DNA replication may be determined in part by the interaction between replication factors, DNA polymerase and the DNA template in the replication complex.

Contribution of 3' leads to 5' exonuclease activity of DNA polymerase III holoenzyme from Escherichia coli to specificity

The effects of deoxynucleoside monophosphates on the 3' leads to 5' exonuclease activity of DNA polymerase III holoenzyme have been correlated with their effects on the fidelity of DNA replication. In particular, dGMP inhibits the proofreading activity of the enzyme and decreases the fidelity in those cases where a "following nucleotide effect" is also noted. This is strong evidence for proofreading. However, the absence of the effects of proofreading inhibitors or following nucleotides need not be evidence against the occurrence of proofreading: a theoretical analysis shows that these effects may not be observed even though there is active proofreading. This is suggested to be the case with the phage T4 enzyme system. The proofreading activity of Pol III appears to be directed primarily towards removing purine x pyrimidine-mediated rather than purine x purine-mediated misincorporations. recA protein inhibits the proofreading activity of Pol III on synthetic templates containing mismatched 3' termini. This is paralleled by a decrease in the fidelity of DNA replication in vitro. The inhibition is increased in the presence of dGMP or dAMP but there is no further increase in the infidelity of replication. The presence of both dNMPs and recA protein does not enable Pol III to copy past pyrimidine photodimers.

Fidelity of replication of bacteriophage phi X174 DNA in vitro and in vivo

Seven different revertants of bacteriophage phi X174am16 (AB5276G leads to T) have been isolated and the nature of the reversions determined by sequencing their DNA. The revertants each differ from am16 by just a single base substitution. These may be distinguished with varying degrees of ease by characteristic temperature sensitivities of growth. This has facilitated the determination of the frequency at which DNA polymerase III catalyses different types of substitution mutations in copying phi X174 DNA in vitro and in vivo. During the replicative form (RF) leads to single-stranded (SS) stage of replication in vitro, four different revertants may be readily produced according to well-defined rate laws on biasing the concentrations of dNTPs. Transversion mutations are found to be formed predominantly by purine x purine mismatching, whilst transitions are formed predominantly by G x T mismatching. The substitutions via G x T and G x A mismatches are estimated to occur at similar frequencies in vivo. The two most common revertants isolated in vivo, however, are not those readily produced during the RF leads to SS stage in vitro but are those produced on purine x purine mismatching in the SS leads to RF stage. The accuracy of the DNA polymerase in vitro appears to be similar to that in this stage in vivo. However, the overall accuracy of the RF leads to SS replication in vivo is more accurate than predicted from the measurements of the accuracy in vitro.

Accepted mutations in a gene family: evolutionary diversification of duplicated DNA

We report and compare the DNA sequences of 14 silkmoth (Antheraea polyphemus) chorion genes, derived from either cDNA or chromosomal DNA clones. Seven of these genes are members of the A multigene family, and seven are members of the B family. Where available, the previously reported (Jones and Kafatos 1980) intronic and extragenic flanking DNA sequences are also considered. Closely related sequences are compared, revealing the types of spontaneous mutations that were fixed during paralogous evolution. Segmental mutations (i.e. mutations other than substitutions) are nearly always interpretable as small duplications or deletions, related to small direct repeats. Segmental mutations are strongly constrained in the coding regions, although they do occur. Nucleotide substitutions also appear to be under selective constraints: relatively few substitutions leading to amino acid replacements are accepted, silent substitutions leading to some codons (especially purine-terminated ones) are disfavored, and different compositional biases are maintained in different parts of the sequences. Other sequence differences can be interpreted as indicative of neutral drift, including most differences in non-coding regions and most T/C transitions in third-base positions. In the non-coding regions, which are thought to be only loosely constrained by selection, transitions are observed more frequently than might be expected: they account for 52% of all substitutions, and they appear to be favored two to threefold over transversions when allowance is made for the skewed base composition of these regions.