A multiprocessor parallel approach to bit-parallel approximate string matching

The purpose of this project is to present with empirical results that a parallel design with the use of multiple processors can be successfully applied along with bit-parallel approximate string matching algorithms to solve practical bioinformatics problems. It will demonstrate that nearly optimal speedup can be achieved with a cluster of between two and eight workstations using MPI (Message Passing Interface), directly decreasing the total latency required to perform a string matching problem

Drosophila Spastin Regulates Synaptic Microtubule Networks and Is Required for Normal Motor Function

Nina Tang Sherwood is with California Institute of Technology, Qi Sun is with California Institute of Technology, Mingshan Xue is with UT Austin, Bing Zhang is with UT Austin, Kai Zinn is with California Institute of Technology.The most common form of human autosomal dominant hereditary spastic paraplegia (AD-HSP) is caused by mutations in the SPG4 (spastin) gene, which encodes an AAA ATPase closely related in sequence to the microtubule-severing protein Katanin. Patients with AD-HSP exhibit degeneration of the distal regions of the longest axons in the spinal cord. Loss-of-function mutations in the Drosophila spastin gene produce larval neuromuscular junction (NMJ) phenotypes. NMJ synaptic boutons in spastin mutants are more numerous and more clustered than in wild-type, and transmitter release is impaired. spastin-null adult flies have severe movement defects. They do not fly or jump, they climb poorly, and they have short lifespans. spastin hypomorphs have weaker behavioral phenotypes. Overexpression of Spastin erases the muscle microtubule network. This gain-of-function phenotype is consistent with the hypothesis that Spastin has microtubule-severing activity, and implies that spastin loss-of-function mutants should have an increased number of microtubules. Surprisingly, however, we observed the opposite phenotype: in spastin-null mutants, there are fewer microtubule bundles within the NMJ, especially in its distal boutons. The Drosophila NMJ is a glutamatergic synapse that resembles excitatory synapses in the mammalian spinal cord, so the reduction of organized presynaptic microtubules that we observe in spastin mutants may be relevant to an understanding of human Spastin's role in maintenance of axon terminals in the spinal cord.Biological Sciences, School o

Lipid ion channels and the role of proteins

Synthetic lipid membranes in the absence of proteins can display quantized conduction events for ions that are virtually indistinguishable from those of protein channel. By indistinguishable we mean that one cannot decide based on the current trace alone whether conductance events originate from a membrane, which does or does not contain channel proteins. Additional evidence is required to distinguish between the two cases, and it is not always certain that such evidence can be provided. The phenomenological similarities are striking and span a wide range of phenomena: The typical conductances are of equal order and both lifetime distributions and current histograms are similar. One finds conduction bursts, flickering, and multistep-conductance. Lipid channels can be gated by voltage, and can be blocked by drugs. They respond to changes in lateral membrane tension and temperature. Thus, they behave like voltage-gated, temperature-gated and mechano-sensitive protein channels, or like receptors. Lipid channels are remarkably under-appreciated. However, the similarity between lipid and protein channels poses an eminent problem for the interpretation of protein channel data. For instance, the Hodgkin-Huxley theory for nerve pulse conduction requires a selective mechanism for the conduction of sodium and potassium ions. To this end, the lipid membrane must act both as a capacitor and as an insulator. Non-selective ion conductance by mechanisms other than the gated protein-channels challenges the proposed mechanism for pulse propagation. ... Some important questions arise: Are lipid and protein channels similar due a common mechanism, or are these similarities fortuitous? Is it possible that both phenomena are different aspects of the same phenomenon? Are lipid and protein channels different at all? ... (abbreviated)Comment: 10 pages, 10 figures - accepted by 'Accounts of Chemical Research

The draft genome of the transgenic tropical fruit tree papaya ( Carica papaya Linnaeus)

Papaya, a fruit crop cultivated in tropical and subtropical regions, is known for its nutritional benefits and medicinal applications. Here we report a 3 draft genome sequence of \u27SunUp\u27 papaya, the first commercial virus-resistant transgenic fruit tree1 to be sequenced. The papaya genome is three times the size of the Arabidopsis genome, but contains fewer genes, including significantly fewer disease-resistance gene analogues. Comparison of the five sequenced genomes suggests a minimal angiosperm gene set of 13,311. A lack of recent genome duplication, atypical of other angiosperm genomes sequenced so far may account for the smaller papaya gene number in most functional groups. Nonetheless, striking amplifications in gene number within particular functional groups suggest roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical daylengths. Transgenesis at three locations is closely associated with chloroplast insertions into the nuclear genome, and with topoisomerase I recognition sites. Papaya offers numerous advantages as a system for fruit-tree functional genomics, and this draft genome sequence provides the foundation for revealing the basis of Carica\u27s distinguishing morpho-physiological, medicinal and nutritional properties

A computational investigation of kinetoplastid trans-splicing

Trans-splicing is an unusual process in which two separate RNA strands are spliced together to yield a mature mRNA. We present a novel computational approach which has an overall accuracy of 82% and can predict 92% of known trans-splicing sites. We have applied our method to chromosomes 1 and 3 of Leishmania major, with high-confidence predictions for 85% and 88% of annotated genes respectively. We suggest some extensions of our method to other systems

Frustration in Biomolecules

Biomolecules are the prime information processing elements of living matter. Most of these inanimate systems are polymers that compute their structures and dynamics using as input seemingly random character strings of their sequence, following which they coalesce and perform integrated cellular functions. In large computational systems with a finite interaction-codes, the appearance of conflicting goals is inevitable. Simple conflicting forces can lead to quite complex structures and behaviors, leading to the concept of "frustration" in condensed matter. We present here some basic ideas about frustration in biomolecules and how the frustration concept leads to a better appreciation of many aspects of the architecture of biomolecules, and how structure connects to function. These ideas are simultaneously both seductively simple and perilously subtle to grasp completely. The energy landscape theory of protein folding provides a framework for quantifying frustration in large systems and has been implemented at many levels of description. We first review the notion of frustration from the areas of abstract logic and its uses in simple condensed matter systems. We discuss then how the frustration concept applies specifically to heteropolymers, testing folding landscape theory in computer simulations of protein models and in experimentally accessible systems. Studying the aspects of frustration averaged over many proteins provides ways to infer energy functions useful for reliable structure prediction. We discuss how frustration affects folding, how a large part of the biological functions of proteins are related to subtle local frustration effects and how frustration influences the appearance of metastable states, the nature of binding processes, catalysis and allosteric transitions. We hope to illustrate how Frustration is a fundamental concept in relating function to structural biology.Comment: 97 pages, 30 figure

The analysis of pendolino (peo) mutants reveals differences in the fusigenic potential among Drosophila telomeres

Drosophila telomeres are sequence-independent structures that are maintained by transposition to chromosome ends of three specialized retroelements (HeT-A, TART and TAHRE; collectively designated as HTT) rather than telomerase activity. Fly telomeres are protected by the terminin complex (HOAP-HipHop-Moi-Ver) that localizes and functions exclusively at telomeres and by non-terminin proteins that do not serve telomere-specific functions. Although all Drosophila telomeres terminate with HTT arrays and are capped by terminin, they differ in the type of subtelomeric chromatin; the Y, XR, and 4L HTT are juxtaposed to constitutive heterochromatin, while the XL, 2L, 2R, 3L and 3R HTT are linked to the TAS repetitive sequences; the 4R HTT is associated with a chromatin that has features common to both euchromatin and heterochromatin. Here we show that mutations in pendolino (peo) cause telomeric fusions (TFs). The analysis of several peo mutant combinations showed that these TFs preferentially involve the Y, XR and 4th chromosome telomeres, a TF pattern never observed in the other 10 telomere-capping mutants so far characterized. peo encodes a non-terminin protein homologous to the E2 variant ubiquitin-conjugating enzymes. The Peo protein directly interacts with the terminin components, but peo mutations do not affect telomeric localization of HOAP, Moi, Ver and HP1a, suggesting that the peodependent telomere fusion phenotype is not due to loss of terminin from chromosome ends. peo mutants are also defective in DNA replication and PCNA recruitment. However, our results suggest that general defects in DNA replication are unable to induce TFs in Drosophila cells. We thus hypothesize that DNA replication in Peodepleted cells results in specific fusigenic lesions concentrated in heterochromatinassociated telomeres. Alternatively it is possible that Peo plays a dual function being independently required for DNA replication and telomere capping

METHODS FOR HIGH-THROUGHPUT COMPARATIVE GENOMICS AND DISTRIBUTED SEQUENCE ANALYSIS

High-throughput sequencing has accelerated applications of genomics throughout the world. The increased production and decentralization of sequencing has also created bottlenecks in computational analysis. In this dissertation, I provide novel computational methods to improve analysis throughput in three areas: whole genome multiple alignment, pan-genome annotation, and bioinformatics workflows. To aid in the study of populations, tools are needed that can quickly compare multiple genome sequences, millions of nucleotides in length. I present a new multiple alignment tool for whole genomes, named Mugsy, that implements a novel method for identifying syntenic regions. Mugsy is computationally efficient, does not require a reference genome, and is robust in identifying a rich complement of genetic variation including duplications, rearrangements, and large-scale gain and loss of sequence in mixtures of draft and completed genome data. Mugsy is evaluated on the alignment of several dozen bacterial chromosomes on a single computer and was the fastest program evaluated for the alignment of assembled human chromosome sequences from four individuals. A distributed version of the algorithm is also described and provides increased processing throughput using multiple CPUs. Numerous individual genomes are sequenced to study diversity, evolution and classify pan-genomes. Pan-genome annotations contain inconsistencies and errors that hinder comparative analysis, even within a single species. I introduce a new tool, Mugsy-Annotator, that identifies orthologs and anomalous gene structure across a pan-genome using whole genome multiple alignments. Identified anomalies include inconsistently located translation initiation sites and disrupted genes due to draft genome sequencing or pseudogenes. An evaluation of pan-genomes indicates that such anomalies are common and alternative annotations suggested by the tool can improve annotation consistency and quality. Finally, I describe the Cloud Virtual Resource, CloVR, a desktop application for automated sequence analysis that improves usability and accessibility of bioinformatics software and cloud computing resources. CloVR is installed on a personal computer as a virtual machine and requires minimal installation, addressing challenges in deploying bioinformatics workflows. CloVR also seamlessly accesses remote cloud computing resources for improved processing throughput. In a case study, I demonstrate the portability and scalability of CloVR and evaluate the costs and resources for microbial sequence analysis