Transitions,transversions, and the molecular evolutionary clock

Summary Nucleotide substitutions in the form of transitions (purine-purine or pyrimidine-pyrimidine interchanges) and transversions (purine-pyrimidine interchanges) occur during evolution and may be complied by aligning the sequences of homologous genes. Referring to the genetic code tables, silent transitions take place in third positions of codons in family boxes and two-codon sets. Silent transversions in third positions occur only in family boxes, except for AC transversions between AGR and CGR arginine codons (R=A or G). Comparisons of several protein genes have been made, and various subclasses of transitional and transversional nucleotide substitutions have been compiled. Considerable variations occur among the relative proportions of transitions and transversions. Such variations could possibly be caused by mutator genes, favoring either transitions or, conversely, transversions, during DNA replication. At earlier stages of evolutionary divergence, transitions are usually more frequent, but there are exceptions. No indication was found that transversions usually originate from multiple substitutions in transitions.     

Maxicircle DNA and edited mRNA sequences of closely related trypanosome species: implications of kRNA editing for evolution of maxicircle genomes.

kRNA editing produces functional mRNAs by uridine insertion and deletion. We analyzed portions of the apocytochrome b and NADH dehydrogenase subunits 7 and 8 (ND7 and 8) genes and their edited mRNAs in Trypanosoma congolense and compared these to the corresponding sequences in T.brucei. We find that these genes are highly diverged between the two species, especially in the positions of thymidines and in nucleotide transitions. Editing eliminates differences in encoded uridines producing edited mRNAs that are identical except for the nucleotide substitutions. The resulting predicted proteins are identical since all nucleotide substitutions are silent. A T.congolense minicircle-encoded gRNA which can specify editing of ND8 mRNA was identified. This gRNA can basepair with both T.congolense and T.brucei ND8 mRNA despite nucleotide transitions due to the flexibility of G:U base-pairing. These results illustrate how editing affects the characteristics of maxicircle sequence divergence and allows protein sequence conservation despite a level of DNA sequence divergence which would be predicted to be intolerable in the absence of editing.     

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

Human Mitochondrial DNA Variation and Evolution: Analysis of Nucleotide Sequences from Seven Individuals

We have analyzed nucleotide sequence variation in an approximately 900-base pair region of the human mitochondrial DNA molecule encompassing the heavy strand origin of replication and the D-loop. Our analysis has focused on nucleotide sequences available from seven humans. Average nucleotide diversity among the sequences is 1.7%, several-fold higher than estimates from restriction endonuclease site variation in mtDNA from these individuals and previously reported for other humans. This disparity is consistent with the rapidly evolving nature of this noncoding region. However, several instances of convergent or parallel gain and loss of restriction sites due to multiple substitutions were observed. In addition, other results suggest that restriction site (as well as pairwise sequence) comparisons may underestimate the total number of substitutions that have occurred since the divergence of two mtDNA sequences from a common ancestral sequence, even at low levels of divergence. This emphasizes the importance of recognizing the large standard errors associated with estimates of sequence variability, particularly when constructing phylogenies among closely related sequences. Analysis of the observed number and direction of substitutions revealed several significant biases, most notably a strand dependence of substitution type and a 32-fold bias favoring transitions over transversions. The results also revealed a significantly nonrandom distribution of nucleotide substitutions and sequence length variation. Significantly more multiple substitutions were observed than expected for these closely related sequences under the assumption of uniform rates of substitution. The bias for transitions has resulted in predominantly convergent or parallel changes among the observed multiple substitutions. There is no convincing evidence that recombination has contributed to the mtDNA sequence diversity we have observed.     

Estimation of DNA Sequence Context-dependent Mutation Rates Using Primate Genomic Sequences

It is understood that DNA and amino acid substitution rates are highly sequence context-dependent, e.g., C --> T substitutions in vertebrates may occur much more frequently at CpG sites and that cysteine substitution rates may depend on support of the context for participation in a disulfide bond. Furthermore, many applications rely on quantitative models of nucleotide or amino acid substitution, including phylogenetic inference and identification of amino acid sequence positions involved in functional specificity. We describe quantification of the context dependence of nucleotide substitution rates using baboon, chimpanzee, and human genomic sequence data generated by the NISC Comparative Sequencing Program. Relative mutation rates are reported for the 96 classes of mutations of the form 5' alphabetagamma 3' --> 5' alphadeltagamma 3', where alpha, beta, gamma, and delta are nucleotides and beta not equal delta, based on maximum likelihood calculations. Our results confirm that C --> T substitutions are enhanced at CpG sites compared with other transitions, relatively independent of the identity of the preceding nucleotide. While, as expected, transitions generally occur more frequently than transversions, we find that the most frequent transversions involve the C at CpG sites (CpG transversions) and that their rate is comparable to the rate of transitions at non-CpG sites. A four-class model of the rates of context-dependent evolution of primate DNA sequences, CpG transitions > non-CpG transitions approximately CpG transversions > non-CpG transversions, captures qualitative features of the mutation spectrum. We find that despite qualitative similarity of mutation rates among different genomic regions, there are statistically significant differences.     

Use of atp6 in fungal phylogenetics: an example from the boletales

Complete nucleotide sequences have been determined for atp6 from Suillus luteus and cox3 from Suillus sinuspaulianus (Boletales, Hymenomycetes, Basidiomycota), which code for ATPase subunit 6 and cytochrome oxidase subunit 3, respectively. These sequences were used to design PCR primers for the amplification of partial atp6 and cox3 sequences from other members of the Boletales and outgroup taxa. In atp6 and cox3 from Russula rosacea, one of the outgroup taxa, we observed a number of in-frame TGA(trp) codons, which imply a Neurospora crassa-type mitochondrial code in R. rosacea and possibly in basidiomycetes in general. Interestingly, however, most basidiomycetes other than R. rosacea appear to strongly prefer the TGG(trp) codon, which is unusual, given the strong A + T bias in fungal mitochondrial genomes. Pairwise comparisons were performed between atp6 sequences from increasingly divergent fungal lineages, and results show that all three codon positions become saturated in substitutions after an estimated divergence time of approx 300 Ma. This means that atp6 is likely to provide phylogenetic resolution within fungal classes but not at higher taxonomic levels. Also, because of the strong A + T bias in fungal mitochondrial genomes, A/T transversions were found to be more common than any other type of substitution, resulting in transversions being about two to three times more common in most pairwise sequence comparisons. Finally, atp6 sequences were used to infer phylogenetic relationships between 27 taxa from the Boletales and 4 outgroup taxa. Analyses were performed (i) on nucleotide sequence data using parsimony (successive approximation) as well as maximum likelihood methods and (ii) on deduced amino acid sequences using distance methods based on empirical substitution probabilities. Results from the various analyses are largely concordant with each other as well as with prior analyses of partial mitochondrial large-subunit rDNA (mtLSU rDNA). Analysis of the combined atp6 and mtLSU rDNA sequences results in increased bootstrap support for several key branches. Relationships that have been resolved for the first time in the current analysis are discussed.     

Mode and Tempo of Molecular Evolution in the Nematode Caenorhabditis: Cytochrome Oxidase II and Calmodulin Sequences

Through direct sequencing methods, the mitochondrial gene for cytochrome oxidase subunit two (CO II) and the single-copy nuclear gene for calmodulin were compared among strains of Caenorhabidits elegans and two other Caenorhabditis species (C. remanei and C. briggsae). In addition the CO II sequence was determined from a distantly related nematode, Steinernema intermedii. Among the 11 strains of C. elegans tested, there are four types of CO II gene, arising from two major lineages. Levels of intraspecific difference in the CO II gene are low (less than 2.0%) compared to the extraordinary divergence between congeneric species, which is about 50% when corrected for multiple hits. Concordant with the increase in divergence between taxa is a change in the pattern of substitution from a strong transition bias (24 transitions compared to two transversions) within species to a substitution pattern that appears to reflect the base composition of the mitochondrial genome when more divergent nematodes are compared. The base composition of the Caenorhabditis CO II gene is strongly biased toward A + T at all three positions of codons and appears to constrain the amino acid composition of the protein. Both the CO II and calmodulin genes show extreme conservation of amino acid sequences. When the accumulation of changes at silent sites in the two genes is compared among strains, it becomes evident that the mitochondrial gene is changing faster than the nuclear gene.     

Patterns of nucleotide change in mitochondrial ribosomal RNA genes and the phylogeny of piranhas

The patterns and rates of nucleotide substitution in mitochondrial ribosomal RNA genes are described and applied in a phylogenetic analysis of fishes of the subfamily Serrasalminae (Teleostei, Characiformes, Characidae). Fragments of 345 bp of the 12S and 535 bp of the 16S genes were sequenced for 37 taxa representing all but three genera in the subfamily. Secondary-structure models based on comparative sequence analysis were derived to characterize the pattern of change among paired and unpaired nucleotides, forming stem and loop regions, respectively. Base compositional biases were in the direction of A-rich loops and G-rich stems. Ninety-five percent of substitutions in stem regions were compensatory mutations, suggesting that selection for maintenance of base pairing is strong and that independence among characters cannot be assumed in phylogenetic analyses of stem characters. The relative rate of nucleotide substitution was similar in both fragments sequenced but higher in loop than in stem regions. In both genes, C-T transitions were the most common type of change, and overall transitions outnumbered transversions by a factor of two in 16S and four in 12S. Phylogenetic analysis of the mitochondrial DNA sequences suggests that a clade formed by the generaPiaractus, Colossoma, andMylossoma is the sister group to all other serrasalmins and that the generaMyleus, Serrasalmus, andPristobrycon are paraphyletic. A previous hypothesis concerning relationships for the serrasalmins, based on morphological evidence, is not supported by the molecular data. However, phylogenetic analysis of host-specific helminth parasites and cytogenetic data support the phylogeny of the Serrasalminae obtained in this study and provide evidence for coevolution between helminth parasites and their fish hosts.     

Methylglyoxal induces G:C to C:G and G:C to T:A transversions in the supF gene on a shuttle vector plasmid replicated in mammalian cells

We previously reported that the majority of base-pair substitutions induced by an endogenous mutagen, methylglyoxal, were G:C-->T:A transversions and G:C-->A:T transitions in wild-type and nucleotide excision repair (NER)-deficient (uvrA or uvrC) Escherichia coli strains. To investigate the mutation spectrum of methylglyoxal in mammalian cells and to compare the spectrum with those detected in other experimental systems, we analyzed mutations in a bacterial suppressor tRNA (supF) gene in the shuttle vector plasmid pMY189. We treated pMY189 with methylglyoxal and immediately transfected it into simian COS-7 cells. The cytotoxicity and the mutation frequency (MF) increased according to the dose of methylglyoxal. In the mutants induced by methylglyoxal, multi-base deletions were predominant (50%), followed by base-pair substitutions (35%), in which 89% of the substitutions occurred at G:C sites. Among them, G:C-->C:G and G:C-->T:A transversions were predominant. The overall distribution of methylglyoxal-induced mutations detected in the supF gene was different from that for the spontaneous mutations. These results suggest that methylglyoxal may take part in causing G:C-->C:G and G:C-->T:A transversions in vivo.     

Rates of molecular evolution and the fraction of nucleotide positions free to vary

Summary Selective constraints on DNA sequence change were incorporated into a model of DNA divergence by restricting substitutions to a subset of nucleotide positions. A simple model showed that both mutation rate and the fraction of nucleotide positions free to vary are strong determinants of DNA divergence over time.When divergence between two species approaches the fraction of positions free to vary, standard methods that correct for multiple mutations yield severe underestimates of the number of substitutions per site. A modified method appropriate for use with DNA sequence, restriction site, or thermal renaturation data is derived taking this fraction into account. The model also showed that the ratio of divergence in two gene classes (e.g., nuclear and mitochondrial) may vary widely over time even if the ratio of mutation rates remains constant.DNA sequence divergence data are used increasingly to detect differences in rates of molecular evolution. Often, variation in divergence rate is assumed to represent variation in mutation rate. The present model suggests that differing divergence rates among comparisons (either among gene classes or taxa) should be interpreted cautiously. Differences in the fraction of nucleotide positions free to vary can serve as an important alternative hypothesis to explain differences in DNA divergence rates.     

Mitochondrial DNA sequences of primates: Tempo and mode of evolution

Summary We cloned and sequenced a segment of mitochondrial DNA from human, chimpanzee, gorilla, orangutan, and gibbon. This segment is 896 bp in length, contains the genes for three transfer RNAs and parts of two proteins, and is homologous in all 5 primates. The 5 sequences differ from one another by base substitutions at 283 positions and by a deletion of one base pair. The sequence differences range from 9 to 19% among species, in agreement with estimates from cleavage map comparisons, thus confirming that the rate of mtDNA evolution in primates is 5 to 10 times higher than in nuclear DNA. The most striking new finding to emerge from these comparisons is that transitions greatly outnumber transversions. Ninety-two percent of the differences among the most closely related species (human, chimpanzee, and gorilla) are transitions. For pairs of species with longer divergence times, the observed percentage of transitions falls until, in the case of comparisons between primates and non-primates, it reaches a value of 45. The time dependence is probably due to obliteration of the record of transitions by multiple substitutions at the same nucleotide site. This finding illustrates the importance of choosing closely related species for analysis of the evolutionary process. The remarkable bias toward transitions in mtDNA evolution necessitates the revision of equations that correct for multiple substitutions at the same site. With revised equations, we calculated the incidence of silent and replacement substitutions in the two protein-coding genes. The silent substitution rate is 4 to 6 times higher than the replacement rate, indicating strong functional constraints at replacement sites. Moreover, the silent rate for these two genes is about 10% per million years, a value 10 times higher than the silent rate for the nuclear genes studied so far. In addition, the mean substitution rate in the three mitochondrial tRNA genes is at least 100 times higher than in nuclear tRNA genes. Finally, genealogical analysis of the sequence differences supports the view that the human lineage branched off only slightly before the gorilla and chimpanzee lineages diverged and strengthens the hypothesis that humans are more related to gorillas and chimpanzees than is the orangutan.Abbreviations mtDNA mitochondrial DNA - bp base pair - URF unidentified reading frame     

Evolutionary genetics of the suiformes as reconstructed using mtDNA sequencing

We have amplified and sequnced the entire mitochondrial DNA cytochromeb gene from four species of Suidae: babirusa, warthog, bearded pig, and some specimens belonging to different subspecies and populations of wild and domestic pigs (Sus scrofa). These sequences were aligned with additional mammalian sequences retrieved from the literature and were used to obtain phylogenetic trees of the Suiformes (Artiodactyla). Several species of Carnivora, Perissodactyla. Cetacea, and other Artiodactyla were used as outgroups. Molecular phylogenetic relationships among the Suiformes reflect their current taxonomy: Hippopotamidae, Tayassuidae, and Suidae are separated by deep genetic gaps, and the division of the Suidae into the subfamilies Babyrousinae., Phacochoerinae, and Suinae has strong genetic correlates. Cytochromeb sequences show differences among Asian and Western populations ofSus scrofa, agreeing with other genetic information (karyotypes blood groups, and protein variability). The two Italian subspecies of wild boar have unique mtDNA cytochromeb haplotypes. The evolutionary rates of cytochromeb sequences are different at transitions versus transversions as well as at first, second, and third positions of codons. Therefore, these classes of substitutions reached different levels of mutational saturation. Only transversions and the conservative first and second position substitutions are linearly related to genetic distances among the Suiformes. Therefore, divergence times were computed using unsaturated conserved nucleotide substitutions and calibrated using paleontological divergence times between some Artiodactyla. Transversions apparently evolve at remarkably regular rates in ungulate taxa which have accumulated less than 20% estimated sequence divergence, corresponding to about 40–45 million years of independent evolution. Molecular, information suggests that Hippopotamidae and Tayassuidae are not closely related (as stated by Pickford, 1986, 1989, 1993) and that the origin of babirusa and warthog (about 10–19 and 5–15 million years ago, respectively) is more recent than supported by current evolutionary reconstructions. The inferred origin of bearded pig is about 2.1 million years old, and genetic divergence among differentSus scrofa populations is probably a Pleistocene event. The addition of new sequences of Suiformes does not help in resolving the phylogenetic position ofHippopotamus amphibius, which shows weak but recurrent linkages with the cetacean evolutionary lineage.To whom correspondence should be addressed.     

The Mitochondrial Genome of the Honeybee Apis Mellifera: Complete Sequence and Genome Organization

The complete sequence of honeybee (Apis mellifera) mitochondrial DNA is reported being 16,343 bp long in the strain sequenced. Relative to their positions in the Drosophila map, 11 of the tRNA genes are in altered positions, but the other genes and regions are in the same relative positions. Comparisons of the predicted protein sequences indicate that the honeybee mitochondrial genetic code is the same as that for Drosophila; but the anticodons of two tRNAs differ between these two insects. The base composition shows extreme bias, being 84.9% AT (cf. 78.6% in Drosophila yakuba). In protein-encoding genes, the AT bias is strongest at the third codon positions (which in some cases lack guanines altogether), and least in second codon positions. Multiple stepwise regression analysis of the predicted products of the protein-encoding genes shows a significant association between the numbers of occurrences of amino acids and %T in codon family, but not with the number of codons per codon family or other parameters associated with codon family base composition. Differences in amino acid abundances are apparent between the predicted Apis and Drosophila proteins, with a relative abundance in the Apis proteins of lysine and a relative deficiency of alanine. Drosophila alanine residues are as often replaced by serine as conserved in Apis. The differences in abundances between Drosophila and Apis are associated with %AT in the codon families, and the degree of divergence in amino acid composition between proteins correlates with the divergence in %AT at the second codon positions. Overall, transversions are about twice as abundant as transitions when comparing Drosophila and Apis protein-encoding genes, but this ratio varies between codon positions. Marked excesses of transitions over chance expectation are seen for the third positions of protein-coding genes and for the gene for the small subunit of ribosomal RNA. For the third codon positions the excess of transitions is adequately explained as due to the restriction of observable substitutions to transitions for conserved amino acids with two-codon families; the excess of transitions over expectation for the small ribosomal subunit suggests that the conservation of nucleotide size is favored by selection.     

Sequence Evolution of Drosophila Mitochondrial DNA

We have compared nucleotide sequences of corresponding segments of the mitochondrial DNA (mtDNA) molecules of Drosophila yakuba and Drosophila melanogaster, which contain the genes for six proteins and seven tRNAs. The overall frequency of substitution between the nucleotide sequences of these protein genes is 7.2%. As was found for mtDNAs from closely related mammals, most substitutions (86%) in Drosophila mitochondrial protein genes do not result in an amino acid replacement. However, the frequencies of transitions and transversions are approximately equal in Drosophila mtDNAs, which is in contrast to the vast excess of transitions over transversions in mammalian mtDNAs. In Drosophila mtDNAs the frequency of C----T substitutions per codon in the third position is 2.5 times greater among codons of two-codon families than among codons of four-codon families; this is contrary to the hypothesis that third position silent substitutions are neutral in regard to selection. In the third position of codons of four-codon families transversions are 4.6 times more frequent than transitions and A----T substitutions account for 86% of all transversions. Ninety-four percent of all codons in the Drosophila mtDNA segments analyzed end in A or T. However, as this alone cannot account for the observed high frequency of A----T substitutions there must be either a disproportionately high rate of A----T mutation in Drosophila mtDNA or selection bias for the products of A----T mutation. --Consideration of the frequencies of interchange of AGA and AGT codons in the corresponding D. yakuba and D. melanogaster mitochondrial protein genes provides strong support for the view that AGA specifies serine in the Drosophila mitochondrial genetic code.     

Evidence from mutation spectra that the UV hypermutability of xeroderma pigmentosum variant cells reflects abnormal, error-prone replication on a template containing photoproducts.

Xeroderma pigmentosum (XP) variant patients are genetically predisposed to sunlight-induced skin cancer. Fibroblasts derived from these patients are extremely sensitive to the mutagenic effect of UV radiation and are abnormally slow in replicating DNA containing UV-induced photoproducts. However, unlike cells from the majority of XP patients, XP variant cells have a normal or nearly normal rate of nucleotide excision repair of such damage. To determine whether their UV hypermutability reflected a slower rate of excision of photoproducts specifically during early S phase when the target gene for mutations, i.e., the hypoxanthine (guanine) phosphoribosyltransferase gene (HPRT), is replicated, we synchronized diploid populations of normal and XP variant fibroblasts, irradiated them in early S phase, and compared the rate of loss of cyclobutane pyrimidine dimers and 6-4 pyrimidine-pyrimidones from DNA during S phase. There was no difference. Both removed 94% of the 6-4 pyrimidine-pyrimidones within 8 h and 40% of the dimers within 11 h. There was also no difference between the two cell lines in the rate of repair during G1 phase. To determine whether the hypermutability resulted from abnormal error-prone replication of DNA containing photoproducts, we determined the spectra of mutations induced in the coding region of the HPRT gene of XP variant cells irradiated in early S and G1 phases and compared with those found in normal cells. The majority of the mutations in both types of cells were base substitutions, but the two types of cells differed significantly from each other in the kinds of substitutions, but the two types differed significantly from each other in the kinds of substitutions observed either in mutants from S phase (P < 0.01) or from G1 phase (P = 0.03). In the variant cells, the substitutions were mainly transversions (58% in S, 73% in G1). In the normal cells irradiated in S, the majority of the substitutions were G.C --> A.T, and most involved CC photoproducts in the transcribed strand. In the variant cells irradiated in S, substitutions involving cytosine in the transcribed strand were G.C --> T.A transversions exclusively. G.C --> A.T transitions made up a much smaller fraction of the substitutions than in normal cells (P < 0.02), and all of them involved photoproducts located in the nontranscribed strand. The data strongly suggest that XP variant cells are much less likely than normal cells to incorporate either dAMP or dGMP opposite the pyrimidines involved in photoproducts. This would account for their significantly higher frequency of mutants and might explain their abnormal delay in replicating a UV-damaged template.     

Sequence evolution and phylogenetic signal in control-region and cytochrome b sequences of rainbow fishes (Melanotaeniidae)

The nucleotide sequences of segments of the cytochrome b gene (351 bp), the tRNA(Pro) gene (49 bp), and the control region (approximately 313 bp) of mitochondrial DNA were obtained from 26 fish representing different populations and species of Melanotaenia and one species of Glossolepis, freshwater rainbow fishes confined to Australia and New Guinea. The purpose was to investigate relative rates and patterns of sequence evolution. Overall levels of divergence were similar for the cytochrome b and tRNA control-region sequences, both ranging from < 1% within subspecies to 15%-19% between genera. However, the patterns of sequence evolution differed. For the cytochrome b gene, transitions consistently exceeded transversions, the bias ranging from 4.2:1 to 2:1, depending on the level of sequence divergence. However, in the control-region sequence, a bias toward transitions (2:1) was observed only in comparisons between very similar sequences, and transversions outnumbered transitions in comparisons of divergent sequences. Graphic comparisons suggested that the control region was saturated for transitions at relatively low levels of sequence divergence but accumulated transversions at a greater rate than did the cytochrome b sequence. These distinct patterns of base substitution are associated with differences in A+T content, which is 70% for the tRNA control- region segment versus 50% for cytochrome b. A test for skewness in the distribution of lengths of random trees indicated that both segments contained phylogenetic signal. Parsimony analyses of the data from the two regions, with or without weighting schemes appropriate to the respective patterns of sequence evolution, identified the same five groupings of sequences, but the relationships among the groups differed. However, in most cases the branches uniting different combinations of groups were poorly supported, and the differences among topologies were insignificant. Considering the observed patterns of base substitution and the results of the phylogenetic analyses, we deduce that both the control region and cytochrome b are appropriate for population genetic studies but that the control region is less effective than cytochrome b for resolving relationships among divergent lineages of rainbow fishes.      

Compensatory substitutions and the evolution of the mitochondrial 12S rRNA gene in mammals

12S ribosomal RNA (rRNA) gene sequences from a suite of mammalian taxa (13 placentals, 4 marsupials, 1 monotreme), for which phylogenetic relationships are well established based on independent criteria, were employed to study the evolution of this gene. Phylogenetic analysis of 12S sequences produces a phylogeny that agrees with expectations. Base composition provides evidence for directional symmetrical substitution pressure in loops; in stems, base composition is much more even. Rates of nucleotide substitution are lower in stems than loops. Patterns of nucleotide substitution show an overall preference for transitions over transversions, with this difference more profound in stems than loops. Among different transversion pathways, there is a wide range of transformation frequencies. An analysis of compensatory substitutions shows that there is strong evidence for their occurrence and that a weighting factor of 0.61 should be applied in phylogenetic analyses to account for the dependence of mutations at stem positions relative to positions where changes are independent. Among stem variables (i.e., stem length, interaction distance, substitution rates, G+C content, and the percentage of bases that are paired), several significant correlations were discovered, but stem length and interaction distance are uncorrelated with other variables.      

Deficient nucleotide excision repair increases base-pair substitutions but decreases TGGC frameshifts induced by methylglyoxal in Escherichia coli.

To investigate the mutation spectrum of a well-known mutagen, methylglyoxal, and the influence of nucleotide excision repair (NER) on methylglyoxal-induced mutations, we treated wild-type and NER-deficient (uvrA or uvrC) Escherichia coli strains with methylglyoxal, and analyzed mutations in the chromosomal lacI gene. In the three strains, the cell death and the mutation frequency increased according to the dose of methylglyoxal added to the culture medium. The frequencies of methylglyoxal-induced base-pair substitutions were higher in the NER-deficient strains than in the wild-type strain, in the presence and absence of mucAB gene. Paradoxically, the frequency of methylglyoxal-induced TGGC frameshifts was higher in the wild-type strain than in the NER-deficient strains. When the methylglyoxal-induced mutation spectra in the presence and absence of mucAB gene are compared, the ratios of base-pair substitutions to frameshifts were increased by the effects of mucAB gene. In the three strains, more than 75% of the base-pair substitutions occurred at G:C sites, independent of the mucAB gene. When the mucAB gene was present, G:C-->T:A transversions were predominant, followed by G:C-->A:T transitions. When the mucAB gene was absent, the predominant mutations differed in the three strains: in the wild-type and uvrC strains, G:C-->A:T transitions were predominant, followed by G:C-->T:A transversions, while in the uvrA strains, G:C-->T:A transversions were predominant, followed by G:C-->A:T transitions. These results suggest that NER may be involved in both the repair and the fixation of methylglyoxal-induced mutations.     

Rapid Rate of Control-Region Evolution in Pacific Butterflyfishes (Chaetodontidae)

Sequence differences in the tRNA-proline (tRNA^pro) end of the mitochondrial control-region of three species of Pacific butterflyfishes accumulated 33–43 times more rapidly than did changes within the mitochondrial cytochrome b gene (cytb). Rapid evolution in this region was accompanied by strong transition/transversion bias and large variation in the probability of a DNA substitution among sites. These substitution constraints placed an absolute ceiling on the magnitude of sequence divergence that could be detected between individuals. This divergence ``ceiling' was reached rapidly and led to a decay in the relative rate of control-region/cytb b evolution. A high rate of evolution in this section of the control-region of butterflyfishes stands in marked contrast to the patterns reported in some other fish lineages. Although the mechanism underlying rate variation remains unclear, all taxa with rapid evolution in the 5′-end of the control-region showed extreme transition biases. By contrast, in taxa with slower control-region evolution, transitions accumulated at nearly the same rate as transversions. More information is needed to understand the relationship between nucleotide bias and the rate of evolution in the 5′-end of the control-region. Despite strong constraints on sequence change, phylogenetic information was preserved in the group of recently differentiated species and supported the clustering of sequences into three major mtDNA groupings. Within these groups, very similar control-region sequences were widely distributed across the Pacific Ocean and were shared between recognized species, indicating a lack of mitochondrial sequence monophyly among species. Received: 30 June 1996 / Accepted: 15 May 1997     

Novel mutational spectrum induced by N-methyl-N'-nitro-N-nitrosoguanidine in the coding region of the hypoxanthine (guanine) phosphoribosyltransferase gene in diploid human fibroblasts

The kinds and locations of mutations in the coding region of the hypoxanthine (guanine) phosphoribosyltransferase (hprt) gene of 75 independent mutants, derived from N-methyl-N'-nitro-N-nitrosoguanidine (MNNG)-treated normal human fibroblasts, were characterized by direct sequencing of mRNA-polymerase chain reaction (mRNA-PCR)-amplified cDNA. Treatment of human cells with low (6 or 8 microM) or high (10 or 12 microM) doses of MNNG resulted in 35-fold or 150-fold average increases in mutation frequency, respectively. A high frequency of mutants lacking a complete exon was observed in both groups. Further characterization of half of these mutants by DNA-PCR amplification of intron-exon boundaries showed that they contained base substitutions. The kinds of base substitutions differed distinctly between these two groups. In the low dose group, a broad mutational spectrum was observed: ten out of the 31 base substitutions were A.T to G.C transitions, six contained G.C to A.T transitions, and the other 15 exhibited transversions. In contrast, the majority (84%) of base substitutions among the high dose group were G.C to A.T transitions; the others (16%) were transversions. All of the 32 G.C to A.T transitions were located on the non-transcribed strand, assuming that the causative premutational lesion was O6-methylguanine. These results indicate preferential repair of lesions located on the transcribed strand. In addition, G.C to A.T and A.T to G.C transitions preferentially occurred at positions with guanine and thymine at the adjacent 5' position, respectively.