DNA序列进化过程中核苷酸替代的非独立性研究

本文评述了DNA序列间核苷酸替代数的估计方法,并通过对七个物种中组蛋白基因的比较对DNA进化的模型进行了考察。发现H2A基因第三位点上的碱基组成在物种间变异很大,并且跟H2A基因第一位点、H4基因第一、三位点及H2A上游,下游序列中的碱基组成有强正相关,提示DNA序列进化过程中存在着物种特异的区域性约束力。可能的原因是高等真核生物中GC含量升高,或者是染色体重组使这些同源序列位于不同的等质区段,从而受到不同的选择突变压。密码内各位点上核苷酸替代的相关性分析表明不同位点的替代是非独立的,其原因可能是一次替代事件引起多个位点的变化。文中讨论了这些结果对进化树推断的意义。     

Conserved nucleotide differences and subfamily structure of porcine short interspersed elements

Interspersed elements are ubiquitous in the genomes of higher eukaryotes and account for over a third of the genomic DNA (Smit 1996). In swine the short interspersed elements, SINEs or PREs (porcine repetitive elements), have been found in a number of introns and 3' untranslated regions of different genes. However, compared to human Alu repeats the number of available PRE DNA sequences is still limited. In this study we have compared 85 PREs selected from DNA sequence database entries. The PREs were aligned and for each nucleotide position the relative frequencies of the four bases were calculated. A consensus sequence was derived from the first base usage. Similar to studies of SINEs in other species, the analysis showed that most mutations in PREs occur at CpG dinucleotide hot spots. The position variability for the two most frequent bases shows a bimodal distribution. The analysis suggests that the porcine SINEs can be divided into three major subfamilies sharing conserved nucleotide similarities.     

Sequence requirements for protein-primed initiation and elongation of phage O29 DNA replication

The double-stranded linear DNA of Bacillus subtilis phage O29 is replicated by a mechanism in which a terminal protein (TP) acts as a primer. The second 3'-terminal nucleotide of the template directs the incorporation of the 5'-terminal nucleotide into the TP, giving rise to the initiation complex TP-dAMP. Elongation then proceeds by a sliding-back mechanism in which the dAMP covalently linked to the TP pairs to the 3'-terminal nucleotide of the template strand to recover full-length DNA. We have studied the sequence requirements for efficient initiation of replication using mutated TP-free double-stranded DNA fragments. Efficient initiation only requires the terminal repetition 5'-AA. The 3'-terminal T, although not used as template, increases the affinity of DNA polymerase for the initiator nucleotide; in addition, although to a minor extent, the third 3'-terminal position also directs the formation of the initiation complex and modulates the initiation rate at the second position. Efficient elongation requires a previous sliding-back, demanding again a repetition of two nucleotides at the 3' end; if the sliding-back is prevented, a residual elongation can proceed directly from the second position or after jumping back from the third to the first position.     

Molecular evolution of the hepatitis B virus genome

The hepatitis B virus (HBV) has a circular DNA genome of about 3,200 base pairs. Economical use of the genome with overlapping reading frames may have led to severe constraints on nucleotide substitutions along the genome and to highly variable rates of substitution among nucleotide sites. Nucleotide sequences from 13 complete HBV genomes were compared to examine such variability of substitution rates among sites and to examine the phylogenetic relationships among the HBV variants. The maximum likelihood method was employed to fit models of DNA sequence evolution that can account for the complexity of the pattern of nucleotide substitution. Comparison of the models suggests that the rates of substitution are different in different genes and codon positions; for example, the third codon position changes at a rate over ten times higher than the second position. Furthermore, substantial variation of substitution rates was detected even after the effects of genes and codon positions were corrected; that is, rates are different at different sites of the same gene or at the same codon position. Such rates after the correction were also found to be positively correlated at adjacent sites, which indicated the existence of conserved and variable domains in the proteins encoded by the viral genome. A multiparameter model validates the earlier finding that the variation in nucleotide conservation is not random around the HBV genome. The test for the existence of a molecular clock suggests that substitution rates are more or less constant among lineages. The phylogenetic relationships among the viral variants were examined. Although the data do not seem to contain sufficient information to resolve the details of the phylogeny, it appears quite certain that the serotypes of the viral variants do not reflect their genetic relatedness. Correspondence to: Z. Yang     

Transition and transversion rate in the evolution of animal mitochondrial DNA

We present a further application of the stochastic model previously described (Lanave et al., 1984, 1985) for measuring the nucleotide substitution rate in the mammalian evolution of the mitochondrial DNA (mtDNA). The applicability of this method depends on the validity of "stationarity conditions" (equal nucleotide frequencies at first, second and third silent codon positions in homologous protein coding genes). In the comparison of homologous sequences satisfying the stationarity condition at the silent sites, only the four codon families (quartets) for which both transitions and transversions are silent at the third position are considered here. This has allowed us to estimate the transition and transversion rates for any pair of species. We have analyzed the third silent codon position of the triplet rat-mouse-cow, of a series of slightly divergent primates and of two Drosophila species. In terms of two external dating input we have then determined the phylogenetic trees for rat, mouse, and cow as well as for a number of primates including man. The phylogenetic tree that we have derived for the triplet rat, mouse and cow agrees with that we had previously determined by analyzing the first, second and third silent codon positions (in both duets and quartets) of mt genes (Lanave et al., 1985). For primates our method leads to the following branching order from the oldest to the most recent: Gibbon, Orangutan, Gorilla, Chimpanzee and Man. In absolute time, fixing the distance Chimpanzee-Man as 5 million years (Myr) we estimate the dating of the divergence nodes as: Gorilla 7 Myr; Orangutan 16 Myr; Gibbon 20 Myr. In all cases analyzed, the transition rate has been found to be substantially higher than the transversion rate. Moreover we have found that the transition/transversion ratio is different in the various lineages. We suggest that this fact is probably related to the nucleotide frequencies at the third silent codon position.     

Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

Nonsynonymous site heteroplasmy in fish mitochondrial DNA

Heteroplasmic nucleotide polymorphisms are rarely observed in wild animal mitochondrial DNA. The occurrence of such site heteroplasmy is expected to be extremely rare at nonsynonymous sites where the number of nucleotide substitutions per site is low due to functional constraints. This report deals with nonsynonymous mitochondrial heteroplasmy from two wild fish species, chum salmon and Japanese flounder. We detected an A/C nonsynonymous heteroplasmic site corresponding to putative amino acids, Ile or Met, in NADH dehydrogenase subunit-5 (ND5) region of chum salmon. The heteroplasmic site was at the 3rd position of 58th codon. As for Japanese flounder we detected a C/T nonsynonymous heteroplasmic site corresponding to putative amino acids, Leu or Pro, in ND4 region. The heteroplasmic site was at the 2nd position of 450th codon. We also verified heteroplasmy at these sites by sequencing cloned fragments.     

Directed termination of the polymerase chain reaction: kinetics and applications in mutation detection.

We describe a PCR-based method (DT-PCR) that integrates both DNA amplification and directed chain termination into a single-step process. This method exploits unbalanced nucleotide concentrations to induce the polymerase chain reaction to terminate at specific nucleotide sites, leading to the generation of two sets of nested termination fragments from genomic DNA. The kinetic mechanism underlying the termination process is outlined and the application of this method to the detection and characterization of mutations in fragments as long as 1 kb is described. The method is effective for the analysis of both haploid and diploid genomes, and not only allows the recognition of both indels (insertions and deletions) and nucleotide substitutions, but also enables the determination of their position in a single-step fashion.     

Determination of netropsin-DNA binding constants from footprinting data

A theory for deriving drug-DNA site binding constants from footprinting data is presented. Plots of oligonucleotide concentration, as a function of drug concentration, for various cutting positions on DNA are required. It is assumed that the rate of cleavage at each nucleotide position is proportional to the concentration of enzyme at that nucleotide and to the probability that the nucleotide is not blocked by drug. The probability of a nucleotide position not being blocked is calculated by assuming a conventional binding equilibrium for each binding site with exclusions for overlapping sites. The theory has been used to evaluate individual site binding constants for the antiviral agent netropsin toward a 139 base pair restriction fragment of pBR-322 DNA. Drug binding constants, evaluated from footprinting data in the presence of calf thymus DNA and poly(dGdC) as carrier and in the absence of carrier DNA, were determined by obtaining the best fit between calculated and experimental footprinting data. Although the strong sites on the fragment were all of the type (T.A)4, the value of the binding constant was strongly sequence dependent. Sites containing the dinucleotide sequence 5'-TA-3' were found to have significantly lower binding constants than those without this sequence, suggesting that an adenine-adenine clash produces a DNA structural alteration in the minor groove which discourages netropsin binding to DNA. The errors, scope, and limitations associated with the method are presented and discussed.     

Expression,purification and biochemical characterization of Methanocaldococcus jannaschii DNA ligase

We describe the biochemical characterization of Methanocaldococcus jannaschii (M. jannaschii) DNA ligase and its potential application in single nucleotide polymorphism (SNP) genotyping. The recombinant M. jannaschii DNA ligase is an ATP-dependent ligase. The ligase activity was dependent on metal ions of Mg²⁺ and Mn²⁺. The optimal concentrations of ATP cofactor and Mg²⁺ ion were 0.01–2 and 10 mM, respectively. The optimal pH value for DNA ligation was 8.5. High concentrations of NaCl inhibited DNA ligation. The effects of mismatches on joining short oligonucleotides by M. jannaschii DNA ligase were fully characterized. The mismatches at the first position 5′ to the nick inhibited ligation more than those at the first position 3′ to the nick. The mismatches at other positions 5′ to the nick (3rd to 7th sites) exhibited less inhibition on ligation. However, the introduction of a C/C mismatch at the third position 5′ to the nick could completely inhibit the ligation of the terminal-mismatched nick of an oligonucleotide duplex by M. jannaschii DNA ligase. Therefore, introducing an additional mismatch at the third position 5′ to the SNP site is a more effective approach in genotyping by M. jannaschii DNA ligase.     

The implications of using mutagenic primers in combination with Taq polymerase having proofreading activity.

Polymerases with proofreading activity provide high fidelity PCR amplifications. In this study we examined the consequences of using a Taq polymerase with proofreading activity, such as Optimase Taq polymerase, in combination with 4 different mutagenic reverse primers for the amplification of a 345-bp FII PCR product. The amplifications were performed with Optimase Taq polymerase (Transgenomic), and Taq DNA polymerase-recombinant (Invitrogen), without proofreading activity. Mutation screening was carried out by DHPLC and restriction fragment analysis. The usage of Optimase Taq polymerase results in complete reversion of the first and second mutated nucleotide introduced at the 3' end of the mutagenic reverse primer. It also partially reverses the missense nucleotide introduced in the third position of the mutagenic primer and leads to misleading DHPLC and restriction fragment analysis patterns. Nevertheless it cannot perform such an activity when an abnormal nucleotide is introduced in the fourth position.     

Sense codon in Escherichia coli are translated in context

The nucleotide frequencies 5' and 3' to the sense codons in highly and weakly expressed genes have been investigated by the chi-squares method. A comparison between the experimental and computer-generated random nucleotide sequences (in which each codon is substituted by a random synonymous one) was made. It was shown that the choice of a particular codon among the synonymous ones in a given position of the gene depends on the three nucleotides 3' and 5' adjacent to the codon in highly expressed genes (the triplet 3' and a single nucleotide 5' to the codons in weakly expressed genes). Concrete patterns for the preferable choice of synonymous codons depending on their contexts are presented. It is suggested that these constraints are related to the efficiency of messenger translation. The constraints on the amino acid sequences of encoded proteins also lead to statistically significant bases in nucleotide frequencies around the sense codons. The biological role of these constraints is discussed.     

Description of the cytochrome c oxidase subunit II gene in some genera of New World monkeys (primates,Platyrrhini)

Nucleotide sequence variation at the mitochondrial cytochrome c oxidase subunit II gene (COII) was analyzed in 27 New World monkey specimens, nine newly reported herein. The study involved comparisons among platyrrhines and also between platyrrhines and catarrhines. The analysis of the frequencies of transitions and transversions at each codon position showed transitional saturation at third codon position. Neighbor-Joining trees obtained from genetic distances estimated by means of Kimura's (1980) two-parameter model showed poor resolution of phylogenetic relationships among platyrrhine genera. Rates of nucleotide substitutions were largely homogeneous except in the genus Saimiri that showed low numbers of unique substitutions suggesting the maintenance of ancestral or plesiomorphic states of platyrrhine mt DNA nucleotide characters probably due to its large population sizes as compared to other platyrrhines.     

Radioprobing of DNA: distribution of DNA breaks produced by decay of 125I incorporated into a triplex-forming oligonucleotide correlates with geometry of the triplex.

The distribution of breaks produced in both strands of a DNA duplex by the decay of 125I carried by a triplex-forming DNA oligonucleotide was studied at single nucleotide resolution. The 125I atom was located in the C5 position of a single cytosine residue of an oligonucleotide designed to form a triple helix with the target sequence duplex. The majority of the breaks (90%) are located within 10 bp around the decay site. The addition of the free radical scavenger DMSO produces an insignificant effect on the yield and distribution of the breaks. These results suggest that the majority of these breaks are produced by the direct action of radiation and are not mediated by diffusible free radicals. The frequency of breaks in the purine strand was two times higher that in the pyrimidine strand. This asymmetry in the yield of breaks correlates with the geometry of this type of triplex; the C5 of the cytosine in the third strand is closer to the sugar-phosphate backbone of the purine strand. Moreover, study of molecular models shows that the yield of breaks at individual bases correlates with distance from the 125I decay site. We suggest the possible use of 125I decay as a probe for the structure of nucleic acids and nucleoprotein complexes.     

Site-specific excision repair of 1-nitrosopyrene-induced DNA adducts at the nucleotide level in the HPRT gene of human fibroblasts: effect of adduct conformation on the pattern of site-specific repair.

Studies showing that different types of DNA adducts are repaired in human cells at different rates suggest that DNA adduct conformation is the major determinant of the rate of nucleotide excision repair. However, recent studies of repair of cyclobutane pyrimidine dimers or benzo[a]pyrene diol epoxide (BPDE)-induced adducts at the nucleotide level in DNA of normal human fibroblasts indicate that the rate of repair of the same adduct at different nucleotide positions can vary up to 10-fold, suggesting an important role for local DNA conformation. To see if site-specific DNA repair is a common phenomenon for bulky DNA adducts, we determined the rate of repair of 1-nitrosopyrene (1-NOP)-induced adducts in exon 3 of the hypoxanthine phosphoribosyltransferase gene at the nucleotide level using ligation-mediated PCR. To distinguish between the contributions of adduct conformation and local DNA conformation to the rate of repair, we compared the results obtained with 1-NOP with those we obtained previously using BPDE. The principal DNA adduct formed by either agent involves guanine. We found that rates of repair of 1-NOP-induced adducts also varied significantly at the nucleotide level, but the pattern of site-specific repair differed from that of BPDE-induced adducts at the same guanine positions in the same region of DNA. The average rate of excision repair of 1-NOP adducts in exon 3 was two to three times faster than that of BPDE adducts, but at particular nucleotides the rate was slower or faster than that of BPDE adducts or, in some cases, equal to that of BPDE adducts. These results indicate that the contribution of the local DNA conformation to the rate of repair at a particular nucleotide position depends upon the specific DNA adduct involved. However, the data also indicate that the conformation of the DNA adduct is not the only factor contributing to the rate of repair at different nucleotide positions. Instead, the rate of repair at a particular nucleotide position depends on the interaction between the specific adduct conformation and the local DNA conformation at that nucleotide.     

Complete sequence of bovine mitochondrial DNA conserved features of the mammalian mitochondrial genome

We present here the complete 16,338 nucleotide DNA sequence of the bovine mitochondrial genome. This sequence is homologous to that of the human mitochondrial genome (Anderson et al., 1981) and the genes are organized in virtually identical fashion. The bovine mitochondrial protein genes are 63 to 79% homologous to their human counterparts, and most of the nucleotide differences occur in the third positions of codons. The minimum rate of base substitution that accounts for the nucleotide differences in the codon third positions is very high: at least 6 × 10^?9 changes per position per year. The bovine and human mitochondrial transfer RNA genes exhibit more interspecies variation than do their cytoplasmic counterparts, with the “TΨC” loop being the most variable part of the molecule. The bovine 12 S and 16 S ribosomal RNA genes, when compared with those from human mitochondrial DNA, show conserved features that are consistent with proposed secondary structure models for the ribosomal RNAs. Unlike the pattern of moderate-to-high homology between the bovine and human mitochondrial DNAs found over most of the genome, the DNA sequence in the bovine D-loop region is only slightly homologous to the corresponding region in the human mitochondrial genome. This region is also quite variable in length, and accounts for the bulk of the size difference between the human and bovine mitochondrial DNAs.     

Characteristics of Nucleotide Substitution in the Hepatitis C Virus Genome: Constraints on Sequence Change in Coding Regions at Both Ends of the Genome

Comparison of complete genome sequences for different variants of hepatitis C virus (HCV) reveals several different constraints on sequence change. Synonymous changes are suppressed in coding regions at both 5′ and 3′ ends of the genome. No evidence was found for the existence of alternative reading frames or for a lower mutation frequency in these regions. Instead, suppression may be due to constraints imposed by RNA secondary structures identified within the core and NS5b genes. Nonsynonymous substitutions are less frequent than synonymous ones except in the hypervariable region of E2 and, to a lesser extent, in E1, NS2, and NS5b. Transitions are more frequent than transversions, particularly at the third position of codons where the bias is 16:1. In addition, nucleotide substitutions may not occur symmetrically since there is a bias toward G or C at the third position of codons, while T ↔ C transitions were twice as frequent as A ↔ G transitions. These different biases do not affect the phylogenetic analysis of HCV variants but need to be taken into account in interpreting sequence change in longitudinal studies. Received: 9 September 1996 / Accepted: 20 April 1997     

Positive selection on the Plasmodium falciparum sporozoite threonine-asparagine-rich protein: analysis of isolates mainly from low endemic areas

The sporozoite threonine-asparagine-rich protein (STARP) of Plasmodium falciparum is an attractive target for a pre-erythrocytic stage malaria vaccine because both naturally acquired and experimentally induced anti-STARP antibodies can block sporozoite invasion of hepatocytes. To explore the extent of sequence variation, we surveyed nucleotide polymorphism across the entire gene, encompassing 2 exons and an intron, of 124 P. falciparum-infected blood samples from Thailand and 10 from 4 other endemic areas. In total 24 haplotypes were identified despite low-level nucleotide diversity at this locus. The mean number of nonsynonymous substitutions per nonsynonymous site (d(N)) significantly exceeded that of synonymous substitutions per synonymous site (d(S)), suggesting that the STARP gene has evolved under positive selection, probably from host immune pressure. The preponderance of conservative amino acid exchanges and a strongly biased T-nucleotide toward the third position of codons in repeat arrays have reflected simultaneous constraints on this molecule, probably from its respective unknown function and nucleotide composition. Sequence conservation in the STARP locus among clinical isolates from different disease endemic areas would not compromise vaccine incorporation.