DOCSLIB.ORG

A Flexible Pipeline for Complete Analysis of Bacterial

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
A Flexible Pipeline for Complete Analysis of Bacterial bioRxiv preprint doi: https://doi.org/10.1101/2020.02.28.969394; this version posted March 2, 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-NC 4.0 International license. 1 Bactopia: a flexible pipeline for complete 2 analysis of bacterial genomes 3 4 5 Robert A. Petit III1 and Timothy D. Read1 6 1Division of Infectious Diseases, Department of Medicine, Emory University School of 7 Medicine. 8 9 Corresponding Author: 10 Timothy Read1 11 12 Email address: [email protected] 13 14 Abstract 15 Sequencing of bacterial genomes using Illumina technology has become such a standard 16 procedure that often data are generated faster than can be conveniently analyzed. We created 17 a new series of pipelines called Bactopia, built using Nextflow workflow software, to provide 18 efficient comparative genomic analyses for bacterial species or genera. Bactopia consists of a 19 dataset setup step (Bactopia Datasets; BaDs) where a series of customizable datasets are 20 created for the species of interest; the Bactopia Analysis Pipeline (BaAP), which performs 21 quality control, genome assembly and several other functions based on the available datasets bioRxiv preprint doi: https://doi.org/10.1101/2020.02.28.969394; this version posted March 2, 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-NC 4.0 International license. 22 and outputs the processed data to a structured directory format; and a series of Bactopia Tools 23 (BaTs) that perform specific post-processing on some or all of the processed data. BaTs include 24 pan-genome analysis, computing average nucleotide identity between samples, extracting and 25 profiling the 16S genes and taxonomic classification using highly conserved genes. It is 26 expected that the number of BaTs will increase to fill specific applications in the future. As a 27 demonstration, we performed an analysis of 1,664 public Lactobacillus genomes, focusing on L. 28 crispatus, a species that is a common part of the human vaginal microbiome. Bactopia is an 29 open source system that can scale from projects as small as one bacterial genome to 30 thousands that allows for great flexibility in choosing comparison datasets and options for 31 downstream analysis. Bactopia code can be accessed at 32 https://www.github.com/bactopia/bactopia. 33 34 Introduction 35 Sequencing a bacterial genome, an activity that once required the infrastructure of a dedicated 36 genome center, is now a routine task that even a small laboratory can undertake. A large 37 number of open source software tools have been created to handle various parts of the process 38 of using raw read data for functions such as SNP calling and de novo assembly. As a result of 39 dedicated community efforts, it has recently become much easier to locally install these 40 bioinformatic tools through package managers (Bioconda (1), Brew) or through the use of 41 software containers (Docker, Singularity). Despite these advances, producers of bacterial 42 sequence data face a bewildering array of choices when considering how to perform analysis, 43 particularly when large numbers of genomes are involved and processing efficiency and 44 scalability become major factors. 45 bioRxiv preprint doi: https://doi.org/10.1101/2020.02.28.969394; this version posted March 2, 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-NC 4.0 International license. 46 Efficient bacterial multi-genome analysis has been hampered by three missing functionalities. 47 First, is the need to have “workflows of workflows'' that can integrate analyses and provide a 48 simplified way to start with a collection of raw genome data, remove low quality sequences and 49 perform the basic analytic steps of de novo assembly, mapping to reference sequence and 50 taxonomic assignment. Second, is the desire to incorporate user-specific knowledge of the 51 species into the input of the main genome analysis pipeline. While many microbiologists are not 52 expert bioinformaticians, they are experts in the organisms they study. Third, is the need to 53 create an output format from the main pipeline that could be used for future customized 54 downstream analysis such as pan-genome analysis and basic visualization of phylogenies. 55 56 Here we introduce Bactopia, an integrated suite of workflows for flexible analysis of Illumina 57 genome sequencing projects of bacteria from the same taxon. Bactopia is based on the 58 Nextflow workflow software (2), and is designed to be scalable, allowing projects as small as a 59 single genome to be run on a local desktop, or many thousands of genomes to be run as a 60 batch on a cloud infrastructure. Running multiple tasks on a single platform standardizes the 61 underlying data quality used for gene and variant calling between projects run in different 62 laboratories. This structure also simplifies the user experience. In Bactopia, complex multi- 63 genome analysis can be run in a small number of commands. However, there are myriad 64 options for fine-tuning datasets used for analysis and the functions of the system. The 65 underlying Nextflow structure ensures reproducibility. To illustrate the functionality of the system 66 we performed a Bactopia analysis of 1,664 public genome projects of the Lactobacillus genus, 67 an important component of the microbiome of humans and animals. 68 bioRxiv preprint doi: https://doi.org/10.1101/2020.02.28.969394; this version posted March 2, 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-NC 4.0 International license. 69 Design and implementation 70 Bactopia links together open source bioinformatics software, available from Bioconda (1), using 71 Nextflow (2). Nextflow was chosen for its flexibility: Bactopia can be run locally, on clusters, or 72 on cloud platforms with simple parameter changes. It also manages the parallel execution of 73 tasks and creates checkpoints allowing users to resume jobs. Nextflow automates installation of 74 the component software of the workflow through integration with Bioconda. For ease of 75 deployment, Bactopia can either be installed through Bioconda, a Docker container, or a 76 Singularity container. All the software programs used by Bactopia (version v1.3.0) described in 77 this manuscript are listed in Table 1 with their individual version numbers. 78 79 There are three main components of Bactopia (Figure 1): 1) Bactopia Datasets (BaDs), a 80 framework for formatting organism-specific datasets to be used by the downstream analysis 81 pipeline; 2) Bactopia Analysis Pipeline (BaAP), a customizable workflow for the analysis of 82 individual bacterial genome projects that is an extension and generalization of the previously 83 published Staphylococcus aureus-specific “Staphopia Analysis Pipeline” (StAP) (3). The inputs 84 to BaAP are FASTQ files from bacterial Illumina sequencing projects, either imported from the 85 National Centers for Biotechnology Information (NCBI) Short Read Archive (SRA) database or 86 provided locally, and any reference data in the BaDs; 3) Bactopia Tools (BaTs), a set of 87 workflows that use the output files from a BaAP project to run genomic analysis on multiple 88 genomes. For this project we used BaTs to a) summarize the results of running multiple 89 bacterial genomes through BaAP, b) extract 16S sequences and create a phylogeny, c) assign 90 taxonomic classifications with the Genome Taxonomy Database (4), d) subset Lactobacillus 91 crispatus samples by average nucleotide identity (ANI) with FastANI (5), e) run pan-genome 92 analysis for L. crispatus using Roary (6) and create a core-genome phylogeny. bioRxiv preprint doi: https://doi.org/10.1101/2020.02.28.969394; this version posted March 2, 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-NC 4.0 International license. 93 Figure 1 - Bactopia overview 94 95 96 Bactopia Datasets 97 The Bactopia pipeline can be run without downloading and formatting Bactopia Datasets 98 (BaDs). However, providing them enriches the downstream analysis. Bactopia can import 99 specific existing public datasets, as well as accessible user-provided datasets in the appropriatete 100 format. A subcommand (bactopia datasets) was created to automate downloading, building, andnd 101 (or) configuring these datasets for Bactopia. 102 103 BaDs can be grouped into those that are general and those that are user-supplied. General 104 datasets include: a Mash (7) sketch of the NCBI RefSeq (8) and PLSDB (9) databases and a 105 Sourmash (10) signature of microbial genomes (includes viral and fungal) from the NCBI 106 GenBank (11) database. Ariba (12), a software for detecting genes in raw read (FASTQ) files, 107 uses a number of default reference databases for virulence and antibiotic resistance. The bioRxiv preprint doi: https://doi.org/10.1101/2020.02.28.969394; this version posted March 2, 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-NC 4.0 International license.
Load more

Recommended publications
  • Computer Applications Making Rapid Advances in High Throughput Microbial Proteomics (HTMP) Balakrishna Anandkumar1,§, Steve W
    Send Orders for Reprints to [email protected] Combinatorial Chemistry & High Throughput Screening, 2014, 17, 173-182 173 Computer Applications Making Rapid Advances in High Throughput Microbial Proteomics (HTMP) Balakrishna Anandkumar1,§, Steve W. Haga2,§ and Hui-Fen Wu*,3,4,5,6 1Department of Biochemistry and Biotechnology, Sourashtra College, Madurai 625004, India 2Department of Computer Science and Engineering, National Sun Yat-sen University, Kaohsiung, 804, Taiwan 3Department of Chemistry, National Sun Yat Sen University, Kaohsiung, 804, Taiwan 4School of Pharmacy, College of Pharmacy, Kaohsiung Medical University, Kaohsiung, 800, Taiwan 5Center for Nanoscience and Nanotechnology, National Sun Yat-Sen University, Kaohsiung, 804, Taiwan 6Doctoral Degree Program in Marine Biotechnology, National Sun Yat-sen University, Kaohsiung, 804, Taiwan Abstract: The last few decades have seen the rise of widely-available proteomics tools. From new data acquisition devices, such as MALDI-MS and 2DE to new database searching softwares, these new products have paved the way for high throughput microbial proteomics (HTMP). These tools are enabling researchers to gain new insights into microbial metabolism, and are opening up new areas of study, such as protein-protein interactions (interactomics) discovery. Computer software is a key part of these emerging fields. This current review considers: 1) software tools for identifying the proteome, such as MASCOT or PDQuest, 2) online databases of proteomes, such as SWISS-PROT, Proteome Web, or the Proteomics Facility of the Pathogen Functional Genomics Resource Center, and 3) software tools for applying proteomic data, such as PSI-BLAST or VESPA. These tools allow for research in network biology, protein identification, functional annotation, target identification/validation, protein expression, protein structural analysis, metabolic pathway engineering and drug discovery.
    [Show full text]
  • For Immediate Release. Oct. 8Th, 2010 Microbesonline.Org Enhanced for Metagenomics and Metabolism
    For Immediate Release. Oct. 8th, 2010 MicrobesOnline.org Enhanced for Metagenomics and Metabolism & The Arkin lab (genomics.lbl.gov) computational biology and bioinformatics team, as part of the Joint BioEnergy Institute (JBEI, www.jbei.org) and the Virtual Institute for Microbial Stress and Survival (VIMSS, vimss.lbl.gov), has added significant new functionality to the MicrobesOnline.org web resource for comparative and functional analysis of microbial genomes. This work, led by Dylan Chivian with John Bates, Keith Keller, Morgan Price, Paramvir Dehal, under the guidance of Adam Arkin, offers the scientific community new capabilities for metagenomic and metabolic analyses. metaMicrobesOnline Figure 1. metaMicrobesOnline permits the user to select metagenomes and isolates they want to study and then perform easy searches for keywords or gene families (e.g. “amylase” or “PFAM00128”). Scanning electron microscopy image of microbial compost community courtesy Bernhard Knierim and Manfred Auer. The new meta.MicrobesOnline.org web resource, which currently contains over 1600 genomes from bacterial, archaeal, and microeukaryotic isolates, offers combined phylogenetic gene tree analysis of millions of genes from over 150 ecological and organismal metagenomes. These trees are built using our FastTree [1] program, which offers rapid highly accurate tree building, even for very large trees. Such combined analysis is superior to BLAST-based homology approaches in that trees offers the ability to place genes from environmental samples into an evolutionary context and permits more precise functional grouping within a gene family, yielding information about the key genes for a given environment. Additionally, comparison with isolate genomes gives researchers clues for which additional genes to look for to complete the components of systems, or may possess phylogenetic markers to aid in assigning the species for environmental sequence fragments, permitting the determination of which community members are responsible for which roles.
    [Show full text]
  • Comparing Genomes: Databases and Computational Tools for Comparative Analysis of Prokaryotic Genomes DOI: 10.3395/Reciis.V1i2.Sup.105En
    [www.reciis.cict.fiocruz.br] ISSN 1981-6286 SUPPLEMENT – BIOINFORMATICS AND HEALTH Review Articles Comparing genomes: databases and computational tools for comparative analysis of prokaryotic genomes DOI: 10.3395/reciis.v1i2.Sup.105en Marcos Antonio Basílio Catanho de Miranda Laboratório de Genômica Laboratório de Genômica Funcional e Bioinformática Funcional e Bioinformática do Instituto Oswaldo Cruz do Instituto Oswaldo Cruz da Fundação Oswaldo da Fundação Oswaldo Cruz, Cruz, Rio de Janeiro, Brazil Rio de Janeiro, Brazil [email protected] [email protected] Wim Degrave Laboratório de Genômica Funcional e Bioinformática do Instituto Oswaldo Cruz da Fundação Oswaldo Cruz, Rio de Janeiro, Brazil [email protected] Abstract Since the 1990’s, the complete genetic code of more than 600 living organisms has been deciphered, such as bacte- ria, yeasts, protozoan parasites, invertebrates and vertebrates, including Homo sapiens, and plants. More than 2,000 other genome projects representing medical, commercial, environmental and industrial interests, or comprising mo- del organisms, important for the development of the scientific research, are currently in progress. The achievement of complete genome sequences of numerous species combined with the tremendous progress in computation that occurred in the last few decades allowed the use of new holistic approaches in the study of genome structure, orga- nization and evolution, as well as in the field of gene prediction and functional classification. Numerous public or proprietary databases and computational tools have been created attempting to optimize the access to this information through the web. In this review, we present the main resources available through the web for comparative analysis of prokaryotic genomes.
    [Show full text]
  • OMA Browser—Exploring Orthologous Relations Across 352 Complete Genomes Adrian Schneider*, Christophe Dessimoz and Gaston H
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by RERO DOC Digital Library Vol. 23 no. 16 2007, pages 2180–2182 BIOINFORMATICS APPLICATIONS NOTE doi:10.1093/bioinformatics/btm295 Genome analysis OMA Browser—Exploring orthologous relations across 352 complete genomes Adrian Schneider*, Christophe Dessimoz and Gaston H. Gonnet Institute of Computational Science, ETH Zurich, Switzerland Received on February 8, 2007; revised and accepted on May 25, 2007 Advance Access publication June 1, 2007 Associate Editor: Alfonso Valencia ABSTRACT employ more basic algorithms, but distinguish themselves by a Motivation: Inference of the evolutionary relation between proteins, large number of analyzed sequences. in particular the identification of orthologs, is a central problem in The OMA project (Dessimoz et al., 2005) is a massive comparative genomics. Several large-scale efforts with various cross-comparison of complete genomes to identify methodologies and scope tackle this problem, including OMA the evolutionary relation between any pair of proteins. (the Orthologous MAtrix project). The main features of OMA are the large number of genomes Results: Based on the results of the OMA project, we introduce here from all kingdoms of life, the strict verification of orthology the OMA Browser, a web-based tool allowing the exploration of assignments and the determination of the phylogenetic orthologous relations over 352 complete genomes. Orthologs can be relationship between any two proteins. These major differences viewed as groups across species, but also at the level of sequence to other orthologs projects will be explained in the following pairs, allowing the distinction among one-to-one, one-to-many and subsections.
    [Show full text]
  • GTL PI Meeting 2006 Abstracts Molecular Interactions
    Milestone 1 Section 4 Molecular Interactions 32 Exploring the Genome and Proteome of Desulfitobacterium hafniense DCB2 for its Protein Complexes Involved in Metal Reduction and Dechlorination James M. Tiedje1* ([email protected]), John Davis2, Sang-Hoon Kim1, David Dewitt1, Christina Harzman1, Robin Goodwin1, Curtis Wilkerson1, Joan Broderick1,3, and Terence L. Marsh1 1Michigan State University, East Lansing, MI; 2Columbus State University, Columbus, GA; and 3Montana State University, Bozeman, MT Desulfitobacterium hafniense is an anaerobic, low GC, Gram-positive spore-forming bacterium that shows considerable promise as a bioremediative competent population of sediments. Of particular relevance to its remediation capabilities are metal reduction, which changes the mobility and toxicity of metals, and chlororespiration, the ability to dechlorinate organic xenobiotics. Our work focuses on these two capabilities in D. hafniense DCB2 whose genome was recently sequenced by JGI and annotated by ORNL. We are determining the complete metal reducing repertoire of D. hafniense and the genes required for respiratory and non-respiratory metal reduction through DNA microarrays produced by ORNL and proteomics. In addition, genomic analysis identified seven putative reduc- tive dehalogenases (RDases) and we are investigating these with genetic and biochemical tools. The status of these efforts is presented below. Physiology D. hafniense is capable of reducing iron as well as uranium, selenium, copper and cobalt. With the exception of selenium, the metals were reduced under conditions where the metal was the only avail- able electron acceptor. Selenium was unable to serve in this capacity but was reduced when grown fermentatively. Grown in the presence of Se, the cellular morphology viewed with light microscopy and confirmed with SEM & TEM was elongated with small polyp-like spheres present on the outer surface and in the medium.
    [Show full text]
  • Comparative Genomics and Functional Analysis of Rhamnose Catabolic Pathways and Regulons in Bacteria Irina A
    Hope College Digital Commons @ Hope College Faculty Publications 12-23-2013 Comparative Genomics And Functional Analysis Of Rhamnose Catabolic Pathways And Regulons In Bacteria Irina A. Rodionova Sanford-Burnham Medical Research Institute Xiaoqing Li Sanford-Burnham Medical Research Institute Vera Thiel Pennsylvania State University Sergey Stolyar Pacific oN rthwest National Laboratory Krista Stanton Hope College See next page for additional authors Follow this and additional works at: http://digitalcommons.hope.edu/faculty_publications Part of the Microbiology Commons Recommended Citation Rodionova IA, Li X, Thiel V, Stolyar S, Stanton K, Fredrickson JK, Bryant DA, Osterman AL, Best AA and Rodionov DA (2013) Comparative genomics and functional analysis of rhamnose catabolic pathways and regulons in bacteria. Front. Microbiol. 4:407. doi: 10.3389/fmicb.2013.00407 This Article is brought to you for free and open access by Digital Commons @ Hope College. It has been accepted for inclusion in Faculty Publications by an authorized administrator of Digital Commons @ Hope College. For more information, please contact [email protected]. Authors Irina A. Rodionova, Xiaoqing Li, Vera Thiel, Sergey Stolyar, Krista Stanton, James K. Fredrickson, Donald A. Bryant, Andrei L. Osterman, Aaron A. Best, and Dmitry A. Rodionov This article is available at Digital Commons @ Hope College: http://digitalcommons.hope.edu/faculty_publications/1116 ORIGINAL RESEARCH ARTICLE published: 23 December 2013 doi: 10.3389/fmicb.2013.00407 Comparative genomics and functional analysis of rhamnose catabolic pathways and regulons in bacteria Irina A. Rodionova 1, Xiaoqing Li 1, Vera Thiel 2, Sergey Stolyar 3†, Krista Stanton 4, James K. Fredrickson 3, Donald A. Bryant 2,5, Andrei L.
    [Show full text]
  • Environmental Conditions Shape the Nature of a Minimal Bacterial Genome
    ARTICLE https://doi.org/10.1038/s41467-019-10837-2 OPEN Environmental conditions shape the nature of a minimal bacterial genome Magdalena Antczak 1, Martin Michaelis 1 & Mark N. Wass 1 Of the 473 genes in the genome of the bacterium with the smallest genome generated to date, 149 genes have unknown function, emphasising a universal problem; less than 1% of proteins have experimentally determined annotations. Here, we combine the results from 1234567890():,; state-of-the-art in silico methods for functional annotation and assign functions to 66 of the 149 proteins. Proteins that are still not annotated lack orthologues, lack protein domains, and/ or are membrane proteins. Twenty-four likely transporter proteins are identified indi- cating the importance of nutrient uptake into and waste disposal out of the minimal bacterial cell in a nutrient-rich environment after removal of metabolic enzymes. Hence, the envir- onment shapes the nature of a minimal genome. Our findings also show that the combination of multiple different state-of-the-art in silico methods for annotating proteins is able to predict functions, even for difficult to characterise proteins and identify crucial gaps for further development. 1 Industrial Biotechnology Centre, School of Biosciences, University of Kent, Canterbury, Kent CT2 7NJ, UK. Correspondence and requests for materials should be addressed to M.M. (email: [email protected]) or to M.N.W. (email: [email protected]) NATURE COMMUNICATIONS | (2019) 10:3100 | https://doi.org/10.1038/s41467-019-10837-2 | www.nature.com/naturecommunications 1 ARTICLE NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-019-10837-2 long-term goal of synthetic biology has been the identi- Unknown class of proteins (7%) have related sequences in Afication of the minimal genome, i.e., the smallest set of eukaryotes or archaea (6%) while just over half (55%) have genes required to support a living organism.
    [Show full text]
  • VIMSS Computational Microbiology Core Research on Comparative and Functional Genomics
    Genomics:GTL Program Projects Lawrence Berkeley National Laboratory 4 VIMSS Computational Microbiology Core Research on Comparative and Functional Genomics Adam Arkin*1,2,3 ([email protected]), Eric Alm1, Inna Dubchak1, Mikhail Gelfand, Katherine Huang1, Vijaya Natarajan1, Morgan Price1, and Yue Wang2 1Lawrence Berkeley National Laboratory, Berkeley, CA; 2University of California, Berkeley, CA; and 3Howard Hughes Medical Institute, Chevy Chase, MD Background. The VIMSS Computational Core group is tasked with data management, statisti- cal analysis, and comparative and evolutionary genomics for the larger VIMSS effort. In the early years of this project, we focused on genome sequence analysis including development of an operon prediction algorithm which has been validated across a number of phylogenetically diverse species. Recently, the Computational Core group has expanded its efforts, integrating large amounts of func- tional genomic data from several species into its comparative genomic framework. Operon Prediction. To understand how bacteria work from genome sequences, before considering experimental data, we developed methods for identifying groups of functionally related genes. Many bacterial genes are organized in linear groups called operons. The problem of identifying operons had been well studied in model organisms such as E. coli, but we wished to predict operons in less studied bacteria such as D. vulgaris, where data to train the prediction method is not available. We used comparisons across dozens of genomes to identify likely conserved operons, and used these conserved operons instead of training data. The predicted operons from this approach show good agreement with known operons in model organisms and with gene expression data from diverse bacteria. Statistical Modeling of Functional Genomics Experiments.
    [Show full text]
  • Pangenome Analysis Toolkit
    PATO: Pangenome Analysis Toolkit Miguel D. Fernández-de-Bobadilla & Val F. Lanza 2021-01-29 Contents Introduction 2 Installation 2 Linux/Unix . .2 MacOS ...................................................3 Functions and Classes 3 Classes . .3 mash class...............................................4 mmseq class ..............................................4 accnet class..............................................4 nr_list class . .4 gff_list class . .4 core_genome and core_snps_genome ...............................5 accnet_enr ..............................................5 Functions . .5 Main Functions . .5 Analysis . 11 Quality Control . 13 Visualization . 14 Examples 15 Quality Control and DataSet manipulation . 15 Species and Outliers . 15 Redundancies . 16 Population Structure / Pangenome Description . 17 Core Genome analysis . 18 Accessory Genome Analysis . 21 Population structure . 25 Annotation . 30 Outbreak / Transmission / Description . 34 Mash Similarity. 36 Core-genome + snp_matrix (roary-like) . 38 Core_snp_genome + snp_matrix (snippy-like) . 41 Snps-pairwaise (more computationally expensive but more accurate) . 44 Pangenomes Analysis . 46 Benchmarking 50 Session Info . 52 1 Introduction PATO is an R package designed to analyze pangenomes (set of genomes) from the same or different species (intra/inter species). PATO allows the analysis of population structure, phylogenetics and horizontal gene transfer using in each case the core-genome (set of genes common to all genomes), the accessory genome (set of “dispensable” genes) or the whole genome. PATO uses external software in its core software: MASH, MMSeq2, Minimap2. This software can handle thousands of genomes using conventional computers without the need of using HPC facilities. We have designed PATO to work with the most common files such as GFF (General Feature Format), FNA (FASTA Nucleotide ), FFN (FASTA Feature Nucleotide) or FAA (FASTA Aminoacid). PATO is able to perform common genomic analysis of phylogenetic, annotation or population structure.
    [Show full text]
  • Myoglobin Primary Structure Reveals Multiple Convergent Transitions To
    RESEARCH ARTICLE Myoglobin primary structure reveals multiple convergent transitions to semi- aquatic life in the world’s smallest mammalian divers Kai He1,2,3,4*, Triston G Eastman1, Hannah Czolacz5, Shuhao Li1, Akio Shinohara6, Shin-ichiro Kawada7, Mark S Springer8, Michael Berenbrink5*, Kevin L Campbell1* 1Department of Biological Sciences, University of Manitoba, Winnipeg, Canada; 2Department of Biochemistry and Molecular Biology, School of Basic Medical Sciences, Southern Medical University, Guangzhou, China; 3State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China; 4Guangdong Provincial Key Laboratory of Single Cell Technology and Application, Southern Medical University, Guangzhou, China; 5Department of Evolution, Ecology and Behaviour, University of Liverpool, Liverpool, United Kingdom; 6Department of Bio-resources, Division of Biotechnology, Frontier Science Research Center, University of Miyazaki, Miyazaki, Japan; 7Department of Zoology, Division of Vertebrates, National Museum of Nature and Science, Tokyo, Japan; 8Department of Evolution, Ecology and Organismal Biology, University of California, Riverside, Riverside, United States *For correspondence: Abstract The speciose mammalian order Eulipotyphla (moles, shrews, hedgehogs, solenodons) [email protected] (KH); combines an unusual diversity of semi-aquatic, semi-fossorial, and fossorial forms that arose from [email protected]. terrestrial forbearers. However, our understanding of
    [Show full text]
  • GTL PI Meeting 2007 Abstracts
    Develop the Knowledgebase, Computational Meth- GTL Milestone 3 ods, and Capabilities to Advance Understanding and Prediction of Complex Biological Systems Section 1 Computing Infrastructure, Bioinformatics, and Data Management GTL 112 Center for Computational Biology at the University of California, Merced Michael Colvin1* ([email protected]), Arnold Kim,1 Masa Watanabe,1* and Felice Lightstone2 1School of Natural Sciences, University of California, Merced, California and 2Biosciences Directorate, Lawrence Livermore National Laboratory, Livermore, California Project Goals: The goals of the UCM-CCB are 1) to help train a new generation of biologists who bridge the gap between the computational and life sciences and to implement a new biology curriculum that can both influence and be adopted by other universities, 2) and to facilitate the development of multidisciplinary research programs in computational and mathematical biology. The Center for Computational Biology (UCM-CCB) was established at the newest campus of the University of California in fall, 2004. The UCM-CCB is sponsoring multidisciplinary scientific projects in which biological understanding is guided by mathematical and computational modeling. The center is also facilitating the development and dissemination of undergraduate and graduate course materials based on the latest research in computational biology. This project is a multi-insti- tutional collaboration including the new University of California campus at Merced, Rice University, Rensselaer Polytechnic Institute,
    [Show full text]
  • The Impact of Long-Distance Horizontal Gene Transfer on Prokaryotic Genome Size
    The impact of long-distance horizontal gene transfer on prokaryotic genome size Otto X. Cordero1 and Paulien Hogeweg Theoretical Biology and Bioinformatics, University of Utrecht, Padualaan 8 3584 CH, Utrecht, The Netherlands Edited by James M. Tiedje, Michigan State University, East Lansing, MI, and approved October 30, 2009 (received for review July 11, 2009) Horizontal gene transfer (HGT) is one of the most dominant forces of the processes underlying complexification in microbial ge- molding prokaryotic gene repertoires. These repertoires can be as nomes and highlight the interdependence between prokaryotic small as Ϸ200 genes in intracellular organisms or as large as Ϸ9,000 genome organization and environmental complexity. genes in large, free-living bacteria. In this article we ask what is the impact of HGT from phylogenetically distant sources, relative to Results and Discussion the size of the gene repertoire. Using different approaches for HGT Cumulative Impact of dHGT Increases Nonlinearly in Large Genomes. detection and focusing on both cumulative and recent evolution- Our first approach to infer the contribution of dHGT to the ary histories, we find a surprising pattern of nonlinear enrichment size of bacterial gene repertoires builds on the reconciliation of long-distance transfers in large genomes. Moreover, we find a of the presence–absence distributions of gene families with the strong positive correlation between the sizes of the donor and tree of life (TOL) (13). This method can detect transfers along recipient genomes. Our results also show that distant horizontal the whole history of descent, as long as they explain patchy transfers are biased toward those functional groups that are presence–absence distributions.
    [Show full text]