Resequencing the Mycobacterium avium subsp. paratuberculosis K10 Genome: Improved Annotation and Revised Genome Sequence

Abstract

We report the resequencing and revised annotation of the Mycobacterium avium subsp. paratuberculosis K10 genome. A total of 90 single-nucleotide errors and a 51-bp indel in the original K10 genome were corrected, and the whole genome annotation was revised. Correction of these sequencing errors resulted in 28 frameshift alterations. The amended genome sequence is accessible via the supplemental section of study SRR060191 in the NCBI Sequence Read Archive and will serve as a valuable reference genome for future studies.The American bovine isolate K10 remains the only Mycobacterium avium subsp. paratuberculosis genome to be fully sequenced and published to date (1). Although this 4.8-Mbp genome likely contains some assembly errors (3), it has provided, and will continue to provide, an invaluable resource for Mycobacterium research. The assembly errors were identified through optical mapping of related M. avium subsp. paratuberculosis strain ATCC 19698, which revealed a 648-kb inversion around the origin of replication and two additional copies of the insertion sequences IS1311 and IS_MAP03 (3). These findings were subsequently validated via PCR, Southern blotting, and (Sanger) sequence analysis in ATCC 19698 and were also confirmed to be present in K10 (3). We designate this interim corrected genome M. avium subsp. paratuberculosis K10Ј. To further improve this resource, we undertook a resequencing project of the original M. avium subsp. paratuberculosis K10 genome.Whole-genome sequencing was performed on the Illumina GAIIx platform using one flow cell lane with 36-cycle pairedend chemistry. Reads were variably trimmed at the 3Ј end based on the Illumina Read Segment Quality Indicator (Illumina manual), and read pairs containing ambiguous bases were removed. Read mapping onto the K10Ј genome sequence was performed using SHRiMP (ver. 1.3.2) (2), and singlenucleotide polymorphisms and indels (deletion and insertion polymorphisms [DIPs]) were called using Nesoni (ver. 0.29; Monash University Victorian Bioinformatics Consortium) with default parameters. Read mapping determined that the data set comprised an average sequence coverage of 72.6 across the K10Ј genome. This high sequence coverage allowed differences between K10{K10Ј and the resequenced version of the genome, designated K10Љ, to be identified with high confidence.Ninety single-nucleotide differences and one 51-bp indel were identified in the K10Љ genome. As confirmation that these differences are likely to represent errors in the original genome sequence, we have also detected these polymorphisms in two additional bovine M. avium subsp. paratuberculosis genomes recently sequenced and assembled within our laboratory (data not shown). Seven of the 90 differences and the 51-bp indel were subjected to PCR and Sanger sequencing for verification. All of the polymorphisms were confirmed to be present in K10Љ compared to the original genome sequence. Thirty-six single-nucleotide deletions and four nucleotide insertions were identified in K10Љ compared to the refere...

Key Role for the Alternative Sigma Factor, SigH, in the Intracellular Life of Mycobacterium avium subsp. paratuberculosis during Macrophage Stress

“…The second circle shows the location of the open reading frames in the M. avium subsp. paratuberculosis genome (35). Genes (magenta) on the forward strand are shown outside the baseline; genes (olive green) on the reverse strand are shown inside the baseline.…”

COV2HTML: A Visualization and Analysis Tool of Bacterial Next Generation Sequencing (NGS) Data for Postgenomics Life Scientists

“…Nowadays, they allow large-scale comparative and evolutionary studies and understanding regulatory networks. The broadest application of NGS (Shendure and Ji, 2008) has incredibly increased the number of new genomes (Monot et al, 2009;Sekizuka et al, 2011) and rapidly re-sequenced existing genomes previously obtained by the automated Sanger sequencing method (Wynne et al, 2010). The regulation of gene expression being a fundamental process within every cell type, the deep transcriptome sequencing (RNA-seq) has recently emerged as a method enabling the study of RNA-based regulatory mechanisms in a genome-wide manner (Wang et al, 2009).…”

Viable Mycobacterium avium ssp. paratuberculosis isolated from calf milk replacer

“…Interestingly, the percentages of nonsynonymous SNP were much higher than synonymous SNP, and only a few SNP were present in the intergenic regions ( Figure 2B). As expected, a phylogenetic tree of SNP predicted from each genome ( Figure 2C) showed a clear distinction from the MAP_ S397 strain originally isolated from a sheep (Bannantine et al, 2012), but was closely related to the K10 isolate from a cow (Wynne et al, 2010). In addition, genomes from all CMR isolates clustered separately from genomes recently isolated from elk circulating in California (David Press, Wildlife Ecologist, Point Reyes National Seashore, Inverness, CA, personal communications).…”