A Study on the Phylogenetics of Gene Transfer: from Pathways to Kingdoms

Total Page:16

File Type:pdf, Size:1020Kb

A Study on the Phylogenetics of Gene Transfer: from Pathways to Kingdoms A STUDY ON THE PHYLOGENETICS OF GENE TRANSFER: FROM PATHWAYS TO KINGDOMS A THESIS SUBMITTED TO THE GRADUATE DIVISION OF UNIVERSITY OF HAWAI'I IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE IN MOLECULAR BIOSCIENCES AND BIOENGINEERING May 2008 By Aren Ewing Thesis Committee: Gemot Presting, Chairperson Anne Alvarez Daniel Rubinoff We certify that we have read this dissertation and that, in our opinion, it is satisfactory in scope and quality as a dissertation for the degree of Masters of Science in Molecular Biosciences and Bioengineering. Thesis Committee 2 Acknowledgments With much gratitude I would like to thank my advisor Dr. Gemot Presting. He provided and outstanding workplace and lab, with the freedom to explore our interests and the facility to follow them. During my time in the Presting lab, I also met some outstanding people and scientists. I would like to thank Thomas Wolfgruber, Anupma Sharma, Kevin Schneider, Jeffery Lai, and Normon Wang. They have provided great support and friends over the years. I would also like to thank my committee who has stuck in there with me for the completion of this degree. It was a struggle to take what was once a PhD project and convert to a MS and while making it as excellent as possible. I would also like to thank my girlfriend Hideko Kasahara for her support and encouragement during the writing of this thesis. 3 Table of Contents Acknowledgments ........................................................................................................................... 3 Table of Contents ............................................................................................................................ 4 List of Tables .................................................................................................................................. 5 List of Figures ................................................................................................................................. 6 Abstract ........................................................................................................................................... 8 Chapter I: Introduction ................................................................................................................... 9 Objectives ................................................................................................................................. 17 Chapter 2 ....................................................................................................................................... 20 Phylogenetic profiling of the Arabidopsis thaliana proteome .................................................. 20 Introduction .................................................................................................. ,........................ 20 Methods ................................................................................................................................. 24 Results ................................................................................................................................... 28 Discussion ............................................................................................................................. 48 Conclusion ............................................................................................................................ 63 Chapter 3 ....................................................................................................................................... 64 Elucidation ofXanthomonas species relationships using comparative genomics .................... 64 Introduction ........................................................................................................................... 64 Methods ................................................................................................................................. 73 Results ................................................................................................................................... 77 Discussion ............................................................................................................................. 89 Conclusion ............................................................................................................................ 96 Thesis Conclusion ......................................................................................................................... 98 1: Pathways evolve as a group as demonstrated with phylogenetic profiling ...................... 98 2: There is a dominant phylogenetic tree which can identifY consistent relationships between closely related species in the genus Xanthomonas . ................................................ 99 3: Lateral gene transfer occurred from a plant to bacteria and possibly confers a selective advantage for the infection of host plants ........................................................................... 100 Appendix ..................................................................................................................................... 101 Evaluating functionality of a putative horizontally transferred gene in Xanthomonas axonopodis pv. citri ........... ........ ............... ........ ....... ........ ...... ...................... ............... ............ 10 1 Introduction ..... ........ ....... ......... ...................... ....... ...................... ........ .............. ........ ........... 10 1 Summary ........................ ....................... ....................... ....... ....... ......................................... 125 References. ............................................. ....... .............................. ......................... ........ ....... ........ 126 4 List of Tables Table 1.1: Distribution of matches with E values from Ie-03 to le-18 ....................................... 32 Table 2.1: Gene counts summed by trees and generation method ................................................ 78 Table 2.2: Agreement on topology by model ............................................................................... 81 Table 2.3: Annotation group and method agreement for each tree topology ............................... 83 Table 2.4: List of 44 Xanthomonas genes shared by 14 organisms in iteration 3 ........................ 88 Table 3.1: Presence of insert sequenced and corresponding match to known LTG! sequence. 110 Table 3.2: List of 16 Xanthomonads ........................................................................................... 112 Table 3.3: LTG 1 presence vs host range .................................................................................... 122 Table 3.4: Infection Experiments ................................................................................................ 123 5 List of Figures Figure 1.1: Distribution of Arabidopsis gene matches in 170 non-plant genomes ..................... 30 Figure 1.2: Greedy outlier removal of pathway members .......................................................... 35 Figure 1.3: Phylogenetic profile clusters of (A) 11 isoprenoid, (B) 24 tryptophan and (C) 25 heme/chlorophyll biosynthesis genes ................................................................................... 40 Figure 1.4: Amino acid biosynthesis genes cluster together ....................................................... 44 Figure 1.5: The heme/chlorophyll biosynthesis pathway ........................................................... 46 Figure 2.1: Three monophyletic Xanthomonas trees ................................................................... 80 Figure 2.2: Consensus Tree for 229 Genes .................................................................................. 85 Figure 2.3: tRNA methyltransferase expanded 14 organism tree ............................................... 87 Figure 3.1: Xanthomonas region comparison ............................................................................ 108 Figure 3.2: Alignment of DUF239 to 28 known plant homologs .............................................. 109 Figure 3.3: PCR using Genomic DNA ...................................................................................... 112 Figure 3.4: Genomic Southern blot ofXanthomonas isolates ................................................... 113 Figure 3.5: Verification of Southern blot for integrity of Genomic DNA re-probed 16s rDNA 113 Figure 3.6: Knockout plasmid pAC3 ......................................................................................... 116 Figure 3.7: Deletion suicide plasmid pUFR080 ......................................................................... 117 Figure 3.8: Insertion plasmid pAC3 ........................................................................................... 119 Figure 3.9: Insertion suicide plasmid pUFR080 ........................................................................ 120 6 Abbreviations Alignment matrix: A conceptual matrix created when comparing two sequences, with 1 sequences on each axix. BLAST: Basic Local Alignment Search Tool BLOSSUM: Stands for (BLOcks of Amino Acid SUbstitution Matrix) primarily used to score alignments of evolutionary divergent protein sequences. Blossum62 is a matrix calculated for proteins with a 62% sequence identity. Gene Set: (Chapter 2) the group of genes that were most similar between all genomes compared based on the query genome. For the first phase of project this consists of 4 genes 3 Xanthomonas
Recommended publications
  • "Phylogenetic Profiling"
    Introductory Review Phylogenetic p rofiling Matteo Pellegrini, Todd O. Yeates and David Eisenberg Howard Hughes Medical Institute, University of California at Los Angeles, Los Angeles, CA, USA Sorel T. Fitz-Gibbon IGPP Center for Astrobiology, University of California at Los Angeles, Los Angeles, CA, USA 1. Introduction Biology has been profoundly changed by the development of techniques to se- quence DNA. The advent of rapid sequencing in conjunction with the capability to assemble sequence fragments into complete genome sequences enables researchers to read and analyze entire genomes of organisms. Parallel progress has been made in algorithms to study the evolutionary history of proteins. The techniques rely on the ability to measure the similarity of protein sequences in order to determine the likelihood that different proteins are descended from a common ancestor. It is therefore possible to reconstruct families of proteins that share a common ancestor. Combining these two capabilities, we can now not only determine which proteins are coded within an organism’s genome but we can also discover the evolutionary relationships between the proteins of multiple organisms. Phylogenetic profiling is the study of which protein types are found in which organisms. In order to perform phylogenetic profiling, one must first establish a classification of proteins into families. An example of such a classification scheme across a broad range of fully sequenced organisms is the Clusters of Orthologous Groups (Tatusov, 1997), where an attempt is made to group together proteins that perform a similar function. Next, each organism is described in terms of which protein families are coded or not coded in its genome.
    [Show full text]
  • Pellegrini M. Using Phylogenetic Profiles To
    Chapter 9 Using Phylogenetic Profiles to Predict Functional Relationships Matteo Pellegrini Abstract Phylogenetic profiling involves the comparison of phylogenetic data across gene families. It is possible to construct phylogenetic trees, or related data structures, for specific gene families using a wide variety of tools and approaches. Phylogenetic profiling involves the comparison of this data to determine which families have correlated or coupled evolution. The underlying assumption is that in certain cases these couplings may allow us to infer that the two families are functionally related: that is their function in the cell is coupled. Although this technique can be applied to noncoding genes, it is more commonly used to assess the function of protein coding genes. Examples of proteins that are functionally related include subunits of protein complexes, or enzymes that perform consecutive steps along biochemical pathways. We hypothesize the deletion of one of the families from a genome would then indirectly affect the function of the other. Dozens of different implementations of the phylogenetic profiling technique have been developed over the past decade. These range from the first simple approaches that describe phylogenetic profiles as binary vectors to the most complex ones that attempt to model to the coevolution of protein families on a phylogenetic tree. We discuss a set of these implementations and present the software and databases that are available to perform phylogenetic profiling. Key words: Phylogenetic profiles, Coevolution, Functional associations, Comparative genomics, Coevolving proteins 1. Introduction The remarkable improvements in sequencing technology that have occurred over the past few decades have made the sequenc- ing of genomes an ever more routine task.
    [Show full text]
  • Integrating Genomic Data to Predict Transcription Factor Binding
    Integrating Genomic Data to Predict Transcription Factor Binding Dustin T. Holloway 1 Mark Kon 2 Charles DeLisi 3 [email protected] [email protected] [email protected] 1 Molecular Biology Cell Biology and Biochemistry, 3 Boston University, Boston, MA 02215, U.S.A. Bioinformatics and Systems Biology, 2 Boston University, Boston, MA 02215, U.S.A. Department of Mathematics and Statistics, Boston University, Boston, MA 02215 , U.S.A. Abstract Transcription factor binding sites (TFBS) in gene promoter regions are often predicted by using position specific scoring matrices (PSSMs), which summarize sequence patterns of experimentally determined TF binding sites. Although PSSMs are more reliable than simple consensus string matching in predicting a true binding site, they generally result in high numbers of false positive hits. This study attempts to reduce the number of false positive matches and generate new predictions by integrating various types of genomic data by two methods: a Bayesian allocation procedure, and support vector machine classification. Several methods will be explored to strengthen the prediction of a true TFBS in the Saccharomyces cerevisiae genome: binding site degeneracy, binding site conservation, phylogenetic profiling, TF binding site clustering, gene expression profiles, GO functional annotation, and k-mer counts in promoter regions. Binding site degeneracy (or redundancy) refers to the number of times a particular transcription factor’s binding motif is discovered in the upstream region of a gene. Phylogenetic conservation takes into account the number of orthologous upstream regions in other genomes that contain a particular binding site. Phylogenetic profiling refers to the presence or absence of a gene across a large set of genomes.
    [Show full text]
  • Scalable Phylogenetic Profiling Using Minhash Uncovers Likely Eukaryotic Sexual Reproduction Genes
    bioRxiv preprint doi: https://doi.org/10.1101/852491; this version posted November 22, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. 1 Scalable Phylogenetic Profiling using MinHash Uncovers Likely Eukaryotic Sexual Reproduction Genes 1,2,3,* 1,2,3 4,5 1,2,3,6,7,* David Moi ​ ,​ Laurent Kilchoer ,​ Pablo S. Aguilar ​ and Christophe Dessimoz ​ ​ ​ ​ 1 2 Department​ of Computational Biology, University of Lausanne, Switzerland; Center​ for Integrative Genomics, ​ 3 4 University of Lausanne, Switzerland; SIB​ Swiss Institute of Bioinformatics, Lausanne, Switzerland; Instituto​ de ​ ​ 5 Investigaciones Biotecnologicas (IIBIO), Universidad Nacional de San Martín Buenos Aires, Argentina; Instituto​ de ​ 6 Fisiología, Biología Molecular y Neurociencias (IFIBYNE-CONICET), Department​ of Genetics, Evolution, and ​ 7 Environment, University College London, UK; Department​ of Computer Science, University College London, UK. ​ *Corresponding authors: [email protected] and [email protected] Abstract Phylogenetic profiling is a computational method to predict genes involved in the same biological process by identifying protein families which tend to be jointly lost or retained across the tree of life. Phylogenetic profiling has customarily been more widely used with prokaryotes than eukaryotes, because the method is thought to require many diverse genomes. There are now many eukaryotic genomes available, but these are considerably larger, and typical phylogenetic profiling methods require quadratic time or worse in the number of genes. We introduce a fast, scalable phylogenetic profiling approach entitled HogProf, which leverages hierarchical orthologous groups for the construction of large profiles and locality-sensitive hashing for efficient retrieval of similar profiles.
    [Show full text]
  • Review Article Protein-Protein Interaction Detection: Methods and Analysis
    Hindawi Publishing Corporation International Journal of Proteomics Volume 2014, Article ID 147648, 12 pages http://dx.doi.org/10.1155/2014/147648 Review Article Protein-Protein Interaction Detection: Methods and Analysis V. Srinivasa Rao,1 K. Srinivas,1 G. N. Sujini,2 and G. N. Sunand Kumar1 1 Department of CSE, VR Siddhartha Engineering College, Vijayawada 520007, India 2 Department of CSE, Mahatma Gandhi Institute of Technology, Hyderabad 500075, India Correspondence should be addressed to V. Srinivasa Rao; [email protected] Received 26 October 2013; Revised 5 December 2013; Accepted 20 December 2013; Published 17 February 2014 Academic Editor: Yaoqi Zhou Copyright © 2014 V. Srinivasa Rao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Protein-protein interaction plays key role in predicting the protein function of target protein and drug ability of molecules. The majority of genes and proteins realize resulting phenotype functions as a set of interactions. The in vitro and in vivo methods like affinity purification, Y2H (yeast 2 hybrid), TAP (tandem affinity purification), and so forth have their own limitations like cost, time, and so forth, and the resultant data sets are noisy and have more false positives to annotate the function of drug molecules. Thus, in silico methods which include sequence-based approaches, structure-based approaches, chromosome proximity, gene fusion, in silico 2 hybrid, phylogenetic tree, phylogenetic profile, and gene expression-based approaches were developed. Elucidation of protein interaction networks also contributes greatly to the analysis of signal transduction pathways.
    [Show full text]
  • Mapping Global and Local Co-Evolution Across 600 Species to Identify Novel Homologous Recombination Repair Genes
    Downloaded from genome.cshlp.org on October 9, 2021 - Published by Cold Spring Harbor Laboratory Press Mapping global and local co-evolution across 600 species to identify novel homologous recombination repair genes Dana Sherill-Rofe1*, Dolev Rahat1,2*, Steven Findlay3,4,#, Anna Mellul1,#, Irene Guberman1, Maya Braun1, Idit Bloch1, Alon Lalezari1, Arash Samiei 3,4, Ruslan Sadreyev5,6,7, Michal 8 3,4,9,10 2,† 1,†, Goldberg , Alexandre Orthwein , Aviad Zick and Yuval Tabach 1Department of Developmental Biology and Cancer Research, Institute for Medical Research- Israel-Canada, Hebrew University of Jerusalem, Jerusalem, Israel. 2Sharett Institute of Oncology, Hadassah Medical Center, Ein-Kerem, Jerusalem, Israel. 3Lady Davis Institute for Medical Research, Segal Cancer Centre, Jewish General Hospital, Montreal, Quebec, Canada. 4Division of Experimental Medicine, McGill University, Montreal, Quebec, Canada. 5Department of Molecular Biology, Massachusetts General Hospital, Boston, MA USA. 6Department of Genetics, Harvard Medical School, Boston, MA USA. 7Department of Pathology, Massachusetts General Hospital and Harvard Medical School, Boston, MA USA. 8Department of Genetics, Alexander Silberman Institute of Life Sciences, Hebrew University of Jerusalem, Jerusalem, Israel. 9Department of Microbiology and Immunology, McGill University, Montreal, Quebec, Canada. 10Gerald Bronfman Department of Oncology, McGill University, Montreal, Quebec, Canada. *These authors contributed equally to this work. # These authors contributed equally to this work.
    [Show full text]
  • Phylogenetic Profiling of the Arabidopsis Thaliana Proteome
    Open Access Research2004Gutiérrezet al.Volume 5, Issue 8, Article R53 Phylogenetic profiling of the Arabidopsis thaliana proteome: what comment proteins distinguish plants from other organisms? Rodrigo A Gutiérrez*†§, Pamela J Green‡, Kenneth Keegstra*† and John B Ohlrogge† Addresses: *Department of Energy Plant Research Laboratory, Michigan State University, East Lansing, MI 48824-1312, USA. †Department of Plant Biology, Michigan State University, East Lansing, MI 48824-1312, USA. ‡Delaware Biotechnology Institute, University of Delaware, 15 § Innovation Way, Newark, DE 19711, USA. Current address: Department of Biology, New York University, 100 Washington Square East, New reviews York, NY 10003, USA. Correspondence: Rodrigo A Gutiérrez. E-mail: [email protected] Published: 15 July 2004 Received: 16 March 2004 Revised: 10 May 2004 Genome Biology 2004, 5:R53 Accepted: 7 June 2004 The electronic version of this article is the complete one and can be found online at http://genomebiology.com/2004/5/8/R53 reports © 2004 Gutiérrez et al.; licensee BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL Phylogenetic<p>Theopportunitycodingof life</p> genes availability to profilingof decipher <it>Arabidopsis of theof the the complete genetic Arabidopsis thaliana genomefactors </it>onthaliana that sequence define the proteome: ofbasis plant <it>Arab of form whattheir idopsisand proteinspattern function. thaliana of distingu sequence To </it>together beginish similarityplants this task, from wi to thwe other organismsthose have organisms? of classified other across organisms the the nuc threlear proe domains protein-vides an deposited research Abstract Background: The availability of the complete genome sequence of Arabidopsis thaliana together with those of other organisms provides an opportunity to decipher the genetic factors that define plant form and function.
    [Show full text]
  • Phylogenetic Profiling of the Arabidopsis Genome
    Figure 1 Phylogenetic Profiling of the Arabidopsis Genome Aren Ewing and Gernot Presting Department of Molecular Biosciences and Bioengineering, University of Hawaii, Honolulu, HI 96822, USA Conclusion: Introduction: Results: The goal of this project is to assign genes to Member proteins of several biosynthetic pathways pathways based on co-evolution. Arabidopsis Clustered pathway proteins were found in the hierarchical tree using the greedy outlier removal clustered based on their phylogenetic profiles using several thaliana is the first fully sequenced plant, and program. Tryptophan and isoprenoid biosynthesis pathways were examined in detail. Most of different hierarchical trees and datasets. Within these many (60%) of its genes are functionally the proteins in the pathways do cluster, those that don’t either catalyze the first enzymatic step compact clusters we found new metabolic pathway annotated. This makes it a good candidate for or are misannotated / mispredicted (figure 3). New pathway members were found near known members that were not annotated in Aracyc. This study phlyogenetic studies. Utilizing existing pathway members in the hierarchical tree (figure 4). demonstrates the utility of a novel tool for assigning new information about pathways of Arabidopsis Figure 1 (left and poster border): All of the Arabidopsis genes with at least one match other members to pathways through the use of a high throughput (AraCyc), a tool was created to find novel than Arabidopsis (19324 rows) clustered in a binary tree. Clusters of unique profiles are pathway/gene list search. genes that may be associated with known Figure 1 scattered throughout the tree. Red color signifies presence of that Arabidopsis gene in target pathways.
    [Show full text]
  • Prediction of Protein-Protein Interactions Using Protein Signature Profiling
    Article Prediction of Protein-Protein Interactions Using Protein Signature Profiling Mahmood A. Mahdavi and Yen-Han Lin* Department of Chemical Engineering, University of Saskatchewan, Saskatoon, SK S7N 5A9, Canada. Protein domains are conserved and functionally independent structures that play an important role in interactions among related proteins. Domain-domain inter- actions have been recently used to predict protein-protein interactions (PPI). In general, the interaction probability of a pair of domains is scored using a trained scoring function. Satisfying a threshold, the protein pairs carrying those domains are regarded as “interacting”. In this study, the signature contents of proteins were utilized to predict PPI pairs in Saccharomyces cerevisiae, Caenorhabditis ele- gans, and Homo sapiens. Similarity between protein signature patterns was scored and PPI predictions were drawn based on the binary similarity scoring function. Results show that the true positive rate of prediction by the proposed approach is approximately 32% higher than that using the maximum likelihood estimation method when compared with a test set, resulting in 22% increase in the area un- der the receiver operating characteristic (ROC) curve. When proteins containing one or two signatures were removed, the sensitivity of the predicted PPI pairs in- creased significantly. The predicted PPI pairs are on average 11 times more likely to interact than the random selection at a confidence level of 0.95, and on aver- age 4 times better than those predicted by either phylogenetic profiling or gene expression profiling. Key words: protein-protein interaction, protein signature, ROC curve Introduction Protein-protein interaction (PPI) is the key element utilized as input to predict more accurate interac- of any biological process in a living cell.
    [Show full text]
  • Genomics of Energy & Environment User Meeting Abstracts
    DEPARTMENT OF ENERGY JOINT GENOME INSTITUTE (DOE JGI) Genomics of Energy & Environment User Meeting Abstracts March 26–28, 2013 / Walnut Creek Marriott / Walnut Creek, CA An electronic version of this document is available at: http://go.usa.gov/g3P All information current as of March 14, 2013 JGI Contact: Marsha Fenner DOE Joint Genome Institute [email protected] The Joint Genome Institute is a user facility of the Department of Energy Office of Science DOE Joint Genome Institute: www.jgi.doe.gov DOE Office of Science: science.doe.gov Eighth Annual DOE Joint Genome Institute User Meeting Sponsored By U.S. Department of Energy Office of Science March 26-28, 2013 Walnut Creek Marriott Walnut Creek, California Contents Speaker Presentations .................................................................................................................................. 1 Poster Presentations ..................................................................................................................................... 8 Attendees .................................................................................................................................................... 58 Author Index ............................................................................................................................................... 64 Speaker Presentations Abstracts alphabetical by speaker Assembly-driven metagenomics of a hypersaline microbial ecosystem Eric E. Allen ([email protected]) Marine Biology Research Division, Scripps Institution
    [Show full text]
  • Phylogene Server for Identification and Visualization of Co-Evolving Proteins Using Normalized Phylogenetic Profiles
    PhyloGene server for identification and visualization of co-evolving proteins using normalized phylogenetic profiles The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Sadreyev, Ilyas R., Fei Ji, Emiliano Cohen, Gary Ruvkun, and Yuval Tabach. 2015. “PhyloGene server for identification and visualization of co-evolving proteins using normalized phylogenetic profiles.” Nucleic Acids Research 43 (Web Server issue): W154-W159. doi:10.1093/nar/gkv452. http://dx.doi.org/10.1093/nar/gkv452. Published Version doi:10.1093/nar/gkv452 Citable link http://nrs.harvard.edu/urn-3:HUL.InstRepos:17820916 Terms of Use This article was downloaded from Harvard University’s DASH repository, and is made available under the terms and conditions applicable to Other Posted Material, as set forth at http:// nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of- use#LAA W154–W159 Nucleic Acids Research, 2015, Vol. 43, Web Server issue Published online 09 May 2015 doi: 10.1093/nar/gkv452 PhyloGene server for identification and visualization of co-evolving proteins using normalized phylogenetic profiles Ilyas R. Sadreyev1, Fei Ji1,2, Emiliano Cohen3, Gary Ruvkun1,2 and Yuval Tabach3,* 1Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA, 2Department of Genetics, Harvard Medical School, Boston, Boston, MA, USA and 3Department of Developmental Biology and Cancer Research, The Institute For Medical Research-Israel-Canada, The Hebrew University-Hadassah Medical School, Jerusalem, Israel Received February 19, 2015; Revised April 16, 2015; Accepted April 24, 2015 ABSTRACT tion together, for example members of the same pathway or protein complex, often show similar pattern of conserva- Proteins that function in the same pathways, pro- tion across phylogenetic clades.
    [Show full text]
  • Phylo Dcor: Distance Correlation As a Novel Metric for Phylogenetic Profiling Gabriella Sferra, Federica Fratini, Marta Ponzi and Elisabetta Pizzi*
    Sferra et al. BMC Bioinformatics (2017) 18:396 DOI 10.1186/s12859-017-1815-5 SOFTWARE Open Access Phylo_dCor: distance correlation as a novel metric for phylogenetic profiling Gabriella Sferra, Federica Fratini, Marta Ponzi and Elisabetta Pizzi* Abstract Background: Elaboration of powerful methods to predict functional and/or physical protein-protein interactions from genome sequence is one of the main tasks in the post-genomic era. Phylogenetic profiling allows the prediction of protein-protein interactions at a whole genome level in both Prokaryotes and Eukaryotes. For this reason it is considered one of the most promising methods. Results: Here, we propose an improvement of phylogenetic profiling that enables handling of large genomic datasets and infer global protein-protein interactions. This method uses the distance correlation as a new measure of phylogenetic profile similarity. We constructed robust reference sets and developed Phylo-dCor, a parallelized version of the algorithm for calculating the distance correlation that makes it applicable to large genomic data. Using Saccharomyces cerevisiae and Escherichia coli genome datasets, we showed that Phylo-dCor outperforms phylogenetic profiling methods previously described based on the mutual information and Pearson’s correlation as measures of profile similarity. Conclusions: In this work, we constructed and assessed robust reference sets and propose the distance correlation as a measure for comparing phylogenetic profiles. To make it applicable to large genomic data, we developed Phylo-dCor, a parallelized version of the algorithm for calculating the distance correlation. Two R scripts that can be run on a wide range of machines are available upon request. Keywords: Phylogenetic profiling, Distance correlation, Protein-protein interaction Background to gain insights on the function of uncharacterized pro- In the last two decades, several computational ap- teins [see for example [6–8].
    [Show full text]