Advances in the analysis of DNA sequence variations using oligonucleotide microchip technology

The analysis of DNA variation (polymorphisms and mutations) on a genome-wide scale is becoming both increasingly important and technically challenging. An integration of a growing number of molecular biological methods of DNA-sequence analysis with the high-throughput feature of oligonucleotide microarray-based technologies is one of the most promising current directions of research and development.     

Advances in sedimentation velocity analysis.

On February 20, 1996, a workshop titled "Advances in Sedimentation Velocity Analysis" was held at the Biophysical Society meeting in Baltimore, Maryland, in honor of Professor David Yphantis's 65th birthday. Although he is known more for his work with sedimentation equilibrium, David's work on instrumentation and data analysis is the foundation for many of the recent advances in both equilibrium and velocity sedimentation. Over the years he has trained numerous graduate students, most of whom have gone on to emphasize the use of analytical ultracentrifugation to answer biochemical questions involving macromolecular assembly. His laboratory was one of very few that continued to use and develop analytical ultracentrifugation during its nadir in the 1970s and early 1980s. The rebirth and resurgence of analytical ultracentrifugation owe a great deal to his persistence and enthusiasm. These efforts have borne fruit. In the last five years, through his work at the National Analytical Ultracentrifugation Facility, he has helped train nearly 100 individuals in the delicate art of nonlinear least-squares analysis of equilibrium sedimentation data. Furthermore, the number of researchers using the ultracentrifuge and the number of papers published has skyrocketed in the last few years. This workshop, then, was a way to thank David for his years of devotion to analytical ultracentrifugation.     

Automated carboxy-terminal sequence analysis of peptides.

Proteins and peptides can be sequenced from the carboxy-terminus with isothiocyanate reagents to produce amino acid thiohydantoin derivatives. Previous studies in our laboratory have focused on solution phase conditions for formation of the peptidylthiohydantoins with trimethylsilylisothiocyanate (TMS-ITC) and for hydrolysis of these peptidylthiohydantoins into an amino acid thiohydantoin derivative and a new shortened peptide capable of continued degradation (Bailey, J. M. & Shively, J. E., 1990, Biochemistry 29, 3145-3156). The current study is a continuation of this work and describes the construction of an instrument for automated C-terminal sequencing, the application of the thiocyanate chemistry to peptides covalently coupled to a novel polyethylene solid support (Shenoy, N. R., Bailey, J. M., & Shively, J. E., 1992, Protein Sci. I, 58-67), the use of sodium trimethylsilanolate as a novel reagent for the specific cleavage of the derivatized C-terminal amino acid, and the development of methodology to sequence through the difficult amino acid, aspartate. Automated programs are described for the C-terminal sequencing of peptides covalently attached to carboxylic acid-modified polyethylene. The chemistry involves activation with acetic anhydride, derivatization with TMS-ITC, and cleavage of the derivatized C-terminal amino acid with sodium trimethylsilanolate. The thiohydantoin amino acid is identified by on-line high performance liquid chromatography using a Phenomenex Ultracarb 5 ODS(30) column and a triethylamine/phosphoric acid buffer system containing pentanesulfonic acid. The generality of our automated C-terminal sequencing methodology was examined by sequencing model peptides containing all 20 of the common amino acids. All of the amino acids were found to sequence in high yield (90% or greater) except for asparagine and aspartate, which could be only partially removed, and proline, which was found not be capable of derivatization. In spite of these current limitations, the methodology should be a valuable new tool for the C-terminal sequence analysis of peptides.     

New tools for protein sequence analysis.

Analysis of protein sequence is an important tool in studies of both native and recombinant proteins. Novel techniques and instrumentation which facilitate determination of protein primary structure have recently been developed.     

Automated genome sequence analysis and annotation.

MOTIVATION: Large-scale genome projects generate a rapidly increasing number of sequences, most of them biochemically uncharacterized. Research in bioinformatics contributes to the development of methods for the computational characterization of these sequences. However, the installation and application of these methods require experience and are time consuming. RESULTS: We present here an automatic system for preliminary functional annotation of protein sequences that has been applied to the analysis of sets of sequences from complete genomes, both to refine overall performance and to make new discoveries comparable to those made by human experts. The GeneQuiz system includes a Web-based browser that allows examination of the evidence leading to an automatic annotation and offers additional information, views of the results, and links to biological databases that complement the automatic analysis. System structure and operating principles concerning the use of multiple sequence databases, underlying sequence analysis tools, lexical analyses of database annotations and decision criteria for functional assignments are detailed. The system makes automatic quality assessments of results based on prior experience with the underlying sequence analysis tools; overall error rates in functional assignment are estimated at 2.5-5% for cases annotated with highest reliability ('clear' cases). Sources of over-interpretation of results are discussed with proposals for improvement. A conservative definition for reporting 'new findings' that takes account of database maturity is presented along with examples of possible kinds of discoveries (new function, family and superfamily) made by the system. System performance in relation to sequence database coverage, database dynamics and database search methods is analysed, demonstrating the inherent advantages of an integrated automatic approach using multiple databases and search methods applied in an objective and repeatable manner. AVAILABILITY: The GeneQuiz system is publicly available for analysis of protein sequences through a Web server at http://www.sander.ebi.ac. uk/gqsrv/submit     

Microcomputer programs for DNA sequence analysis.

Computer programs are described which allow (a) analysis of DNA sequences to be performed on a laboratory microcomputer or (b) transfer of DNA sequences between a laboratory microcomputer and another computer system, such as a DNA library. The sequence analysis programs are interactive, do not require prior experience with computers and in many other respects resemble programs which have been written for larger computer systems (1-7). The user enters sequence data into a text file, accesses this file with the programs, and is then able to (a) search for restriction enzyme sites or other specified sequences, (b) translate in one or more reading frames in one or both directions in order to find open reading frames, or (c) determine codon usage in the sequence in one or more given reading frames. The results are given in table format and a restriction map is generated. The modem program permits collection of large amounts of data from a sequence library into a permanent file on the microcomputer disc system, or transfer of laboratory data in the reverse direction to a remote computer system.     

Minimal-risk scoring matrices for sequence analysis.

We introduce a minimal-risk method for estimating the frequencies of amino acids at conserved positions in a protein family. Our method, called minimal-risk estimation, finds the optimal weighting between a set of observed amino acid counts and a set of pseudofrequencies, which represent prior information about the frequencies. We compute the optimal weighting by minimizing the expected distance between the estimated frequencies and the true population frequencies, measured by either a squared-error or a relative-entropy metric. Our method accounts for the source of the pseudofrequencies, which arise either from the background distribution of amino acids or from applying a substitution matrix to the observed data. Our frequency estimates therefore depend on the size and composition of the observed data as well as the source of the pseudofrequencies. We convert our frequency estimates into minimal-risk scoring matrices for sequence analysis. A large-scale cross-validation study, involving 48 variants of seven methods, shows that the best performing method is minimal-risk estimation using the squared-error metric. Our method is implemented in the package EMATRIX, which is available on the Internet at http://motif.stanford.edu/ematrix.     

RNA sequence analysis using covariance models.

We describe a general approach to several RNA sequence analysis problems using probabilistic models that flexibly describe the secondary structure and primary sequence consensus of an RNA sequence family. We call these models 'covariance models'. A covariance model of tRNA sequences is an extremely sensitive and discriminative tool for searching for additional tRNAs and tRNA-related sequences in sequence databases. A model can be built automatically from an existing sequence alignment. We also describe an algorithm for learning a model and hence a consensus secondary structure from initially unaligned example sequences and no prior structural information. Models trained on unaligned tRNA examples correctly predict tRNA secondary structure and produce high-quality multiple alignments. The approach may be applied to any family of small RNA sequences.     

Lipoamide dehydrogenase of Staphylococcus aureus: nucleotide sequence and sequence analysis.

A complex of four proteins was previously isolated from Staphylococcus aureus. The complex had a strong interaction with membrane bound ribosomes, which suggested that it may be involved in protein secretion. However, the complex was identified as pyruvate dehydrogenase (PDH), which disproved the direct role of the complex in protein secretion. Here we report the nucleotide sequence of the last gene of the S. aureus pyruvate dehydrogenase operon, pdhD, which encodes lipoamide dehydrogenase (LPD). The pdhD gene encodes a protein of 468 amino acids, with a molecular mass of 49.5 kDa. The protein is closely related to other lipoamide dehydrogenases from bacteria and eukaryotes. The possible role of membrane bound lipoamide dehydrogenase is briefly discussed.     

Cloning and sequence analysis of myostatin promoter in sheep.

To better understand the structure and function of the myostatin's gene promoter region in sheep, we cloned and sequenced a 1.517 kb fragment containing the 5'-regulatory region of the sheep myostatin gene (GenBank accession number is AY918121). The promoter sequence consists of three TATA boxes, one CAAT box, and eight putative E-boxes. Some putative muscle growth response elements for Octamer-binding factor 1(Octamer), Activator protein 1(AP1), Growth factor independence 1 zinc finger protein (Gfi-1B), Myocyte enhancer factor 2 (MEF2), Muscle-specific Mt binding site (MTBF), Glucocorticoid response elements (GRE) and Progesterone receptor binding site (PRE) were detected. Some of the motifs are conserved as compared to with that in the goat, bovine and porcine myostatin promoters. However, some differences were also found.     

ESPript: analysis of multiple sequence alignments in PostScript.

MOTIVATION: The program ESPript (Easy Sequencing in PostScript) allows the rapid visualization, via PostScript output, of sequences aligned with popular programs such as CLUSTAL-W or GCG PILEUP. It can read secondary structure files (such as that created by the program DSSP) to produce a synthesis of both sequence and structural information. RESULTS: ESPript can be run via a command file or a friendly html-based user interface. The program calculates an homology score by columns of residues and can sort this calculation by groups of sequences. It offers a palette of markers to highlight important regions in the alignment. ESPript can also paste information on residue conservation into coordinate files, for subsequent visualization with a graphics program. AVAILABILITY: ESPript can be accessed on its Web site at http://www.ipbs.fr/ESPript. Sources and helpfiles can be downloaded via anonymous ftp from ftp.ipbs.fr. A tar file is held in the directory pub/ESPript.     

A subcloning strategy for DNA sequence analysis.

We describe here a new strategy of fragment preparation for sequencing procedures using endlabelled DNA fragments as substrates (2,3) which is directly applicable to DNA fragments cloned into the Pst I site of pBR322, or in modified form, to inserts into the BamH I or Sal I site of the same plasmid. Ordered sets of subclones of predetermined overlap are are generated. These can be sequenced directly without further strand- or fragment separation steps.     

Rat sperm protamine. Isolation and sequence analysis

Rat protamine was extracted from S-(pyridylethylated) epididymal sperm cell nuclei with dilute hydrochloric acid. The final purification was achieved by reversed-phase high-performance liquid chromatography. The primary structure was determined by N-terminal sequencing of the total S-(pyridylethylated) protein, and of endoproteinase Lys-C- and thermolysin-derived fragments. Rat protamine consists of 50 amino-acid residues. It is a typical type 1 protamine and differs in two and ten positions from the corresponding mouse and rabbit protamine, respectively. Only 26 positions are invariant in all type 1 mammalian protamines.     

DNA sequence analysis by MALDI mass spectrometry.

Conventional DNA sequencing is based on gel electrophoretic separation of the sequencing products. Gel casting and electrophoresis are the time limiting steps, and the gel separation is occasionally imperfect due to aberrant mobility of certain fragments, leading to erroneous sequence determination. Furthermore, illegitimately terminated products frequently cannot be distinguished from correctly terminated ones, a phenomenon that also obscures data interpretation. In the present work the use of MALDI mass spectrometry for sequencing of DNA amplified from clinical samples is implemented. The unambiguous and fast identification of deletions and substitutions in DNA amplified from heterozygous carriers realistically suggest MALDI mass spectrometry as a future alternative to conventional sequencing procedures for high throughput screening for mutations. Unique features of the method are demonstrated by sequencing a DNA fragment that could not be sequenced conventionally because of gel electrophoretic band compression and the presence of multiple non-specific termination products. Taking advantage of the accurate mass information provided by MALDI mass spectrometry, the sequence was deduced, and the nature of the non-specific termination could be determined. The method described here increases the fidelity in DNA sequencing, is fast, compatible with standard DNA sequencing procedures, and amenable to automation.