11068219
An approximation algorithm for haplotype inference by maximum parsimony
This paper studies haplotype inference by maximum parsimony using population data. We define the optimal haplotype inference (OHI) problem as given a set of genotypes and a set of related haplotypes, find a minimum subset of haplotypes that can resolve all the genotypes. We prove that OHI is NP-hard and can be formulated as an integer quadratic programming (IQP) problem. To solve the IQP problem, we propose an iterative semi-definite programming based approximation algorithm, (called SDPHapInfer). We show that this algorithm finds a solution within a factor of O(logn) of the optimal solution, where n is the number of genotypes. This algorithm has been implemented and tested on a variety of simulated and biological data. In comparison with three other methods: HAPAR, HAPLOTYPER, and PHASE, the experimental results indicate that SDPHapInfer and HAPLOTYPER have similar error rates. In addition, the results generated by PHASE have lower error rates on some data but higher error rates on others. The error rates of HAPAR are higher than the others on biological data. In terms of efficiency, SDPHapInfer, HAPLOTYPER, and PHASE output a solution in a stable and consistent way, and they run much faster than HAPAR when the number of genotypes becomes large.
2057857
Gene-finding via tandem mass spectrometry
We propose a new gene-finding methodology that combines high performance liquid chromatograph (HPLC)-tandem mass spectrometry experiments with a fast computer algorithm to locate coding regions and introns. Proteins are first extracted from cells and digested by enzymes, and then the resulting peptides are separated and analyzed by HPLC-tandem mass spectrometry. We designed an algorithm to find DNA coding sequences, corresponding to open reading frames (ORF), in the genome such that their translated amino acid sequences are optimally correlated with these tandem mass spectra. In this algorithm, we also allow one gap, corresponding to an intron, between two DNA coding sequences, such that their concatenation becomes one coding sequence. Finally, the algorithm assembles these candidate coding sequences and introns into gene structures. Our algorithm was implemented to predict genes on 4 contigs with a total of 123 kbps using two sets of simulated digestion- HPLC-tandem mass spectrometry data of 2523 Caenorhabditis elegans Chromosome IV proteins, digested by trypsin and Asp-N respectively. Among 15 annotated genes in the forward strand, all 98 exons are hit by the predicted no-gap coding sequences, and 60 out of 83 introns are correctly predicted. We also tested gene structure prediction in a contig containing 3 genes. Combining splicing site predictions with predicted coding sequences and introns, we found all 3 gene structures.
2822926
Sorting with Fixed-length Reversals
Abstract A popular puzzle TOP-SPIN consists of a permutation of 20 numbered disks on an oval track, with a turnstile capable of reversing a string of 4 consecutive disks. The goal is to sort the disks into identity permutation using reversals. We consider the more general case of sorting n element permutations with a turnstile of size k. The problem of computing the reversal distances is of considerable importance in reconstructing the evolutionary history of the genome. The minimum reversal sequence sorting the one genome to another corresponds to the most likely evolutionary path between them. We give a complete solution for all n and k, of the number of equivalence classes of n-permutations under k-reversal, for both permutations and circular permutations. We also prove the upper and lower bounds on the reversal distances between two permutations.
31567164
Prediction of protein function using protein-protein interaction data
Assigning functions to novel proteins is one of the most important problems in the post-genomic era. We develop a novel approach that applies the theory of Markov random fields to infer a protein's functions using protein-protein interaction data and the functional annotations of its interaction protein partners. For each function of interest and a protein, we predict the probability that the protein has that function using Bayesian approaches. Unlike in other available approaches for protein annotation where a protein has or does not have a function of interest, we give a probability for having the function. This probability indicates how confident we are about the prediction. We apply our method to predict cellular functions (43 categories including a category "others") for yeast proteins defined in the Yeast Proteome Database, using the protein-protein interaction data from the Munich Information Center for Protein Sequences. We show that our approach outperforms other available methods for function prediction based on protein interaction data.
10712046
Algorithms for identifying protein cross-links via tandem mass spectrometry
Cross-linking technology combined with tandem mass spectrometry (MS-MS) is a powerful method that provides a rapid solution to the discovery of protein-protein interactions and protein structures. We studied the problem of detecting cross-linked peptides and cross-linked amino acids from tandem mass spectral data. Our method consists of two steps: the first step finds two protein subsequences whose mass sum equals a given mass measured from the mass spectrometry; and the second step finds the best cross-linked amino acids in these two peptide sequences that are optimally correlated to a given tandem mass spectrum. We designed fast and space-efficient algorithms for these two steps, and implemented and tested them on experimental data of cross-linked Hemoglobin proteins.
8315237
Local Lazy Regression: Making Use of the Neighborhood to Improve QSAR Predictions
Traditional quantitative structure-activity relationship (QSAR) models aim to capture global structure-activity trends present in a data set. In many situations, there may be groups of molecules which exhibit a specific set of features which relate to their activity or inactivity. Such a group of features can be said to represent a local structure-activity relationship. Traditional QSAR models may not recognize such local relationships. In this work, we investigate the use of local lazy regression (LLR), which obtains a prediction for a query molecule using its local neighborhood, rather than considering the whole data set. This modeling approach is especially useful for very large data sets because no a priori model need be built. We applied the technique to three biological data sets. In the first case, the root-mean-square error (RMSE) for an external prediction set was 0.94 log units versus 0.92 log units for the global model. However, LLR was able to characterize a specific group of anomalous molecules with much better accuracy (0.64 log units versus 0.70 log units for the global model). For the second data set, the LLR technique resulted in a decrease in RMSE from 0.36 log units to 0.31 log units for the external prediction set. In the third case, we obtained an RMSE of 2.01 log units versus 2.16 log units for the global model. In all cases, LLR led to a few observations being poorly predicted compared to the global model. We present an analysis of why this was observed and possible improvements to the local regression approach.
44900435
Two test data selection strategies towards testing of Boolean specifications
Test data selection for software represented as Boolean formulae has not received much attention until the last decade. Most of the approaches are expression driven because generation of test cases is based on the Boolean formulae being tested. The paper discusses test data generation using a fault based approach in the sense that generation of test cases is based on particular types of faults occurring in the Boolean expression. Two special types of faults are considered. Moreover, two strategies of test case selection that guarantee the detection of these two types of faults are proposed.
418775
Progressive image watermarking
Progressive transmission of images is very useful in many applications, especially in image transmission over the Internet. To view an image, people would want to see part of the image while the image transmission is in progress, rather than having to wait until the end of the image transmission. On the other hand, the ease of transmission and copying of images creates the need to use digital watermarking to embed the copyright information seamlessly into the media. We propose a progressive image watermarking scheme. In this scheme, the watermark is embedded in such a way that we can retrieve part of it even when the watermarked image is still being transmitted. As transmission progresses, the retrieved watermark has a decreasing bit error rate. Our proposed methods not only transmit the watermarked image progressively, but also intelligently select watermark embedding locations robust to various attacks.
1816267
From Region Encoding To Extended Dewey: On Efficient Processing of XML Twig Pattern Matching
Finding all the occurrences of a twig pattern in an XML database is a core operation for efficient evaluation of XML queries. A number of algorithms have been proposed to process a twig query based on region encoding labeling scheme. While region encoding supports efficient determination of structural relationship between two elements, we observe that the information within a single label is very limited. In this paper, we propose a new labeling scheme, called extended Dewey. This is a powerful labeling scheme, since from the label of an element alone, we can derive all the elements names along the path from the root to the element. Based on extended Dewey, we design a novel holistic twig join algorithm, called TJFast. Unlike all previous algorithms based on region encoding, to answer a twig query, TJFast only needs to access the labels of the leaf query nodes. Through this, not only do we reduce disk access, but we also support the efficient evaluation of queries with wildcards in branching nodes, which is very difficult to be answered by algorithms based on region encoding. Finally, we report our experimental results to show that our algorithms are superior to previous approaches in terms of the number of elements scanned, the size of intermediate results and query performance.
43913114
Second-generation error concealment for video transport over error-prone channels
Video transport over error-prone channels may result in loss or erroneous decoding of the video data. Error concealment is an effective mechanism to reconstruct the video data. In this paper, we review error-concealment methods and introduce a new framework, which we refer to as second-generation error concealment. All the error-concealment methods reconstruct the lost video data by making use of certain a priori knowledge about the video content. First-generation error concealment builds such a priori in a heuristic manner. The proposed second-generation error concealment builds the a priori by modeling the statistics of the video content explicitly, typically in the region of interest (ROI). Context-based models are trained with the correctly received video data and then used to replenish the lost video data. Trained models capture the statistics of the video content and thus reconstruct the lost video data better than reconstruction by heuristics. A new dynamic model ‘updating principal components’ (UPC) is proposed as a model for second-generation error concealment. UPC can be applied to pixel values to conceal loss of pixel data. In addition, UPC can be applied to motion vectors, which results in ‘updating eigenflows’ (U-Eigenflow), to conceal loss of motion vectors. With UPC applied to both pixel values and motion vectors, hybrid temporal/spatial error concealment can be achieved. The proposed second-generation error-concealment method provides superior performances to first-generation error-concealment methods. Copyright © 2002 John Wiley & Sons, Ltd.
7383373
Scale-Driven Iterative Optimization for Brain Extraction and Registration
We present a novel framework to automatically separate brain region from other non-brain regions in head images. The idea of the proposed method is to estimate larger scale patterns in an image and then correct the boundaries iteratively. The scale estimation is based on the recently proposed total variation (TV) regularized L^1 functional. An iterative optimization method is used to refine non-convex and acute angle boundaries. The final algorithm is able to extract large scale patterns with arbitrary shapes, which is particularly suitable for brain extraction. In order to reduce the computation overhead in 3D data, a multi-level technique is proposed to exponentially improve the speed of the brain extraction process. Based on accurate results of brain extraction, a non-rigid brain registration algorithm is proposed to improve accuracy and consistency of existing registration methods. Experimental results on real 3D brain MR images demonstrate that the proposed methods outperform existing solutions. In addition, results are provided to show that the proposed algorithm can also be used to segment large scale patterns in general images.