Genome sequencing of ovine isolates of Mycobacterium avium subspecies paratuberculosis offers insights into host association
© Bannantine et al; licensee BioMed Central Ltd. 2012
Received: 7 October 2011
Accepted: 12 March 2012
Published: 12 March 2012
The genome of Mycobacterium avium subspecies paratuberculosis (MAP) is remarkably homogeneous among the genomes of bovine, human and wildlife isolates. However, previous work in our laboratories with the bovine K-10 strain has revealed substantial differences compared to sheep isolates. To systematically characterize all genomic differences that may be associated with the specific hosts, we sequenced the genomes of three U.S. sheep isolates and also obtained an optical map.
Our analysis of one of the isolates, MAP S397, revealed a genome 4.8 Mb in size with 4,700 open reading frames (ORFs). Comparative analysis of the MAP S397 isolate showed it acquired approximately 10 large sequence regions that are shared with the human M. avium subsp. hominissuis strain 104 and lost 2 large regions that are present in the bovine strain. In addition, optical mapping defined the presence of 7 large inversions between the bovine and ovine genomes (~ 2.36 Mb). Whole-genome sequencing of 2 additional sheep strains of MAP (JTC1074 and JTC7565) further confirmed genomic homogeneity of the sheep isolates despite the presence of polymorphisms on the nucleotide level.
Comparative sequence analysis employed here provided a better understanding of the host association, evolution of members of the M. avium complex and could help in deciphering the phenotypic differences observed among sheep and cattle strains of MAP. A similar approach based on whole-genome sequencing combined with optical mapping could be employed to examine closely related pathogens. We propose an evolutionary scenario for M. avium complex strains based on these genome sequences.
KeywordsM. paratuberculosis Evolution Johne's disease Genome Optical mapping
Mycobacterium avium subspecies paratuberculosis (MAP) causes Johne's disease in sheep, cattle, goats and other ruminant animals. This disease is chronic in nature with multiple years separating the initial infection from clinical stages of disease . The details of the pathogenic mechanisms occurring during this long incubation period still need further study, but it has been demonstrated that MAP colonizes the small intestine through invasion of both M cells and epithelial cells . The disease is of considerable economic significance to livestock industries, particularly the dairy industry.
Generally, MAP is a genetically homogenous subspecies, especially among bovine, human and wildlife isolates [3–5]. However, three lineages of MAP have emerged following extensive molecular strain typing and comparative genomic studies-type I and type III strains (ovine) and type II (bovine) strains. The type III strains were originally called intermediate strains and are highly similar genetically, and thus, difficult to distinguish from type I strains. Early on, the type I (MAP-S) and type II (MAP-C) strains were distinguished based on their molecular fingerprints using IS1311 polymorphism , representational difference analysis , MLSSR typing [8–10] and hsp65 sequencing . On the other hand, type III (a sub-lineage of the MAP-S strains) was genotyped based on gyrA and gyrB genes .
In addition to these recently published genotypic distinctions between "S" and "C" strains of MAP, phenotypic differences have been noted since the middle of the last century . More recently, Motiwala et al.  have shown transcriptional changes in human macrophages infected with MAP-C, human and bison isolates induce an anti-inflammatory gene expression pattern, while the MAP-S isolates showed expression of pro-inflammatory cytokines. Furthermore, some of the ovine strains are pigmented . The ovine and bovine strains likewise are distinct in their growth characteristics. The MAP-S strains are more fastidious and slower in their growth rate than the MAP-C counterpart. In contrast to MAP-C strains, the MAP-S strains do not grow readily on Herold's egg yolk media or Middlebrook 7H9 media that is not supplemented with egg yolk . Nutrient limitation will kill MAP-S strains but it is only bacteriostatic for MAP-C strains . On the transcriptional level, RNA extracted in low iron and heat stressed environments is divergent between MAP-S and MAP-C strains . Recently, iron storage in low iron conditions was only observed in the MAP-C strains but not MAP-S strains . Because of these well-documented phenotypic differences, we hypothesized that sequencing of the genomes of ovine isolates and comparing them to other genomes in the MAC group could provide some clues for these host-specific variations.
The MAP-C strain K-10 was sequenced in 2005 to obtain a complete genome 4.8 Mb in size . It was subsequently found to possess an inversion due to misalignment that was resolved by optical mapping . Very recently, draft sequences of ten MAP isolates have been reported with the presence of two large duplications, especially among human isolates . Finally, another M. avium subspecies (strain 104) has also been sequenced but not published as yet. This genome of subspecies hominissuis is 5.4 Mb in size and greater than 95% homologous to the MAP K-10 genome [3, 5, 22]. Both of these genomes have served as reference genomes in the current project to assist in assembly, open reading frame (ORF) predictions, and annotation. With the help of next-generation sequencing and optical mapping, we were able to assemble a draft of the standard sheep strain of MAP S397 and compare its sequence to other clinical isolates from sheep or the K-10 strain. Interestingly, several inversion regions and single nucleotide polymorphisms distinguished the MAP-S strains from their MAP-C counterpart. Insights into the evolution of MAP strains have been gained through this analysis.
Genome general features
A summary of the genomic features of M. avium subspecies isolates from different hosts
Genome size (bp)
G + C (%)
Protein coding (%)
Total protein coding genes (PCG)
PCG without function prediction
PCG connected to KEGG pathways
Reference genome assembly of clinical ovine isolates using simulated MAP S397 genome
%Homologya to S397
%Homology to K-10
Non-specific matches read countc
Paired read distance distribution
No. of SNPs
The IS elements usually play a role in the genomic diversity among strains of mycobacteria  and could act as a good target for molecular diagnostics . Similar to K-10, the S397 genome has all the well-studied insertion sequences (e.g. IS900, IS1311 and IS_map02). IS900 is generally considered a MAP specific element that was originally discovered in 1989 [26, 27]. A total of 17 copies of IS900 were found in the S397 genome, which is identical to the K-10 strain. Another element, IS_map02, is a MAP specific insertion sequence that was discovered by sequencing the K-10 genome. A total of 6 copies of IS_map02 are present in both S397 and K-10. Likewise, IS1311 is present 7 times in each genome. No IS elements were found to be unique to one or the other genome.
Organization of the MAP S397 genome
Boundaries and flanking ORFs of aligned segments between MAP K-10 and MAP S397 genomes
K-10 coordinates by Wynne
K-10 coordinates by Li
Approx. size (Kb)
Alignment between K-10 and S397
IS_MAP03 and IS1311
IS_MAP03 and IS1311
No known IS element
No known IS element
No known IS element
No known IS element
No known IS element
No known IS element
No known IS element
No known IS element
Large sequences present in the three sheep strain genomes but absent in MAP K-10.
Several new or only partially described LSPs common to MAP S397 and MAH 104 strains were also identified. A good example here is the novel LSP found in MAP sheep and MAH 104 genomes is comprised of 14 ORFs (MAPs_17580 - MAPs_17710), predicted to encode proteins involved in the biosynthesis of glycopeptidolipids . This region in MAP S397 revealed the presence of four additional ORFs (hyp, hlpA, dhgA and mtfC) with homology to glycopeptidolipid biosynthesis genes immediately downstream. The additional 4 ORFs were also not present in the MAH 104 sequence. Finally, a putative transcriptional regular labeled as MAPs_44910 is present in MAP S397. The protein encoded by this ORF has homology to the GntR-family of transcriptional regulators, which are widely distributed across bacterial species and regulate a variety of cellular processes [31, 32].
A second subset of sequence polymorphism was represented by 32 ORFs that were present in the MAP K-10 genome but absent from the genome of MAP S397 (Additional file 6: Figure S1). Several of these deletions have already been described earlier. The deletion encompassing MAP1485c-MAP1491 was previously identified by Marsh and coworkers as S strain deletion #1 in an Australian MAP sheep isolate  and by Semret and coworkers as LSPA20 . An additional larger deletion in the MAP S397 included the cluster of ORFs between MAP1728c and MAP1744. This deletion was partially identified by Marsh and coworkers as RDA3 , and later fully described as S deletion #2 .
A novel deletion comprising the ORFs MAP1432-MAP1438c (partial) was identified in the current study as absent from MAP S397. This deletion, termed sΔ-1, was originally discovered by comparative genomic analysis and subsequently confirmed by PCR analysis. This gene cluster is predicted to encode four energy metabolism enzymes as well as a lipase (MAP1438c). MAP1432 encodes a hypothetical protein with homology to the REP13E12, a family of repetitive elements that were originally described in M. tuberculosis and have been shown to be targets of phage integration . There is a homolog to MAP1434 that is present in S397 (MAPs_13210). The region around MAPs_13210 is not near the end of a contig and is nearly identical to an inverted stretch in K-10, thus leading to the conclusion that MAPs_13210 is only a homolog of MAP1434, but that the gene itself is not present in the S397 genome.
Interestingly, MAP2656 was initially identified as absent via microarray analysis  but sequencing of MAP S397 identified a homologue with 100% identity (MAPs_10401 & MAPs_10402). Likewise, MAP2325 was identified as being absent from Australian sheep isolates of MAP . This ORF was not identified as missing from MAP S397 as sequencing confirmed the presence of an ORF (MAPs_34380) with 100% identity to MAP2325. These discrepancies may represent a geographic difference between MAP isolates recovered from sheep in Australia and the United States or it may be an error from the microarray experiment. These were the only observed differences between the microarray and sequence data. Overall, genomic alignments indicated the presence of a significant number of insertions and deletions between ovine and bovine strains of MAP that are suggested to be associated with their respective host.
Evolutionary analysis of the MAP S397 genome
Comparative genomic hybridizations using DNA microarrays have revealed large sequence polymorphisms (LSPs) between MAP-S and MAP-C strains [36, 41]. Two large deletions of an Australian sheep isolate were found by genomic hybridization to the MAP K-10 array . One deletion encompassed 8 ORFs (MAP1485c-MAP1491) and a second deletion encompassed 17 ORFs extending from MAP1728c to MAP1744. These deletions relative to the bovine strains were later observed in U.S. ovine MAP isolates [5, 13]. Construction of a MAP array containing MAH sequences revealed LSPs in the ovine strains that were missing in the bovine K-10 strain [5, 42]. These documented differences formed the basis for whole-genome sequencing of a sheep isolate to enable comprehensive description of all genetic differences from MAP-S and MAP-C strains. We took advantage of next-generation sequencing technology combined with optical mapping  to decipher the complete genome of MAP isolates from sheep flocks raised in the USA. Our analysis confirmed earlier polymorphisms among MAP-S and MAP-C strains and revealed novel regions of difference. Surprisingly, both genome sequencing and optical mapping showed remarkable differences between MAP-S and MAP-C strains despite the overall similarity in the clinical signs of Johne's disease in sheep and cattle. Recently, a study using a large number of MAP isolates provided an example of such a genomic polymorphism including 2 large regions of duplication, termed vGI-17 (containing 63 ORFs) and vGI-18 (containing 109 ORFs), observed in most MAP-C strains but not MAP-S isolates . Both of these duplications were also missing in our sequenced MAP-S genomes as determined by PCR amplification using outward facing primers reported by Wynne et al. (data not shown).
There are 70 genes present in all three ovine isolates that are absent from the K-10 strain, an indication for MAP adaptation to specific hosts (in this case sheep). Analysis of additional ovine and bovine isolates is needed to strengthen any linkage between these genes with host association. Within this subset, we identified a surprising number of genes annotated as hypothetical proteins (N = 30). Six transcriptional regulators were also present among these genes with the remaining genes showing weak homology to sequences in the GenBank database. We hypothesize that these genes could be responsible for the observed phenotypic differences between ovine and bovine strains and warrant future studies to address this hypothesis.
Genome sequencing of MAP-S strains have revealed extensive genome inversions and previously characterized deletions when compared to the K-10 strain. Furthermore, there appears to be a high degree of homology within US MAP-S strains as suggested by the remarkably low number of SNPs present in the three isolates sequenced. Evolutionary analysis based on whole genome sequencing suggests MAH is the progenitor strain, followed by MAP-S, followed by MAP-C strains.
Overall, Next-generation sequencing combined with optical mapping provided us with a high resolution tool to decipher the evolution of important pathogenic mycobacteria. Comparative sequence analysis of the MAP isolates from sheep has improved our understanding of the evolutionary history of members of MAC and provided the foundation for novel insights into the pathogenesis of this important pathogen. Similar approaches can be used to examine other closely related pathogens.
MAP ovine isolates
Isolates were cultured in Middlebrook 7H9 broth (BD Biosciences, San Jose, CA) media supplemented with 10% OADC (2% glucose, 5% bovine serum albumin factor V, and 0.85% NaCl), 0.05% Tween 80 and 2 μg/ml of Mycobactin J at 37°C . The MAP ovine S397 strain was obtained from a Suffolk breed in Iowa. It was isolated from the distal ileum at necropsy in 2004. The other 2 sheep isolates of MAP (JTC1074 and JTC7565) were isolated from the intestine of infected sheep in Texas and obtained from the Johne's Testing Center at the University of Wisconsin-Madison. All isolates were genotyped using the IS1311 restriction endonuclease, which yielded the 2-band pattern typical of ovine strains .
Genomic DNA was extracted as described in detail previously [3, 46]. For the S397 strain, the DNA (1-5 μg) was sequenced using Roche 454 pyrosequencing (GS20 and FLX) at the National Animal Disease Center. A whole-genomic shotgun sequencing library was prepared according to Roche protocols. The library was used with the appropriate emulsion based PCR kits to produce sufficient beads for sequencing using the Roche Standard Chemistry GS-LR 70 sequencing kit. For the JTC1074 and JTC7565, the purified genomic DNA (~5 μg) of each strain was sent to Genomic Resource Center at the University of Maryland for Illumina whole genome sequencing (Multiplexing Sample Preparation oligonucleotide Kit) as outline before . The adapters and indexing oligonucleotides were purchased from Illumina (5 Paired End Cluster Generation Kits-v4). The CLC Genomic Workbench software (version 4.0.3) was used to perform reference and de novo assembly on all sequenced genomes.
The S397 sequence was annotated using the Integrated Microbial Genomes Expert Review (IMG-ER) pipeline . The sequences of the JTC isolates were annotated based on S397. Genes were each designated by the locus tag "MAPs" to distinguish it as a MAP sheep strain gene. This locus tag is followed by a five digit unique identifier, which incrementally increases by ten (i.e. MAPs_45660... MAPs_45670... MAPs_45680...). With this numbering configuration, additional genes can easily be added as they are discovered or when remaining gaps are closed.
The genome data for MAP K-10 (accession no. GenBank: NC_002944.2) and M. avium subsp. hominissuis (MAH) strain 104 (GenBank: NC_008595.1) were used in alignments in the Artemis and Artemis Comparison Tool (ACT) programs or Mauve 2.3.1 . BLASTP analysis was used for similarity searches and protein sequence analysis. In addition, Mauve algorithm was used to align two or more genomes . For detecting single nucleotide polymorphisms (SNPs) among sheep isolates, the CLC Genomic Workbench was used. The coverage range setting for each strain was at 10-55 reads, and the frequency of the mutation was at least in 50% of the reads.
Shotgun optical mapping, as previously described [20, 51–55], was used to construct a physical restriction map for the S397 genome. Genomic DNA, in agarose inserts , was electroeluted into a solution containing a lambda DNA sizing standard (30 pg/μl), and then were mounted on cleaned, derivitized glass surfaces using a microfluidic device  followed by polymerization of a thin layer of polyacrylamide (3.3% containing 0.02% Triton X-100). Mounted DNA was digested with 20 units of Bsi WI (NEB, Ipswich, MA) for 1 to 2 hrs at 37°C. Fluorochrome-stained DNA fragments were imaged by fluorescence microscopy with a 63 × objective lens (Carl Zeiss, Thornwood, NY) and a high-resolution digital camera (Princeton Instruments, Trenton, NJ). Images were acquired and processed using "ChannelCollect" and "Pathfinder" -custom software  that converts captured images into map data sets. Bayesian inference and an efficient dynamic programming algorithm were also being used to fine-tune the parameters including standard deviation, digestion rate, false cut, and false match probability etc. [54, 58, 59]. The final circular optical map contig was built using an iterative assembly process  including rounds of pair-wise alignment (single molecule maps vs. seed maps; provisional assemblies) and assembly [52, 54]. Due to the high G + C content of MAP, which skews fragment sizing by integrated fluorescence intensity measurement, the final maps were globally scaled (0.95) to correct this problem [20, 61]. A laboratory software implementation of an optical map alignment algorithm  was used to align between optical fragments generated from MAP S397 and the in silico restriction maps of MAP K-10, which provided a whole-genome rearrangement comparison between the two genomes. This restriction framework was used to generate a temporary rearranged genome as the reference sequence to guide the assembly of MAP S397 de novo contigs with the function "move contigs" in Mauve 2.3.1 .
PCR amplification of inversion breakpoints and deletions
PCR reactions were performed in 25-μl reaction mixture containing 1 M betaine (Sigma-Aldrich, St. Louis, MO), 50 mM potassium glutamate (Sigma-Aldrich), 10 mM Tris-HCl pH 8.8, 0.1% Triton X-100, 2 mM magnesium chloride, 0.2 mM dNTPs, 0.5 μM each primer, 0.5 U of GoTaq® Flexi DNA Polymerase (Promega, Madison, WI) and 25 ng of genomic DNA. The amplification thermocycle started with an initial step of 94°C for 5 minutes followed by 5 cycles of 94°C for 30 s, 62°C for 30 s with 1°C decrease for each cycle and 72°C for 3.5 min, and followed by 30 cycles of 94°C for 30 s, 57°C for 30 s and 72°C for 3.5 min. PCR primers used for examining the breakpoints included control F: AAGCATCACCTGCATGAGC, control R: CGGGAATTTATCCGTTTCAG, F1: GGGATCGATCTTGACCACAT, R1: GTGCCTGGACTCGATTTTGT, F2: AAGAGGTCGGAGGTTCGAGT and R2: CGGTGAGAGATTTCGTCACA. Primers used to demonstrate the S397 sΔ-1 deletion included F18: CGTCTTCCCCGTCGTCGTTC, B24: CGATGAGAGTCCGTGCGTGG, F15: CGGCGGGCGGTCAGGGTTTG, B17: GCAGGTTGGGGTTCGGCTTG, F7: GGTGGTCGGCGTCCTCGTAG, B9: CGTCGTCACAGCGAAAACGG, F3: CCACCCGCCTCACACCACTC, B4: AGGACGCCGACCACCAAACG. Conditions for the amplifications are essentially as described immediately above except that Advantage GC Genomic LA PCR Polymerase kit (Clontech) was used for each reaction.
Nucleotide sequence accession number
This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession AFIF00000000.
The authors would like to thank members of the Genomic Resource Center at the University of Maryland-Baltimore for Illumina sequencing and Janis K. Hansen (USDA-ARS) for technical assistance. This work was supported by the USDA-Agricultural Research Service (JPB, MLP, DPA and DOB), NRI 2007-35204-18400 and JDIP -Q6286224301 grants from the USDA and US-Egypt Joint Scientific Baord#1937 to AMT.
- Wu CW, Livesey M, Schmoller SK, Manning EJ, Steinberg H, Davis WC, et al: Invasion and persistence of Mycobacterium avium subsp. paratuberculosis during early stages of Johne's disease in calves. Infect Immun. 2007, 75: 2110-2119. 10.1128/IAI.01739-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Bermudez LE, Petrofsky M, Sommer S, Barletta RG: Peyer's patch-deficient mice demonstrate that Mycobacterium avium subsp. paratuberculosis translocates across the mucosal barrier via both M cells and enterocytes but has inefficient dissemination. Infect Immun. 2010, 78: 3570-3577. 10.1128/IAI.01411-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu CW, Glasner J, Collins M, Naser S, Talaat AM: Whole-genome plasticity among Mycobacterium avium subspecies: insights from comparative genomic hybridizations. J Bacteriol. 2006, 188: 711-723. 10.1128/JB.188.2.711-723.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Taylor AW: Varieties of Mycobacterium johnei isolated from sheep. J Pathol Bacteriol. 1951, 63: 333-336. 10.1002/path.1700630217.View ArticlePubMedGoogle Scholar
- Paustian ML, Zhu X, Sreevatsan S, Robbe-Austerman S, Kapur V, Bannantine JP: Comparative genomic analysis of Mycobacterium avium subspecies obtained from multiple host species. BMC Genomics. 2008, 9: 135-149. 10.1186/1471-2164-9-135.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh I, Whittington R, Cousins D: PCR-restriction endonuclease analysis for identification and strain typing of Mycobacterium avium subsp. paratuberculosis and Mycobacterium avium subsp. avium based on polymorphisms in IS1311. Mol Cell Probes. 1999, 13: 115-126. 10.1006/mcpr.1999.0227.View ArticlePubMedGoogle Scholar
- Dohmann K, Strommenger B, Stevenson K, de JL, Stratmann J, Kapur V, et al: Characterization of genetic differences between Mycobacterium avium subsp. paratuberculosis type I and type II isolates. J Clin Microbiol. 2003, 41: 5215-5223. 10.1128/JCM.41.11.5215-5223.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Amonsin A, Li LL, Zhang Q, Bannantine JP, Motiwala AS, Sreevatsan S, et al: Multilocus short sequence repeat sequencing approach for differentiating among Mycobacterium avium subsp. paratuberculosis strains. J Clin Microbiol. 2004, 42: 1694-1702. 10.1128/JCM.42.4.1694-1702.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Sevilla I, Li L, Amonsin A, Garrido JM, Geijo MV, Kapur V, et al: Comparative analysis of Mycobacterium avium subsp. paratuberculosis isolates from cattle, sheep and goats by short sequence repeat and pulsed-field gel electrophoresis typing. BMC Microbiol. 2008, 8: 204-10.1186/1471-2180-8-204.PubMed CentralView ArticlePubMedGoogle Scholar
- Thibault VC, Grayon M, Boschiroli ML, Willery E, Allix-Beguec C, Stevenson K, et al: Combined multilocus short sequence repeat and mycobacterial interspersed repetitive unit- variable-number tandem repeat typing of Mycobacterium avium subsp. paratuberculosis isolates. J Clin Microbiol. 2008, 46: 4091-4094. 10.1128/JCM.01349-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Turenne CY, Semret M, Cousins DV, Collins DM, Behr MA: Sequencing of hsp6 distinguishes among subsets of the Mycobacterium aviu complex. J Clin Microbiol. 2006, 44: 433-440. 10.1128/JCM.44.2.433-440.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Castellanos E, Juan Ld, Domínguez L, Aranaz A: Progress in molecular typing of Mycobacterium avium subspecies paratuberculosis. Res Vet Sci. 2011, doi:10.1016/j.rvsc.2011.05.017Google Scholar
- Motiwala AS, Janagama HK, Paustian ML, Zhu X, Bannantine JP, Kapur V, et al: Comparative transcriptional analysis of human macrophages exposed to animal and human isolates of Mycobacterium avius subspecies paratuberculosis with diverse genotypes. Infect Immun. 2006, 74: 6046-6056. 10.1128/IAI.00326-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Stevenson K, Hughes VM, de JL, Inglis NF, Wright F, Sharp JM: Molecular characterization of pigmented and nonpigmented isolates of Mycobacterium avium subsp. paratuberculosis. J Clin Microbiol. 2002, 40: 1798-1804. 10.1128/JCM.40.5.1798-1804.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Whittington RJ, Marsh IB, Saunders V, Grant IR, Juste R, Sevilla IA, et al: Culture Phenotypes of Genomically and Geographically Diverse Mycobacterium avium subsp. paratuberculosis Isolates from Different Hosts. J Clin Microbiol. 2011, 49: 1822-1830. 10.1128/JCM.00210-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Gumber S, Taylor DL, Marsh IB, Whittington RJ: Growth pattern and partial proteome of Mycobacterium avium subsp. paratuberculosis during the stress response to hypoxia and nutrient starvation. Vet Microbiol. 2009, 133: 344-357. 10.1016/j.vetmic.2008.07.021.View ArticlePubMedGoogle Scholar
- Gumber S, Whittington RJ: Analysis of the growth pattern, survival and proteome of Mycobacterium avium subsp. paratuberculosis following exposure to heat. Vet Microbiol. 2009, 136: 82-90. 10.1016/j.vetmic.2008.10.003.View ArticlePubMedGoogle Scholar
- Janagama HK, Senthilkumar TM, Bannantine JP, Rodriguez GM, Smith I, Paustian ML, et al: Identification and functional characterization of the iron-dependent regulator (IdeR) of Mycobacterium avium subsp. paratuberculosis. Microbiology. 2009, 155: 3683-3690. 10.1099/mic.0.031948-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Li L, Bannantine JP, Zhang Q, Amonsin A, May BJ, Alt D, et al: The complete genome sequence of Mycobacterium avium subspecies paratuberculosis. Proc Natl Acad Sci USA. 2005, 102: 12344-12349. 10.1073/pnas.0505662102.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu CW, Schramm TM, Zhou S, Schwartz DC, Talaat AM: Optical mapping of the Mycobacterium avium subspecies paratuberculosis genome. BMC Genomics. 2009, 10: 25-10.1186/1471-2164-10-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Wynne JW, Bull TJ, Seemann T, Bulach DM, Wagner J, Kirkwood CD, et al: Exploring the zoonotic potential of Mycobacterium avium subspecies paratuberculosis through comparative genomics. PLoS One. 2011, 6: e22171-10.1371/journal.pone.0022171.PubMed CentralView ArticlePubMedGoogle Scholar
- Bannantine JP, Zhang Q, Li LL, Kapur V: Genomic homogeneity between Mycobacterium avium subsp. avium and Mycobacterium avium subsp. paratuberculosis belies their divergent growth rates. BMC Microbiol. 2003, 3: 10-10.1186/1471-2180-3-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R, et al: PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics. 2010, 26: 1608-1615. 10.1093/bioinformatics/btq249.PubMed CentralView ArticlePubMedGoogle Scholar
- Sreevatsan S, Pan X, Stockbauer KE, Connell ND, Kreiswirth BN, Whittam TS, et al: Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci USA. 1997, 94: 9869-9874. 10.1073/pnas.94.18.9869.PubMed CentralView ArticlePubMedGoogle Scholar
- Nodieva A, Jansone I, Broka L, Pole I, Skenders G, Baumanis V: Recent nosocomial transmission and genotypes of multidrug-resistant Mycobacterium tuberculosis. Int J Tuberc Lung Dis. 2010, 14: 427-433.PubMedGoogle Scholar
- Collins DM, Gabric DM, de Lisle GW: Identification of a repetitive DNA sequence specific to Mycobacterium paratuberculosis. FEMS Microbiol Lett. 1989, 60: 175-178. 10.1111/j.1574-6968.1989.tb03440.x.View ArticleGoogle Scholar
- Green EP, Tizard ML, Moss MT, Thompson J, Winterbourne DJ, McFadden JJ, et al: Sequence and characteristics of IS900, an insertion element identified in a human Crohn's disease isolate of Mycobacterium paratuberculosis. Nucleic Acids Res. 1989, 17: 9063-9073. 10.1093/nar/17.22.9063.PubMed CentralView ArticlePubMedGoogle Scholar
- Wynne JW, Seemann T, Bulach D, Coutts SA, Talaat AM, Michalski WP: Re-sequencing the Mycobacterium avium subsp. paratuberculosis K10 genome: improved annotation and revised genome sequence. J Bacteriol. 2010, 192: 6319-6320. 10.1128/JB.00972-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Semret M, Turenne CY, de HP, Collins DM, Behr MA: Differentiating host-associated variants of Mycobacterium avium by PCR for detection of large sequence polymorphisms. J Clin Microbiol. 2006, 44: 881-887. 10.1128/JCM.44.3.881-887.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Eckstein TM, Belisle JT, Inamine JM: Proposed pathway for the biosynthesis of serovar-specific glycopeptidolipids in Mycobacterium avium serovar 2. Microbiology. 2003, 149: 2797-2807. 10.1099/mic.0.26528-0.View ArticlePubMedGoogle Scholar
- Haydon DJ, Guest JR: A new family of bacterial regulatory proteins. FEMS Microbiol Lett. 1991, 63: 291-295.View ArticlePubMedGoogle Scholar
- Vindal V, Suma K, Ranjan A: GntR family of regulators in Mycobacterium smegmatis: a sequence and structure based characterization. BMC Genomics. 2007, 8: 289-10.1186/1471-2164-8-289.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh IB, Bannantine JP, Paustian ML, Tizard ML, Kapur V, Whittington RJ: Genomic comparison of Mycobacterium avium subsp. paratuberculosis sheep and cattle strains by microarray hybridization. J Bacteriol. 2006, 188: 2290-2293. 10.1128/JB.188.6.2290-2293.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh IB, Whittington RJ: Deletion of an mmp gene and multiple associated genes from the genome of the S strain of Mycobacterium avium subsp. paratuberculosis identified by representational difference analysis and in silic analysis. Mol Cell Probes. 2005, 19: 371-384. 10.1016/j.mcp.2005.06.005.View ArticlePubMedGoogle Scholar
- Gordon SV, Brosch R, Billault A, Garnier T, Eiglmeier K, Cole ST: Identification of variable regions in the genomes of tubercle bacilli using bacterial artificial chromosome arrays. Mol Microbiol. 1999, 32: 643-655. 10.1046/j.1365-2958.1999.01383.x.View ArticlePubMedGoogle Scholar
- Alexander DC, Turenne CY, Behr MA: Insertion and deletion events that define the pathogen Mycobacterium avium subsp. paratuberculosis. J Bacteriol. 2009, 191: 1018-1025. 10.1128/JB.01340-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Ghosh P, Hsu C-Y, Alyamani E, Shehata MM, Al-Dubaib MA, Al-Naeem A, et al: Genome-wide Analysis of the Emerging Infection with Mycobacterium avium subspecies paratuberculosis in the Arabian Camels (Camelus dromedarius). PLoS ONE. 2012,Google Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Tamura K, Kumar S: Evolutionary distance estimation under heterogeneous substitution pattern among lineages. Mol Biol Evol. 2002, 19: 1727-1736. 10.1093/oxfordjournals.molbev.a003995.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralView ArticlePubMedGoogle Scholar
- Semret M, Alexander DC, Turenne CY, de HP, Overduin P, van Soolingen D, et al: Genomic polymorphisms for Mycobacterium avium subsp. paratuberculosis diagnostics. J Clin Microbiol. 2005, 43: 3704-3712. 10.1128/JCM.43.8.3704-3712.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Castellanos E, Aranaz A, Gould KA, Linedale R, Stevenson K, Alvarez J, et al: Discovery of Stable and Variable Differences in the Mycobacterium avium subsp paratuberculosi Type I, II, and III Genomes by Pan-Genome Microarray Analysis. Appl Environ Microbiol. 2009, 75: 676-686. 10.1128/AEM.01683-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Sohal JS, Singh SV, Singh PK, Singh AV: On the evolution of 'Indian Bison type' strains of Mycobacterium avium subspecies paratuberculosis. Microbiol Res. 2010, 165: 163-171. 10.1016/j.micres.2009.03.007.View ArticlePubMedGoogle Scholar
- Darling AC, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14: 1394-1403. 10.1101/gr.2289704.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu CW, Schmoller SK, Shin SJ, Talaat AM: Defining the stressome of Mycobacterium avium subsp. paratuberculosis in vitro and in naturally infected cows. J Bacteriol. 2007, 189: 7877-7886. 10.1128/JB.00780-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Bannantine JP, Baechler E, Zhang Q, Li L, Kapur V: Genome scale comparison of Mycobacterium avium subsp. paratuberculosis with Mycobacterium avium subsp. avium reveals potential diagnostic sequences. J Clin Microbiol. 2002, 40: 1303-1310. 10.1128/JCM.40.4.1303-1310.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Hegedus Z, Zakrzewska A, Agoston VC, Ordas A, Racz P, Mink M, et al: Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol Immunol. 2009, 46: 2918-2930. 10.1016/j.molimm.2009.07.002.View ArticlePubMedGoogle Scholar
- Markowitz VM, Mavromatis K, Ivanova NN, Chen IM, Chu K, Kyrpides NC: IMG ER: a system for microbial genome annotation expert review and curation. Bioinformatics. 2009, 25: 2271-2278. 10.1093/bioinformatics/btp393.View ArticlePubMedGoogle Scholar
- Darling AE, Mau B, Perna NT: progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010, 5: e11147-10.1371/journal.pone.0011147.PubMed CentralView ArticlePubMedGoogle Scholar
- Perna NT, Mayhew GF, Posfai G, Elliott S, Donnenberg MS, Kaper JB, et al: Molecular evolution of a pathogenicity island from enterohemorrhagic Escherichia coli O157:H7. Infect Immun. 1998, 66: 3810-3817.PubMed CentralPubMedGoogle Scholar
- Lin J, Qi R, Aston C, Jing J, Anantharaman TS, Mishra B, et al: Whole-genome shotgun optical mapping of Deinococcus radiodurans. Science. 1999, 285: 1558-1562. 10.1126/science.285.5433.1558.View ArticlePubMedGoogle Scholar
- Zhou S, Bechner MC, Place M, Churas CP, Pape L, Leong SA, et al: Validation of rice genome sequence by optical mapping. BMC Genomics. 2007, 8: 278-295. 10.1186/1471-2164-8-278.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou S, Deng W, Anantharaman TS, Lim A, Dimalanta ET, Wang J, et al: A whole-genome shotgun optical map of Yersinia pestis strain KIM. Appl Environ Microbiol. 2002, 68: 6321-6331. 10.1128/AEM.68.12.6321-6331.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou S, Kile A, Kvikstad E, Bechner M, Severin J, Forrest D, et al: Shotgun optical mapping of the entire Leishmania major Friedlin genome. Mol Biochem Parasitol. 2004, 138: 97-106. 10.1016/j.molbiopara.2004.08.002.View ArticlePubMedGoogle Scholar
- Zhou S, Wei F, Nguyen J, Bechner M, Potamousis K, Goldstein S, et al: A single molecule scaffold for the maize genome. PLoS Genet. 2009, 5: e1000711-10.1371/journal.pgen.1000711.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwartz DC, Cantor CR: Separation of yeast chromosome-sized DNAs by pulsed field gradient gel electrophoresis. Cell. 1984, 37: 67-75. 10.1016/0092-8674(84)90301-5.View ArticlePubMedGoogle Scholar
- Dimalanta ET, Lim A, Runnheim R, Lamers C, Churas C, Forrest DK, et al: A microfluidic system for large DNA molecule arrays. Anal Chem. 2004, 76: 5293-5301. 10.1021/ac0496401.View ArticlePubMedGoogle Scholar
- Lim A, Dimalanta ET, Potamousis KD, Yen G, Apodoca J, Tao C, et al: Shotgun optical maps of the whole Escherichia coli O157:H7 genome. Genome Res. 2001, 11: 1584-1593. 10.1101/gr.172101.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou S, Kvikstad E, Kile A, Severin J, Forrest D, Runnheim R, et al: Whole-genome shotgun optical mapping of Rhodobacter sphaeroides strain 2.4.1 and its use for whole-genome shotgun sequence assembly. Genome Res. 2003, 13: 2142-2151. 10.1101/gr.1128803.PubMed CentralView ArticlePubMedGoogle Scholar
- Teague B, Waterman MS, Goldstein S, Potamousis K, Zhou S, Reslewic S, et al: High-resolution genome structure by single molecule analysis. Proc Natl Acad Sci USA. 2010, 107: 10848-10853. 10.1073/pnas.0914638107.PubMed CentralView ArticlePubMedGoogle Scholar
- Lai Z, Jing J, Aston C, Clarke V, Apodaca J, Dimalanta ET, et al: A shotgun optical map of the entire Plasmodium falciparum genome. Nat Genet. 1999, 23: 309-313. 10.1038/15484.View ArticlePubMedGoogle Scholar
- Valouev A, Schwartz DC, Zhou S, Waterman MS: An algorithm for assembly of ordered restriction maps from single DNA molecules. Proc Natl Acad Sci USA. 2006, 103: 15770-15775. 10.1073/pnas.0604040103.PubMed CentralView ArticlePubMedGoogle Scholar