An algorithmic approach to multiple complete digest mapping.

Multiple Complete Digest (MCD) mapping is a method of determining the locations of restriction sites along a target DNA molecule. The resulting restriction map has many potential applications in DNA sequencing and genetics. In this work, we present a heuristic algorithm for fragment identification, a key step in the process of constructing an MCD map. Given measurements of the restriction fragment sizes from one or more complete digestions of each clone in a clone library covering the molecule to be mapped, the algorithm identifies groups of restriction fragments on different clones that correspond to the same region of the target DNA. Once these groups are correctly determined the desired map can be constructed by solving a system of simple linear inequalities. We demonstrate the effectiveness of our algorithm on real data provided by the Genome Center at the University of Washington.     

A rapid method to identify cosmids containing rare restriction sites.

A procedure for identifying specific cosmid clones containing recognition sites for "rare cutting" restriction enzymes has been developed. Cosmid clones containing human inserts were selected by hybridisation to human repetitive DNA. An oligonucleotide corresponding to the NotI recognition site, eight bases long, was labelled and used to probe DNA samples from one hundred cosmids. By optimising the difference in melting characteristics between eight-base perfect match and six-base match/two base mismatch hybrids, we were able to detect the cosmids containing either NotI (8 bp match) or XmaIII/EagI (6 bp match) sites. The generation of a map for rare cutter sites along a human chromosome, or a chromosome region, should be simplified using this approach, which will enable the identification of a set of "milestones" at intervals of several hundred kilobases (kb) along the DNA sequence.     

Algorithms for optical mapping.

Optical mapping is a novel technique for determining the restriction sites on a DNA molecule by directly observing a number of partially digested copies of the molecule under a light microscope. The problem is complicated by uncertainty as to the orientation of the molecules and by erroneous detection of cuts. In this paper we study the problem of constructing a restriction map based on optical mapping data. We give several variants of a polynomial reconstruction algorithm, as well as an algorithm that is exponential in the number of cut sites, and hence is appropriate only for small number of cut sites. We give a simple probabilistic model for data generation and for the errors and prove probabilistic upper and lower bounds on the number of molecules needed by each algorithm in order to obtain a correct map, expressed as a function of the number of cut sites and the error parameters. To the best of our knowledge, this is the first probabilistic analysis of algorithms for the problem. We also provide experimental results confirming that our algorithms are highly effective on simulated data.     

'DNA Strider': a 'C' program for the fast analysis of DNA and protein sequences on the Apple Macintosh family of computers

DNA Strider is a new integrated DNA and Protein sequence analysis program written with the C language for the Macintosh Plus, SE and II computers. It has been designed as an easy to learn and use program as well as a fast and efficient tool for the day-to-day sequence analysis work. The program consists of a multi-window sequence editor and of various DNA and Protein analysis functions. The editor may use 4 different types of sequences (DNA, degenerate DNA, RNA and one-letter coded protein) and can handle simultaneously 6 sequences of any type up to 32.5 kB each. Negative numbering of the bases is allowed for DNA sequences. All classical restriction and translation analysis functions are present and can be performed in any order on any open sequence or part of a sequence. The main feature of the program is that the same analysis function can be repeated several times on different sequences, thus generating multiple windows on the screen. Many graphic capabilities have been incorporated such as graphic restriction map, hydrophobicity profile and the CAI plot- codon adaptation index according to Sharp and Li. The restriction sites search uses a newly designed fast hexamer look-ahead algorithm. Typical runtime for the search of all sites with a library of 130 restriction endonucleases is 1 second per 10,000 bases. The circular graphic restriction map of the pBR322 plasmid can be therefore computed from its sequence and displayed on the Macintosh Plus screen within 2 seconds and its multiline restriction map obtained in a scrolling window within 5 seconds.     

Baboon endogenous virus genome. I. Restriction enzyme map of the unintegrated DNA genome of a primate retrovirus

A detailed restriction map was deduced for the genome of an endogenous retrovirus of a higher primate, that of baboon. The cleavage sites for 12 restriction enzymes were mapped. The unintegrated linear viral DNA intermediate that is produced by infection of permissive cells with baboon endogenous virus was isolated. Hybridization with a strong-stop complementary DNA probe demonstrated presence of a terminal repetition in the linear viral DNA. The positions of restriction sites for two particular enzymes, SmaI and XhoI, near each end were consistent with this result and indicated that the length of the repetition is 0.55 +/- 0.01 kilobase. The linear viral DNA had a unique restriction map indicating that it is not a set of random circular permutations of the RNA genome. From hybridization with a 3'-specific probe, the DNA restriction map was aligned relative to the 5'-to-3' orientation of the viral RNA. We observed a minor heterogeneity in a BamHI recognition site 1.95 kilobases from the right end of the linear map.     

Derivation of a restriction map of bacteriophage T3 DNA and comparison with the map of bacteriophage T7 DNA.

The DNA of bacteriophage T3 was characterized by cleavage with seven restriction endonucleases. AvaI, XbaI, BglII, and HindIII each cut T3 DNA at 1 site, KpnI cleaved it at 2 sites, MboI cleaved it at 9 sites, and HpaI cleaved it at 17 sites. The sizes of the fragments produced by digestion with these enzymes were determined by using restriction fragments of T7 DNA as molecular weight standards. As a result of this analysis, the size of T3 DNA was estimated to be 38.74 kilobases. The fragments were ordered with respect to each other and to the genetic map to produce a restriction map of T3 DNA. The location and occurrence of the restriction sites in T3 DNA are compared with those in the DNA of the closely related bacteriophage T7.     

Restriction endonuclease mapping of linear unintegrated proviral DNA of bovine leukemia virus.

A detailed restriction map was deduced for the genome of the exogenous bovine leukemia virus. The cleavage sites for nine restriction enzymes were mapped. The unintegrated linear viral DNA intermediate that is produced by infection of permissive cells with bovine leukemia virus was isolated. The linear viral DNA had a unique restriction map, indicating that it is not a set of random circular permutations of the RNA genome. From hybridization with a 3'-enriched probe, the DNA restriction map was aligned relative to the 5'-to-3' orientation of the viral RNA. Restriction enzyme analysis of integrated bovine leukemia virus information present in animals with enzootic bovine leukosis provided evidence for the existence of genetic variants of the virus.     

The distribution of restriction enzyme sites in Escherichia coli.

A statistical analysis of physical map data for eight restriction enzymes covering nearly the entire genome of E. coli is presented. The methods of analysis are based on a top-down modeling approach which requires no knowledge of the statistical properties of the base sequence. For most enzymes, the distribution of mapped sites is found to be fairly homogeneous. Some heterogeneity in the distribution of sites is observed for the enzymes Pstl and HindIII. In addition, BamHI sites are found to be more evenly dispersed than we would expect for random placement and we speculate on a possible mechanism. A consistent departure from a uniform distribution, observed for each of the eight enzymes, is found to be due to a lack of closely spaced sites. We conclude from our analysis that this departure can be accounted for by deficiencies in the physical map data rather than non-random placement of actual restriction sites. Estimates of the numbers of sites missing from the map are given, based both on the map data itself and on the site frequencies in a sample of sequenced E. coli DNA. We conclude that 5 to 15% of the mapped sites represent multiple sites in the DNA sequence.     

DNA restriction dependent on two recognition sites: activities of the SfiI restriction-modification system in Escherichia coli

In contrast to many type II restriction enzymes, dimeric proteins that cleave DNA at individual recognition sites 4-6 bp long, the SfiI endonuclease is a tetrameric protein that binds to two copies of an elongated sequence before cutting the DNA at both sites. The mode of action of the SfiI endonuclease thus seems more appropriate for DNA rearrangements than for restriction. To elucidate its biological function, strains of Escherichia coli expressing the SfiI restriction-modification system were transformed with plasmids carrying SfiI sites. The SfiI system often failed to restrict the survival of a plasmid with one SfiI site, but plasmids with two or more sites were restricted efficiently. Plasmids containing methylated SfI sites were not restricted. No rearrangements of the plasmids carrying SfiI sites were detected among the transformants. Hence, provided the target DNA contains at least two recognition sites, SfiI displays all of the hallmarks of a restriction-modification system as opposed to a recombination system in E. coli cells. The properties of the system in vivo match those of the enzyme in vitro. For both restriction in vivo and DNA cleavage in vitro, SfiI operates best with two recognition sites on the same DNA.     

Complete physical map of the Bacillus subtilis 168 chromosome constructed by a gene-directed mutagenesis method

All the SfiI sites and most of the NotI sites were located precisely on the chromosome of Bacillus subtilis 168 by a novel method, termed gene-directed mutagenesis. The stepwise elimination of these restriction sites by this method allowed not only the physical connection of the restriction fragments but also the accurate determination of the position of the restriction sites themselves. The resulting physical map of the 4165 x 10(3) base-pair B. subtilis chromosome has been correlated with the genetic map by determination of the exact location of known genes. The complete physical map provides a rapid and accurate way for mapping of new genes as well as analysis of large DNA rearrangements on the chromosome. The novel strategy is, in principle, applicable to the analysis of the genome of other organisms.     

Restriction endonuclease map of Euglena gracilis chloroplast DNA.

A physical map of the Euglena gracilis chloroplast genome has been constructed, based on cleavage sites of Euglena gracilis chloroplast DNA treated with bacterial restriction endonucleases. Covalently close, circular chloroplast DNA is cleaved by restriction endonuclease SalI into three fragments and by restriction endonuclease BamHI into six fragments. These nine cleavage sites have been ordered by fragment molecular weight analysis, double digestions, partial digestions, and by digestion studies of isolated DNA fragments. A fragment pattern of the products of EcoRI restriction endonuclease digestion of Euglena chloroplast DNA is also described. One of these fragments has been located on the cleavage site map.     

Cleavage of adeno-associated virus DNA with Sali,Psti and Haeii restriction endonucleases.

Duplex AAV-2 DNA was digested with SalI, PstI or HaeII restriction endonucleases and the cleavage sites were mapped. SalI cleaves AAV DNA at 0.310 map units, PstI at 0.106, 0.422 and 0.914 and the five HaeII sites were mapped at 0.110. 0.156, 0.181, 0.536 and 0.600 map units. These cleavage products will be useful for the isolation of specific regions from the AAV DNA, located outside of the stably transcribed region of the genome, and will also help to map more complex restriction enzyme cleavages.     

AFLP: a new technique for DNA fingerprinting.

A novel DNA fingerprinting technique called AFLP is described. The AFLP technique is based on the selective PCR amplification of restriction fragments from a total digest of genomic DNA. The technique involves three steps: (i) restriction of the DNA and ligation of oligonucleotide adapters, (ii) selective amplification of sets of restriction fragments, and (iii) gel analysis of the amplified fragments. PCR amplification of restriction fragments is achieved by using the adapter and restriction site sequence as target sites for primer annealing. The selective amplification is achieved by the use of primers that extend into the restriction fragments, amplifying only those fragments in which the primer extensions match the nucleotides flanking the restriction sites. Using this method, sets of restriction fragments may be visualized by PCR without knowledge of nucleotide sequence. The method allows the specific co-amplification of high numbers of restriction fragments. The number of fragments that can be analyzed simultaneously, however, is dependent on the resolution of the detection system. Typically 50-100 restriction fragments are amplified and detected on denaturing polyacrylamide gels. The AFLP technique provides a novel and very powerful DNA fingerprinting technique for DNAs of any origin or complexity.     

Restriction map of the temperate phage E105 from the polylysogenic culture of Erwinia carotovora 268

DNA structure of the temperate bacteriophage E105 from polylysogenic culture of Erwinia carotovora 268 has been studied. The viral 29.12-29.17 MD DNA has been shown to be linear and nonpermuted. The complete restriction map of the viral DNA has been constructed for MvaI and HpaI and partial for Eco31 restriction endonucleases based on the pair hydrolysis of the native DNA as well as its fragments. Altogether, 19 sites for restriction endonucleases have been localized on bacteriophage DNA.     

Restriction site map of the Chlorella virus PBCV-1 genome

The virus PBCV-1, which replicates in a Chlorella-like green alga, has a dsDNA genome. The DNA was mapped for BamHI, HindIII, and PstI restriction sites. The resulting map has a size of 333 kbp and is circular—indicating either covalently closed circular DNA or circularly permuted linear DNA. Several regions of repetitive DNA were also identified and located on the restriction map.     

Size and physical map of the chromosome of Haemophilus influenzae.

A variation of pulse-field electrophoresis, field-inversion gel electrophoresis, was used to determine the size and physical map of the chromosome of Haemophilus influenzae. The DNA of H. influenzae had a low G + C content (39%) and no restriction sites for the enzymes NotI or SfiI. However, a number of restriction enzymes (SmaI, ApaI, NaeI, and SacII) that recognized 6-base-pair sequences containing only G and C nucleotides were found to generate a reasonable number of DNA fragments that were separable in agarose gels by field-inversion gel electrophoresis. The sizes of the DNA fragments were calibrated with a lambda DNA ladder and lambda DNA restriction fragments. The sum of fragment sizes obtained with restriction digests yielded a value for the chromosome of 1,980 kilobase pairs. Hybridization of a labeled fragment with two or more fragments from a digest with a different restriction enzyme provided the information needed to construct a circular map of the H. influenzae chromosome.     

Physical map of PM2 DNA.

The cleavage sites of the restriction nucleases HpaI, HpaII, HindII, HindIII, and PstI have been mapped on the DNA of the bacteriophage PM2. This map has been used for the localization of two strong binding sites of Escherichia coli RNA polymerase on PM2 DNA.     

鳞翅目基因组IIB型内切酶位点的统计分析

IIB型限制内切酶能够识别并切割特异酶切位点两端特定距离的DNA,形成粘性末端的30 bp左右的等长DNA片段。利用其特性与限制性酶切位点关联测序技术(RAD)相结合发展出2b-RAD简化基因组测序技术,应用于遗传图谱构建、种群遗传结构分析、性状定位以及细菌分型等多种研究领域。构建2b-RAD测序文库之前,需要对基因组中的IIB型限制内切酶位点进行预测与统计分析,制定有效的测序文库构建方案。本文利用Python语言构建分析基因组中IIB型限制内切酶位点的流程,预测并统计6个鳞翅目代表物种基因组含有的8个商业化IIB型限制内切酶的酶切位点,比较了各个基因组与IIB型限制内切酶之间含有的酶切位点总量、重复序列数量以及酶切间隔长度的关系,为在昆虫基因组中进一步试行2b-RAD研究提供了参考。     

Alignment of Sfi I sites with the Not I restriction map of Schizosaccharomyces pombe genome.

A Sfi I restriction map of the fission yeast Schizosaccharomyces pombe genome was aligned with the Not I restriction map. There are 16 Sfi I sites in the S. pombe genome. Three Sfi I sites are on chromosome III which is devoid of Not I sites. The sizes of the entire genome and individual chromosomes, calculated from the Sfi I fragment sizes, are consistent with that calculated from the Not I fragment sizes. The Sfi I map provides greater physical characterization of the S. pombe genome and further validates the use of S. pombe chromosomal DNA as size standard. These maps have allowed detection of polymorphism on all three chromosomes.     

Purification of BsuE methyltransferase and its application in genome mapping.

We have used a combination of BsuE methyltransferase (M-BsuE) and NotI restriction enzyme to cut genomic DNA at a subset of NotI sites. The usefulness of this system is shown in a re-examination of the restriction map of the human MHC. Combinations of methylases and restriction enzymes can be used to generate cuts at different frequencies in genomic DNA, such that they generate ends complementary to NotI ends, and can be used in conjunction with NotI linking clones in chromosome jumping experiments. These enzyme combinations have the potential to produce cutting sites in genomic DNA spaced at intervals favorable for extensive mapping, fragment enrichment, and cloning efforts.