PubMed Health⌕ Search

Biomedical subjects

Michael Spitzer

Publications and source records attributed to Michael Spitzer.

5 recordsLinked to original sources

IsoSVM--distinguishing isoforms and paralogs on the protein level.

BACKGROUND: Recent progress in cDNA and EST sequencing is yielding a deluge of sequence data. Like database search results and proteome databases, this data gives rise to inferred protein sequences without ready access to the underlying genomic data. Analysis of this information (e.g. for EST clustering or phylogenetic reconstruction from proteome data) is hampered because it is not known if two protein sequences are isoforms (splice variants) or not (i.e. paralogs/orthologs). However, even without knowing the intron/exon structure, visual analysis of the pattern of similarity across the alignment of the two protein sequences is usually helpful since paralogs and orthologs feature substitutions with respect to each other, as opposed to isoforms, which do not. RESULTS: The IsoSVM tool introduces an automated approach to identifying isoforms on the protein level using a support vector machine (SVM) classifier. Based on three specific features used as input of the SVM classifier, it is possible to automatically identify isoforms with little effort and with an accuracy of more than 97%. We show that the SVM is superior to a radial basis function network and to a linear classifier. As an example application we use IsoSVM to estimate that a set of Xenopus laevis EST clusters consists of approximately 81% cases where sequences are each other's paralogs and 19% cases where sequences are each other's isoforms. The number of isoforms and paralogs in this allotetraploid species is of interest in the study of evolution. CONCLUSION: We developed an SVM classifier that can be used to distinguish isoforms from paralogs with high accuracy and without access to the genomic data. It can be used to analyze, for example, EST data and database search results. Our software is freely available on the Web, under the name IsoSVM.

Algorithms↗

Correspondence of function and phylogeny of ABC proteins based on an automated analysis of 20 model protein data sets.

Using our BLAST-based procedure RiPE (Retrieval-induced Phylogeny Environment), which automates the evolutionary analysis of a protein family, we assembled a set of 1138 ABC protein components [adenosine triphosphate (ATP)-binding cassette and transmembrane domain] from the protein data sets of 20 model organisms and subjected them to phylogenetic and functional analysis. For maximum speed, we based the alignment directly on a homology search with a profile of all known human ABC proteins and used neighbor-joining tree estimation. All but 11 sequences from Homo sapiens, Arabidopsis thaliana, Drosophila melanogaster, and Saccharomyces cerevisiae were placed into the correct subtree/subfamily, reproducing published classifications of the individual organisms. By following a simple "function transfer rule", our comparative phylogenetic analysis successfully predicted the known function of human ABC proteins in 19 of 22 cases. Three functional predictions did not correspond, and 10 were novel. Predictions based on BLAST alone were inferior in five cases and superior in two. Bacterial sequences were placed close to the root of most subtrees. This placement coincides with domain architecture, suggesting an early diversification of the ABC family before the kingdoms split apart. Our approach can, in principle, be used to annotate any protein family of any organism included in the study.

ATP-Binding Cassette Transporters↗

Cloning, cellular localization, genomic organization, and tissue-specific expression of the TGFbeta1-inducible SMAP-5 gene.

SMAP-5 is a member of the five-pass transmembrane protein family localizing in the Golgi apparatus and the endoplasmic reticulum. These proteins have been implicated in intracellular trafficking, in secretion and in vesicular transport. Phylogenetic analyses revealed that SMAP-5 is a member of a small Rab GTPase interacting factor protein family. The human SMAP-5 gene spans about 12.5 kb and comprises 6 exons on chromosomal locus 5q32. The proximal 5'-flanking region of the gene lacks a TATA box and is highly GC rich. Consistent with this, the SMAP-5 gene is expressed in all tissues. The highest level of expression was found in coronary smooth muscle cells, in which expression of the SMAP-5 gene was induced by transforming growth factor beta1, thus indicating that this protein may play an important role in inflammation.

Alternative Splicing↗

VisCoSe: visualization and comparison of consensus sequences.

We introduce visualization and comparison of consensus sequences (VisCoSe) as a WWW service and a stand-alone command line Perl script for visualizing and comparing consensus sequences of protein and nucleotide sequences. VisCoSe is the only interface available that simultaneously calculates consensus sequences of multiple data sets and automatically compares these consensus sequences. Furthermore, VisCoSe allows visualization of chemical properties of amino acids.

Algorithms↗

BLASTing proteomes, yielding phylogenies.

We develop a procedure called RiPE (Retrieval-induced Phylogeny Environment) that automatically performs an evolutionary analysis of a protein (sub)family, (i) by retrieving the relevant sequences via a homology search, (ii) by using the search report to construct the alignment using only homologous subsequences (taking into account their neighborhood with a low chance of homology), (iii) by realigning, and (iv) by generating phylogenetic trees based on the alignment. In a first implementation of our scheme, we start with the available proteome data of model organisms, perform a PSI-BLAST search, use MView to convert hits into a multiple alignment, and perform realignment and tree building. As a test case, we have investigated the human ABC transporters of the subfamily G, starting with the five known human ABCG transporters. Our method retrieved homologous sequences not previously analyzed, generating a tree that is more plausible and better supported than previously published trees. The RiPE 0.1 prototype is available at the RiPE website, http://ifg-izkf.uni-muenster.de/fuellen/RiPE/ripe.html.

ATP Binding Cassette Transporter, Subfamily G, Mem↗