Metagenomic Binning Through Low-Density Hashing

Total Page:16

File Type:pdf, Size:1020Kb

Metagenomic Binning Through Low-Density Hashing Metagenomic binning through low-density hashing The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Lou, Yunan, et al. "Metagenomic binning through low-density hashing." Bioinformatics 35, 2, (January 2019): 219–226 © 2018 The Author(s) As Published http://dx.doi.org/10.1093/bioinformatics/bty611 Publisher Oxford University Press (OUP) Version Final published version Citable link https://hdl.handle.net/1721.1/122806 Terms of Use Creative Commons Attribution NonCommercial License 4.0 Detailed Terms https://creativecommons.org/licenses/by-nc/4.0/ Bioinformatics, 35(2), 2019, 219–226 doi: 10.1093/bioinformatics/bty611 Advance Access Publication Date: 13 July 2018 Original Paper Genome analysis Metagenomic binning through low-density Downloaded from https://academic.oup.com/bioinformatics/article-abstract/35/2/219/5053315 by MIT Libraries user on 07 November 2019 hashing Yunan Luo1,†, Yun William Yu2,4,†, Jianyang Zeng3, Bonnie Berger4,* and Jian Peng1,* 1Department of Computer Science, University of Illinois at Urbana-Champaign, Champaign, IL 61801, USA, 2Department of Biomedical Informatics, Harvard Medical School, Boston, MA 02115, USA, 3Machine Learning and Computational Biology Group, Institute for Interdisciplinary Information Sciences, Tsinghua University, Beijing 100084, China and 4Department of Mathematics and Computer Science and AI Laboratory, Massachusetts Institute of Technology, Cambridge, MA 02139, USA *To whom correspondence should be addressed. †The authors wish it to be known that, in their opinion, the first two authors should be regarded as Joint First Authors. Associate Editor: Inanc Birol Received on April 8, 2018; revised on June 18, 2018; editorial decision on July 4, 2018; accepted on July 10, 2018 Abstract Motivation: Vastly greater quantities of microbial genome data are being generated where envir- onmental samples mix together the DNA from many different species. Here, we present Opal for metagenomic binning, the task of identifying the origin species of DNA sequencing reads. We introduce ‘low-density’ locality sensitive hashing to bioinformatics, with the addition of Gallager codes for even coverage, enabling quick and accurate metagenomic binning. Results: On public benchmarks, Opal halves the error on precision/recall (F1-score) as compared with both alignment-based and alignment-free methods for species classification. We demonstrate even more marked improvement at higher taxonomic levels, allowing for the discovery of novel lineages. Furthermore, the innovation of low-density, even-coverage hashing should itself prove an essential methodological advance as it enables the application of machine learning to other bio- informatic challenges. Availability and implementation: Full source code and datasets are available at http://opal.csail. mit.edu and https://github.com/yunwilliamyu/opal. Contact: [email protected] or [email protected] Supplementary information: Supplementary data are available at Bioinformatics online. 1 Introduction these questions are vastly more complex than standard genomic datasets due to the mixing of unknown amounts of different Metagenomics is the study of the microbiome—the many genomes genomes present. These data result in major identification challenges (bacterial, fungal and even viral) that make up a particular environ- for certain bacterial, as well as viral, species, strains and genera ment. The microbiome has been linked to human health: soil sam- (Janda and Abbott, 2007; Tu et al., 2014). ples can lead to the discovery of new antibiotics (Forsberg et al., We focus on whole-genome metagenomic DNA sequencing, 2012), and the human gut microbiome has been linked to Crohn’s since cheaper Amplicon-based sequencing methods, which concen- Disease (Erickson et al., 2012), obesity (Turnbaugh and Gordon, trate on the diversity of given marker genes (e.g. the 16S rRNA 2009) and even autism spectrum disorder (MacFabe, 2012). gene) and only analyze protein-coding regions, are limited in their Metagenomics fundamentally asks what organisms are present in a ability to provide microbial functions from the samples (1000 genomic sample with the goal of gaining insight into function. Genomes Project Consortium, 2012; Altschul et al., 1990; Berlin However, the sequencing datasets required to shine any light on et al., 2015). Unfortunately, the mixing of DNA from many VC The Author(s) 2018. Published by Oxford University Press. 219 This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected] 220 Y.Luo et al. different, sometimes related organisms in varying quantities poses et al., 2015; Keich et al., 2004; Ma et al., 2002). Spaced-seeds are an substantial computational and statistical challenges to metagenomic example of LSH, but are generally high-density because they have binning, the process of grouping reads and assigning them to an ori- historically been designed to optimize for the probability of at least gin organism. This important first step occurs before downstream one hit over the expected number of hits in a pair of homologous data analysis can be applied to elucidate the structure of microbial but non-identical regions. The original PatternHunter, which intro- populations and assign functional annotations (1000 Genomes duced spaced-seeds to bioinformatics used 11 locations in a 19-mer Project Consortium, 2012; Bromberg and Rost, 2007). Existing se- (Keich et al., 2004; Ma et al., 2002). Brinda et al. (2015) used quence alignment tools, such as BWA-MEM (Li, 2013), Bowtie 2 spaced-seeds in Kraken for metagenomic binning, but optimized for (Langmead and Salzberg, 2012) or BLAST (Altschul et al., 1990), hit count and coverage, using e.g. 16 locations in a 28-mer, or 24 Downloaded from https://academic.oup.com/bioinformatics/article-abstract/35/2/219/5053315 by MIT Libraries user on 07 November 2019 can readily be used and usually provide high-resolution alignments locations in a 36-mer. Employing a ‘high density’ of over half the and accurate results by simply finding the highest scoring matching locations in a k-mer can decrease spurious matches, necessary when genome; they have the added advantage of tolerance to small num- using a discriminative k-mer framework like that of Kraken, which bers of mismatches or gaps. However, the computational cost of can only distinguish between two species when a k-mer in the read alignment-based methods becomes prohibitive as metagenomic matches one species but not another. However, this approach still datasets continue to grow (Wood and Salzberg, 2014; Yu et al., runs into the same exponential space problem of large k-mer sizes 2015a). Similarly, contig-assembly based methods, while effective, when used in a ML framework because most of the positions in the are also computationally expensive (Alneberg et al., 2014). read still contribute to the dimensionality (Section 2). Thus, using Alternatively, the field has turned to alignment-free metagenomic lower density hashing is necessary, but it comes at the cost of not binning (also known as compositional binning) (Nawy, 2015), which capturing long discriminative k-mers. Fortunately, discriminative k- assigns sequence fragments to their taxonomic origins according to mers are not needed in SVM-based models such as the one we use, specific patterns of their constituent k-mers. State-of-the-art tools such as well as other ML-based models, so even using lengths as short as as Kraken (Wood and Salzberg, 2014)andCLARK(Ounit et al., 6–12 positions can be sufficient. Thus, low-density hashing is 2015) use exact occurrences of uniquely discriminating k-mers in enabled by SVM-based models, and furthermore, allows those mod- reads (i.e. k-mers that appear in species A but not species B) and are els to access longer-range correlations that cannot be captured using very efficient, but as a result are limited in both their sensitivity and contiguous short k-mers. To the best of our knowledge, low-density ability to detect unknown organisms. Alternately, Kaiju (Menzel hashing—using significantly fewer than half the locations in a k-mer et al.,2016) uses matches of reads to a protein database instead of ref- (e.g. 25% by default in the method we introduce)—has not previ- erence genomes for the same purpose. Other approaches rely on ously been used for bioinformatic analysis. supervised machine learning (ML) classifiers, such as Naive Bayes or Here we newly overcome these bottlenecks in handling long k- support vector machines (SVMs), trained on a set of reference genome mers by developing a novel compositional metagenomic binning algo- sequences to classify the origins of metagenomic fragments (Wang rithm, Opal, which efficiently encodes long k-mers using low- et al.,2007; Patil et al.,2011; Vervier et al.,2016; Brady and dimensional profiles generated using even-coverage, low-density hash- Salzberg, 2009) using the relative k-mer frequency vector of a read. ing. We take inspiration from low-density parity check (LDPC) error More recently, latent strain analysis performs covariance analysis of correcting codes (also known as Gallager codes)—not previously used k-mers to partition reads for low-abundance strain assembly and de- in bioinformatics—to
Recommended publications
  • Efficient Discretization Approaches for Machine Learning Techniques To
    Advances in Science, Technology and Engineering Systems Journal Vol. 5, No. 3, 547-556 (2020) ASTES Journal www.astesj.com ISSN: 2415-6698 Special Issue on Multidisciplinary Innovation in Engineering Science & Technology Efficient Discretization Approaches for Machine Learning Techniques to Improve Disease Classification on Gut Microbiome Composition Data Hai Thanh Nguyen*;1, Nhi Yen Kim Phan1, Huong Hoang Luong2 , Trung Phuoc Le2, Nghi Cong Tran3 1College of Information and Communication Technology, Can Tho University, Can Tho city, 900100, Vietnam. 2Department of Information Technology, FPT University, Can Tho city, 900000, Vietnam. 3National Central University, Taoyuan, 320317, Taiwan, R.O.C. ARTICLEINFOABSTRACT Article history: The human gut environment can contain hundreds to thousands bacterial species which Received: 10 May, 2020 are proven that they are associated with various diseases. Although Machine learning has Accepted: 18 June, 2020 been supporting and developing metagenomic researches to obtain great achievements in Online: 28 June, 2020 personalized medicine approaches to improve human health, we still face overfitting issues in Bioinformatics tasks related to metagenomic data classification where the performance in Keywords: the training phase is rather high while we get low performance in testing. In this study, we Personalized medicine present discretization methods on metagenomic data which include Microbial Compositions Bacterial composition to obtain better results in disease prediction tasks. Data types used in the experiments consist Classic machine learning of species abundance and read counts on various taxonomic ranks such as Genus, Family, Discretization Order, etc. The proposed data discretization approaches for metagenomic data in this work Species abundance are unsupervised binning approaches including binning with equal width bins, considering Read counts the frequency of values and data distribution.
    [Show full text]
  • Metagenomics Approaches for the Detection and Surveillance of Emerging and Recurrent Plant Pathogens
    microorganisms Review Metagenomics Approaches for the Detection and Surveillance of Emerging and Recurrent Plant Pathogens Edoardo Piombo 1,2 , Ahmed Abdelfattah 3,4 , Samir Droby 5, Michael Wisniewski 6,7, Davide Spadaro 1,8,* and Leonardo Schena 9 1 Department of Agricultural, Forest and Food Sciences (DISAFA), University of Torino, 10095 Grugliasco, Italy; [email protected] 2 Department of Forest Mycology and Plant Pathology, Uppsala Biocenter, Swedish University of Agricultural Sciences, P.O. Box 7026, 75007 Uppsala, Sweden 3 Institute of Environmental Biotechnology, Graz University of Technology, Petersgasse 12, 8010 Graz, Austria; [email protected] 4 Department of Ecology, Environment and Plant Sciences, University of Stockholm, Svante Arrhenius väg 20A, 11418 Stockholm, Sweden 5 Department of Postharvest Science, Agricultural Research Organization (ARO), The Volcani Center, Rishon LeZion 7505101, Israel; [email protected] 6 U.S. Department of Agriculture—Agricultural Research Service (USDA-ARS), Kearneysville, WV 25430, USA; [email protected] 7 Department of Biological Sciences, Virginia Technical University, Blacksburg, VA 24061, USA 8 AGROINNOVA—Centre of Competence for the Innovation in the Agroenvironmental Sector, University of Torino, 10095 Grugliasco, Italy 9 Department of Agriculture, Università Mediterranea, 89122 Reggio Calabria, Italy; [email protected] * Correspondence: [email protected]; Tel.: +39-0116708942 Abstract: Globalization has a dramatic effect on the trade and movement of seeds, fruits and vegeta- bles, with a corresponding increase in economic losses caused by the introduction of transboundary Citation: Piombo, E.; Abdelfattah, A.; plant pathogens. Current diagnostic techniques provide a useful and precise tool to enact surveillance Droby, S.; Wisniewski, M.; Spadaro, protocols regarding specific organisms, but this approach is strictly targeted, while metabarcoding D.; Schena, L.
    [Show full text]
  • Nitritation/Anammox System by an Iterative Hybrid Assembly Method
    Liu et al. Microbiome (2020) 8:155 https://doi.org/10.1186/s40168-020-00937-3 RESEARCH Open Access High-quality bacterial genomes of a partial- nitritation/anammox system by an iterative hybrid assembly method Lei Liu1,2,3, Yulin Wang1, You Che1, Yiqiang Chen1, Yu Xia1,2,3, Ruibang Luo4, Suk Hang Cheng5, Chunmiao Zheng2,3* and Tong Zhang1,3* Abstract Background: Genome-centric approaches are widely used to investigate microbial compositions, dynamics, ecology, and interactions within various environmental systems. Hundreds or even thousands of genomes could be retrieved in a single study contributed by the cost-effective short-read sequencing and developed assembly/ binning pipelines. However, conventional binning methods usually yield highly fragmented draft genomes that limit our ability to comprehensively understand these microbial communities. Thus, to leverage advantage of both the long and short reads to retrieve more complete genomes from environmental samples is a must-do task to move this direction forward. Results: Here, we used an iterative hybrid assembly (IHA) approach to reconstruct 49 metagenome-assembled genomes (MAGs), including 27 high-quality (HQ) and high-contiguity (HC) genomes with contig number ≤ 5, eight of which were circular finished genomes from a partial-nitritation anammox (PNA) reactor. These 49 recovered MAGs (43 MAGs encoding full-length rRNA, average N50 of 2.2 Mbp), represented the majority (92.3%) of the bacterial community. Moreover, the workflow retrieved HQ and HC MAGs even with an extremely low coverage (relative abundance < 0.1%). Among them, 34 MAGs could not be assigned to the genus level, indicating the novelty of the genomes retrieved using the IHA method proposed in this study.
    [Show full text]
  • Metabcc-LR: Metagenomics Binning by Coverage and Composition for Long Reads Anuradha Wickramarachchi 1,∗, Vijini Mallawaarachchi 1, Vaibhav Rajan 2 and Yu Lin 1,∗
    ✐ ✐ “output” — 2020/5/25 — 7:49 — page 1 — #1 ✐ ✐ Bioinformatics doi.10.1093/bioinformatics/xxxxxx Advance Access Publication Date: Day Month Year Manuscript Category Subject Section MetaBCC-LR: Metagenomics Binning by Coverage and Composition for Long Reads Anuradha Wickramarachchi 1,∗, Vijini Mallawaarachchi 1, Vaibhav Rajan 2 and Yu Lin 1,∗ 1Research School of Computer Science, Australian National University, Canberra, ACT 0200, Australia and 2Department of Information Systems and Analytics, School of Computing, National University of Singapore, Singapore 117417. ∗To whom correspondence should be addressed. Associate Editor: XXXXXXX Received on XXXXX; revised on XXXXX; accepted on XXXXX Abstract Motivation: Metagenomics studies have provided key insights into the composition and structure of microbial communities found in different environments. Among the techniques used to analyse metagenomic data, binning is considered a crucial step to characterise the different species of microorganisms present. The use of short-read data in most binning tools poses several limitations, such as insufficient species-specific signal, and the emergence of long-read sequencing technologies offers us opportunities to surmount them. However, most current metagenomic binning tools have been developed for short reads. The few tools that can process long reads either do not scale with increasing input size or require a database with reference genomes that are often unknown. In this paper, we present MetaBCC- LR, a scalable reference-free binning method which clusters long reads directly based on their k -mer coverage histograms and oligonucleotide composition. Results: We evaluate MetaBCC-LR on multiple simulated and real metagenomic long-read datasets with varying coverages and error rates. Our experiments demonstrate that MetaBCC-LR substantially outperforms state-of-the-art reference-free binning tools, achieving ∼13% improvement in F1-score and ∼30% improvement in ARI compared to the best previous tools.
    [Show full text]
  • Download From
    Rotimi et al. BMC Bioinformatics (2018) 19:309 https://doi.org/10.1186/s12859-018-2320-1 SOFTWARE Open Access Selection of marker genes for genetic barcoding of microorganisms and binning of metagenomic reads by Barcoder software tools Adeola M. Rotimi1 , Rian Pierneef1,2 and Oleg N. Reva1* Abstract Background: Metagenomic approaches have revealed the complexity of environmental microbiomes with the advancement in whole genome sequencing displaying a significant level of genetic heterogeneity on the species level. It has become apparent that patterns of superior bioactivity of bacteria applicable in biotechnology as well as the enhanced virulence of pathogens often requires distinguishing between closely related species or sub-species. Current methods for binning of metagenomic reads usually do not allow for identification below the genus level and generally stops at the family level. Results: In this work, an attempt was made to improve metagenomic binning resolution by creating genome specific barcodes based on the core and accessory genomes. This protocol was implemented in novel software tools available for use and download from http://bargene.bi.up.ac.za/. The most abundant barcode genes from the core genomes were found to encode for ribosomal proteins, certain central metabolic genes and ABC transporters. Performance of metabarcode sequences created by this package was evaluated using artificially generated and publically available metagenomic datasets. Furthermore, a program (Barcoding 2.0) was developed to align reads against barcode sequences and thereafter calculate various parameters to score the alignments and the individual barcodes. Taxonomic units were identified in metagenomic samples by comparison of the calculated barcode scores to set cut-off values.
    [Show full text]
  • Genomes from Uncultivated Prokaryotes: a Comparison of Metagenome-Assembled and Single-Amplified Genomes Johannes Alneberg1†, Christofer M
    Alneberg et al. Microbiome (2018) 6:173 https://doi.org/10.1186/s40168-018-0550-0 RESEARCH Open Access Genomes from uncultivated prokaryotes: a comparison of metagenome-assembled and single-amplified genomes Johannes Alneberg1†, Christofer M. G. Karlsson2†, Anna-Maria Divne3, Claudia Bergin3, Felix Homa3, Markus V. Lindh2,5, Luisa W. Hugerth1,4, Thijs J. G. Ettema3, Stefan Bertilsson6, Anders F. Andersson1* and Jarone Pinhassi2* Abstract Background: Prokaryotes dominate the biosphere and regulate biogeochemical processes essential to all life. Yet, our knowledge about their biology is for the most part limited to the minority that has been successfully cultured. Molecular techniques now allow for obtaining genome sequences of uncultivated prokaryotic taxa, facilitating in- depth analyses that may ultimately improve our understanding of these key organisms. Results: We compared results from two culture-independent strategies for recovering bacterial genomes: single- amplified genomes and metagenome-assembled genomes. Single-amplified genomes were obtained from samples collected at an offshore station in the Baltic Sea Proper and compared to previously obtained metagenome-assembled genomes from a time series at the same station. Among 16 single-amplified genomes analyzed, seven were found to match metagenome-assembled genomes, affiliated with a diverse set of taxa. Notably, genome pairs between the two approaches were nearly identical (average 99.51% sequence identity; range 98.77–99.84%) across overlapping regions (30–80% of each genome). Within matching pairs, the single-amplified genomes were consistently smaller and less complete, whereas the genetic functional profiles were maintained. For the metagenome-assembled genomes, only on average 3.6% of the bases were estimated to be missing from the genomes due to wrongly binned contigs.
    [Show full text]
  • Compact and Evenly Distributed K-Mer Binning for Genomic Sequences
    bioRxiv preprint doi: https://doi.org/10.1101/2020.10.12.335364; this version posted November 3, 2020. 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. Compact and evenly distributed k-mer binning for genomic sequences Johan Nystr¨om-Persson1,2, Gabriel Keeble-Gagn`ere3, Niamat Zawad2 1 JNP Solutions, Tokyo, Japan 2 Lifematics Inc, Tokyo, Japan 3 Agriculture Victoria, Victoria, Australia Abstract The processing of k-mers (subsequences of length k) is at the foundation of many sequence processing algorithms in bioinformatics, including k-mer counting for genome size estimation, genome assembly, and taxonomic classification for metagenomics. Minimizers - ordered m-mers where m < k - are often used to group k-mers into bins as a first step in such processing. However, minimizers are known to generate bins of very different sizes, which can pose challenges for distributed and parallel processing, as well as generally increase memory requirements. Furthermore, although various minimizer orderings have been proposed, their practical value for improving tool efficiency has not yet been fully explored. Here we present Discount, a distributed k-mer counting tool based on Apache Spark, which we use to investigate the behaviour of various minimizer orderings in practice when applied to metagenomics data. Using this tool, we then introduce the universal frequency ordering, a new combination of frequency counted minimizers and universal k-mer hitting sets, which yields both evenly distributed binning and small bin sizes.
    [Show full text]
  • MEGAN-LR: New Algorithms Allow Accurate Binning and Easy Interactive Exploration of Metagenomic Long Reads and Contigs Daniel H
    Huson et al. Biology Direct (2018) 13:6 https://doi.org/10.1186/s13062-018-0208-7 RESEARCH Open Access MEGAN-LR: new algorithms allow accurate binning and easy interactive exploration of metagenomic long reads and contigs Daniel H. Huson1,2* , Benjamin Albrecht1,CanerBagcı˘ 1,5, Irina Bessarab4, Anna Górska1,5, Dino Jolic3,5 and Rohan B. H. Williams4 Abstract Background: There are numerous computational tools for taxonomic or functional analysis of microbiome samples, optimized to run on hundreds of millions of short, high quality sequencing reads. Programs such as MEGAN allow the user to interactively navigate these large datasets. Long read sequencing technologies continue to improve and produce increasing numbers of longer reads (of varying lengths in the range of 10k-1M bps, say), but of low quality. There is an increasing interest in using long reads in microbiome sequencing, and there is a need to adapt short read tools to long read datasets. Methods: We describe a new LCA-based algorithm for taxonomic binning, and an interval-tree based algorithm for functional binning, that are explicitly designed for long reads and assembled contigs. We provide a new interactive tool for investigating the alignment of long reads against reference sequences. For taxonomic and functional binning, we propose to use LAST to compare long reads against the NCBI-nr protein reference database so as to obtain frame-shift aware alignments, and then to process the results using our new methods. Results: All presented methods are implemented in the open source edition of MEGAN, and we refer to this new extension as MEGAN-LR (MEGAN long read).
    [Show full text]
  • Genomic and Evolutionary Exploration of Asgard Archaea
    Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 1861 Genomic and evolutionary exploration of Asgard archaea EVA F. CACERES ACTA UNIVERSITATIS UPSALIENSIS ISSN 1651-6214 ISBN 978-91-513-0761-9 UPPSALA urn:nbn:se:uu:diva-393710 2019 Dissertation presented at Uppsala University to be publicly examined in B22, Biomedical Center (BMC), Husargatan 3, Uppsala, Thursday, 14 November 2019 at 09:15 for the degree of Doctor of Philosophy. The examination will be conducted in English. Faculty examiner: Professor Simonetta Gribaldo (Institut Pasteur, Department of Microbiology). Abstract Caceres, E. F. 2019. Genomic and evolutionary exploration of Asgard archaea. Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 1861. 88 pp. Uppsala: Acta Universitatis Upsaliensis. ISBN 978-91-513-0761-9. Current evolutionary theories postulate that eukaryotes emerged from the symbiosis of an archaeal host with, at least, one bacterial symbiont. However, our limited grasp of microbial diversity hampers insights into the features of the prokaryotic ancestors of eukaryotes. This thesis focuses on the study of a group of uncultured archaea to better understand both existing archaeal diversity and the origin of eukaryotes. In a first study, we used short-read metagenomic approaches to obtain eight genomes of Lokiarchaeum relatives. Using these data we described the Asgard superphylum, comprised of at least four different phyla: Lokiarchaeota, Odinarchaeota, Thorarchaeota and Heimdallarchaoeta. Phylogenetic analyses suggested that eukaryotes affiliate with the Asgard group, albeit the exact position of eukaryotes with respect to Asgard archaea members remained inconclusive. Comparative genomics showed that Asgard archaea genomes encoded homologs of numerous eukaryotic signature proteins (ESPs), which had never been observed in Archaea before.
    [Show full text]
  • Yet Another Deep-Learning Approach to Characterize Bacterial Genome Sequences
    bioRxiv preprint doi: https://doi.org/10.1101/2021.08.09.455756; this version posted August 10, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Genomic style: yet another deep-learning approach to characterize bacterial genome sequences Yuka Yoshimura1, Akifumi Hamada1, Yohann Augey1, Manato Akiyama1, and Yasubumi Sakakibara1* 1 Department of Biosciences and Informatics, Keio University, Yokohama, Kanagawa 223- 8522, Japan Corresponding Author * Correspondence to Yasubumi Sakakibara ([email protected]) Abstract Motivation: Biological sequence classification is the most fundamental task in bioinformatics analysis. For example, in metagenome analysis, binning is a typical type of DNA sequence classification. In order to classify sequences, it is necessary to define sequence features. The k- mer frequency, base composition, and alignment-based metrics are commonly used. In contrast, in the field of image recognition using machine learning, image classification is broadly divided into those based on shape and those based on style. A style matrix was introduced as a method of expressing the style of an image (e.g., color usage and texture). Results: We propose a novel sequence feature, called genomic style, inspired by image classification approaches, for classifying and clustering DNA sequences. As with the style of images, the DNA sequence is considered to have a genomic style unique to the bacterial species, and the style matrix concept is applied to the DNA sequence. Our main aim is to introduce the genomics style as yet another basic sequence feature for metagenome binning problem in replace of the most commonly used sequence feature k-mer frequency.
    [Show full text]
  • Environmental DNA Metabarcoding of Animal and Plant Communities
    Received: 22 July 2016 | Revised: 31 August 2017 | Accepted: 5 September 2017 DOI: 10.1111/mec.14350 INVITED REVIEWS AND SYNTHESES Environmental DNA metabarcoding: Transforming how we survey animal and plant communities Kristy Deiner1 | Holly M. Bik2 | Elvira Machler€ 3,4 | Mathew Seymour5 | Ana€ıs Lacoursiere-Roussel 6 | Florian Altermatt3,4 | Simon Creer5 | Iliana Bista5,7 | David M. Lodge1 | Natasha de Vere8,9 | Michael E. Pfrender10 | Louis Bernatchez6 1Atkinson Center for a Sustainable Future, Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA 2Department of Nematology, University of California, Riverside, CA, USA 3Eawag, Swiss Federal Institute of Aquatic Science and Technology, Department of Aquatic Ecology, Dubendorf,€ Switzerland 4Department of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich,€ Switzerland 5Molecular Ecology and Fisheries Genetics Laboratory, School of Biological Sciences, Environment Centre Wales Building, Bangor University, Bangor, Gwynedd, UK 6IBIS (Institut de Biologie Integrative et des Systemes), Universite Laval, Quebec, QC, Canada 7Wellcome Trust Sanger Institute, Hinxton, Cambridgeshire, UK 8Conservation and Research Department, National Botanic Garden of Wales, Llanarthne, Carmarthenshire, UK 9Institute of Biological, Environmental and Rural Sciences, Aberystwyth University, Aberystwyth, UK 10Department of Biological Sciences and Environmental Change Initiative, University of Notre Dame, Notre Dame, IN, USA Correspondence Kristy Deiner, Department of Ecology and Abstract Evolutionary Biology, Cornell University, The genomic revolution has fundamentally changed how we survey biodiversity on Ithaca, NY, USA. Email: [email protected] earth. High-throughput sequencing (“HTS”) platforms now enable the rapid sequenc- ing of DNA from diverse kinds of environmental samples (termed “environmental Funding information U.S. Department of Defense, Grant/Award DNA” or “eDNA”).
    [Show full text]
  • Genomic Style: Yet Another Deep-Learning Approach to Characterize Bacterial Genome Sequences
    bioRxiv preprint doi: https://doi.org/10.1101/2021.08.09.455756; this version posted August 10, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Genomic style: yet another deep-learning approach to characterize bacterial genome sequences Yuka Yoshimura1, Akifumi Hamada1, Yohann Augey1, Manato Akiyama1, and Yasubumi Sakakibara1* 1 Department of Biosciences and Informatics, Keio University, Yokohama, Kanagawa 223- 8522, Japan Corresponding Author * Correspondence to Yasubumi Sakakibara ([email protected]) Abstract Motivation: Biological sequence classification is the most fundamental task in bioinformatics analysis. For example, in metagenome analysis, binning is a typical type of DNA sequence classification. In order to classify sequences, it is necessary to define sequence features. The k- mer frequency, base composition, and alignment-based metrics are commonly used. In contrast, in the field of image recognition using machine learning, image classification is broadly divided into those based on shape and those based on style. A style matrix was introduced as a method of expressing the style of an image (e.g., color usage and texture). Results: We propose a novel sequence feature, called genomic style, inspired by image classification approaches, for classifying and clustering DNA sequences. As with the style of images, the DNA sequence is considered to have a genomic style unique to the bacterial species, and the style matrix concept is applied to the DNA sequence. Our main aim is to introduce the genomics style as yet another basic sequence feature for metagenome binning problem in replace of the most commonly used sequence feature k-mer frequency.
    [Show full text]