'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.     

Plotting genetic maps on a microcomputer

Maps of genetic linkage and restriction enzyme cleavage sites can be quickly prepared on an IBM PC microcomputer with the commercially available program Lotus 1-2-3. Data can be entered on the keyboard or imported from other programs. The maps can be displayed on the screen or with a printer or plotter. These procedures should be useful in the research laboratory, in preparing figures for publication and in teaching.     

Plasmid map: a microcomputer program for display and storage of plasmid data

We describe a plasmid map program which runs on an IBM PC microcomputer and facilitates the drawing of circular plasmid maps. The user enters information from the keyboard in the form of restriction enzyme sites, genes and their locations, and other plasmid markers such as promoters, origins, or other sites. This information can then be stored in a file for future reference. The plasmid map can be displayed on the screen, printed on a dot-matrix printer, or plotted on a Hewlett Packard HP7475A plotter.     

A computer program for comparative analysis of nucleic acid sequences.

The programs offer the possibility of comparing pairs of homologous sequences in order to find out percentage of homology, number of identical and deviating nucleotides, of transitions and transversions and, derived from these, KNUC-values according to Kimura (1) and the corresponding standard error sigmaK. The sequences can be printed in pairs underneath each other, homologies are indicated by asterisks between the identical nucleotides. Out of a set of homologous sequences stored on a disk any number of sequences can be compared in pairs in this way, and a matrix containing either the percentage of homology values, the number of deviating nucleotides or the KNUC-values together with the corresponding standard errors can be sent to screen, printer or disk. A program will be available soon which creates a dendrogram representing the similarity between the sequences by use of an average linkage clustering method deduced from this matrix. The programs are written for Apple II computers using UCSD-PASCAL and for Sirius I/Victor 9000 computers using TURBO-PASCAL.     

CLONING: a microcomputer program for cloning simulations

A comprehensive computational tool is presented that performs cloning simulations using IBM PC/XT/AT or compatible microcomputers. The CLONING program contains a specific data base for restriction sites, gene markers, fragment sources and reference comments. It draws complete linear or circular maps either on the screen or employing conventional dot-matrix printers. The design of new recombinant molecules is a totally interactive process.     

A flexible new computer program for handling DNA sequence data.

A compact new computer program for handling nucleic acid sequence data is presented. It consists of a number of different subsets, which may be used according to a given code system. The program is designed for the determination of restriction enzyme and other recognition sites in correlation with translation patterns, and allows tabulation of codon frequencies and protein molecular weights within specified gene boundaries. The program is especially designed for detection of overlapping genes. The language, is FORTRAN and thus the program may be used on small computers; it may also be used without any prior computer experience. Copies are available on request.     

Autoscreening of Restriction Endonucleases for PCR-Restriction Fragment Length Polymorphism Identification of Fungal Species, with Pleurotus spp. as an Example

A molecular method based on PCR-restriction fragment length polymorphism (RFLP) analysis of internal transcribed spacer (ITS) ribosomal DNA sequences was designed to rapidly identify fungal species, with members of the genus Pleurotus as an example. Based on the results of phylogenetic analysis of ITS sequences from Pleurotus, a PCR-RFLP endonuclease autoscreening (PRE Auto) program was developed to screen restriction endonucleases for discriminating multiple sequences from different species. The PRE Auto program analyzes the endonuclease recognition sites and calculates the sizes of the fragments in the sequences that are imported into the program in groups according to species recognition. Every restriction endonuclease is scored through the calculation of the average coefficient for the sequence groups and the average coefficient for the sequences within a group, and then virtual electrophoresis maps for the selected restriction enzymes, based on the results of the scoring system, are displayed for the rapid determination of the candidate endonucleases. A total of 85 haplotypes representing 151 ITS sequences were used for the analysis, and 2,992 restriction endonucleases were screened to find the candidates for the identification of species. This method was verified by an experiment with 28 samples representing 12 species of Pleurotus. The results of the digestion by the restriction enzymes showed the same patterns of DNA fragments anticipated by the PRE Auto program, apart from those for four misidentified samples. ITS sequences from 14 samples (of which nine sequences were obtained in this study), including four originally misidentified samples, confirmed the species identities revealed by the PCR-RFLP analysis. The method developed here can be used for the identification of species of other living microorganisms.     

Comparison of the transposon-like structures encoding clindamycin resistance in Bacteroides R-plasmids

The R-plasmids pBF4, pBFTM10, and pBI136 encode transmissible clindamycin resistance (Ccr) in Bacteroides spp. These plasmids are distinct replicons but the regions implicated in Ccr share some homology and appear to have a transposon-like structure. To better understand the mechanism of dissemination and to locate the Ccr determinant(s), the genetic and structural properties of the Ccr regions of each plasmid were compared and contrasted. For this work a single EcoRI restriction fragment containing the Ccr region from each plasmid was cloned into pBR322 in Escherichia coli. Results of restriction mapping and heteroduplex experiments showed that the pBF4 EcoRI-D and pBFTM10 EcoRI-B fragments shared more than 90% base sequence homology but that the EcoRI-C fragment of pBI136 had diverged significantly. The pBI136 fragment also did not confer tetracycline resistance in E. coli as shown for the pBF4 EcoRI-D fragment (D.G. Guiney, P. Hasegawa, and C. E. Davis, 1984, Plasmid 11, 248-252). Heteroduplex experiments showed that the pBI136 EcoRI-C and pBF4 EcoRI-D fragments shared a 1.2-kb region of homology attributed to a directly repeated sequence which bounds the Ccr region. Southern hybridization studies indicated that an additional 0.85 kb of the pBI136 EcoRI-C fragment was homologous to the EcoRI-D fragment of pBF4. This region was characterized by its sequential restriction endonuclease sites for HindIII, AvaII, and DdeI, and it is proposed that the Ccr gene(s) resides in this area.     

Predict7, a program for protein structure prediction

We describe a program for protein sequence analysis which runs in IBM PC computers. Protein sequences are loaded from files in Mount-Conrad and Lipman-Pearson format. Seven features are analyzed: hydrophilicity, hydropathy, surface probability, side chain flexibility, antigenicity, secondary structure and N-glycosylation sites. Numeric results can be shown, printed or stored in files exportable to other programs. Graphics of up to four predictions can be displayed on the screen, printed out or plotted, with several definable options. This program has been designed to be fast, user-friendly and to be shared with the scientific community.     

alpha-Tubulin genes of Drosophila

Four Drosophila alpha-tubulin genes have been isolated on recombinant DNA molecules. The identity of two of these genes (T alpha 1 and T alpha 2) was established by isolating complementary mRNAs and then examining the in vitro translation products of the mRNAs. The one- and two-dimensional gel patterns and the peptide maps of the in vitro products were indistinguishable from those of embryonic alpha-tubulin. In turn, the embryonic tubulin was identified by determining its amino-terminal sequence. We identified two other cloned alpha-tubulin genes (T alpha 3 and T alpha 4) by their complementarity to T alpha 1 and T alpha 2. Maps of restriction endonuclease sites indicate that the four genes are different. DNA hybridization studies demonstrated, however, that three of them have extensive sequence homology with each other and slight homology with the fourth, T alpha 4. Hybridization to genomic DNA fragments indicated that the four clones genes account for all of the different alpha-tubulin genes of Drosophila melanogaster. Three of them are present only once in the haploid genome; the other, T alpha 1, is present in either one or two copies. Each of the four genes hybridizes in situ to a different site on the third chromosome.     

Patterns of digestion of human chromosomes by restriction endonucleases demonstrated by in situ nick translation

Summary A restriction enzyme-nick translation procedure has been developed for localizing sites of restriction endonuclease action on chromosomes. This method involves digestion of fixed chromosome preparations with a restriction enzyme, nick translation with DNA polymerase I in the presence of biotinylated-dUTP, detection of the incorporated biotin label with streptavidinalkaline phosphatase, and finally staining for alkaline phosphatase. Results obtained on human chromosomes using a wide variety of restriction enzymes are described, and compared with results of Giemsa and Feulgen staining after restriction enzyme digestion. Results of nick translation are not in general the opposite of those obtained with Giemsa staining, as might have been expected. Although the nick translation procedure is believed to give a more accurate picture of the distribution of restriction enzyme recognition sites on chromosomes than Giemsa staining, it is clear that the results of the nick translation experiments are affected by accessibility to the enzymes of the chromosomal DNA, as well as by the extractability of the DNA.     

Cloning and sequence analysis of the human brain beta-adrenergic receptor. Evolutionary relationship to rodent and avian beta-receptors and porcine muscarinic receptors

Two cDNA clones, lambda-CLFV-108 and lambda-CLFV-119, encoding for the beta-adrenergic receptor, have been isolated from a human brain stem cDNA library. One human genomic clone, LCV-517 (20 kb), was characterized by restriction mapping and partial sequencing. The human brain beta-receptor consists of 413 amino acids with a calculated Mr of 46480. The gene contains three potential glucocorticoid receptor-binding sites. The beta-receptor expressed in human brain was homology with rodent (88%) and avian (52%) beta-receptors and with porcine muscarinic cholinergic receptors (31%), supporting our proposal [(1984) Proc. Natl. Acad. Sci. USA 81, 272 276] that adrenergic and muscarinic cholinergic receptors are structurally related. This represents the first cloning of a neurotransmitter receptor gene from human brain.     

Computer-aided gene design.

A computer program, which runs on MS-DOS personal computers, is described that assists in the design of synthetic genes coding for proteins. The goal of the program is the design of a gene which (i) contains as many unique restriction sites as possible and (ii) uses a specific codon usage. The gene designed according to the criteria above is (i) suitable for 'modular mutagenesis' experiments and (ii) optimized for expression. The program 'reverse-translates' protein sequences into degenerated DNA sequences, generates a map of potential restriction sites and locates sequence positions where unique restriction sites can be accommodated. The nucleic acid sequence is then 'refined' according to a specific codon usage to remove any degeneration. Unique restriction sites, if potentially present, can be 'forced' into the degenerated nucleic acid sequence by using 'priority codes' assigned to different restriction sequences.     

Automatic construction of restriction site maps.

A computer program is described which constructs maps of restriction endonuclease cleavage sites in DNA molecules, given only the fragment lengths. The program utilizes fragment length data from single and double restriction enzyme digests to generate maps for linear or circular molecules. The search for a map can be limited to the unknown (insert) region of a recombinant phage or plasmid. Typical restriction maps with four or five enzymes which cut at three to five unknown sites can be calculated in a few minutes.     

MiCA: A Web-Based Tool for the Analysis of Microbial Communities Based on Terminal-Restriction Fragment Length Polymorphisms of 16S and 18S rRNA Genes

A web-based resource, Microbial Community Analysis (MiCA), has been developed to facilitate studies on microbial community ecology that use analyses of terminal-restriction fragment length polymorphisms (T-RFLP) of 16S and 18S rRNA genes. MiCA provides an intuitive web interface to access two specialized programs and a specially formatted database of 16S ribosomal RNA sequences. The first program performs virtual polymerase chain reaction (PCR) amplification of rRNA genes and restriction of the amplicons using primer sequences and restriction enzymes chosen by the user. This program, in silico PCR and Restriction (ISPaR), uses a binary encoding of DNA sequences to rapidly scan large numbers of sequences in databases searching for primer annealing and restriction sites while permitting the user to specify the number of mismatches in primer sequences. ISPaR supports multiple digests with up to three enzymes. The number of base pairs between the 5′ and 3′ primers and the proximal restriction sites can be reported, printed, or exported in various formats. The second program, APLAUS, infers a plausible community structure(s) based on T-RFLP data supplied by a user. APLAUS estimates the relative abundances of populations and reports a listing of phylotypes that are consistent with the empirical data. MiCA is accessible at .     

The phi 80 and P22 attachment sites. Primary structure and interaction with Escherichia coli integration host factor

Although the lambdoid bacteriophage phi 80 and P22 possess site-specific recombination systems analogous to bacteriophage lambda, they have different attachment (att) site specificities. We have identified and determined the nucleotide sequences of the att sites of phi 80 and P22 and have examined the interaction of these sites with purified Escherichia coli integration host factor (IHF). The sizes of the homologous core regions of the att sites vary greatly: P22 has a 46-base pair core, while phi 80 and lambda have 17- and 15-base pair cores, respectively. The core sequences of the three phage show no significant homology, although dispersed regions of homology in arm sequences indicate that the three phage att sites are related. All three att sites have a high A + T composition, and restriction fragments carrying these sites migrate anomalously upon polyacrylamide gel electrophoresis. IHF binds to a site to the left of the common core in the phi 80 and P22 phage att sites (attP) and to a site to the right of the core in P22 attP and attB (the bacterial att site). In the lambda system, IHF interacts with three regions on attP (designated H1, H2, and H') and none on attB (Craig N., and Nash, H.A. (1984) Cell 39, 707-716). Alignment of the IHF sites of all three phage results in a consensus sequence for IHF binding, Pyr-AANNNNTTGATAT. Among the three phage, the number of IHF sites differs; however, the location and orientation of the binding sites in relation to the respective core regions are well conserved. An IHF site analogous to lambda H2 is present in both phi 80 and P22 attP, while a site analogous to lambda H' is present in P22 attP. This conservation suggests that IHF plays a very similar role in the site-specific recombination pathways of all three phage, and that the flanking arm sequences are necessary for phi 80 and P22 attP function, as is the case for lambda attP function. These structural similarities presumably reflect a conservation of the mechanism of site-specific recombination for the three phage.     

A novel bacterial gene-finding system with improved accuracy in locating start codons.

Although a number of bacterial gene-finding programs have been developed, there is still room for improvement especially in the area of correctly detecting translation start sites. We developed a novel bacterial gene-finding program named GeneHacker Plus. Like many others, it is based on a hidden Markov model (HMM) with duration. However, it is a 'local' model in the sense that the model starts from the translation control region and ends at the stop codon of a coding region. Multiple coding regions are identified as partial paths, like local alignments in the Smith-Waterman algorithm, regardless of how they overlap. Moreover, our semiautomatic procedure for constructing the model of the translation control region allows the inclusion of an additional conserved element as well as the ribosome-binding site. We confirmed that GeneHacker Plus is one of the most accurate programs in terms of both finding potential coding regions and precisely locating translation start sites. GeneHacker Plus is also equipped with an option where the results from database homology searches are directly embedded in the HMM. Although this option does not raise the overall predictability, labeled similarity information can be of practical use. GeneHacker Plus can be accessed freely at http://elmo.ims.u-tokyo.ac.jp/GH/.     

Complete sequences of the ribosome recognition sites in vesicular stomatitis virus mRNAs: recognition by the 40S and 80S complexes.

Nucleotide sequences of the ribosome-protected translation initiation sites from the vesicular stomatitis virus (VSV) M and L protein mRNAs have been determined, completing the sequences of the sites from all the VSV mRNAs. A low level of protection at two internal AUG-containing sites in the N mRNA is also described. Small homologies are evident among some of the sites, but there are no obvious features common to all the sites other than a single AUG codon. In contrast, a large homology between the VSV M mRNA site and the alfalfa mosaic virus coat mRNA site (Koper-Zwarthoff et al., 1977) is noted. This homology suggests the existence of a common ancestral gene for these two apparently unrelated viruses. For each VSV mRNA species, the smallest sites protected in either the 40S or 80S initiation complexes are identical. These sites always contained the initiation codon, but only contained the capped 5' end in those mRNAs having the 5' end near the initiation site. If 40S ribosomes bind to the capped 5' end, either they do not protect it from nuclease digestion or the protection is only transitory in some VSV mRNAs. Consideration of the structures of the ribosome binding sites suggests that the differential effects of hypertonic shock on translation (Nuss and Koch, 1976) may be related to the distance between the 5' end of the mRNA and the initiation codon.     

A new type of illegitimate recombination is dependent on restriction and homologous interaction.

Illegitimate (nonhomologous) recombination requires little or no sequence homology between recombining DNAs and has been regarded as being a process distinct from homologous recombination, which requires a long stretch of homology between recombining DNAs. Under special conditions in Escherichia coli, we have found a new type of illegitimate recombination that requires an interaction between homologous DNA sequences. It was detected when a plasmid that carried 2-kb-long inverted repeats was subjected to type II restriction in vitro and type I (EcoKI) restriction in vivo within a delta rac recBC recG ruvC strain. Removal of one of the repeats or its replacement with heterologous DNA resulted in a reduction in the level of recombination. The recombining sites themselves shared, at most, a few base pairs of homology. Many of the recombination events joined a site in one of the repeats with a site in another repeat. In two of the products, one of the recombining sites was at the end of one of the repeats. Removal of one of the EcoKI sites resulted in decreased recombination. We discuss the possibility that some structure made by homologous interaction between the long repeats is used by the EcoKI restriction enzyme to promote illegitimate recombination. The possible roles and consequences of this type of homologous interaction are discussed.     

In silico analysis of complete bacterial genomes: PCR, AFLP-PCR and endonuclease restriction

We have developed a website, www.in-silico.com, which runs a software program that performs three basic tasks in completely sequenced bacterial genomes by in silico analysis: PCR amplification, amplified fragment length polymorphism (AFLP-PCR) and endonuclease restriction. For PCR, after selection of the genome and introduction of primers, fragment size, DNA sequence and corresponding open reading frame (ORF) identity of the resulting PCR product is computed. Plasmids of sequenced species may be included in the analysis. Theoretical AFLP-PCR analyzes similar parameters, and includes a suggestion tool providing a list of commercial restriction enzyme pairs yielding up to 50 amplicons in the selected genome. Endonuclease restriction analysis of complete genomes and plasmids calculates the number of restriction sites for endonucleases in a given genome. If the number of fragments is 50 or fewer, pulsed field gel electrophoresis image and restriction maps are illustrated. Other tools that have been included in this site are ORF search by name and DNA to protein translation as well as restriction digestion of user-defined DNA sequences. AVAILABILITY: This is a new molecular biology resource freely available over the Internet at http://www.in-silico.com