Biases of the estimates of DNA divergence obtained by the restriction enzyme technique

Summary A mathematical formula for the relationship between the average number of nucleotide substitutions per site and the proportion of shared restriction sites between two homologous nucleons is developed by taking into account the unequal rates of substitution among different pairs of nucleotides. Using this formula, the possible amount of bias of the estimate of the number of nucleotide substitutions obtained by the Upholt-Nei-Li formula for restriction site data is investigated. The results obtained indicate that the bias depends upon the nucleotides in the recognition sequence of the restriction enzyme used, the unequal rates of substitution among different nucleotides, and the unequal nucleotide frequencies, but the primary factor is the unequal rates of nucleotide substitution. The amount of bias is generally larger for four-base enzymes than for six-base enzymes. However, when many restriction enzymes are used for the study of DNA divergence, the bias is unlikely to be very large unless the rate of substitution greatly varies from nucleotide to nucleotide.     

Restriction endonuclease recognition and Southern hybridization of bromodeoxyuridine-substituted genomic DNA

Abstract. Human glioma cell lines exposed to various concentrations of bromode-oxyuridine (BrdUrd) were studied to determine the effect of BrdUrd substitution on restriction endonuclease recognition and Southern hybridization of genomic DNA. BrdUrd substitution had no effect on the recognition of restriction endonucleases. When the exposure to BrdUrd was 2 h or less and the BrdUrd substitution rate was less than 40%, there was no difference in the density of hybridized bands after Southern hybridization using human non-recombinant complementary DNA as a probe. Hybridization was suppressed significantly by exposures longer than 24 h or BrdUrd substitution rates greater than 40%. These results suggest that the BrdUrd substitution rate and the exposure time to BrdUrd influence the hybridization reaction by a DNA probe. Brief exposure (up to 2 h) to BrdUrd does not influence restriction endonuclease recognition or Southern hybridization of genomic DNA.     

Evolutionary change of restriction sites under unequal rates of nucleotide substitution among the three positions of codons

Summary Under the assumption of unequal rates of nucleotide substitution among the three positions of codons, mathematical formulas are derived for the probability that a restriction site observed at time 0 in a DNA sequence will be present at time t, the probability that a new restriction site will emerge at a particular place in two descendant sequences, and the proportion of identical restriction sites between two such sequences. All three quantitites are shown to be larger for the case of unequal rates than for the case of equal rates. As a consequence, estimates of nucleotide divergence (2t) between two sequences based on restriction-site data tend to be lower than the actual values of the assumption of equal rates is made but actually does not hold. However, the degree of underestimation is slight if 2t is 0.10 or smaller. The underestimation can be serious when 2t becomes 0.20 or larger, particularly if the substitution rates at the three codon positions are very different or if transitional substitution occurs more frequently than transversional substitution. The underestimation is more serious for four-base enzymes than for sixbase enzymes.     

Bovine mitochondrial DNA polymorphism in restriction endonuclease cleavage patterns and the location of the polymorphic sites

Cleavage patterns of mitochondrial DNA (mtDNA) by restriction endonuclease analysis were examined in four Japanese Black cows, three Japanese Shorthorn cows, and six Holstein cows. Seventeen restriction enzymes which recognize six base pairs and two restriction enzymes which recognize four base pairs were used in this study. Polymorphism was observed with three restriction enzymes, HindIII, TaqI, and MspI, and was detected within the breeds. Nucleotide substitution was determined in the HindIII polymorphic site by DNA cloning and sequencing; this is C----T at position 10126 of the URF-3 region. Furthermore, the MspI and TaqI polymorphic sites were located on the physical map.     

Sequence-specific responses of restriction endonucleases to bromodeoxyuridine substitution in mammalian DNA

Substitution of BrdU for dT in mammalian DNA alters the rates of DNA cleavage by restriction endonucleases in a manner that can be related to the specificity of cleavage. A formula is proposed that describes inhibitory and stimulatory contributions arising from the substitution of a Br atom for the CH3 group on T. The larger Br atom is postulated to sterically hinder the nuclease from binding to adjacent groups in the DNA cleavage site, while allowing a tighter binding to itself. The inhibition caused by steric hindrance is predicted to vary inversely with distance from the point of cleavage, whereas the stimulation caused by tighter binding is predicted to be independent of distance. The resultant formula gives a good fit to the data obtained for thirteen different restriction nucleases of known specificity. The parameters in the formula appear to be simple functions of ionic strength. This formula can be used to predict the effect of BrdU substitution on any endonuclease whose specificity of cleavage is known.     

An ApaLI restriction site polymorphism is associated with the MB19 polymorphism in apolipoprotein B

An apolipoprotein (apo) B-specific monoclonal antibody, MB19, detects a commonly occurring two-allele genetic polymorphism in human apoB (Young, S. G., S. J. Bertics, L. K. Curtiss, D. C. Casal, and J. L. Witztum. 1986. Proc. Natl. Acad. Sci. USA. 83: 1101-1105). Antibody MB19 binds to two different allotypes of apoB, MB19(1) and MB19(2), with high and low affinity, respectively. The epitope for antibody MB19 is located within apoB-100 thrombolytic fragment T4 (apoB-100 amino acid residues 1-1297). In this study, we examined the relationship of the MB19 polymorphism to a C----T nucleotide substitution at apoB cDNA nucleotide 421. This nucleotide substitution results in a Thr----Ile substitution at apoB-100 amino acid 71, and it changes an ApaLI restriction endonuclease site in the apoB gene. The nucleotide substitution was easily detectable by ApaLI digestion of a 141-base pair fragment of the apoB gene obtained by enzymatic amplification of genomic DNA. In 62 subjects, the MB19 phenotype, as determined by radioimmunoassays, correlated perfectly with the ApaLI restriction site polymorphism in the amplified DNA. The apoB allotype MB19(1) is associated with an Ile at residue 71 and the absence of the ApaLI site, whereas the apoB allotype MB19(2) is associated with a Thr at residue 71 and the presence of the ApaLI site. We conclude that the amino acid substitution at residue 71 probably accounts for the MB19 polymorphism in apoB.     

Influence of bromodeoxyuridine substitution of thymidine on sister-chromatid exchanges and chromosomal aberrations induced by restriction endonucleases

Different levels of replacement of thymidine by 5-bromodeoxyuridine in mammalian DNA have been used to analyze restriction endonuclease-dependent induction of sister-chromatid exchanges and chromosomal aberrations. Data regarding enzyme action in whole cells and in isolated nuclei are presented and discussed. The results indicate a lack of correlation between enzyme effectiveness and the degree of 5-bromodeoxyuridine substitution in the target sequences, specific to the tested restriction endonucleases.     

Identification of a novel transthyretin variant (Val30----Leu) associated with familial amyloidotic polyneuropathy.

A novel variant transthyretin which contains a leucine-for-valine substitution at position 30 was isolated and identified in the serum of a patient with familial amyloidotic polyneuropathy (FAP). The amino acid substitution was proven to result from a guanine-to-cytosine change at the first base of codon 30 located in exon 2 in the mutated transthyretin gene by restriction fragment length analysis on the amplified transthyretin gene using Cfr13 I. The study indicates that the point mutation of the transthyretin gene is a cause of the disorder.     

Generation and characterization of a stable full-length ecotropic murine leukemia virus molecular clone that produces novel phenotypes to Fv1 restriction

Retrovirus tropism can be restricted by host cell factors such as Fv1, TRIM5alpha, and Lv1 that inhibit infection by targeting the incoming viral capsid. The Fv1 gene inhibits murine leukemia virus infection in mice, but the precise mechanism of Fv1-mediated restriction is poorly understood. Our previous studies had demonstrated that Fv1-mediated viral tropism can be determined within the capsid protein at position 114 (Jung and Kozak. 2000. J. Virol. 74: 5385-7). To study the interaction between Fv1 and CA, we introduced amino acid substitution and deletion at this site in the N-tropic AKV capsid gene. The mutated two-LTR proviral DNAs were introduced into SC-1 cells by transfection. After transfection, cell supernatants collected from transfected cells were tested for host range susceptibility. The result indicated that substitution of amino acids did not alter tropism, but the deletion of 114His produced a virus with unusual tropism. The novel phenotype produced here failed to replicate in Fv1-expressing cells. This mutant virus showing such an extreme restriction pattern would be useful for studying the mechanism of Fv1- mediated restriction.     

New polymorphism of FASN gene in chicken

Sequencing, Polymerase Chain Reaction (PCR) and Restriction Fragment Length Polymorphism (RFLP) were carried out to detect polymorphism in the last intron of the FASN gene of the Campero broiler line. The analysis of the sequences presents a G to A substitution located at base 459 of the PCR product (GenBank accession number J02839, located at base Nr. 1222), resulting in a site recognized by restriction enzymes Hae III and Ava II. Also, the sequence presents a G to A substitution (located at base 603 of the PCR product and Nr. 1366 of the J02839 GenBank accession) resulting in a site recognized by restriction enzyme Pst I. Alleles and genotype frequencies were calculated for endonucleases Hae III, Ava II and Pst I for 44 broilers.     

Restriction endonuclease mapping of bacteriophage phi105 and closely related temperate Bacillus subtilis bacteriophages rho10 and rho14.

Cleavage maps of the three similar Bacillus subtilis temperate bacteriophages, phi105, rho10, and rho14, were constructed by partial digestion analysis utilizing the restriction endonuclease EcoRI. Comparison of the topography of these maps indicates that all phage DNAs posses cohesive ends and a number of EcoRI restriction sites; the fragments are conserved, and the estimated base substitution/nucleotide divergence between these phages is 0.03 to 0.07 based on conserved fragments or between 0.03 and 0.11 based on conserved cleavage sites. These lines of evidence indicate that phi105, rho10, and rho14 are closely related. Double-enzyme digestion analysis reveals that rho14 DNA has unique SalGI and BglII restriction sites and phi105 DNA has a unique SalGI restriction site, making these phages possible cloning vectors for B. subtilis.     

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     

An analysis of the action of nonspecific and restriction nucleases on metaphase chromosomes in situ

Patterns of differential staining of Drosophila, mouse, rat, cattle and pig chromosomes were examined after the treatment with nucleases (DNAase I, DNAase II) and restriction enzymes (AluI, HpaII, MspI, BpE, EcoRI). The above effects depend on the species used, on the enzymes and substitution of thymine for bromodeoxyuridine in the chromosomal DNA. It is supposed that such a phenomenon may not only result from the irregular distribution of specific restriction sites along chromosomes but also depend on the specificity of supramolecular organization of the chromosomal DNA.     

Genetic polymorphism of human plasma apolipoprotein A-IV is due to nucleotide substitutions in the apolipoprotein A-IV gene

Genetic polymorphism of human plasma apolipoprotein A-IV has been detected by isoelectric focusing techniques followed by immunoblotting. The molecular basis for this apoA-IV polymorphism has been elucidated. Analysis of the protein coding sequences of the apoA-IV alleles 1 and 2 revealed a single G to T substitution in the apoA-IV-2 allele. The point mutation, occurring in a region highly conserved among the mouse, rat, and human A-IV apolipoproteins, converts the glutamine at position 360 of the mature protein to a histidine. This amino acid substitution adds one positive charge unit to the apoA-IV-1 isoprotein (pI 4.97) thus creating the more basic apoA-IV-2 isoprotein (pI 5.02). Computer analysis of the apoA-IV-2 allele revealed that the single G to T substitution results in the loss of a BbvI and a Fnu4HI restriction enzyme site and in the formation of a new restriction site for the enzyme SfaNI. Protein primary and secondary structure predictions were largely unaffected by this amino acid exchange. These results on the structure of the apoA-IV-1 and apoA-IV-2 alleles suggest that the three other rare isoproteins (apoA-IV-0, apoA-IV-3, and apoA-IV-4) are also due to nucleotide and subsequent amino acid substitutions in the apoA-IV sequence.     

A method for inducing the point C----T transitions in DNA sequences corresponding to the ends of restriction sites

A new method for inducing of C----T substitutions into cytosine-containing restriction sites is developed. The method, based on the selective modification of cytosine residues in DNA sticky ends by sodium bisulfite, was illustrated by induction of a base substitution (C----T at the BamHI site of pBR322 plasmid DNA.     

Mammalian cell processing of a unique uracil residue in simian virus 40 DNA.

The processing of a unique uracil in DNA has been studied in mammalian cells. A synthetic oligodeoxyribonucleotide carrying a potential Bgl II restriction site, where one base has been substituted with a uracil, was inserted in the early intron of SV40 genome. Various heteroduplexes were constructed in such a manner that the restitution of an active Bgl II restriction site corresponds in each case to the specific substitution of the uracil by one of the four bases normally present in the DNA. DNA cuts by this restriction enzyme in one or several constructed heteroduplexes immediately determine the type of base pair substitution produced at the site of the U residue. When the uracil is inserted opposite a purine it is fully repaired; when facing a guanine it is replaced by a cytosine and opposite an adenine it is replaced by a thymine. These results indicate the error-free repair of uracil when it appears in the cell with the usual mechanisms such as cytosine deamination or incorporation of dUTP in place of dTTP during replication. When the uracil is inserted opposite a pyrimidine no error free repair at all is detected for U:C or U:T mismatches. It appears, moreover, that in approximately 18% of the cases U:T mismatch leads to a C:G base pairing. In the majority of the U:pyrimidine mismatches, mutations occur in the vicinity of the uracil, including base substitutions and frameshifts by addition of one or several bases.     

Comparison between Pyrococcus horikoshii and Pyrococcus abyssi genome sequences reveals linkage of restriction-modification genes with large genome polymorphisms

Recent work suggests that restriction-modification gene complexes are mobile genetic elements that insert themselves into the genome and cause various genome rearrangements. In the present work, the complete genome sequences of Pyrococcus horikoshii and Pyrococcus abyssi, two species in a genus of hyperthermophilic archaeon (archaebacterium), were compared to detect large genome polymorphisms linked with restriction-modification gene homologs. Sequence alignments, GC content analysis, and codon usage analysis demonstrated the diversity of these homologs and revealed a possible case of relatively recent acquisition (horizontal transfer). In two cases out of the six large polymorphisms identified, there was insertion of a DNA segment with a modification gene homolog, accompanied by target deletion (simple substitution). In two other cases, homologous DNA segments carrying a modification gene homolog were present at different locations in the two genomes (transposition). In both cases, substitution (insertion/deletion) in one of the two loci was accompanied by inversion of adjacent chromosomal segment. In the fifth case, substitution by a DNA segment carrying type I restriction, modification, and specificity gene homologs was likewise accompanied by adjacent inversion. In the last case, two homologous DNA segments, were found at different loci in the two genomes (transposition), but only one of them had insertion of a modification homolog and an unknown ORF. The possible relationship of these polymorphisms to attack by restriction enzymes on the chromosome will be discussed.     

Relative efficiencies of the maximum-parsimony and distance-matrix methods of phylogeny construction for restriction data.

The relative efficiencies of the maximum-parsimony (MP), UPGMA, and neighbor-joining (NJ) methods in obtaining the correct tree (topology) for restriction-site and restriction-fragment data were studied by computer simulation. In this simulation, six DNA sequences of 16,000 nucleotides were assumed to evolve following a given model tree. The recognition sequences of 20 different six-base restriction enzymes were used to identify the restriction sites of the DNA sequences generated. The restriction-site data and restriction-fragment data thus obtained were used to reconstruct a phylogenetic tree, and the tree obtained was compared with the model tree. This process was repeated 300 times. The results obtained indicate that when the rate of nucleotide substitution is constant the probability of obtaining the correct tree (Pc) is generally higher in the NJ method than in the MP method. However, if we use the average topological deviation from the model tree (dT) as the criterion of comparison, the NJ and MP methods are nearly equally efficient. When the rate of nucleotide substitution varies with evolutionary lineage, the NJ method is better than the MP method, whether Pc or dT is used as the criterion of comparison. With 500 nucleotides and when the number of nucleotide substitutions per site was very small, restriction-site data were, contrary to our expectation, more useful than sequence data. Restriction-fragment data were less useful than restriction-site data, except when the sequence divergence was very small. UPGMA seems to be useful only when the rate of nucleotide substitution is constant and sequence divergence is high.     

Genotyping of apolipoprotein A-IV by digestion of amplified DNA with restriction endonuclease Fnu4HI: use of a tailored primer to abolish additional recognition sites during the gene amplification.

Apolipoprotein A-IV (apoA-IV) is involved in the metabolism of chylomicrons and high density lipoproteins. It displays genetic polymorphism due to two co-dominant alleles apoA-IV1 and apoA-IV2. The mutation that causes the polymorphism is a G to T substitution in the third base of codon 360 in the apoA-IV2 allele which results in a glutamine (Gln) to histidine (His) change of amino acid 360. This substitution leads to the abolition of a recognition site for the restriction enzyme Fnu4HI. Part of the third exon of the apoA-IV gene is amplified by the polymerase chain reaction (PCR) using a tailored primer, which abolishes the downstream recognition sites during the DNA amplification. The PCR products are digested with the restriction enzyme Fnu4HI and electrophoresed on a polyacrylamide gel. The apoA-IV genotypes are determined after staining with either ethidium bromide or silver. To validate the method, we determined the inheritance of the apoA-IV alleles in a three-generation kindred of 8 subjects and analyzed amplified DNA of 32 subjects of different apoA-IV phenotypes with this method. The results were compared to those obtained from isoelectric focusing and immunoblotting. In all cases studied, the two methods gave concordant results.     

Enzymatic Cleavage of Type II Restriction Endonucleases on the 2′-O-Methyl Nucleotide and Phosphorothioate Substituted DNA

The effects of nucleotide analogue substitution on the cleavage efficiencies of type II restriction endonucleases have been investigated. Six restriction endonucleases (EcoRV, SpeI, XbaI, XhoI, PstI and SphI) were investigated respectively regarding their cleavage when substrates were substituted by 2′-O-methyl nucleotide (2′-OMeN) and phosphorothioate (PS). Substitutions were made in the recognition sequence and the two nucleotides flanking the recognition sequence for each endonuclease. The endonuclease cleavage efficiencies were determined using FRET-based assay. Results demonstrated a position-dependent inhibitory effect of substitution on the cleavage efficiency for all the six endonucleases. In general, the 2′-OMeN substitutions had greater impact than the PS substitutions on the enzymatic activities. Nucleotides of optimal substitutions for protection against RE cleavage were identified. Experimental results and conclusions in this study facilitate our insight into the DNA-protein interactions and the enzymatic cleavage mechanism, particularly for those whose detailed structure information is not available. In addition, the information could benefit the development of bioengineering and synthetic biology.