riboSeed: leveraging prokaryotic genomic architecture to assemble across ribosomal regions

Nicholas R. Waters, Florence Abram, Fiona Brennan, Ashleigh Holmes, Leighton Pritchard

Research output: Contribution to journalArticle

1 Citation (Scopus)

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.

Original languageEnglish
Article numbere68
Number of pages10
JournalNucleic Acids Research
Volume46
Issue number11
Early online date28 Mar 2018
DOIs
Publication statusPublished - 20 Jun 2018

    Fingerprint

Keywords

  • bacterial genome sequencing
  • ribosomal operons (rDNAs)
  • riboSeed

Cite this