Rodger Staden

The staden sequence analysis package

I describe the current version of the sequence analysis package developed at the MRC Laboratory of Molecular Biology, which has come to be known as the “Staden Package.” The package covers most of the standard sequence analysis tasks such as restriction site searching, translation, pattern searching, comparison, gene finding, and secondary structure prediction, and provides powerful tools for DNA sequence determination. Currently the programs are only available for computers running the UNIX operating system. Detailed information about the package is available from our WWW site: http://www.mrc-lmb.cam.ac.uk/pubseq/.

Computer methods to locate signals in nucleic acid sequences.

This paper describes computer methods for locating signals in nucleic acid sequences. The signals include ribosome binding sites, promoter sequences and splice junctions. The methods are of use both to those trying to interpret the function of newly determined sequences and to those studying the molecular mechanisms involved in the recognition of these special signal sequences.

Codon preference and its use in identifying protein coding regions in long DNA sequences

This paper describes a computer method that uses codon preference to help find protein coding regions in long DNA sequences. The method can distinguish between introns and exons and can help to detect sequencing errors.

The current status and portability of our sequence handling software

I describe the current status of our sequence analysis software. The package contains a comprehensive suite of programs for managing large shotgun sequencing projects, a program containing 61 functions for analysing single sequences and a program for comparing pairs of sequences for similarity. The programs that have been described before have been improved by the addition of new functions and by being made very much easier to use. The major interactive programs have 125 pages of online help available from within them. Several new programs are described including screen editing of aligned gel readings for shotgun sequencing projects; a method to highlight errors in aligned gel readings, new methods for searching for putative signals in sequences. We use the programs on a VAX computer but the whole package has been rewritten to make it easy to transport it to other machines. I believe the programs will now run on any machine with a FORTRAN77 compiler and sufficient memory. We are currently putting the programs onto an IBM PC XT/AT and another micro running under UNIX.

Graphic methods to determine the function of nucleic acid sequences

We describe an interactive computer program (ANALYSEQ) that is used from a simple graphics terminal. The main purpose of the program is to determine the function of nucleic acid sequences but it also offers the simpler listing, searching and counting options. It contains methods to locate genes by looking for the effects that coding for a protein has on the coding sequence, to locate tRNA genes by looking for secondary structure and conserved bases, and methods to locate signals such as promoters. Techniques to identify unusual regions of sequence and to search for potential Z DNA-forming regions are also included. Most of the routines produce graphical output which gives ease of interpretation and allows superposition of several independent forms of analysis.

Measurements of the effects that coding for a protein has on a DNA sequence and their use for finding genes

Protein genes can be found either by searching the DNA sequence for signals such as ribosome binding sites or by looking for the effects that coding for a protein has on the coding sequence. This paper describes how these coding effects can be measured and used to detect protein coding regions.

Methods for calculating the probabilities of finding patterns in sequences

This paper describes the use of probability-generating functions for calculating the probabilities of finding motifs in nucleic acid and protein sequences. Equations and algorithms are given for calculating the probabilities associated with nine different ways of defining motifs. Comparisons are made with searches of random sequences. A higher level structure--the pattern--is defined as a list of motifs. A pattern also specifies the permitted ranges of spacing allowed between its constituent motifs. Equations for calculating the expected numbers of matches to patterns are given.

Software for genome mapping by fingerprinting techniques

A genome mapping package has been developed for reading and assembling data from clones analysed by restriction enzyme fragmentation and polyacrylamide gel electrophoresis. The package comprises: data entry; matching; assembly; statistical analysis; modelling. Data entry can be either manual or by a semiautomatic system based on a scanning densitometer. The primary emphasis in the analytical routines is on flexibility and interactive convenience, so that the operator has full knowledge of and control over the growing map, but a variety of automatic options are included. The package continually grows to meet the needs of the Caenorhabditis project.

Deviations from Expected Frequencies of CpG Dinucleotides in Herpesvirus DNAs May Be Diagnostic of Differences in the States of Their Latent Genomes

The DNA sequences of genomes from G + C-rich and A + T-rich lymphotropic herpesviruses [i.e. gammaherpesviruses; Epstein-Barr virus and herpesvirus saimiri (HVS)] are deficient in CpG dinucleotides and contain an excess of TpG and CpA dinucleotides relative to frequencies predicted from their mononucleotide compositions. In contrast, for sequences from genomes of G + C-rich and A + T-rich neurotropic herpesviruses (i.e. alphaherpesviruses; herpes simplex virus and varicellazoster virus) and human cytomegalovirus (HCMV; a betaherpesvirus) the mean observed frequencies of these dinucleotides are close to those expected from their mononucleotide compositions. Comparisons between DNA sequences that encode proteins conserved in all these viruses also show that sequences of these lymphotropic viruses are CpG-deficient whereas the homologous genes from the neurotropic viruses and the HCMV are not. Analyses of local variations in dinucleotide frequencies reveal some occurrences of clustered CpG dinucleotides in generally deficient genomes (e.g. upstream of the thymidylate synthase gene of HVS) and locally CpG-deficient regions within some generally non-deficient genomes (e.g. the major immediate early genes of human, simian and murine CMVs). A relative deficiency in CpG and an excess of TpG and CpA dinucleotides is a diagnostic feature of higher eukaryotic DNA sequences that have been subjected to methylation of cytosine residues in CpG doublets with the resulting increase in mutations to give TpG (and thereby its complement, CpA). The available evidence implicates the latent genome as the site of methylation of these herpesviruses. We conclude that in the neurotropic herpesviruses the normal latent precursors to infectious progeny are not methylated whereas there is local methylation of the immediate early locus in the latent genomes of CMVs, and the latent genomes of these lymphotropic herpesviruses are extensively methylated.

Methods for discovering novel motifs in nucleic acid sequences.

We describe a computer tool to aid the discovery of new motifs in nucleic acid sequences. A typical use would be to analyse a set of upstream regions from a family of related genes in order to find possible control sequences. The heart of the method is the creation of dictionaries of related subsequences. These dictionaries can then be analysed to look for the commonest or best-defined subsequences, those that occur in the highest number of different sequences, or for those in equivalent positions within the family. We show the application of the method to a set of E. coli promoter sequences.