In Silico Analysis of the Restriction Fragment Length Distribution in the Human Genome

The Restriction On Computer (ROC) program (freely available at http://www.mcb.harvard.edu/ gilbert/ROC) was developed and used to analyze the restriction fragment length distribution in the human genome. In contrast to other programs searching for restriction sites, ROC simultaneously analyzes several long nucleotide sequences, such as the entire genomes, and in essence simulates electrophoretic analysis of DNA restriction fragments. In addition, this program extracts and analyzes DNA repeats that account for peaks in the restriction fragment length distribution. The ROC analysis data are consistent with the experimental data obtained via in vitro restriction enzyme analysis (DNA taxonoprint). A difference between the in vitro and in silico results is explained by underrepresentation of tandem DNA repeats in genomic databases. The ROC analysis of individual genome fragments elucidated the nature of several DNA markers, which were earlier revealed by DNA taxonoprint, and showed that L1 and Alurepeats are nonrandomly distributed in various chromosomes. Another advantage is that the ROC procedure makes it possible to analyze the nonrandom character of a genomic distribution of short DNA sequences. The ROC analysis showed that a low poly(G) frequency is characteristic of the entire human genome, rather than of only coding sequences. The method was proposed for a more complex in silico analysis of the genome. For instance, it is possible to simulate DNA restriction together with blot hybridization and then to analyze the nature of markers revealed.     

In silico analysis of the restriction fragments length distribution in the human genome.]

The Restriction On Computer (ROC) program (freely available at http://www.mcb.harvard.edu/gilbert/ROC) was developed and used to analyze the restriction fragment length distribution in the human genome. In contrast to other programs searching for restriction sites, ROC simultaneously analyzes several long nucleotide sequences, such as the entire genomes, and in essence simulates electrophoretic analysis of DNA restriction fragments. In addition, this program extracts and analyzes DNA repeats that account for peaks in the restriction fragment length distribution. The ROC analysis data are consistent with the experimental data obtained via in vitro restriction enzyme analysis (taxonomic printing). A difference between the in vitro and in silico results is explained by underrepresentation of tandem DNA repeats in genomic databases. The ROC analysis of individual genome fragments elucidated the nature of several DNA markers, which were earlier revealed by taxonomic printing, and showed that L1 and Alu repeats are nonrandomly distributed in various chromosomes. Another advantage is that the ROC procedure makes it possible to analyze the nonrandom character of a genomic distribution of short DNA sequences. The ROC analysis showed that a low poly(G) frequency is characteristic of the entire human genome, rather than of only coding sequences. The method was proposed for a more complex in silico analysis of the genome. For instance, it is possible to simulate DNA restriction together with blot hybridization and then to analyze the nature of markers revealed.     

ELISA detection of restriction site polymorphisms in the pig ryanodine receptor locus

We compare two strategies for ELISA detection of restriction site polymorphisms (EDRSP) that are suitable for high-throughput genotyping of the pig ryanodine receptor point mutation (RYR1 ^hal ). In both procedures, target DNA is amplified by PCR with one primer that is 5′ biotinylated and a second primer that is 5′ fluoresceinylated. PCR products are captured in duplicate wells on a streptavidin-coated, 96-well plate. The duplicates may be treated in two ways. In a single restriction enzyme assay, one duplicate is exposed to a restriction enzyme that cuts one allele specifically, and the second duplicate is exposed to no restriction enzyme. In a dual restriction enzyme assay, the second replicate is exposed to a second restriction enzyme that cuts the alternate allele specifically. Thereafter, the two procedures are similar; anti-fluorescein antibodies conjugated to peroxidase are allowed to bind to the fluoresceinylated ends, the plate is washed, and a substrate is converted to a colored end product. The ratio of the absorbances in the two wells is used to classify subjects by genotype. When the dual restriction enzyme assay is run, three genotype groups are easily distinguishable. When the single restriction enzyme assay is run, heterozygotes generate values that may overlap with those of the homozygotes that are not cut by the restriction enzyme. Dual restriction enzyme assays are more accurate than single restriction enzyme assays; however, single restriction enzyme assays are sufficient for identifying pigs that carry RYR1 ^hal . Received: 30 December 1997 / Accepted: 20 April 1998     

Mapping genes within a YAC by computer-assisted interpretation of partial restriction digestions.

Partial restriction digestion is used to map restriction sites and the location of genes within yeast artificial chromosomes (YACs). Locus-specific probes are hybridised to the partially digested YAC DNA and the fragments to which they hybridise are compared with the pattern of partial digestion products that include each map region. A least squares criterion is presented which allows for error in fragment length determination. This rapidly defines the most likely location of a marker within the restriction map and permits the combination of results from digestions with different restriction enzymes. Approximate confidence intervals may be assigned to gene locations, and tests of goodness-of-fit of the data may be performed. Since the number of erroneously matched fragments increases in proportion to the square of the number of sites, denser maps are not necessarily more informative. Simulations indicate that the optimal number of internal restriction sites given typical experimental error (1% of YAC length) is about five sites; the associated broad support interval (on average one third of YAC length) may be reduced by combining results from different enzyme digestions. Application of a computer implementation of this model to experimental data showed that the model fitted well, and estimates of location were found to be consistent with other evidence.     

Estimating the fraction of clonable genomic DNA

Tang and Waterman (1990, Bulletin of Mathematical Biology, 52, 455–475) give an estimate of the fraction of clonable genomic DNA by complete digestion with a single restriction enzyme. However, their argument neglects the ‘end effects’ of restriction enzyme digestion. This effectively assumes that the DNA is of infinite length. However, neglecting the end effects is not necessary to obtain useful estimates of the fraction of clonable DNA. An alternative approach is taken to estimating the fraction of clonable DNA that accounts for the finite nature of the DNA strand. This approach yields the Tang and Waterman estimate by allowing the genome length to go to infinity. Comparisons based on simulation studies show the finite approach to be more accurate than Tang and Waterman’s estimate, especially for shorter genomes. Sequence-based studies also show that the proposed method is often superior to the Tang-Waterman estimate.     

Complexes formed between the restriction endonuclease EcoK and heteroduplex DNA

The complexes between the Escherichia coli K restriction endonuclease and heteroduplex DNA (one strand methylated and one unmethylated) have been characterized and shown to have different properties from those formed with unmodified DNA. The nature of the heteroduplex complex appears to commit the enzyme to its methylase mode.     

Restriction fragment length polymorphisms at the human parathyroid hormone gene locus

Summary Two common Pst I and Taq I restriction enzyme fragment length polymorphisms (RFLPs) were detected at the human parathyroid hormone (PTH) gene locus. The allele frequencies in a Northern German population were 0.578/0.422 (Pst I) and 0.628/0.372 (Taq I). The allele distributions follow Hardy-Weinberg expectations of equilibrium in the population. The Mendelian nature of the polymorphisms were confirmed in family studies.     

Designation by restriction fragment length polymorphism of major histocompatibility complex class IV haplotypes in meat-type chickens

Major histocompatiblity complex (MHC) class IV haplotypes were identified in a population of meat-type chickens by restriction fragment length polymorphism (RFLP) analysis. Fourteen different haplotypes were designated on the basis of restriction patterns obtained from Southern blots of PvuII- or BglII-digested DNA, hybridized with the MHC class IV cDNA probe bg32.1. Digestion with each restriction enzyme yielded the same level of polymorphism among individuals. For each haplotype, 4–10 restriction fragments ranging from 0–8 to 8 kb were observed. Such a designation of meat-type chicken MHC class IV haplotypes enables a rapid recognition of previously defined haplotypes, is readily adjustable to additional, newly found restriction patterns and could prove useful in practical breeding programmes.     

Morphological and biochemical effects of endonucleases on isolated mammalian chromosomes in vitro

Endonuclease digestion of isolated and unfixed mammalian metaphase chromosomes in vitro was examined as a means to study the higher-order regional organization of chromosomes related to banding patterns and the mechanisms of endonuclease-induced banding. Isolated mouse LM cell chromosomes, digested with the restriction enzymes AluI, HaeIII, EcoRI, BstNI, AvaII, or Sau96I, demonstrated reproducible G- and/or C-banding at the cytological level depending on the enzyme and digestion conditions. At the molecular level, specific DNA alterations were induced that correlated with the banding patterns produced. The results indicate that: (1) chromatin extraction is intimately involved in the mechanism of endonuclease induced chromosome banding. (2) The extracted DNA fragments are variable in size, ranging from 200 bp to more than 4 kb in length. (3) For HaeIII, there appears to be variation in the rate of restriction site cleavage in G- and R-bands; HaeIII sites appear to be more rapidly cleaved in R-bands than in G-bands. (4) AluI and HaeIII ultimately produce banding patterns that reflect regional differences in the distribution of restriction sites along the chromosome. (5) BstNI restriction sites in the satellite DNA of constitutive heterochromatin are not cleaved intrachromosomally, probably reflecting an inaccessibility of the BstNI sites to enzyme due to the condensed nature of this chromatin or specific DNA-protein interactions. This implies that some enzymes may induce banding related to regional differences in the accessibility of restriction sites along the chromosome. (6) Several specific nonhistone protein differences were noted in the extracted and residual chromatin following an AluI digestion. Of these, some nonhistones were primarily detected in the extracted chromatin while others were apparently resistant to extraction and located principally in the residual chromatin. (7) The chromatin in constitutive heterochromatin is transiently resistant to cleavage by micrococcal nuclease.     

Mouse genomic DNA sequences homologous to sea urchin TU elements are genetically stable polydispersed repeats useful for analysis of multiple RFLPs.

Segments of the murine genome that hybridize to the inverted repeat regions of the transposable TU elements of sea urchins include tandem repeats of a sequence (CTCC) that encodes the recognition site for the restriction enzyme Mnl1, as do the analogous polypurine/polypyrimidine (pPu/pPy) stretches in humans. The Mnl1-sensitive repeats, which exist as a microsatellite sequence 200-300 bp in length, lack the terminal dyad symmetry characteristic of the TU elements and are structurally and functionally distinct from these elements. DNA fragments containing these repeat units that are isolated from different generations of isogenic (or congenic) mice or from different tissues of genetically identical individuals are indistinguishable by RFLP analysis; however, they show restriction fragment length polymorphism in different strains. This polymorphism appears to reflect DNA sequence changes occurring at sites flanking the repeats rather than variability in the number of repeats. Their genetic stability and occurrence in a wide variety of animal species make the Mnl1 repeats useful in studying genetic variation that has occurred over an evolutionary time scale of greater duration than can be examined conveniently by VNTR analysis.     

Development of a fluorophore-ribosomal DNA restriction typing method for monitoring structural shifts of microbial communities

DNA restriction fragment polymorphism technologies such as amplified ribosomal DNA restriction analysis (ARDRA) and terminal restriction fragment length polymorphism (T-RFLP) have been widely used in investigating microbial community structures. However, these methods are limited due to either the low resolution or sensitivity. In this study, a fluorophore-ribosomal DNA restriction typing (f-DRT) approach is developed for structural profiling of microbial communities. 16S rRNA genes are amplified from the community DNA and digested by a single restriction enzyme Msp I. All restriction fragments are end-labeled with a fluorescent nucleotide Cy5-dCTP via a one-step extension reaction and detected with an automated DNA sequencer. All 50 predicted restriction fragments between 100 and 600 bp were detected when twelve single 16S rRNA gene sequences were analyzed using f-DRT approach; 92% of these fragments were determined with accuracy of ±2 bp. In the defined model communities containing five components with different ratios, relative abundance of each component was correctly revealed by this method. The f-DRT analysis also showed structural shifts of intestinal microbiota in carcinogen-treated rats during the formation of precancerous lesions in the colon, as sensitive as multiple digestion-based T-RFLP analysis. This study provides a labor and cost-saving new method for monitoring structural shifts of microbial communities.     

Comparison of the ribosomal RNA genes in four closely relatedCucurbitaceae

Cucurbitaceae are characterized by a high copy number for nuclear ribosomal RNA genes. We have investigated the genomic ribosomal DNA (rDNA) of four closely related species of this family with respect to structure, length heterogeneity, and evolution. InCucumis melo (melon) there are two main length variants of rDNA repeats with 10.7 and 10.55kb.Cucumis sativus (cucumber) shows at least three repeat types with 11.5, 10.5, and 10.2kb.Cucurbita pepo (zucchini) has two different repeat types with 10.0 and 9.3kb. There are also two different repeat types inCucurbita maxima (pumpkin) of about 11.2 and 10.5kb. Restriction enzyme mapping of the genomic rDNA of these four plants and of cloned repeats ofC. sativus shows further heterogeneities which are due to methylation or point mutations. By comparison of the restriction enzyme maps it was possible to trace some evolutionary events in the family ofCucurbitaceae. Some aspects of regulation and function of the middle repetitive rRNA genes (here between 2000 and 10000 copies) are discussed.     

Genomic analysis of a virulent and a less virulent strain of the entomopathogenic fungus Beauveria bassiana, using restriction fragment length polymorphisms.

The genomic DNA of two strains of the entomopathogenic fungus Beauveria bassiana, strain GK2016, a "wild type" (virulent), and strain GK2051, a less virulent mutant to grasshoppers, was digested with 12 restriction endonucleases. Gel electrophoresis conditions were established to show restriction fragment length patterns visually in the digested DNA stained with ethidium bromide. The less virulent mutant was generated by ultraviolet illumination of conidiospores at a 95% lethal dose. Both strains of the fungi were identical in morphology as well as in 16 of 22 API-ZYM kit enzyme assays. Differences in levels of total enzyme activity were observed for esterase, esterase-lipase, beta-galactosidase, chitinase, and protease. A Neurospora crassa beta-tubulin gene (heterologous gene) and two homologous DNA probes (pJK16 and pJK18) hybridized to several specific DNA bands in B. bassiana strain GK2016 but not in strain GK2051. Strain GK2051 gave different restriction fragment length pattern when compared with its parent strain. Taken together, the data show restriction fragment length differences between the genomic DNA of the two strains, including the loss of some DNA sequences from the mutant strain, which may be involved in pathogenicity. Finally, B. bassiana GK2016 contains a beta-tubulin gene with at least partial homology to that of N. crassa.     

Genomic typing of Pseudomonas aeruginosa isolates by comparison of Riboprinting and PFGE: correlation of experimental results with those predicted from the complete genome sequence of isolate PAO1

The whole genomic typing of 21 isolates of Pseudomonas aeruginosa from 15 intensive care unit (ICU) patients was performed by pulsed-field gel electrophoresis (PFGE using SpeI) and Riboprinting (using EcoRI and PvuII), and then the results were compared with predictions made from the whole genome sequence of P. aeruginosa PAO1. The analysis of electronic images from PFGE and Riboprinting by GelComparII demonstrated similar discrimination between PFGE and Riboprinting with PvuII enzyme; however, Riboprinting by EcoRI had reduced banding patterns and was shown to be of lower discrimination than PvuII. When analyzing isolates from patients, both PFGE and Riboprinting using PvuII enzyme gave equivalent results, with the exception of two isolates that were closely related by PvuII Riboprinting and unrelated by PFGE. These discrepancies in typing results can be explained and adjusted for by comparisons with the rrn properties and the SpeI restriction fragments predicted from the whole genome of P. aeruginosa PAO1. Properties of the rrn operon that need to be taken into account include: (i) restriction enzyme sites that produce one or two fragments for each rrn operon; (ii) genomic variability in ISR sequence length; (iii) different enzymes need to be used to determine differences in rrn operon copy number from Riboprints; and (iv) choice of a restriction enzyme that produces riboprinter bands derived from rrn operon regions that are highly variable within the genome and between isolates. This knowledge has ramifications for PFGE and Riboprinter design and analysis so that for each new species to be typed comparisons can be made using the whole genome sequence.     

Length variation in the non-transcribed spacer of Calliphora erythrocephala ribosomal DNA is due to a 350 base-pair repeat

The non-transcribed spacer regions in the ribosomal DNA cistrons of Calliphora erythrocephala vary in length. This length variation is shown to be due to variable numbers of small repeated units found in certain areas of the NTS² regions, as has been found in several other systems. In contrast to the other organisms, however, the length variation in C. erythrocephala leads to the formation of only two major size classes with a length difference of about 1 kb. We have investigated the nature of this length difference by means of electron microscope heteroduplex and restriction enzyme digestion analyses. We demonstrate that the length variation in C. erythrocephala is due to a variation in the number of 350 bp repeating units defined by the presence of sites for three restriction enzymes, HhaI, Sau3A and AluI. Furthermore, the existence of an XbaI site within one 350 bp unit and the lack of a HhaI restriction site within another 350 bp unit reveals that at least some sequence differences exist among the repeating units. These restriction site differences can be viewed as genetic markers and lead to the hypothesis that the longer NTS class originated from the shorter NTS class by an unequal crossover event that served to duplicate the region containing repetitive units. An interesting feature of our observation is that not all NTS length classes are equally represented, suggesting special rules for the unequal recombination events often hypothesized as the basis for such variation.     

Restriction enzyme cutting site distribution regularity for DNA looping technology

The restriction enzyme cutting site distribution regularity and looping conditions were studied systematically. We obtained the restriction enzyme cutting site distributions of 13 commonly used restriction enzymes in 5 model organism genomes through two novel self-compiled software programs. All of the average distances between two adjacent restriction sites fell sharply with increasing statistic intervals, and most fragments were 0–499 bp. A shorter DNA fragment resulted in a lower looping rate, which was also directly proportional to the DNA concentration. When the length was more than 500 bp, the concentration did not affect the looping rate. Therefore, the best known fragment length was longer than 500 bp, and did not contain the restriction enzyme cutting sites which would be used for digestion. In order to make the looping efficiencies reach nearly 100%, 4–5 single cohesive end systems were recommended to digest the genome separately.     

Nonrandom DNA sequencing of exonuclease III-deleted complementary DNA

The nonrandom DNA sequence analysis procedure of Poncz et al. [Proc. Natl. Acad. Sci. USA 79, 4298-4302 (1982)] was extensively modified to permit the determination of complementary DNA (cDNA) sequences containing G-C homopolymer regions. The recombinant cDNA plasmid was cleaved at a unique restriction enzyme site close to the cDNA and treated with Exonuclease III under controlled conditions to generate a set of overlapping fragments having deletions 50-1500 bases in length at the free 3' termini. After removal of single-stranded DNA regions by Bal31 and DNA polymerase I large fragment, the unique restriction enzyme site was recreated by blunt end ligation of synthetic oligonucleotides to the deleted DNA fragments and restriction enzyme digestion. The cDNA fragment was excised from the cloning vector using a second different restriction enzyme having a unique site that flanks the cDNA fragment and subsequently force-cloned into either M13 mp10 or mp11. This method should also be particularly useful for the sequencing of other types of DNA molecules with lengths 1500 bp or smaller.     

Species-specific polymorphisms in transcribed ribosomal DNA of fivePythium species

Twenty-five isolates representing fivePythium species collected from diverse hosts and geographic origins were evaluated using polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP) analysis. DNA regions coding for the small-subunit ribosomal RNA (SrDNA) and the internal transcribed spacer (ITS) were amplified and analyzed by restriction enzyme digestion. The amplified SrDNA was about 1800 bp long and uniform in length among the five species. However, restriction digestion revealed three polymorphic groups. They areP. arrhenomanes andP. graminicola,P. irregulare andP. spinosum, andP. ultimum. The amplified-ITS region showed three different lengths which corresponded to the three polymorphic groups of SrDNA. Each length variant of the ITS showed distinct banding patterns after restriction enzyme digestion. In addition,P. irregulare andP. spinosum each showed distinct banding patterns after digestion with enzymesHinfI andMboI. Physical maps of the restriction sites in the SrDNA and the ITS were determined. Length variation occurred primarily in the spacer between the SrDNA and 5.8 S rDNA; although, it also was detected in the ITS-2 region. Little intraspecific variation was observed in the SrDNA and ITS, and species could be reliably distinguished by RFLP analysis of the amplified rDNA regions. Data presented do not support the maintenance ofP. arrhenomanes andP. graminicola as distinct species. Results indicate that PCR-RFLP can be used as a simple and speedy taxonomical tool for ecological studies ofPythium species.     

HDAPD: a web tool for searching the disease-associated protein structures

Background

PCR-restriction fragment length polymorphism (RFLP) assay is a cost-effective method for SNP genotyping and mutation detection, but the manual mining for restriction enzyme sites is challenging and cumbersome. Three years after we constructed SNP-RFLPing, a freely accessible database and analysis tool for restriction enzyme mining of SNPs, significant improvements over the 2006 version have been made and incorporated into the latest version, SNP-RFLPing 2.

Results

The primary aim of SNP-RFLPing 2 is to provide comprehensive PCR-RFLP information with multiple functionality about SNPs, such as SNP retrieval to multiple species, different polymorphism types (bi-allelic, tri-allelic, tetra-allelic or indels), gene-centric searching, HapMap tagSNPs, gene ontology-based searching, miRNAs, and SNP500Cancer. The RFLP restriction enzymes and the corresponding PCR primers for the natural and mutagenic types of each SNP are simultaneously analyzed. All the RFLP restriction enzyme prices are also provided to aid selection. Furthermore, the previously encountered updating problems for most SNP related databases are resolved by an on-line retrieval system.

Conclusions

The user interfaces for functional SNP analyses have been substantially improved and integrated. SNP-RFLPing 2 offers a new and user-friendly interface for RFLP genotyping that can be used in association studies and is freely available at http://bio.kuas.edu.tw/snp-rflping2.