Ola Brynildsrud1, Lars-Gustav Snipen1, Jon Bohlin1. 1. Section for Biostatistics and Epidemiology, Norwegian University of Life Sciences (NMBU), Oslo, Department of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences (NMBU), Ås and Norwegian Institute of Public Health, Division of Epidemiology, 0403 Oslo, Norway.
Abstract
MOTIVATION: The explosion of whole-genome sequencing (WGS) as a tool in the mapping and understanding of genomes has been accompanied by an equally massive report of tools and pipelines for the analysis of DNA copy number variation (CNV). Most currently available tools are designed specifically for human genomes, with comparatively little literature devoted to CNVs in prokaryotic organisms. However, there are several idiosyncrasies in prokaryotic WGS data. This work proposes a step-by-step approach for detection and quantification of copy number variants specifically aimed at prokaryotes. RESULTS: After aligning WGS reads to a reference genome, we count the individual reads in a sliding window and normalize these counts for bias introduced by differences in GC content. We then investigate the coverage in two fundamentally different ways: (i) Employing a Hidden Markov Model and (ii) by repeated sampling with replacement (bootstrapping) on each individual gene. The latter bypasses the complex problem of breakpoint determination. To demonstrate our method, we apply it to real and simulated WGS data and benchmark it against two popular methods for CNV detection. The proposed methodology will in some cases represent a significant jump in accuracy from other current methods. AVAILABILITY AND IMPLEMENTATION: CNOGpro is written entirely in the R programming language and is available from the CRAN repository (http://cran.r-project.org) under the GNU General Public License.
MOTIVATION: The explosion of whole-genome sequencing (WGS) as a tool in the mapping and understanding of genomes has been accompanied by an equally massive report of tools and pipelines for the analysis of DNA copy number variation (CNV). Most currently available tools are designed specifically for human genomes, with comparatively little literature devoted to CNVs in prokaryotic organisms. However, there are several idiosyncrasies in prokaryotic WGS data. This work proposes a step-by-step approach for detection and quantification of copy number variants specifically aimed at prokaryotes. RESULTS: After aligning WGS reads to a reference genome, we count the individual reads in a sliding window and normalize these counts for bias introduced by differences in GC content. We then investigate the coverage in two fundamentally different ways: (i) Employing a Hidden Markov Model and (ii) by repeated sampling with replacement (bootstrapping) on each individual gene. The latter bypasses the complex problem of breakpoint determination. To demonstrate our method, we apply it to real and simulated WGS data and benchmark it against two popular methods for CNV detection. The proposed methodology will in some cases represent a significant jump in accuracy from other current methods. AVAILABILITY AND IMPLEMENTATION: CNOGpro is written entirely in the R programming language and is available from the CRAN repository (http://cran.r-project.org) under the GNU General Public License.
Authors: Samuel J T Wardell; Attika Rehman; Lois W Martin; Craig Winstanley; Wayne M Patrick; Iain L Lamont Journal: Antimicrob Agents Chemother Date: 2019-09-30 Impact factor: 5.191
Authors: Chandler Roe; Marc Stegger; Berit Lilje; Thor Bech Johannesen; Kim Lee Ng; Raphael N Sieber; Elizabeth Driebe; David M Engelthaler; Paal Skytt Andersen Journal: Microb Genom Date: 2020-07-15
Authors: Jaime E Zlamal; Semen A Leyn; Mallika Iyer; Marinela L Elane; Nicholas A Wong; James W Wamsley; Maarten Vercruysse; Fernando Garcia-Alcalde; Andrei L Osterman Journal: mBio Date: 2021-06-22 Impact factor: 7.867
Authors: James Emmanuel San; Shakuntala Baichoo; Aquillah Kanzi; Yumna Moosa; Richard Lessells; Vagner Fonseca; John Mogaka; Robert Power; Tulio de Oliveira Journal: Front Microbiol Date: 2020-01-30 Impact factor: 5.640