ISMapper: identifying transposase insertion sites in bacterial genomes from short read sequence data.
BACKGROUND
Insertion sequences (IS) are small transposable elements, commonly found in bacterial genomes. Identifying the location of IS in bacterial genomes can be useful for a variety of purposes including epidemiological tracking and predicting antibiotic resistance. However IS are commonly present in multiple copies in a single genome, which complicates genome assembly and the identification of IS insertion sites. Here we present ISMapper, a mapping-based tool for identification of the site and orientation of IS insertions in bacterial genomes, directly from paired-end short read data.


RESULTS
ISMapper was validated using three types of short read data: (i) simulated reads from a variety of species, (ii) Illumina reads from 5 isolates for which finished genome sequences were available for comparison, and (iii) Illumina reads from 7 Acinetobacter baumannii isolates for which predicted IS locations were tested using PCR. A total of 20 genomes, including 13 species and 32 distinct IS, were used for validation. ISMapper correctly identified 97 % of known IS insertions in the analysis of simulated reads, and 98 % in real Illumina reads. Subsampling of real Illumina reads to lower depths indicated ISMapper was able to correctly detect insertions for average genome-wide read depths >20x, although read depths >50x were required to obtain confident calls that were highly-supported by evidence from reads. All ISAba1 insertions identified by ISMapper in the A. baumannii genomes were confirmed by PCR. In each A. baumannii genome, ISMapper successfully identified an IS insertion upstream of the ampC beta-lactamase that could explain phenotypic resistance to third-generation cephalosporins. The utility of ISMapper was further demonstrated by profiling genome-wide IS6110 insertions in 138 publicly available Mycobacterium tuberculosis genomes, revealing lineage-specific insertions and multiple insertion hotspots.


CONCLUSIONS
ISMapper provides a rapid and robust method for identifying IS insertion sites directly from short read data, with a high degree of accuracy demonstrated across a wide range of bacteria.
Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, VIC, 3010, Australia. hawkey.jane@gmail.com.|Faculty of Veterinary and Agricultural Science, The University of Melbourne, Parkville, VIC, 3010, Australia. hawkey.jane@gmail.com.
Jane
J
Hawkey
School of Molecular Bioscience, The University of Sydney, Sydney, 2006, Australia. mohammad.hamidian@sydney.edu.au.
Mohammad
M
Hamidian
Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, VIC, 3010, Australia. rwick@student.unimelb.edu.au.
Ryan R
RR
Wick
Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, VIC, 3010, Australia. d.edwards2@student.unimelb.edu.au.
David J
DJ
Edwards
Faculty of Veterinary and Agricultural Science, The University of Melbourne, Parkville, VIC, 3010, Australia. hbj@unimelb.edu.au.
Helen
H
Billman-Jacobe
School of Molecular Bioscience, The University of Sydney, Sydney, 2006, Australia. ruth.hall@sydney.edu.au.
Ruth M
RM
Hall
Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, VIC, 3010, Australia. kholt@unimelb.edu.au.
Kathryn E
KE
Holt
D040981:Acinetobacter baumannii / Q000235:genetics; D001483:Base Sequence; D003198:Computer Simulation; D004352:Drug Resistance, Microbial / Q000235:genetics; D016680:Genome, Bacterial; D008969:Molecular Sequence Data; D016254:Mutagenesis, Insertional / Q000235:genetics*; D009169:Mycobacterium tuberculosis / Q000235:genetics; D016133:Polymerase Chain Reaction; D015203:Reproducibility of Results; D017422:Sequence Analysis, DNA; D012984:Software; D019895:Transposases / Q000378:metabolism*
D019895:Transposases