A New Tool for De Novo Extraction of Strains from Metagenomes Christopher Quince1* , Tom O

Total Page:16

File Type:pdf, Size:1020Kb

A New Tool for De Novo Extraction of Strains from Metagenomes Christopher Quince1* , Tom O Quince et al. Genome Biology (2017) 18:181 DOI 10.1186/s13059-017-1309-9 METHOD Open Access DESMAN: a new tool for de novo extraction of strains from metagenomes Christopher Quince1* , Tom O. Delmont2, Sébastien Raguideau1, Johannes Alneberg3, Aaron E. Darling4, Gavin Collins5,6 and A. Murat Eren2,7 Abstract We introduce DESMAN for De novo Extraction of Strains from Metagenomes. Large multi-sample metagenomes are being generated but strain variation results in fragmentary co-assemblies. Current algorithms can bin contigs into metagenome-assembled genomes but are unable to resolve strain-level variation. DESMAN identifies variants in core genes and uses co-occurrence across samples to link variants into haplotypes and abundance profiles. These are then searched for against non-core genes to determine the accessory genome of each strain. We validated DESMAN on a complex 50-species 210-genome 96-sample synthetic mock data set and then applied it to the Tara Oceans microbiome. Keywords: Metagenomes, Strain, Niche Background become uncertain and fragment into multiple contigs [4]. Metagenomics, the direct sequencing of DNA extracted Metagenomes contain many conserved regions between from an environment, offers a unique opportunity to strains. These act effectively as inter-genome repeats, and study whole microbial communities in situ. The major- hence produce highly fragmented assemblies. This is par- ity of contemporary metagenomics studies use shotgun ticularly true when multiple samples are co-assembled sequencing, where DNA is fragmented prior to sequenc- together. It is possible to bin these contigs into parti- ing with short reads, of the order of hundreds of base tions that derive from the same species using sequence pairs (bps). To realise the potential of metagenomics fully, composition [5, 6] and more powerfully, the varying methods capable of resolving both the species and the coverage of individual co-assembled contigs over multi- strains present in this data are needed. Reference-based ple samples [7–10]. However, the resulting genome bins, solutions for strain identification have been developed or metagenome-assembled genomes (MAGs), represent [1, 2] but for the vast majority of microbial species, com- aggregates of multiple similar strains. These strains will prehensive strain-level databases do not exist. This situ- vary both in the precise sequence of shared genes, when ation is unlikely to change, particularly given the great that variation is below the resolution of the assembler, but diversity of microbes that are elusive to standard culti- also in gene complement, because not all genes and hence, vation techniques [3]. This motivates de novo strategies contigs will be present in all strains. capable of resolving novel variation at high resolution Modified experimental approaches can be used to sim- directly from metagenomic data. plify the challenge of metagenomics assembly by reduc- It is not usually possible simply to assemble metage- ing individual sample complexity, for example through nomic reads into individual genomes that provide strain- enrichment cultures that preferentially grow organisms level resolution. This is because in the presence of repeats adapted to particular growth conditions [11] or with (identical regions that exceed the read length), assemblies potentially less bias by selecting small subsets of cells using flow cytometry and sequencing with low-input *Correspondence: [email protected] DNA techniques [12]. The latter has been coupled with 1Warwick Medical School, University of Warwick, Gibbet Hill Road, CV4 7AL the sequencing of a standard whole-community metage- Coventry, UK nomics sample in a novel binning pipeline, MetaSort Full list of author information is available at the end of the article [13], which exploits the assembly graph to map the flow © The Author(s). 2017 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Quince et al. Genome Biology (2017) 18:181 Page 2 of 22 cytometry sample sequences onto those from the commu- which strain. From the analysis of core genes, we know nity metagenome and to extract genomes. However, for how many strains are present and their relative abun- the majority of studies that do not perform enrichment or dances across samples. The signature of relative frequen- flow cytometry, improved bioinformatics algorithms will cies across samples associated with each strain will also be required to resolve strain variation from metagenome be observed on the non-core gene variants but, crucially, data sets. not all strains will possess these genes and potentially A number of methods exist that map reads against ref- they may be in multiple copies. The relative strain fre- erence genes or genomes to resolve strain-level variation quencies have to be adjusted, therefore, to reflect these de novo [14–16]. The most straightforward approach is copy numbers. For instance, if a gene is present in just a to take the consensus single-nucleotide polymorphisms single copy in one strain, it can have no variants. In addi- (SNPs) in individual samples to be the haplotypes [16, 17]. tion, the total coverage associated with a gene will also This cannot, however, resolve mixtures and will entirely depend on which strains possess that gene being a sim- miss strains that are not dominant at least somewhere. ple sum of the individual strain coverages. Here, we do These shortcomings can be addressed by using the fre- not address the multi-copy problem, just gene presence or quency of the variants across multiple samples to resolve absence in a strain. We infer these given the observed vari- de novo strain-level variation and abundances. This is the ant base frequencies and gene coverages across samples approach taken in the Lineage algorithm of O’Brien whilst keeping the strain signatures fixed at those com- et al. [14] and ConStrains [15]. However, no method puted from the SCSGs and SCGs. This also provides a has yet been developed that works from assembled con- strategy for inferring non-core gene haplotypes on strains. tigs, avoiding the need for any reference genomes, and, Taken together, these two steps provide a procedure for hence, is applicable to microbial populations that lack cul- resolving both strain haplotypes on the core genome and tured representatives. Here, we show that it is possible their gene complements entirely de novo from short-read to combine this principle with contig-binning algorithms metagenome data. We recommend applying this strat- and resolve the strain-level variation in MAGs, both in egy to genes, but crucially genes called on the assem- terms of nucleotide variation on core genes and variation bled contigs. If contig assignments are preferred, the in gene complement. same methodology could be applied directly to the con- We denote our strategy DESMAN forDenovoExtrac- tigs themselves, or a consensus assignment of genes tion of Strains from Metagenomes. We assume that a on a contig used to determine its presence or absence co-assembly has been performed and the contigs binned in a given strain. The DESMAN pipeline is summarised into MAGs. Any binning algorithm could be used for this, in Fig. 1. but here we applied CONCOCT [9]. We also assume that The advantage of using base frequencies across samples reads have been mapped back onto these contigs as part to resolve strains, rather than existing haplotype resolu- of this process. To resolve strain variation within a MAG tion algorithms that link variants using reads [18], is that or group of MAGs deriving from a single species, we first it enables us to resolve variation that is less divergent than identify core genes that are present in all strains as a sin- the reciprocal of the read length and to link strains across gle copy. In the absence of any reference genomes, these contigs. The intuition behind frequency-based strain will simply be those genes known to be core for all bacteria inference is similar to that of contig binning. The frequen- and archaea (single-copy core genes or SCGs), e.g. the 36 cies of variants associated with a strain fluctuate across clusters of orthologous groups of proteins (COGs) iden- samples with the abundance of that strain. However, in tified in [9]. If reference genomes from the same species this case it is necessary to consider that multiple strains or related taxa are available, then these can be used to may share the same nucleotide at a given variant position. identify further genes that will satisfy the criteria of being To solve this problem, we develop a full Bayesian model, presentinallstrainsinasinglecopy,inwhichcasewe fitted by a Markov chain Monte Carlo (MCMC) Gibbs denote these as single-copy core species genes (SCSGs). sampler, to learn the strain frequencies, their haplotypes Using the read mappings, we calculate the base frequen- and also sequencing error rates. To improve convergence, cies at each position on the SCSGs or SCGs. Next, we we initialise the Gibbs sampler using non-negative matrix determine variant positions using a likelihood ratio test factorisation (NMF), or more properly non-negative ten- applied to the frequencies of each base summed across sor factorisation (NTF), a method from machine learning samples. We then use the base frequencies across samples that is equivalent to the maximum likelihood solution on these variant positions to resolve the number of strains [19]. Our approach is like the Lineage algorithm devel- present, their abundance and their unique sequence or oped by O’Brien et al.
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]
  • 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.
    [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]