Inferring pattern and process: maximum-likelihood implementation of a nonhomogeneous model of DNA sequence evolution for phylogenetic analysis

A nonhomogeneous, nonstationary stochastic model of DNA sequence evolution allowing varying equilibrium G + C contents among lineages is devised in order to deal with sequences of unequal base compositions. A maximum-likelihood implementation of this model for phylogenetic analyses allows handling of a reasonable number of sequences. The relevance of the model and the accuracy of parameter estimates are theoretically and empirically assessed, using real or simulated data sets. Overall, a significant amount of information about past evolutionary modes can be extracted from DNA sequences, suggesting that process (rates of distinct kinds of nucleotide substitutions) and pattern (the evolutionary tree) can be simultaneously inferred. G + C contents at ancestral nodes are quite accurately estimated. The new method appears to be useful for phylogenetic reconstruction when base composition varies among compared sequences. It may also be suitable for molecular evolution studies.      

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.     

Directional mutation pressure,mutator mutations,and dynamics of molecular evolution

Using a general form of the directional mutation theory, this paper analyzes the effect of mutations in mutator genes on the G + C content of DNA, the frequency of substitution mutations, and evolutionary changes (cumulative mutations) under various degrees of selective constraints. Directional mutation theory predicts that when the mutational bias between A/T and G/C nucleotide pairs is equilibrated with the base composition of a neutral set of DNA nucleotides, the mutation frequency per gene will be much lower than the frequency immediately after the mutator mutation takes place. This prediction explains the wide variation of the DNA G + C content among unicellular organisms and possibly also the wide intragenomic heterogeneity of third codon positions for the genes of multicellular eukaryotes. The present analyses lead to several predictions that are not consistent with a number of the frequently held assumptions in the field of molecular evolution, including belief in a constant rate of evolution, symmetric branching of phylogenetic trees, the generality of higher mutation frequency for neutral sets of nucleotides, the notion that mutator mutations are generally deleterious because of their high mutation rates, and teleological explanations of DNA base composition. Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

Nonpolar nucleobase analogs illuminate requirements for site-specific DNA cleavage by vaccinia topoisomerase

Vaccinia DNA topoisomerase forms a covalent DNA-(3'-phosphotyrosyl)-enzyme intermediate at a specific target site 5'-C(+5)C(+4)C(+3)T(+2)T(+1)p downward arrow N(-1) in duplex DNA. Here we study the effects of nonpolar pyrimidine isosteres difluorotoluene (F) and monofluorotoluene (D) and the nonpolar purine analog indole at individual positions of the scissile and nonscissile strands on the rate of single-turnover DNA transesterification and the cleavage-religation equilibrium. Comparison of the effects of nonpolar base substitution to the effects of abasic lesions reported previously allowed us to surmise the relative contributions of base-stacking and polar edge interactions to the DNA transesterification reactions. For example, the deleterious effects of eliminating the +2T base on the scissile strand were rectified by introducing the nonpolar F isostere, whereas the requirement for the +1T base was not elided by F substitution. We impute a role for +1T in recruiting the catalytic residue Lys-167 to the active site. Topoisomerase is especially sensitive to suppression of DNA cleavage upon elimination of the +4G and +3G bases of the nonscissile strand. Indole provided little or no gain of function relative to abasic lesions. Inosine substitutions for +4G and +3G had no effect on transesterification rate, implying that the guanine exocyclic amine is not a critical determinant of DNA cleavage. Prior studies of 2-aminopurine and 7-deazaguanine effects had shown that the O6 and N7 of guanine were also not critical. These findings suggest that either the topoisomerase makes functionally redundant contacts with polar atoms (likely via Tyr-136, a residue important for precleavage active site assembly) or that it relies on contacts to N1 or N3 of the purine ring. The cleavage-religation equilibrium is strongly skewed toward trapping of the covalent intermediate by elimination of the +1A base of the nonscissile strand; the reaction equilibrium is restored by +1 indole, signifying that base stacking flanking the nick is critical for the religation step. Our findings highlight base isosteres as valuable tools for the analysis of proteins that act on DNA in a site-specific manner.     

Directional mutation pressure,selective constraints,and genetic equilibria

Summary Rates of substitution mutations in two directions, v [from an A-T or T-A nucleotide pair (AT-pair) to a G-C or C-G nucleotide pair (GC-pair)] and u [from a GC-pair to an AT-pair], are usually not the same. The net effect, v/(u + v), has previously been defined as directional mutation pressure ( _d ), which explains the wide interspecific variation and narrow intragenomic heterogeneity of DNA G+C content in bacteria. In this article, first, a theory of the evolution of DNA G+C content is presented that is based on the equilibrium among three components: directional mutation pressure, DNA G+C content, and selective constraints. According to this theory, consideration of both u and v as well as selective constraints is essential to explain the molecular evolution of the DNA base composition and sequence. Second, the theory of directional mutation pressure is applied to the analysis of the wide intragenomic heterogeneity of DNA G+C content in multicellular eukaryotes. The theory explains the extensive intragenomic heterogeneity of G+C content of higher eukaryotes primarily as the result of the intragenomic differences of directional mutation pressure and selective constraints rather than the result of positive selections for functional advantages of the DNA G+C content itself.     

Individual nucleotide bases, not base pairs, are critical for triggering site-specific DNA cleavage by vaccinia topoisomerase

Vaccinia DNA topoisomerase forms a covalent DNA-(3'-phosphotyrosyl)-enzyme intermediate at a specific target site 5'-C(+5)C(+4)C(+3)T(+2)T(+1)p downward arrow N(-1) in duplex DNA. Here we study the effects of abasic lesions at individual positions of the scissile and nonscissile strands on the rate of single-turnover DNA transesterification and the cleavage-religation equilibrium. The rate of DNA incision was reduced by factors of 350, 250, 60, and 10 when abasic sites replaced the -1N, +1T, +2T, and +4C bases of the scissile strand, but abasic lesions at +5C and +3C had little or no effect. Abasic lesions in the nonscissile strand in lieu of +4G, +3G, +2A, and +1A reduced the rate of cleavage by factors of 130, 150, 10, and 5, whereas abasic lesions at +5G and -1N had no effect. The striking positional asymmetry of abasic interference on the scissile and nonscissile strands highlights the importance of individual bases, not base pairs, in promoting DNA cleavage. The rate of single-turnover DNA religation by the covalent topoisomerase-DNA complex was insensitive to abasic sites within the CCCTT sequence of the scissile strand, but an abasic lesion at the 5'-OH nucleoside (-1N) of the attacking DNA strand slowed the rate of religation by a factor of 600. Nonscissile strand abasic lesions at +1A and -1N slowed the rate of religation by factors of approximately 140 and 20, respectively, and strongly skewed the cleavage-religation equilibrium toward the covalent complex. Thus, abasic lesions immediately flanking the cleavage site act as topoisomerase poisons.     

Asymmetric directional mutation pressures in bacteria

Background

When there are no strand-specific biases in mutation and selection rates (that is, in the substitution rates) between the two strands of DNA, the average nucleotide composition is theoretically expected to be A = T and G = C within each strand. Deviations from these equalities are therefore evidence for an asymmetry in selection and/or mutation between the two strands. By focusing on weakly selected regions that could be oriented with respect to replication in 43 out of 51 completely sequenced bacterial chromosomes, we have been able to detect asymmetric directional mutation pressures.

Results

Most of the 43 chromosomes were found to be relatively enriched in G over C and T over A, and slightly depleted in G+C, in their weakly selected positions (intergenic regions and third codon positions) in the leading strand compared with the lagging strand. Deviations from A = T and G = C were highly correlated between third codon positions and intergenic regions, with a lower degree of deviation in intergenic regions, and were not correlated with overall genomic G+C content.

Conclusions

During the course of bacterial chromosome evolution, the effects of asymmetric directional mutation pressures are commonly observed in weakly selected positions. The degree of deviation from equality is highly variable among species, and within species is higher in third codon positions than in intergenic regions. The orientation of these effects is almost universal and is compatible in most cases with the hypothesis of an excess of cytosine deamination in the single-stranded state during DNA replication. However, the variation in G+C content between species is influenced by factors other than asymmetric mutation pressure.

     

Major groove interactions of vaccinia Topo I provide specificity by optimally positioning the covalent phosphotyrosine linkage

Vaccinia DNA topoisomerase (vTopo) is a prototypic eukaryotic type I topoisomerase that shows high specificity for nucleophilic substitution at a single phosphodiester linkage in the pentapyrimidine recognition sequence 5'-(C/T)+5 C+4 C+3 T+2 T+1 p / N(-1). This reaction involves reversible transesterification where the active site tyrosine of the enzyme and a 5'-hydroxyl nucleophile of DNA compete for attack at the phosphoryl group. The finite lifetime of the covalent phosphotyrosine adduct allows the enzyme to relax multiple supercoils by rotation of the 5'-OH strand before the DNA backbone is religated. To dissect the nature of the unique sequence specificity, subtle modifications to the major groove of the GGGAA 5'-sequence of the nonscissile strand were introduced and their effects on each step of the catalytic cycle were measured. Although these modifications had no effect on noncovalent DNA binding (K(D)) or the rate of reversible DNA cleavage (k(cl)), significant decreases in the cleavage equilibrium (K(cl) = k(cl)/k(r)) arising from increased rates of 5'-hydroxyl attack (k(r)) at the phosphotyrosine linkage were observed. These data and other findings support a model in which major groove interactions are used to position the phosphotyrosine linkage relative to the mobile 5'-hydroxyl nucleophile. In the absence of native sequence interactions, the phosphotyrosine has a higher probability of encountering the 5'-hydroxyl nucleophile, leading to an enhanced rate of ligation and a diminished equilibrium constant for cleavage. By this unusual specificity mechanism, the enzyme prevents formation of stable covalent adducts at nonconsensus sites in genomic DNA.     

The Role of Context-Dependent Mutations in Generating Compositional and Codon Usage Bias in Grass Chloroplast DNA

Abstract The influence of local base composition on mutations in chloroplast DNA (cpDNA) is studied in detail and the resulting, empirically derived, mutation dynamics are used to analyze both base composition and codon usage bias. A 4 × 4 substitution matrix is generated for each of the 16 possible flanking base combinations (contexts) using 17,253 noncoding sites, 1309 of which are variable, from an alignment of three complete grass chloroplast genome sequences. It is shown that substitution bias at these sites is correlated with flanking base composition and that the A+T content of these flanking sites as well as the number of flanking pyrimidines on the same strand appears to have general influences on substitution properties. The context-dependent equilibrium base frequencies predicted from these matrices are then applied to two analyses. The first examines whether or not context dependency of mutations is sufficient to generate average compositional differences between noncoding cpDNA and silent sites of coding sequences. It is found that these two classes of sites exist, on average, in very different contexts and that the observed mutation dynamics are expected to generate significant differences in overall composition bias that are similar to the differences observed in cpDNA. Context dependency, however, cannot account for all of the observed differences: although silent sites in coding regions appear to be at the equilibrium predicted, noncoding cpDNA has a significantly lower A+T content than expected from its own substitution dynamics, possibly due to the influence of indels. The second study examines the codon usage of low-expression chloroplast genes. When context is accounted for, codon usage is very similar to what is predicted by the substitution dynamics of noncoding cpDNA. However, certain codon groups show significant deviation when followed by a purine in a manner suggesting some form of weak selection other than translation efficiency. Overall, the findings indicate that a full understanding of mutational dynamics is critical to understanding the role selection plays in generating composition bias and sequence structure.     

The role of DNA replication and isochores in generating mutation and silent substitution rate variance in mammals.

It has been suggested that isochores are maintained by mutation biases, and that this leads to variation in the rate of mutation across the genome. A model of DNA replication is presented in which the probabilities of misincorporation and proofreading are affected by the composition and concentration of the free nucleotide pools. The relationship between sequence G+C content and the mutation rate is investigated. It is found that there is very little variation in the mutation rate between sequences of different G+C contents if the total concentration of the free nucleotides remains constant. However, variation in the mutation rate can be arbitrarily large if some mismatches are proofread and the total concentration of free nucleotides varies. Hence the model suggests that the maintenance of isochores by the replication of DNA in free nucleotide pools of biased composition does not lead per se to mutation rate variance. However, it is possible that changes in composition could be accompanied by changes in concentration, thus generating mutation rate variance. Furthermore, there is the possibility that germ-line selection could lead to alterations in the overall free nucleotide concentration through the cell cycle. These findings are discussed with reference to the variance in mammalian silent substitution rates.     

Estimating evolutionary distance from restriction maps of mitochondrial DNA with arbitrary G+C content

Summary We develop a mathematical model for estimating evolutionary distance from restriction enzyme maps, which incorporates non-uniformity of the rate of base substitution into the theory and allows for an arbitrary G+C content at equilibrium. When the G+C content differs significantly from 1/2, the traditional model of base changes can introduce a systematic bias which depends upon the base composition of the restriction site. In addition, the accuracy of estimated evolutionary distance depends heavily upon the choice of restriction enzyme in that the expected number of sites is also affected. Monte Carlo experiments are conducted to check the validity of the present theoretical treatment and from which we draw several cautionary notes on estimation. An application is made to the available data on restriction enzyme maps of human mitochondrial DNA where the G+C content is approximately 1/3.Contribution No. 1372 from the National Institute of Genetics, Mishima, 411 Japan     

Site-specific DNA transesterification by vaccinia topoisomerase: effects of benzo[alpha]pyrene-dA, 8-oxoguanine, 8-oxoadenine and 2-aminopurine modifications

Vaccinia DNA topoisomerase forms a covalent DNA-(3'-phosphotyrosyl)-enzyme intermediate at a specific target site 5'-C+5C+4C+3T+2T+1p downward arrow N-1 in duplex DNA. Here we study the effects of base modifications on the rate and extent of single-turnover DNA transesterification. Chiral trans opened C-10 R and S adducts of benzo[a]pyrene (BP) 7,8-diol 9,10-epoxide were introduced at single N6-deoxyadenosine (dA) positions within the 3'-G+5G+4G+3A+2A+1T-1A-2 sequence of the nonscissile DNA strand. The R and S BPdA adducts intercalate from the major groove on the 5' and 3' sides of the modified base, respectively, and perturb local base stacking. We found that R and S BPdA modifications at +1A reduced the transesterification rate by a factor of 700-1000 without affecting the yield of the covalent topoisomerase-DNA complex. BPdA modifications at +2A reduced the extent of transesterification and elicited rate decrements of 200- and 7000-fold for the S and R diastereomers, respectively. In contrast, BPdA adducts at the -2 position had no effect on the extent of the reaction and relatively little impact on the rate of cleavage. A more subtle probe of major groove contacts entailed substituting each of the purines of the nonscissile strand with its 8-oxo analog. The +3 oxoG modification slowed transesterification 35-fold, whereas other 8-oxo modifications were benign. 8-Oxo substitutions at the -1 position in the scissile strand slowed single-turnover cleavage by a factor of six but had an even greater slowing effect on religation, which resulted in an increase in the cleavage equilibrium constant. 2-Aminopurine at positions +3, +4, or +5 in the nonscissile strand had no effect on transesterification per se but had synergistic effects when combined with 8-oxoA at position -1 in the scissile strand. These findings illuminate the functional interface of vaccinia topoisomerase with the DNA major groove.     

DNA G+C content of the third codon position and codon usage biases of human genes

The human genome, as in other eukaryotes, has a wide heterogeneity in the DNA base composition. The evolutionary basis for this heterogeneity has been unknown. A previous study of the human genome (846 genes analyzed) has shown that, in the major range of the G+C content in the third codon position (0.25-0.75), biases from the Parity Rule 2 (PR2) among the synonymous codons of the four-codon amino acids are similar except in the highest G+C range (Sueoka, N., 1999. Translation-coupled violation of Parity Rule 2 in human genes is not the cause of heterogeneity of the DNA G+C content of third codon position. Gene 238, 53-58.). PR2 is an intra-strand rule where A=T and G=C are expected when there are no biases between the two complementary strands of DNA in mutation and selection rates (substitution rates). In this study, 14,026 human genes were analyzed. In addition, the third codon positions of two-codon amino acids were analyzed. New results show the following: (a) The G+C contents of the third codon position of human genes are scattered in the G+C range of 0.22-0.96 in the third codon position. (b) The PR2 biases are similar in the range of 0.25-0.75, whereas, in the high G+C range (0.75-0.96; 13% of the genes), the PR2-bias fingerprints are different from those of the major range. (c) Unlike the PR2 biases, the G+C contents of the third codon position for both four-codon and two-codon amino acids are all correlated almost perfectly with the G+C content of the third codon position over the total G+C ranges. These results support the notion that the directional mutation pressure, rather than the directional selection pressure, is mainly responsible for the heterogeneity of the G+C content of the third codon position.     

Effect of DNA base composition on the intercalation of proflavine. A kinetic study.

The effect of DNA base composition on the kinetics of the association between DNA and proflavine has been investigated using the temperature jump relaxation method. It is found that, regardless of the G + C base composition the results fit a two step mechanism, the second of which exhibits characteristics of intercalation of proflavine into DNA. However, they two equilibrium constants corresponding to these steps, KI and KII, depend on the nature of the DNAs. The constant KI is found to be an order of magnitude greater for M. lysodeikticus DNA (72% G + C) than for calf thymus DNA (48% G + C). Increasing G-C content thus appears to favor the intermediate non-intercalated complex of proflavine with DNA. Methylation of M. lysodeikticus DNA with dimethyl sulfate, preferentially yielding N7 methyl guanine as the modified base, again leads to an apparent two step mechanism, with the value of KI unchanged with respect to untreated DNA, while the affinity of proflavine for the intercalated complex measured by the value of KII increases for methylated DNA.     

Metabolic rate and directional nucleotide substitution in animal mitochondrial DNA

There is marked heterogeneity of nucleotide composition in mitochondrial DNA across divergent animals. Differences in nucleotide composition presumably reflect differences in directional nucleotide substitution for A+T or G+C nucleotides. In mitochondrial DNA, there is A+T directional nucleotide substitution in most (if not all) animals surveyed, and the magnitude of directional A+T nucleotide substitution differs greatly within and among groups. Differences in directional nucleotide substitution among lineages of mammals can be explained by changes in metabolic physiology. This relationship is thought to be mediated by the effect of oxygen radicals because these toxic compounds are by-products of aerobic metabolism and are known mutagens. Association between metabolism and nucleotide composition provides additional evidence in favor of the hypothesis that rates and patterns of nucleotide substitution in mitochondrial DNA can be influenced by factors that impinge on rates of endogenous DNA damage.      

Fidelity of replication of the leading and the lagging DNA strands opposite N-methyl-N-nitrosourea-induced DNA damage in human cells.

Semi-conservative replication of double-stranded DNA in eukaryotic cells is an asymmetric process involving leading and lagging strand synthesis and different DNA polymerases. We report a study to analyze the effect of these asymmetries when the replication machinery encounters alkylation-induced DNA adducts. The model system is an EBV-derived shuttle vector which replicates in synchrony with the host human cells and carries as marker gene the bacterial gpt gene. A preferential distribution of N-methyl-N-nitrosourea (MNU)-induced mutations in the non transcribed DNA strand of the shuttle vector pF1-EBV was previously reported. The hypermutated strand was the leading strand. To test whether the different fidelity of DNA polymerases synthesizing the leading and the lagging strands might contribute to MNU-induced mutation distribution the mutagenesis study was repeated on the shuttle vector pTF-EBV which contains the gpt gene in the inverted orientation. We show that the base substitution error rates on an alkylated substrate are similar for the replication of the leading and lagging strands. Moreover, we present evidence that the fidelity of replication opposite O6-methylguanine adducts of both the leading and lagging strands is not affected by the 3' flanking base. The preferential targeting of mutations after replication of alkylated DNA is mainly driven by the base at the 5' side of the G residues.