riboSeed : leveraging prokaryotic genomic architecture to assemble across ribosomal regions

Waters, Nicholas R. and Abram, Florence and Brennan, Fiona and Holmes, Ashleigh and Pritchard, Leighton (2018) riboSeed : leveraging prokaryotic genomic architecture to assemble across ribosomal regions. Nucleic Acids Research, 46 (11). e68. ISSN 0305-1048 (https://doi.org/10.1093/nar/gky212)

[thumbnail of Waters-etal-NAR-2018-leveraging-prokaryotic-genomic-architecture-to-assemble-across-ribosomal-regions]
Preview
Text. Filename: Waters_etal_NAR_2018_leveraging_prokaryotic_genomic_architecture_to_assemble_across_ribosomal_regions.pdf
Final Published Version
License: Creative Commons Attribution 4.0 logo

Download (1MB)| Preview

Abstract

The vast majority of bacterial genome sequencing has been performed using Illumina short reads. Because of the inherent difficulty of resolving repeated regions with short reads alone, only ∼10% of sequencing projects have resulted in a closed genome. The most common repeated regions are those coding for ribosomal operons (rDNAs), which occur in a bacterial genome between 1 and 15 times, and are typically used as sequence markers to classify and identify bacteria. Here, we exploit the genomic context in which rDNAs occur across taxa to improve assembly of these regions relative to de novo sequencing by using the conserved nature of rDNAs across taxa and the uniqueness of their flanking regions within a genome. We describe a method to construct targeted pseudocontigs generated by iteratively assembling reads that map to a reference genome's rDNAs. These pseudocontigs are then used to more accurately assemble the newly sequenced chromosome. We show that this method, implemented as riboSeed, correctly bridges across adjacent contigs in bacterial genome assembly and, when used in conjunction with other genome polishing tools, can assist in closure of a genome.