MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

Metal ion-mediated substrate-assisted catalysis in type II restriction endonucleases

The 2.15-Å resolution cocrystal structure of EcoRV endonuclease mutant T93A complexed with DNA and Ca2+ ions reveals two divalent metals bound in one of the active sites. One of these metals is ligated through an inner-sphere water molecule to the phosphate group located 3′ to the scissile phosphate. A second inner-sphere water on this metal is positioned approximately in-line for attack on the scissile phosphate. This structure corroborates the observation that the pro-SP phosphoryl oxygen on the adjacent 3′ phosphate cannot be modified without severe loss of catalytic efficiency. The structural equivalence of key groups, conserved in the active sites of EcoRV, EcoRI, PvuII, and BamHI endonucleases, suggests that ligation of a catalytic divalent metal ion to this phosphate may occur in many type II restriction enzymes. Together with previous cocrystal structures, these data allow construction of a detailed model for the pretransition state configuration in EcoRV. This model features three divalent metal ions per active site and invokes assistance in the bond-making step by a conserved lysine, which stabilizes the attacking hydroxide ion nucleophile.

Structure of FokI has implications for DNA cleavage

FokI is a member an unusual class of restriction enzymes that recognize a specific DNA sequence and cleave nonspecifically a short distance away from that sequence. FokI consists of an N-terminal DNA recognition domain and a C-terminal cleavage domain. The bipartite nature of FokI has led to the development of artificial enzymes with novel specificities. We have solved the structure of FokI to 2.3 Å resolution. The structure reveals a dimer, in which the dimerization interface is mediated by the cleavage domain. Each monomer has an overall conformation similar to that found in the FokI–DNA complex, with the cleavage domain packing alongside the DNA recognition domain. In corroboration with the cleavage data presented in the accompanying paper in this issue of Proceedings, we propose a model for FokI DNA cleavage that requires the dimerization of FokI on DNA to cleave both DNA strands.

Efficient construction of long DNA duplexes with internal non-Watson–Crick motifs and modifications

We have developed a semi-synthetic approach for preparing long stretches of DNA (>100 bp) containing internal chemical modifications and/or non-Watson–Crick structural motifs which relies on splint-free, cell-free DNA ligations and recycling of side-products by non-PCR thermal cycling. A double-stranded DNA PCR fragment containing a polylinker in its middle is digested with two restriction enzymes and a small insert (∼20 bp) containing the modification or non-Watson–Crick motif of interest is introduced into the middle. Incorrect products are recycled to starting materials by digestion with appropriate restriction enzymes, while the correct product is resistant to digestion since it does not contain these restriction sites. This semi-synthetic approach offers several advantages over DNA splint-mediated ligations, including fewer steps, substantially higher yields (∼60% overall yield) and ease of use. This method has numerous potential applications, including the introduction of modifications such as fluorophores and cross-linking agents into DNA, controlling the shape of DNA on a large scale and the study of non-sequence-specific nucleic acid–protein interactions.

DNA methyltransferases of the cyanobacterium Anabaena PCC 7120

From the characterization of enzyme activities and the analysis of genomic sequences, the complement of DNA methyltransferases (MTases) possessed by the cyanobacterium Anabaena PCC 7120 has been deduced. Anabaena has nine DNA MTases. Four are associated with Type II restriction enzymes (AvaI, AvaII, AvaIII and the newly recognized inactive AvaIV), and five are not. Of the latter, four may be classified as solitary MTases, those whose function lies outside of a restriction/modification system. The group is defined here based on biochemical and genetic characteristics. The four solitary MTases, DmtA/M.AvaVI, DmtB/M.AvaVII, DmtC/M.AvaVIII and DmtD/M.AvaIX, methylate at GATC, GGCC, CGATCG and rCCGGy, respectively. DmtB methylates cytosines at the N4 position, but its sequence is more similar to N6-adenine MTases than to cytosine-specific enzymes, indicating that it may have evolved from the former. The solitary MTases, appear to be of ancient origin within cyanobacteria, while the restriction MTases appear to have arrived by recent horizontal transfer as did five now inactive Type I restriction systems. One Mtase, M.AvaV, cannot reliably be classified as either a solitary or restriction MTase. It is structurally unusual and along with a few proteins of prokaryotic and eukaryotic origin defines a structural class of MTases distinct from all previously described.

NOMAD: a versatile strategy for in vitro DNA manipulation applied to promoter analysis and vector design.

Molecular analysis of complex modular structures, such as promoter regions or multi-domain proteins, often requires the creation of families of experimental DNA constructs having altered composition, order, or spacing of individual modules. Generally, creation of every individual construct of such a family uses a specific combination of restriction sites. However, convenient sites are not always available and the alternatives, such as chemical resynthesis of the experimental constructs or engineering of different restriction sites onto the ends of DNA fragments, are costly and time consuming. A general cloning strategy (nucleic acid ordered assembly with directionality, NOMAD; WWW resource locator http:@Lmb1.bios.uic.edu/NOMAD/NOMAD.htm l) is proposed that overcomes these limitations. Use of NOMAD ensures that the production of experimental constructs is no longer the rate-limiting step in applications that require combinatorial rearrangement of DNA fragments. NOMAD manipulates DNA fragments in the form of "modules" having a standardized cohesive end structure. Specially designed "assembly vectors" allow for sequential and directional insertion of any number of modules in an arbitrary predetermined order, using the ability of type IIS restriction enzymes to cut DNA outside of their recognition sequences. Studies of regulatory regions in DNA, such as promoters, replication origins, and RNA processing signals, construction of chimeric proteins, and creation of new cloning vehicles, are among the applications that will benefit from using NOMAD.

Z-DNA-forming sites within the human beta-globin gene cluster.

Agarose-encapsulated, metabolically active, permeabilized nuclei from human hematopoietic cell lines were tested for Z-DNA formation in the beta-globin gene cluster. Biotinylated monoclonal antibodies against Z-DNA were diffused into the nuclei and cross-linked to DNA with a 10-ns laser exposure at 266 nm. Following digestion with restriction enzymes, fragments that had formed Z-DNA were isolated. Seventeen regions with Z-DNA sequence motifs in the 73-kb region were studied by PCR amplification, and five were found in the Z conformation.

Restriction-modification systems as genomic parasites in competition for specific sequences.

Restriction-modification (RM) systems are believed to have evolved to protect cells from foreign DNA. However, this hypothesis may not be sufficient to explain the diversity and specificity in sequence recognition, as well as other properties, of these systems. We report that the EcoRI restriction endonuclease-modification methylase (rm) gene pair stabilizes plasmids that carry it and that this stabilization is blocked by an RM of the same sequence specificity (EcoRI or its isoschizomer, Rsr I) but not by an RM of a different specificity (PaeR7I) on another plasmid. The PaeR7I rm likewise stabilizes plasmids, unless an rm gene pair with identical sequence specificity is present. Our analysis supports the following model for stabilization and incompatibility: the descendants of cells that have lost an rm gene pair expose the recognition sites in their chromosomes to lethal attack by any remaining restriction enzymes unless modification by another RM system of the same specificity protects these sites. Competition for specific sequences among these selfish genes may have generated the great diversity and specificity in sequence recognition among RM systems. Such altruistic suicide strategies, similar to those found in virus-infected cells, may have allowed selfish RM systems to spread by effectively competing with other selfish genes.

Heterogeneity in molecular recognition by restriction endonucleases: osmotic and hydrostatic pressure effects on BamHI, Pvu II, and EcoRV specificity.

The cleavage specificity of the Pvu II and BamHI restriction endonucleases is found to be dramatically reduced at elevated osmotic pressure. Relaxation in specificity of these otherwise highly accurate and specific enzymes, previously termed "star activity," is uniquely correlated with osmotic pressure between 0 and 100 atmospheres. No other colligative solvent property exhibits a uniform correlation with star activity for all of the compounds tested. Application of hydrostatic pressure counteracts the effects of osmotic pressure and restores the natural selectivity of the enzymes for their canonical recognition sequences. These results indicate that water solvation plays an important role in the site-specific recognition of DNA by many restriction enzymes. Osmotic pressure did not induce an analogous effect on the specificity of the EcoRV endonuclease, implying that selective hydration effects do not participate in DNA recognition in this system. Hydrostatic pressure was found to have little effect on the star activity induced by changes in ionic strength, pH, or divalent cation, suggesting that distinct mechanisms may exist for these observed alterations in specificity. Recent evidence has indicated that BamHI and EcoRI share similar structural motifs, while Pvu II and EcoRV belong to a different structural family. Evidently, the use of hydration water to assist in site-specific recognition is a motif neither limited to nor defined by structural families.

Development of a DNA-based method for detection and identification of Phytophthora species

Phytophthora diseases cause major losses to agricultural and horticultural production in Australia and worldwide. Most Phytophthora diseases are soilborne and difficult to control, making disease prevention an important component of many disease management strategies. Detection and identification of the causal agent, therefore, is an essential part of effective disease management. This paper describes the development and validation of a DNA-based diagnostic assay that can detect and identify 27 different Phytophthora species. We have designed PCR primers that are specific to the genus Phytophthora. The resulting amplicon after PCR is subjected to digestion by restriction enzymes to yield a specific restriction pattern or fingerprint unique to each species. The restriction patterns are compared with a key comprising restriction patterns of type specimens or representative isolates of 27 different Phytophthora species. A number of fundamental issues, such as genetic diversity within and among species which underpin the development and validation of DNA-based diagnostic assays, are addressed in this paper.

Binding of small molecules to DNA junctions

Previous results in our laboratory suggest that the (CG) 4 segments whether present in a right-handed or a left-handed conformation form distinctive junctions with adjacent random sequences. These junctions and their associated sequences have unique structural and thermodynamic properties that may be recognized by DNA-binding molecules. This study probes these sequences by using the following small ligands: actinomycin D, 1,4-bis(((di(aminoethyl)amino)ethyl)amino)anthracene-9,10-dione, ametantrone, and tris(phenanthroline)ruthenium (II). These ligands may recognize the distinctive features associated to the (CG)4 segment and its junctions and thus interact preferentially near these sequences. Restriction enzyme inhibition assays were used to determine whether or not binding interactions took place, and to approximate locations of these interactions. These binding studies are first carried out using two small synthetic oligomers BZ-III and BZ-IV. The (5meCG)4 segment present in BZ-III adopts the Z-conformation in the presence of 50 m M Co(NH3)63+. In BZ-IV, the unmethylated (CG)4 segment changes to a non-B conformation in the presence of 50 m M Co(NH3)63+. BZ-IV, containing the (CG)4 segment, was inserted into a clone plasmid then digested with the restriction enzyme Hinf I to produce a larger fragment that contains the (CG)4 segment. The results obtained on the small oligomers and on the larger fragment for restriction enzyme Mbo I indicate that 1,4-bis(((di(aminoethyl)amino)ethyl)amino)anthracene-9,10-dione binds more efficiently at or near the (CG)4 segment. Restriction enzymes EcoRV, Sac I and Not I with cleavage sites upstream and downstream of the (CG)4 insert were used to further localize binding interactions in the vicinity of the (CG)4 insert. RNA polymerase activity was studied in a plasmid which contained the (CG)4 insert downstream from the promoter sites of SP6 and T7 RNA polymerases. Activities of these two polymerases were studied in the presence of each one of the ligands used throughout the study. Only actinomycin D and spider, which bind at or near the (CG)4 segment, alter the activities of SP6 and T7 RNA polymerases. Surprisingly, enhancement of polymerase activity was observed in the presence of very low concentrations of actinomycin D. These results suggest that the conformational features of (CG) segments may serve in regulatory functions of DNA. ^

RestrictionDigest: A powerful Perl module for simulating genomic restriction digests

Background: Reduced-representation sequencing technology iswidely used in genotyping for its economical and efficient features. A popular way to construct the reduced-representation sequencing libraries is to digest the genomic DNA with restriction enzymes. A key factor of this method is to determine the restriction enzyme(s). But there are few computer programs which can evaluate the usability of restriction enzymes in reduced-representation sequencing. SimRAD is an R package which can simulate the digestion of DNA sequence by restriction enzymes and return enzyme loci number as well as fragment number. But for linkage mapping analysis, enzyme loci distribution is also an important factor to evaluate the enzyme. For phylogenetic studies, comparison of the enzyme performance across multiple genomes is important. It is strongly needed to develop a simulation tool to implement these functions. Results: Here, we introduce a Perl module named RestrictionDigest with more functions and improved performance. It can analyze multiple genomes at one run and generate concise comparison of enzyme performance across the genomes. It can simulate single-enzyme digestion, double-enzyme digestion and size selection process and generate comprehensive information of the simulation including enzyme loci number, fragment number, sequences of the fragments, positions of restriction sites on the genome, the coverage of digested fragments on different genome regions and detailed fragment length distribution. Conclusions: RestrictionDigest is an easy-to-use Perl module with flexible parameter settings.With the help of the information produced by the module, researchers can easily determine the most appropriate enzymes to construct the reduced-representation libraries to meet their experimental requirements.

Characterization by restriction fragment length polymorphism and sequence analysis of field and vaccine strains of infectious laryngotracheitis virus involved in severe outbreaks

At the end of 2002 and throughout 2003, there was a severe outbreak of infectious laryngotracheitis (ILT) in an intensive production area of commercial hens in the Sao Paulo State of Brazil. ILT virus was isolated from 28 flocks, and 21 isolates were genotyped by polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP) using four genes and eight restriction enzymes, and by partial sequencing of the infected cell protein 4 (ICP4) and thymidine kinase (TK) genes. Three groups resulted from the combinations of PCR-RFLP patterns: 19 field isolates formed Group I, and the remaining two isolates together with the chicken embryo origin (CEO) vaccine strains formed Group II. Group III comprised the tissue-culture origin (TCO) vaccine strain by itself. The PCR-RFLP results agreed with the sequencing results of two ICP4 gene fragments. The ICP4 gene sequence analysis showed that the 19 field isolates classified into Group I by RFLP-PCR were identical among themselves, but were different to the TCO and CEO vaccines. The two Group II isolates could not be distinguished from one of the CEO vaccines. The nucleotide and amino acid sequence analyses discriminated between the Brazilian and non-Brazilian isolates, as well as between the TCO and CEO vaccines. Sequence analysis of the TK gene enabled classification of the field isolates (Group I) as virulent and non-vaccine. This work shows that the severe ILT outbreak was caused by a highly virulent, non-vaccine strain.

Constraint-based modeling of minimum set covering: application to species differentation

Work presented in the context of the European Master in Computational Logics, as partial requisit for the graduation as Master in Computational Logics

Mapping of a Leishmania major gene/locus that confers pentamidine resistance by deletion and insertion of transposable element

Pentamidine (PEN) is an alternative compound to treat antimony-resistant leishmaniasis patients, which cellular target remains unclear. One approach to the identification of prospective targets is to identify genes able to mediate PEN resistance following overexpression. Starting from a genomic library of transfected parasites bearing a multicopy episomal cosmid vector containing wild-type Leishmania major DNA, we isolated one locus capable to render PEN resistance to wild type cells after DNA transfection. In order to map this Leishmania locus, cosmid insert was deleted by two successive sets of partial digestion with restriction enzymes, followed by transfection into wild type cells, overexpression, induction and functional tests in the presence of PEN. To determine the Leishmania gene related to PEN resistance, nucleotide sequencing experiments were done through insertion of the transposon Mariner element of Drosophila melanogaster (mosK) into the deleted insert to work as primer island. Using general molecular techniques, we described here this method that permits a quickly identification of a functional gene facilitating nucleotide sequence experiments from large DNA fragments. Followed experiments revealed the presence of a P-Glycoprotein gene in this locus which role in Leishmania metabolism has now been analyzed.