Fourfold Degenerate Sites Research Articles

The genome of Campylobacter jejuni: codon and amino acid usage.

The genes from the genome of the AT-rich bacterium Campylobacter jejuni were analysed and characterised with respect to usage and amino acid usage. Codon usage is generally biased for all amino acids having synonymous codons, so that AT-rich synonyms are most frequently used. Markov chain analysis showed that codon bias and over- or underrepresentation of the corresponding tri-letter words are not related. Predicted secondary structure, lipophilicity, codon position within the gene, strand, and position on the (+)-strand were all shown to be determinants of codon usage, and these effects were in part directly explained by compositional phenomena. Codon context and the GC-content at the wobble position of the fourfold degenerate sites exert indirect effects on codon usage. The factors that affect codon usage seem to affect all amino acids, rather than selected amino acids. The usage of amino acids correlates well with the GC-content of genes, i.e. usage of amino acids encoded by GC-rich codons increases with GC-content and vice versa.

Base compositions of genes encoding alpha-actin and lactate dehydrogenase-A from differently adapted vertebrates show no temperature-adaptive variation in G + C content.

There is a long-standing debate in molecular evolution concerning the putative importance of GC content in adapting the thermal stabilities of DNA and RNA. Most studies of this relationship have examined broad-scale compositional patterns, for example, total GC percentages in genomes and occurrence of GC-rich isochores. Few studies have systematically examined the GC contents of individual orthologous genes from differently thermally adapted species. When this has been done, the emphasis has been on comparing large numbers of genes in only a few species. We have approached the GC-adaptation temperature hypothesis in a different manner by examining patterns of base composition of genes encoding lactate dehydrogenase-A (ldh-a) and alpha-actin (alpha-actin) from 51 species of vertebrates whose adaptation temperatures ranged from -1.86 degrees C (Antarctic fishes) to approximately 45 degrees C (desert reptile). No significant positive correlation was found between any index of GC content (GC content of the entire sequence, GC content of the third codon position [GC(3)], and GC content at fourfold degenerate sites [GC(4)]) and any index of adaptation temperature (maximal, mean, or minimal body temperature). For alpha-actin, slopes of regression lines for all comparisons did not differ significantly from zero. For ldh-a, negative correlations between adaptation temperature and total GC content, GC(3), and GC(4) were observed but were shown to be due entirely to phylogenetic influences (as revealed by independent contrast analyses). This comparison of GC content across a wide range of ectothermic ("cold-blooded") and endothermic ("warm-blooded") vertebrates revealed that frogs of the genus Xenopus, which have commonly been used as a representative cold-blooded species, in fact are outliers among ectotherms for the alpha-actin analyses, raising concern about the appropriateness of choosing these amphibians as representative of ectothermic vertebrates in general. Our study indicates that, whereas GC contents of isochores may show variation among different classes of vertebrates, there is no consistent relationship between adaptation temperature and the percentage of thermal stability-enhancing G + C base pairs in protein-coding genes.

Evolution of the A+T-rich region of mitochondrial DNA in the melanogaster species subgroup of Drosophila.

We determined the nucleotide sequences of two regions in the A+T-rich region of mitochondrial DNA (mtDNA) in the siI and siII types of D. simulans, the maII type of D. mauritiana, and D. sechellia. The sequences were aligned with those of the corresponding regions of siIII of D. simulans and maI of D. mauritiana, D. melanogaster, and D. yakuba. The type I and type II elements and the T-stretches were detected in all eight of the mtDNA types compared, indicating that the three elements are essential in the A+T-rich region of this species subgroup. The alignment revealed several short repetitive sequences and relatively large deletions in the central portions of the region. In the highly conserved sequence elements in the type II elements, the substitution rates were not uniform among lineages and acceleration in the substitution rate might have been due to loss of functional constraint in the stem-loop-forming sequences predicted in the type II elements. Patterns of nucleotide substitutions observed in the A+T-rich region were further compared with those in the coding regions and in the intergenic regions of mtDNA. Substitutions between A and T were particularly repressed in the highly conserved sequence elements and in the intergenic regions compared with those in the A+T-rich region excluding the highly conserved sequence elements and in the fourfold degenerate sites in the coding regions. The functional and structural characteristics of the A+T-rich region that might be involved in this substitutional bias are discussed.

Measuring genome divergence in bacteria: a case study using chlamydian data.

We have studied the relative contribution of inversions, transpositions, deletions, and nucleotide substitutions to the evolution of Chlamydia trachomatis and Chlamydia pneumoniae. The minimal number of rearrangement events required for converting the gene order structure of one genome into that of the other was estimated to 59 +/- 6 events, including 13% inversions, 38% short inversions, and 49% transpositions. In contrast to previous findings, no examples of horizontal gene transfer subsequent to species divergence were identified, nor any evidence for an excessive number of tandem gene duplications. A statistical model was used to compare nucleotide frequencies for a set of genes uniquely present in one species to a set of orthologous genes present in both species. The two data sets were not significantly different, which is indicative of a low frequency of horizontal gene transfer events. This is based on the assumption that a foreign gene of different nucleotide content will not have become completely ameliorated, as verified by simulations of the amelioration rate at twofold and fourfold degenerate codon sites. The frequencies of nucleotide substitutions at twofold and fourfold degenerate sites, deletions, inversions, and translocations were estimated to 1.42, 0.62, 0.18, 0.01, and 0.01 per site, respectively.

On the estimation of the rate of nucleotide substitution for the control region of human mitochondrial DNA

To apply molecular clock for studying human evolution, the pattern of nucleotide substitution for the control region of human mtDNA was analyzed in detail. It is well known that the rate of nucleotide substitution for the control region is much higher than that for any other part of mtDNA. In this study, the higher substitution rate was attributed to the higher rate of transition-type substitution between pyrimidines within the D-loop part, whereas the rates of other types of substitution were essentially the same over the entire mtDNA molecule. Even within the control region, the rate and pattern of nucleotide substitution were different between the D-loop part and the rest. The rate and pattern for the non-D-loop part were very similar to those for fourfold-degenerate sites in the protein-coding region. In contrast, the D-loop and non-D-loop parts showed similarities in the base composition, whereas the base composition of fourfold-degenerate sites slightly different from that of the both parts of the control region. It is concluded, therefore, that the nucleotide frequencies of the control region should be used to estimate the number of substitutions (d) between the control region sequences. However, a method to verify the accuracy of the estimation of d by means of the transition/transversion (s/v) ratio was theoretically studied. It was suggested that the s/v ratio becomes constant over a wide range of d values only when the estimation of d is unbiased. On the basis of this result, the estimates of d previously obtained between human sequences were evaluated.

The codon-degeneracy model of molecular evolution.

Mitochondrial genetic codons can be categorized by four patterns of nucleotide-site degeneracy based on varying combinations of twofold- or nondegenerate sites at first codon positions and twofold- or fourfold-degenerate sites at third codon positions. Herein, a model of molecular evolution is introduced that uses these patterns to calculate expected substitution frequencies for each codon position and substitution type relative to overall number of synonymous or nonsynonymous substitutions. Regions of the pocket gopher cytochrome oxidase subunit I (COI) and cytochrome b (cyt-b) genes are analyzed using this model. Chi-square distributions are used to produce relative goodness-of-fit (GF) scores for measuring the difference between substitution frequencies predicted by the codon-degeneracy model (CDM), and frequencies inferred using a well-supported phylogenetic tree of closely related species. The GF scores for expected and observed synonymous (GF(syn) = 0.429, p = 0.807) and nonsynonymous (GF(ns) = 2.309, p = 0.679) substitution frequencies resulted in a failure to reject the CDM as a null hypothesis for the molecular evolution of COI and cyt-b in pocket gophers. Alternative tree topologies and calculations of transition bias for these data result in higher GF scores.

P3B39 - In vivo covalent binding of acetaminophen to canalicular and basolateral membrane proteins in mouse liver

Fusarium graminearum and 21 related species comprising the F. sambucinum species complex lineage 1 (FSAMSC-1) are the most important Fusarium Head Blight pathogens of cereal crops world-wide. FSAMSC-1 species typically produce type B trichothecenes. However, some F. graminearum strains were recently found to produce a novel type A trichothecene (NX-2) resulting from functional variation in the trichothecene biosynthetic enzyme Tri1. We used a PCR-RFLP assay targeting the TRI1 gene to identify the NX-2 allele among a global collection of 2515 F. graminearum. NX-2 isolates were only found in southern Canada and the northern U.S., where they were observed at low frequency (1.8%), but over a broader geographic range and set of cereal hosts than previously recognized. Phylogenetic analyses of TRI1 and adjacent genes produced gene trees that were incongruent with the history of species divergence within FSAMSC-1, indicating trans-species evolution of ancestral polymorphism. In addition, placement of NX-2 strains in the TRI1 gene tree was influenced by the accumulation of nonsynonymous substitutions associated with the evolution of the NX-2 chemotype, and a significant (P < 0.001) change in selection pressure was observed along the NX-2 branch (ω = 1.16) in comparison to other branches (ω = 0.17) in the TRI1 phylogeny. Parameter estimates were consistent with positive selection for specific amino-acid changes during the evolution of NX-2, but direct tests of positive selection were not significant. Phylogenetic analyses of fourfold degenerate sites and intron sequences in TRI1 indicated the NX-2 chemotype had a single evolutionary origin and evolved recently from a type B ancestor. Our results indicate the NX-2 chemotype may be indigenous, and possibly endemic, to southern Canada and the northern U.S. In addition, we demonstrate that the evolution of TRI1 within FSAMSC-1 has been complex, with evidence of trans-species evolution and chemotype-specific shifts in selective constraint.

Molecular evidence from the nuclear genome for the time frame of human evolution.

Evolutionary divergence times can be inferred from molecular distances if a molecular clock can be assumed and if the substitution rate can be estimated. We present new evidence from relative rate tests that the rate of substitution at fourfold degenerate sites of nuclear genome-coding DNA is uniform in primate and rodent lineages. We also review recent relative rate test results showing substitution rate uniformity in the nuclear genome of simian primates. DNA distances between a range of mammalian taxa shows that a molecular clock is inconsistent with many assumed divergence times irrespective of the assumed substitution rate. We find that the substitution rate that implies the best compromise fit with divergence times across the range of taxa is 2.0-2.25 x 10(-9). This range of substitution rates implies a divergence time of humans and chimpanzees of 4.0-3.6 million years ago. This postdates the occurrence of Ardipithecus ramidus and the earliest occurrence of Australopithecus afarensis, suggesting that the common ancestor of humans and chimpanzees was bipedal and that the trait has been lost in chimpanzees rather than gained in humans.

Evolution of eutherian cytochrome c oxidase subunit II: heterogeneous rates of protein evolution and altered interaction with cytochrome c.

Cytochrome c oxidase subunit II (COII), encoded by the mitochondrial genome, exhibits one of the most heterogeneous rates of amino acid replacement among placental mammals. Moreover, it has been demonstrated that cytochrome c oxidase has undergone a structural change in higher primates which has altered its physical interaction with cytochrome c. We collected a large data set of COII sequences from several orders of mammals with emphasis on primates, rodents, and artiodactyls. Using phylogenetic hypotheses based on data independent of the COII gene, we demonstrated that an increased number of amino acid replacements are concentrated among higher primates. Incorporating approximate divergence dates derived from the fossil record, we find that most of the change occurred independently along the New World monkey lineage and in a rapid burst before apes and Old World monkeys diverged. There is some evidence that Old World monkeys have undergone a faster rate of nonsynonymous substitution than have apes. Rates of substitution at four-fold degenerate sites in primates are relatively homogeneous, indicating that the rate heterogeneity is restricted to nondegenerate sites. Excluding the rate acceleration mentioned above, primates, rodents, and artiodactyls have remarkably similar nonsynonymous replacement rates. A different pattern is observed for transversions at four-fold degenerate sites, for which rodents exhibit a higher rate of replacement than do primates and artiodactyls. Finally, we hypothesize specific amino acid replacements which may account for much of the structural difference in cytochrome c oxidase between higher primates and other mammals.

On transition bias in mitochondrial genes of pocket gophers.

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (micros Ps)/(microv Pv), where micros and microv denote mutation rate of transitions and transversions, respectively, and Ps and Pv denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions, Ps/Pv approximately 1 at fourfold degenerate sites, so that the s/v ratio at fourfold degenerate sites is approximately micros/microv, which is a measure of mutational contribution to transition bias. Similarly, the s/v ratio at nondegenerate sites is also an estimate of micros/microv if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochrome b (cyt-b) in pocket gophers, the s/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, the s/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the micros/microv ratio is about 2, the Ps/Pv ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content on s/v ratio of protein-coding genes.

Intraspecific nuclear DNA variation in Drosophila.

We have summarized and analyzed all available nuclear DNA sequence polymorphism studies for three species of Drosophila, D. melanogaster (24 loci), D. simulans (12 loci), and D. pseudoobscura (5 loci). Our major findings are: (1) The average nucleotide heterozygosity ranges from about 0.4% to 2% depending upon species and function of the region, i.e., coding or noncoding. (2) Compared to D. simulans and D. pseudoobscura (which are about equally variable), D. melanogaster displays a low degree of DNA polymorphism. (3) Noncoding introns and 3' and 5' flanking DNA shows less polymorphism than silent sites within coding DNA. (4) X-linked genes are less variable than autosomal genes. (5) Transition (Ts) and transversion (Tv) polymorphisms are about equally frequent in non-coding DNA and at fourfold degenerate sites in coding DNA while Ts polymorphisms outnumber Tv polymorphisms by about 2:1 in total coding DNA. The increased Ts polymorphism in coding regions is likely due to the structure of the genetic code: silent changes are more often Ts's than are replacement substitutions. (6) The proportion of replacement polymorphisms is significantly higher in D. melanogaster than in D. simulans. (7) The level of variation in coding DNA and the adjacent noncoding DNA is significantly correlated indicating regional effects, most notably recombination. (8) Surprisingly, the level of polymorphism at silent coding sites in D. melanogaster is positively correlated with degree of codon usage bias. (9) Three proposed tests of the neutral theory of DNA polymorphisms have been performed on the data: Tajima's test, the HKA test, and the McDonald-Kreitman test. About half of the loci fail to conform to the expectations of neutral theory by one of the tests. We conclude that many variables are affecting levels of DNA polymorphism in Drosophila, from properties of nucleotides to population history and, perhaps, mating structure. No simple, all encompassing explanation satisfactorily accounts for the data.

Tempo and mode of synonymous substitutions in mitochondrial DNA of primates.

Nucleotide substitutions of the four-fold degenerate sites and the total third codon positions of mitochondrial DNA from human, common chimpanzee, bonobo, gorilla, and orangutan were examined in detail by three alternative Markov models; (1) Hasegawa, Kishino, and Yano's (1985) model, (2) Tamura and Nei's (1993) model, and (3) the general reversible Markov model. These sites are expected to be relatively free from constraint, and therefore their tempo and mode in evolution should reflect those of mutation. It turned out that, among the alternative models, the general reversible Markov model best approximates the nucleotide substitutions of the four-fold degenerate sites and the total third codon positions, while the maximum likelihood estimates of the numbers of nucleotide substitutions along each branch do not differ significantly among the three models. It was further shown that the transition rate of these sites during evolution, and therefore transitional mutation rate of mtDNA, are higher in humans than in chimpanzees and gorillas probably by about two times. However, transversional mutation rate and amino acid substitution rate do not differ significantly between humans and the African apes. These and additional observations suggest heterogeneity of the mutation rate as well as of the constraint operating on the mtDNA-encoded proteins among different lineages of Hominoidea.

Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes

Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes

Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes

Three statistics (%GC, GC-skew, and AT-skew) can be used to describe the overall patterns of nucleotide composition in DNA sequences. Fourfold degenerate third codon positions from 16 animal mitochondrial genomes were analyzed. The overall composition, as measured by %GC, varies from 3.6 %GC in the honeybee to 47.2 %GC in human mtDNA. Compositional differences between strands of the mitochondrial genome were quantified using the two skew statistics presented in this paper. Strand-specific distribution of bases varies among animal taxa independently of overall %GC. Compositional patterns reflect the substitution process. Description of these patterns may aid in the formation of hypotheses about substitutional mechanisms.

Disparate rates of molecular evolution in cospeciating hosts and parasites.

DNA sequences for the gene encoding mitochondrial cytochrome oxidase I in a group of rodents (pocket gophers) and their ectoparasites (chewing lice) provide evidence for cospeciation and reveal different rates of molecular evolution in the hosts and their parasites. The overall rate of nucleotide substitution (both silent and replacement changes) is approximately three times higher in lice, and the rate of synonymous substitution (based on analysis of fourfold degenerate sites) is approximately an order of magnitude greater in lice. The difference in synonymous substitution rate between lice and gophers correlates with a difference of similar magnitude in generation times.

Rates of Transition and Transversion in Coding Sequences since the Human-Rodent Divergence

Protein-coding sequences of 337 human genes were compared with those of homologous genes from rodent (mouse or rat). A composite alignment containing 477,189 nucleotide positions was constructed, and 21,570 amino acid replacements were inferred. The rates of transitional and transversional silent substitutions in fourfold degenerate sites are estimated as 1.71 × 10 -9 and 1.22 × 10 -9 site -1 year -1, respectively. Rates of substitutions in replacement sites, subject to selective constraints mediated by the genetic code, are lower, but also reflect a transitional bias. The amino acid exchange rejected least often during evolution is Asp/Glu, which is fixed at 30% the rate of transversions in silent sites. The most mutable amino acids in this survey are threonine and serine; serine coded by AGY is more mutable than serine coded by TCN. A scoring matrix for evaluating amino acid similarity was derived from this study.

DNA mismatch repair and synonymous codon evolution in mammals.

It has been suggested that the differences in synonymous codon use between mammalian genes within a genome are due to differences in the efficiency of DNA mismatch repair. This hypothesis was tested by developing a model of mismatch repair, which was used to predict the expected relationship between the rate of substitution and G+C content at silent sites. It was found that the silent-substitution rate should decline with increasing G+C content over most of the G+C-content range, if it is assumed that mismatch repair is G+C biased, an assumption which is supported by data. This prediction was then tested on a set of 58 primate and artiodactyl genes. There was no evidence of a direct decline in substitution rate with increasing G+C content, for either twofold- or fourfold-degenerate sites. It was therefore concluded that variation in the efficiency of mismatch repair is not responsible for the differences in synonymous codon use between mammalian genes. In support of this conclusion, analysis of the model also showed that the parameter range over which mismatch repair can explain the differences in synonymous codon use between genes is very small.

Compositional heterogeneity and patterns of molecular evolution in the Drosophila genome.

The rates and patterns of molecular evolution in many eukaryotic organisms have been shown to be influenced by the compartmentalization of their genomes into fractions of distinct base composition and mutational properties. We have examined the Drosophila genome to explore relationships between the nucleotide content of large chromosomal segments and the base composition and rate of evolution of genes within those segments. Direct determination of the G + C contents of yeast artificial chromosome clones containing inserts of Drosophila melanogaster DNA ranging from 140-340 kb revealed significant heterogeneity in base composition. The G + C content of the large segments studied ranged from 36.9% G + C for a clone containing the hunchback locus in polytene region 85, to 50.9% G + C for a clone that includes the rosy region in polytene region 87. Unlike other organisms, however, there was no significant correlation between the base composition of large chromosomal regions and the base composition at fourfold degenerate nucleotide sites of genes encompassed within those regions. Despite the situation seen in mammals, there was also no significant association between base composition and rate of nucleotide substitution. These results suggest that nucleotide sequence evolution in Drosophila differs from that of many vertebrates and does not reflect distinct mutational biases, as a function of base composition, in different genomic regions. Significant negative correlations between codon-usage bias and rates of synonymous site divergence, however, provide strong support for an argument that selection among alternative codons may be a major contributor to variability in evolutionary rates within Drosophila genomes.

Three patterns of mitochondrial DNA nucleotide divergence in the meadow vole, Microtus pennsylvanicus.

The DNA sequence was determined for the cytochrome c oxidase II (COII), tRNALys, and ATPase 8 genes from the mitochondrial genome of the meadow vole, Microtus pennsylvanicus. When compared to other rodents, three different patterns of evolutionary divergence were found. Nucleotide variation in tRNALys is concentrated in the T psi C loop. Nucleotide variation in the COII gene in three genera of rodents (Microtus, Mus, Rattus) consists predominantly of transitions in the third base positions of codons. The predicted amino acid sequence in highly conserved (greater than 92% similarity). Analysis of the ATPase 8 gene among four genera (Microtus, Cricetulus, Mus, Rattus) revealed more detectable transversions than transitions, many fixed first and second position mutations, and considerable amino acid divergence. The rate of nucleotide substitution at nonsynonymous sites in the ATPase 8 gene is 10 times the rate in the COII gene. In contrast, the estimated absolute mutation rate as determined by analysis of nucleotide substitutions at fourfold degenerate sites probably is the same for the two genes. The primary sequences of the ATPase 8 and COII peptides are constrained differently, but each peptide is conserved in terms of predicted secondary-level configuration.

Low nucleotide diversity in man.

The nucleotide diversity (pi) in humans is studied by using published cDNA and genomic sequences that have been carefully checked for sequencing accuracy. This measure of genetic variability is defined as the number of nucleotide differences per site between two randomly chosen sequences from a population. A total of more than 75,000 base pairs from 49 loci are compared. The DNA regions studied are the 5' and 3' untranslated regions and the amino acid coding regions. The coding regions are divided into nondegenerate sites (i.e., sites at which all possible changes are nonsynonymous), twofold degenerate sites (i.e., sites at each of which one of the three possible changes is synonymous) and fourfold degenerate sites (i.e., sites at which all three possible changes are synonymous). The pi values estimated are, respectively, 0.03 and 0.04% for the 5' and 3' UT regions, and 0.03, 0.06 and 0.11% for nondegenerate, twofold degenerate and fourfold degenerate sites. Since the highest pi value is only 0.11%, which is about one order of magnitude lower than those in Drosophila populations, the nucleotide diversity in humans is very low. The low diversity is probably due to a relatively small long-term effective population size rather than any severe bottleneck during human evolution.