Polymorphisms in gyrA and gyrB genes among Mycobacterium avium subsp. paratuberculosis type I, II, and III isolates

Genetic and chemical control of tuberculostearic acid production in Mycobacterium avium subspecies paratuberculosis

Tuberculostearic acid (TBSA) is a fatty acid unique to mycobacteria and some corynebacteria and has been studied due to its diagnostic value, biofuel properties, and role in membrane dynamics. In this study, we demonstrate that TBSA production can be abrogated either by addition of pivalic acid to mycobacterial growth cultures or by a bfaA gene knockout encoding a flavin adenine dinucleotide (FAD)-binding oxidoreductase. Mycobacterium avium subspecies paratuberculosis (Map) growth and TBSA production were inhibited in 0.5-mg/mL pivalic acid-supplemented cultures, but higher concentrations were needed to have a similar effect in other mycobacteria, including Mycobacterium smegmatis. While Map C-type strains, isolated from cattle and other ruminants, will produce TBSA in the absence of pivalic acid, the S-type Map strains, typically isolated from sheep, do not produce TBSA in any condition. A SAM-dependent methyltransferase encoded by bfaB and FAD-binding oxidoreductase are both required in the two-step biosynthesis of TBSA. However, S-type strains contain a single-nucleotide polymorphism in the bfaA gene, rendering the oxidoreductase enzyme vestigial. This results in the production of an intermediate, termed 10-methylene stearate, which is detected only in S-type strains. Fatty acid methyl ester analysis of a C-type Map bfaA knockout revealed the loss of TBSA production, but the intermediate was present, similar to the S-type strains. Collectively, these results demonstrate the subtle biochemical differences between two primary genetic lineages of Map and other mycobacteria as well as explain the resulting phenotype at the genetic level. These data also suggest that TBSA should not be used as a diagnostic marker for Map.IMPORTANCEBranched-chain fatty acids are a predominant cell wall component among species belonging to the Mycobacterium genus. One of these is TBSA, which is a long-chain middle-branched fatty acid used as a diagnostic marker for Mycobacterium tuberculosis. This fatty acid is also an excellent biolubricant. Control of its production is important for industrial purposes as well as understanding the biology of mycobacteria. In this study, we discovered that a carboxylic acid compound termed pivalic acid inhibits TBSA production in mycobacteria. Furthermore, Map strains from two separate genetic lineages (C-type and S-type) showed differential production of TBSA. Cattle-type strains of Mycobacterium avium subspecies paratuberculosis produce TBSA, while the sheep-type strains do not. This important phenotypic difference is attributed to a single-nucleotide deletion in sheep-type strains of Map. This work sheds further light on the mechanism used by mycobacteria to produce tuberculostearic acid.

Development and evaluation of genomics informed real-time PCR assays for the detection and strain typing of Mycobacterium avium subsp. paratuberculosis

AIMS

This study aimed to identify specific genomic targets for the detection and strain typing of Map and analyse their sensitivity and specificity, and detect Map directly from faeces.

METHODS AND RESULTS

A comparative genomics approach was used to identify specific genomic targets for the detection and strain typing of Map. A Map specific qPCR using the primer pair 7132 that targets a DNA segregation ATPase protein was able to detect all strains of Map and is more sensitive than the current Johne's disease PCR assays with a sensitivity of 0.0002 fg µl-1. A strain specific qPCR using the Atsa primer pair that targets the arylsulfase gene was able to differentiate between Type S and Type C strains of Map and was more sensitive than the IS1311 PCR and REA with a sensitivity of 40 fg µl-1 and was specific for Type S Map. Both assays successfully detected Map directly from faeces.

CONCLUSION

This study developed and validated two genomics informed qPCR assays, 7132B Map and Atsa Type S and found both assays to be highly specific and sensitive for the detection of Map from culture and directly from faeces. This is the first time that a probe-based qPCR has been designed and developed for Map strain typing, which will greatly improve the response time during outbreak investigations.

Characterisation of IS1311 in Mycobacterium avium subspecies paratuberculosis genomes: Typing, continental clustering, microbial evolution and host adaptation

Johne's disease (JD), caused by Mycobacterium avium subspecies paratuberculosis (MAP) is a global burden for livestock producers and has an association with Crohn's disease in humans. Within MAP there are two major lineages, S/Type I/TypeIII and C/Type II, that vary in phenotype including culturability, host preference and virulence. These lineages have been identified using the IS1311 element, which contains a conserved, single nucleotide polymorphism. IS1311 and the closely related IS1245 element belong to the IS256 family of insertion sequences, are dispersed throughout M. avium taxa but remain poorly characterised. To investigate the distribution and diversity of IS1311 in MAP, 805 MAP genomes were collated from public databases. IS1245 was absent, while IS1311 sequence, copy number and insertion loci were conserved between MAP S lineages and varied within the MAP C lineage. One locus was specific to the S strains, which contained nine IS1311 copies. In contrast, C strains contained either seven or eight IS1311 loci. Most insertion loci were associated with the boundaries of homologous regions that had undergone genome rearrangement between the MAP lineages, suggesting that this sequence may be a driver of recombination. Phylogenomic geographic clustering of MAP subtypes was demonstrated for the first time, at continental scale, and indicated that there may have been recent MAP transmission between Europe and North America, in contrast to Australia where importation of live ruminants is generally prohibited. This investigation confirmed the utility of IS1311 typing in epidemiological studies and resolved anomalies in past studies. The results shed light on potential mechanisms of niche/host adaptation, virulence of MAP and global transmission dynamics.

Improved DNA Amplification of the Hallmark IS900 Element in Mycobacterium avium subsp. paratuberculosis: a Reexamination Based on Whole-Genome Sequence Analysis

Amplification of the IS900 multicopy element is a hallmark nucleic acid-based diagnostic test for Mycobacterium avium subsp. paratuberculosis, which causes Johne's disease in ruminants. This assay is frequently used to determine the presence of the bacterium in feces of infected cattle and sheep. Two IS900 primer sets developed in the 1990s were widely used for decades, and their use has continued in current studies. However, these primers were developed prior to the availability of complete genome sequences. Recent sequence analysis of the binding locations for one primer pair (P90/P91) identified errors and binding inefficiencies that can be easily corrected to further increase detection sensitivity. The P90 primer is missing two nucleotides that should be present near the 3' end, and it does not bind all copies of IS900 due to 5' deletions at some IS900 loci. These IS900 primer pairs, along with newly developed primers, were tested by real-time PCR on purified genomic DNA to determine which primer set performed the best and how primer design errors affect amplification efficiencies. The newly designed PCR primer set (JB5) showed increased sensitivity by two to three quantification cycles using purified genomic DNA and was similar in efficiency to 150C/921. These tests were extended using DNA from feces and tissues of infected cows, which showed similar results. Finally, a 167-bp partial duplication of IS900 was found in type I strains. Although P90 and P91 primers successfully amplify M. avium subsp. paratuberculosis DNA, their use should be discontinued in favor of more efficient primer pairs in future studies. IMPORTANCE This study is an example of how applied genomic analysis can aid diagnostic test improvements. Detection of Mycobacterium avium subsp. paratuberculosis infection of livestock prior to the appearance of clinical disease signs is very difficult but essential for identifying animals shedding the bacterium to prevent transmission of Johne's disease. Total M. avium subsp. paratuberculosis quantity in the feces as determined by real-time PCR (qPCR) using the IS900 target indicates bacterial shedding status and potential for transmission of the pathogen. However, legacy primers designed prior to the availability of complete genome sequences that are used in these tests to detect M. avium subsp. paratuberculosis were based on data from only a single copy of IS900 and not considering all copies collectively as a group. This approach resulted in primer design errors which can be easily corrected to improve test sensitivities. We tested original primers that contain these errors and their corrected versions by qPCR and showed improved sensitivity on purified genomic DNA as well as fecal and tissue samples. These findings may help detect the organism from environmental samples on farms where sensitivity is currently lacking.

Residual enrofloxacin in cattle manure increased persistence and dissemination risk of antibiotic resistance genes during anaerobic digestion

Anaerobic digestion is a common approach to dispose and recycle livestock manures, and the agricultural application of anaerobic digestives represents an important pathway of spreading antibiotic resistance genes (ARGs) from livestock manures to soils. Enrofloxacin is a clinically important fluoroquinolone antibiotic with high residual concentrations in livestock manure, and propagation of fluoroquinolone resistance genes poses a huge risk to public health. Compared with other antibiotics, enrofloxacin is relatively durable in anaerobic digestion system. However, its effect on the persistence of ARGs during anaerobic digestion and its mechanism are not clear. In this study, we investigated effects of 0, 4, and 8 mg/L enrofloxacin on the abundance, persistence, and transferring risk of five plasmid-mediated fluroquinolone ARGs and five typic clinically important non-fluoroquinolone ARGs during cattle manure digestion. The responses of integrons and microbial communities to enrofloxacin were assessed to uncover the underlying mechanisms. All the ten detected ARGs were highly persistent in anaerobic digestion, among them seven ARGs increased over 8.2 times after digestion. Network analysis revealed that the potential hosts of ARGs were critical functional taxa during anaerobic digestion, which can explain the high persistence of ARGs. Residual enrofloxacin significantly increased the abundance of aac(6')-ib-cr, sul1, intI1, and intI2 throughout the digestion, but had no impact on the other ARGs, demonstrating its role in facilitating horizontal gene transfer of the plasmid-mediated aac(6')-ib-cr. The influence of enrofloxacin on microbial communities disappeared at the end of digestion, but the ARG profiles remained distinctive between the enrofloxacin treatments and the control, suggesting the high persistence of enrofloxacin induced ARGs. Our results suggested the high persistence of ARGs in anaerobic digestion system, and highlighted the role of residual enrofloxacin in livestock manure in increasing dissemination risk of fluroquinolone resistance genes.

Strain diversity in Mycobacterium avium subsp. paratuberculosis-positive bovine fecal samples collected in Switzerland

Paratuberculosis or Johne's disease is a chronic intestinal disease in domestic and wild ruminants. It affects global dairy economy and is caused by Mycobacterium avium subsp. paratuberculosis (MAP). The objective of this study was to analyze strain diversity in MAP-positive fecal samples by using a particular single nucleotide polymorphism (SNP) distinguishing between cattle (C-) and sheep (S-) type MAP and analysis of SNPs within gyrA and gyrB genes differentiating between Types I, II, and III. Moreover, mycobacterial interspersed repetitive unit and variable-number tandem repeat (MIRU-VNTR) analysis using eight established loci was performed. A total of 90 fecal samples from diseased animals presenting diarrhea and/or weight loss, originating from 59 bovine herds across 16 cantons of Switzerland were screened by PCR for the MAP-specific F57 and IS900 genes and were further subtyped. 96.7% and 3.3% of the samples contained C- and S-type MAP, respectively. Ten INRA Nouzilly MIRU-VNTR (INMV) profiles, with a discriminatory index of 0.802, calculated based on 65 epidemiological independent genotypes, were detected: INMV 1 (33.8%), INMV 2 (23.1%), INMV 6 (16.9%), INMV 9 (9.2%), INMV 116 (4.6%), INMV 3 (3.1%), INMV 5 (3.1%) and INMV 72 (1.5%), including two novel INMV profiles, namely INMV 253 (3.1%; S-type III) and INMV 252 (1.5%; C-type). INMV 1, INMV 2, and INMV 6 comprised almost 75% of the F57- and IS900-positive samples. Typing data from 11 herds suggest that there are some herds with intra-herd diversity of genotypes. The results of this study indicate a heterogeneity of MAP in Switzerland.

A Comparative Study on the Efficiency of Two Mycobacterium avium subsp. paratuberculosis (MAP)-Derived Lipopeptides of L3P and L5P as Capture Antigens in an In-House Milk ELISA Test

Mycobacterium avium subsp. paratuberculosis (MAP) surface-exposed lipopeptides could be specific capture-antigen molecules targeting antibodies against MAP, in milk, through ELISA. Previous studies have revealed that MAP strains, isolated from sheep (S) or cow (C), could produce specific lipopeptides, L3P or L5P, respectively. In this study, we used L3P and L5P as capture antigens in an in-house milk ELISA (H-MELISA) to assess how these antigens perform, in comparison with other ELISA tests, on well-defined milk samples from MAP-infected sheep. The overall positivity rates of H-MELISA via L3P and L5P varied by the source of milk samples, in which, at bulk-tank-milk (BTM) level, the majority of positive cases (63.83%) reacted more against L5P, whereas a predominant number (69.14%) of milk samples were more responsive against L3P at the individual level. To clarify whether the positivity status of milk samples in H-MELISA L3P/L5P were predictive of MAP strain-types (S/C), strain-typing was carried out using PCR IS1311-restriction enzyme analysis. Although the presence of three MAP strains (S/C/bison types) was detected among the milk samples, the C-type (46.67%) and S-type (75%) MAP strains were detected with higher incidence among BTMs and individual milk samples, respectively. However, further examination on the H-MELISA L3P/L5P-positivity pattern of each C/S-type-MAP sample revealed that some samples had a reverse reactivity against both L3P and L5P. These results could be the consequence of either cross-reactivity between L3P and L5P (due to the similarity in the structures of the two epitopes) or simply a within-herd mixed infection with MAP strains of C and S types. These findings suggest that lipopeptide antigens could contribute a diagnostic test with optimal performance, considering the diversity of MAP strains.

Draft Genome Sequences of Two Bison-Type and Two Sheep-Type Strains of Mycobacterium avium subsp. paratuberculosis

Genome sequences of two type B and two type S strains of Mycobacterium avium subsp. paratuberculosis are presented. These strains were isolated in the United States from sheep, bison, and cattle suffering from Johne’s disease. These genomes will increase our understanding of the minor differences that exist among strains of this genetically stable subspecies.

Molecular characterisation of Mycobacterium avium subsp. paratuberculosis in Australia

BACKGROUND

Mycobacterium avium subsp. paratuberculosis (Map) causes Johne's disease (JD), a chronic enteritis widespread in ruminants, resulting in substantial economic losses, especially to the dairy industry. Understanding the genetic diversity of Map in Australia will assist epidemiological studies for tracking disease transmission and identify subtype characteristics for use in development of improved diagnostic typing methods. Here we investigated the phylogenetic relationships of 351 Map isolates and compared different subtyping methods to assess their suitability for use in diagnostics and accuracy.

RESULTS

SNP-based phylogenetic analysis of 228 Australian isolates and 123 publicly available international isolates grouped Type S and Type C strains into two distinct lineages. Type C strains were highly monomorphic with only 20 SNP differences separating them. Type S strains, when aligned separately to the Telford strain, fell into two distinct clades: The first clade contained seven international isolates while the second clade contained one international isolate from Scotland and all 59 Australian isolates. The Australian Type B strain clustered with US bison strains. IS1311 PCR and Restriction Enzyme Analysis (REA) intermittently generated incorrect results when compared to Long Sequence Polymorphism (LSP) analysis, whole genome SNP-based phylogenetic analysis, IS1311 sequence alignment and average nucleotide identity (ANI). These alternative methods generated consistent Map typing results. A published SNP based assay for genotyping Map was found to be unsuitable for differentiating between Australian and international strain types of Map.

CONCLUSION

This is the first phylogenetic analysis of Australian Map isolates. The Type C lineage was highly monomorphic, and the Type S lineage clustered all Australian isolates into one clade with a single Scottish sheep strain. The Australian isolate classified as Type B by IS1311 PCR and REA is likely to be descended from bison and most closely related to US bison strains. Limitations of the current typing methods were identified in this study.

Complete Genome Sequence of a Type III Ovine Strain of Mycobacterium avium subsp. paratuberculosis

The complete genome sequence of a type III strain of Mycobacterium avium subsp. paratuberculosis was determined. The genome size for this pathogen of sheep is 4,895,755 bp with no plasmid DNA. The chromosome contains 19 copies of the hallmark IS900 element, which is routinely used to identify this subspecies.

The complete genome sequence of a type III strain of Mycobacterium avium subsp. paratuberculosis was determined. The genome size for this pathogen of sheep is 4,895,755 bp with no plasmid DNA. The chromosome contains 19 copies of the hallmark IS900 element, which is routinely used to identify this subspecies.

Comparative Genomics of Mycobacterium avium Subspecies Paratuberculosis Sheep Strains

Mycobacterium avium subspecies paratuberculosis (MAP) is the aetiological agent of Johne's disease (JD), a chronic enteritis that causes major losses to the global livestock industry. Further, it has been associated with human Crohn's disease. Several strains of MAP have been identified, the two major groups being sheep strain MAP, which includes the Type I and Type III sub-lineages, and the cattle strain or Type II MAP lineage, of which bison strains are a sub-grouping. Major genotypic, phenotypic and pathogenic variations have been identified in prior comparisons, but the research has predominately focused on cattle strains of MAP. In countries where the sheep industries are more prevalent, however, such as Australia and New Zealand, ovine JD is a substantial burden. An information gap exists regarding the genomic differences between sheep strain sub-lineages and the relevance of Type I and Type III MAP in terms of epidemiology and/or pathogenicity. We therefore investigated sheep MAP isolates from Australia and New Zealand using whole genome sequencing. For additional context, sheep MAP genome datasets were downloaded from the Sequence Read Archive and GenBank. The final dataset contained 18 Type III and 16 Type I isolates and the K10 cattle strain MAP reference genome. Using a pan-genome approach, an updated global phylogeny for sheep MAP from de novo assemblies was produced. When rooted with the K10 cattle reference strain, two distinct clades representing the lineages were apparent. The Australian and New Zealand isolates formed a distinct sub-clade within the type I lineage, while the European type I isolates formed another less closely related group. Within the type III lineage, isolates appeared more genetically diverse and were from a greater number of continents. Querying of the pan-genome and verification using BLAST analysis revealed lineage-specific variations (n = 13) including genes responsible for metabolism and stress responses. The genetic differences identified may represent important epidemiological and virulence traits specific to sheep MAP. This knowledge will potentially contribute to improved vaccine development and control measures for these strains.

Serological and Molecular Characterization of Mycobacterium avium Subsp. paratuberculosis (MAP) from Sheep, Goats, Cattle and Camels in the Eastern Province, Saudi Arabia

Simple Summary

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne’s disease, affecting small and large ruminants and causing chronic diarrhea and severe emaciation. MAP is prevalent in many countries, including Saudi Arabia. Serological and molecular characterization of MAP and determination of the prevalent strains are essential for the control strategies. The results obtained from 31 herds showed that the sheep type (S-type) was the most prevalent MAP type and the molecular characterization revealed different strain profiles distributed among the sheep, goat, cattle, and camel herds in Eastern Province, Saudi Arabia.

Abstract

The objectives of the present study were to characterize Mycobacterium avium subsp. paratuberculosis (MAP) infection using serological and molecular tools and investigate the distribution and molecular characterization of MAP strains (cattle (C) and sheep (S) types) in sheep, goat, cattle, and camel herds in Eastern Province, Saudi Arabia. Serum and fecal samples were collected from all animals aged >2 years old in 31 herds (sheep = 8, goats = 6, cattle = 8 and camels = 9) from January to December 2019. Serum samples were tested by ELISA for the detection of MAP antibodies. Fecal samples were tested by PCR for the detection of MAP IS900 gene and the identification of MAP strains. MAP antibodies were detected in 19 (61.3%) herds. At the animal level, antibodies against MAP were detected in 43 (19.5%) sheep, 21 (17.1%) goats, 13 (19.7%) cattle and 22 (9.1%) camels. The IS900 gene of MAP was detected in 23 (74.2%) herds and was directly amplified from fecal samples of 59 (26.8%) sheep, 34 (27.6%) goats, 20 (30.3%) cattle and 36 (15.0%) camels. The S-type was the most prevalent MAP type identified in 15 herds, and all were identified as type-I, while the C-type was identified in only 8 herds. The IS900 sequences revealed genetic differences among the MAP isolates recovered from sheep, goats, cattle and camels. Results from the present study show that MAP was prevalent and confirm the distribution of different MAP strains in sheep, goat, cattle and camel herds in Eastern Province, Saudi Arabia.

Complete Genome Sequence of Ovine Mycobacterium avium subsp. paratuberculosis Strain JIII-386 (MAP-S/type III) and Its Comparison to MAP-S/type I, MAP-C, and M. avium Complex Genomes

Mycobacterium avium (M. a.) subsp. paratuberculosis (MAP) is a worldwide-distributed obligate pathogen in ruminants causing Johne’s disease. Due to a lack of complete subtype III genome sequences, there is not yet conclusive information about genetic differences between strains of cattle (MAP-C, type II) and sheep (MAP-S) type, and especially between MAP-S subtypes I, and III. Here we present the complete, circular genome of MAP-S/type III strain JIII-386 (DE) closed by Nanopore-technology and its comparison with MAP-S/type I closed genome of strain Telford (AUS), MAP-S/type III draft genome of strain S397 (U.S.), twelve closed MAP-C strains, and eight closed M.-a.-complex-strains. Structural comparative alignments revealed clearly the mosaic nature of MAP, emphasized differences between the subtypes and the higher diversity of MAP-S genomes. The comparison of various genomic elements including transposases and genomic islands provide new insights in MAP genomics. MAP type specific phenotypic features may be attributed to genes of known large sequence polymorphisms (LSP^Ss) regions I–IV and deletions #1 and #2, confirmed here, but could also result from identified frameshifts or interruptions of various virulence-associated genes (e.g., mbtC in MAP-S). Comprehensive core and pan genome analysis uncovered unique genes (e.g., cytochromes) and genes probably acquired by horizontal gene transfer in different MAP-types and subtypes, but also emphasized the highly conserved and close relationship, and the complex evolution of M.-a.-strains.

Genetic Diversity Among Mycobacterium avium Subspecies Revealed by Analysis of Complete Genome Sequences

Mycobacterium avium comprises four subspecies that contain both human and veterinary pathogens. At the inception of this study, twenty-eight M. avium genomes had been annotated as RefSeq genomes, facilitating direct comparisons. These genomes represent strains from around the world and provided a unique opportunity to examine genome dynamics in this species. Each genome was confirmed to be classified correctly based on SNP genotyping, nucleotide identity and presence/absence of repetitive elements or other typing methods. The Mycobacterium avium subspecies paratuberculosis (Map) genome size and organization was remarkably consistent, averaging 4.8 Mb with a variance of only 29.6 kb among the 13 strains. Comparing recombination events along with the larger genome size and variance observed among Mycobacterium avium subspecies avium (Maa) and Mycobacterium avium subspecies hominissuis (Mah) strains (collectively termed non-Map) suggests horizontal gene transfer occurs in non-Map, but not in Map strains. Overall, M. avium subspecies could be divided into two major sub-divisions, with the Map type II (bovine strains) clustering tightly on one end of a phylogenetic spectrum and Mah strains clustering more loosely together on the other end. The most evolutionarily distinct Map strain was an ovine strain, designated Telford, which had >1,000 SNPs and showed large rearrangements compared to the bovine type II strains. The Telford strain clustered with Maa strains as an intermediate between Map type II and Mah. SNP analysis and genome organization analyses repeatedly demonstrated the conserved nature of Map versus the mosaic nature of non-Map M. avium strains. Finally, core and pangenomes were developed for Map and non-Map strains. A total of 80% Map genes belonged to the Map core genome, while only 40% of non-Map genes belonged to the non-Map core genome. These genomes provide a more complete and detailed comparison of these subspecies strains as well as a blueprint for how genetic diversity originated.

Genotyping methods and molecular epidemiology of Mycobacterium avium subsp. paratuberculosis (MAP)

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne’s disease (JD) which affects mainly ruminants and is characterized by chronic diarrhea and emaciation. Johne’s disease is highly prevalent in many countries around the world and leads to high economic losses associated with decreased production. Genotyping of the involved pathogen could be used in the study of population genetics, pathogenesis and molecular epidemiology including disease surveillance and outbreak investigation. Principally, researchers have first assumed the presence of two different MAP strains that are associated with the animal host species (cattle and sheep). However, nowadays MAP characterization depends mainly upon genetic testing using genetic markers such as insertion elements, repetitive sequences and single nucleotide polymorphisms. This work aims to provide an overview of the advances in molecular biological tools used for MAP typing in the last two decades, discuss how these methods have been used to address interesting epidemiological questions, and explore the future prospects of MAP molecular epidemiology given the ever decreasing costs of the high throughput sequencing technology.

Viable Mycobacterium avium ssp. paratuberculosis isolated from calf milk replacer

When advising farmers on how to control Johne's disease in an infected herd, one of the main recommendations is to avoid feeding waste milk to calves and instead feed calf milk replacer (CMR). This advice is based on the assumption that CMR is free of viable Mycobacterium avium ssp. paratuberculosis (MAP) cells, an assumption that has not previously been challenged. We tested commercial CMR products (n = 83) obtained from dairy farms around the United States by the peptide-mediated magnetic separation (PMS)-phage assay, PMS followed by liquid culture (PMS-culture), and direct IS900 quantitative PCR (qPCR). Conventional microbiological analyses for total mesophilic bacterial counts, coliforms, Salmonella, coagulase-negative staphylococci, streptococci, nonhemolytic Corynebacterium spp., and Bacillus spp. were also performed to assess the overall microbiological quality of the CMR. Twenty-six (31.3%) of the 83 CMR samples showed evidence of the presence of MAP. Seventeen (20.5%) tested positive for viable MAP by the PMS-phage assay, with plaque counts ranging from 6 to 1,212 pfu/50 mL of reconstituted CMR (average 248.5 pfu/50 mL). Twelve (14.5%) CMR samples tested positive for viable MAP by PMS-culture; isolates from all 12 of these samples were subsequently confirmed by whole-genome sequencing to be different cattle strains of MAP. Seven (8.4%) CMR samples tested positive for MAP DNA by IS900 qPCR. Four CMR samples tested positive by both PMS-based tests and 5 CMR samples tested positive by IS900 qPCR plus one or other of the PMS-based tests, but only one CMR sample tested positive by all 3 MAP detection tests applied. All conventional microbiology results were within current standards for whole milk powders. A significant association existed between higher total bacterial counts and presence of viable MAP indicated by either of the PMS-based assays. This represents the first published report of the isolation of viable MAP from CMR. Our findings raise concerns about the potential ability of MAP to survive manufacture of dried milk-based products.

Trends and advances in the diagnosis and control of paratuberculosis in domestic livestock

Paratuberculosis (pTB) is a chronic granulomatous enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP) in a wide variety of domestic and wild animals. Control of pTB is difficult due to the lack of sensitive, efficacious and cost-effective diagnostics and marker vaccines. Microscopy, culture, and PCR have been used for the screening of MAP infection in animals for quite a long time. Besides, giving variable sensitivity and specificity, these tests have not been considered ideal for large-scale screening of domestic livestock. Serological tests like ELISA easily detects anti-MAP antibodies. However, it cannot differentiate between the vaccinated and infected animals. Nanotechnology-based diagnostic tests are underway to improve the sensitivity and specificity. Newer generation diagnostic tests based on recombinant MAP secretory proteins would open new paradigm for the differentiation between infected and vaccinated animals and for early detection of the infection. Due to higher seroreactivity of secretory proteins vis-à-vis cellular proteins, the secretory proteins may be used as marker vaccine, which may aid in the control of pTB infection in animals. Secretory proteins can be potentially used to develop future diagnostics, surveillance and monitoring of the disease progression in animals and the marker vaccine for the control and eradication of pTB.

Phylogenomic exploration of the relationships between strains of Mycobacterium avium subspecies paratuberculosis

Background

Mycobacterium avium subspecies paratuberculosis (Map) is an infectious enteric pathogen that causes Johne’s disease in livestock. Determining genetic diversity is prerequisite to understanding the epidemiology and biology of Map. We performed the first whole genome sequencing (WGS) of 141 global Map isolates that encompass the main molecular strain types currently reported. We investigated the phylogeny of the Map strains, the diversity of the genome and the limitations of commonly used genotyping methods.

Results

Single nucleotide polymorphism (SNP) and phylogenetic analyses confirmed two major lineages concordant with the former Type S and Type C designations. The Type I and Type III strain groups are subtypes of Type S, and Type B strains are a subtype of Type C and not restricted to Bison species.

We found that the genome-wide SNPs detected provided greater resolution between isolates than currently employed genotyping methods. Furthermore, the SNP used for IS1311 typing is not informative, as it is likely to have occurred after Type S and C strains diverged and does not assign all strains to the correct lineage. Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) differentiates Type S from Type C but provides limited resolution between isolates within these lineages and the polymorphisms detected do not necessarily accurately reflect the phylogenetic relationships between strains.

WGS of passaged strains and coalescent analysis of the collection revealed a very high level of genetic stability, with the substitution rate estimated to be less than 0.5 SNPs per genome per year.

Conclusions

This study clarifies the phylogenetic relationships between the previously described Map strain groups, and highlights the limitations of current genotyping techniques. Map isolates exhibit restricted genetic diversity and a substitution rate consistent with a monomorphic pathogen. WGS provides the ultimate level of resolution for differentiation between strains. However, WGS alone will not be sufficient for tracing and tracking Map infections, yet importantly it can provide a phylogenetic context for affirming epidemiological connections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2234-5) contains supplementary material, which is available to authorized users.

Novel Single Nucleotide Polymorphism-Based Assay for Genotyping Mycobacterium avium subsp. paratuberculosis

Typing of Mycobacterium avium subspecies paratuberculosis strains presents a challenge, since they are genetically monomorphic and traditional molecular techniques have limited discriminatory power. The recent advances and availability of whole-genome sequencing have extended possibilities for the characterization of Mycobacterium avium subspecies paratuberculosis, and whole-genome sequencing can provide a phylogenetic context to facilitate global epidemiology studies. In this study, we developed a single nucleotide polymorphism (SNP) assay based on PCR and restriction enzyme digestion or sequencing of the amplified product. The SNP analysis was performed using genome sequence data from 133 Mycobacterium avium subspecies paratuberculosis isolates with different genotypes from 8 different host species and 17 distinct geographic regions around the world. A total of 28,402 SNPs were identified among all of the isolates. The minimum number of SNPs required to distinguish between all of the 133 genomes was 93 and between only the type C isolates was 41. To reduce the number of SNPs and PCRs required, we adopted an approach based on sequential detection of SNPs and a decision tree. By the analysis of 14 SNPs Mycobacterium avium subspecies paratuberculosis isolates can be characterized within 14 phylogenetic groups with a higher discriminatory power than mycobacterial interspersed repetitive unit-variable number tandem repeat assay and other typing methods. Continuous updating of genome sequences is needed in order to better characterize new phylogenetic groups and SNP profiles. The novel SNP assay is a discriminative, simple, reproducible method and requires only basic laboratory equipment for the large-scale global typing of Mycobacterium avium subspecies paratuberculosis isolates.

Genetic diversity of Mycobacterium avium subspecies paratuberculosis and the influence of strain type on infection and pathogenesis: a review

Mycobacterium avium subspecies paratuberculosis (Map) is an important pathogen that causes a chronic, progressive granulomatous enteritis known as Johne's disease or paratuberculosis. The disease is endemic in many parts of the world and responsible for considerable losses to the livestock and associated industries. Diagnosis and control are problematic, due mostly to the long incubation period of the disease when infected animals show no clinical signs and are difficult to detect, and the ability of the organism to survive and persist in the environment. The existence of phenotypically distinct strains of Map has been known since the 1930s but the genetic differentiation of Map strain types has been challenging and only recent technologies have proven sufficiently discriminative for strain comparisons, tracing the sources of infection and epidemiological studies. It is important to understand the differences that exist between Map strains and how they influence both development and transmission of disease. This information is required to develop improved diagnostics and effective vaccines for controlling Johne's disease. Here I review the current classification of Map strain types, the sources of the genetic variability within strains, growth characteristics and epidemiological traits associated with strain type and the influence of strain type on infection and pathogenicity.

SNP genotyping of animal and human derived isolates of Mycobacterium avium subsp. paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) is the aetiological agent of Johne's disease (JD), a chronic granulomatous enteritis that affects ruminants worldwide. While the ability of MAP to cause disease in animals is clear, the role of this bacterium in human inflammatory bowel diseases remains unresolved. Previous whole genome sequencing of MAP isolates derived from human and three animal hosts showed that human isolates were genetically similar and showed a close phylogenetic relationship to one bovine isolate. In contrast, other animal derived isolates were more genetically diverse. The present study aimed to investigate the frequency of this human strain across 52 wild-type MAP isolates, collected predominantly from Australia. A Luminex based SNP genotyping approach was utilised to genotype SNPs that had previously been shown to be specific to the human, bovine or ovine isolate types. Fourteen SNPs were initially evaluated across a reference panel of isolates with known genotypes. A subset of seven SNPs was chosen for analysis within the wild-type collection. Of the seven SNPs, three were found to be unique to paediatric human isolates. No wild-type isolates contain these SNP alleles. Interestingly, and in contrast to the paediatric isolates, three additional adult human isolates (derived from adult Crohn's disease patients) also did not contain these SNP alleles. Furthermore we identified two SNPs, which demonstrate extensive polymorphism within the animal-derived MAP isolates. One of which appears unique to ovine and a single camel isolate. From this study we suggest the existence of genetic heterogeneity between human derived MAP isolates, some of which are highly similar to those derived from bovine hosts, but others of which are more divergent.

Genetic diversity and phylogeny of Mycobacterium avium

Mycobacterium avium, one of the species of the M. avium complex (MAC), includes 4 subspecies, i.e., M. avium subsp. hominissuis (MAH), M. avium subsp. avium (MAA), M. avium subsp. silvaticum (MAS) and M. avium subsp. paratuberculosis (MAP), in turn classified into the S (sheep) and C (cattle) types. These subspecies, although closely related, represent distinct organisms, each endowed with specific pathogenetic and host range characteristics, ranging from environmental opportunistic bacteria that cause infections in swine and immunocompromised patients to pathogens of birds and ruminants. The present review summarizes the basic epidemiological and pathological features of the M. avium subspecies, describes the major genomic events responsible of M. avium subspecies diversity (insertion sequences, sequence variations in specific chromosome loci or genes, deletions, duplications and insertions of large genomic regions) and then reconstructs the phylogenetic relationships among the M. avium subspecies.

How does a Mycobacterium change its spots? Applying molecular tools to track diverse strains of Mycobacterium avium subspecies paratuberculosis

Defining genetic diversity in the wake of the release of several Mycobacterium avium subsp. paratuberculosis (MAP) genome sequences has become a major emphasis in the molecular biology and epidemiology of Johne's disease research. These data can now be used to define the extent of strain diversity on the farm. However, to perform these important tasks, researchers must have a way to distinguish the many MAP isolates/strains that are present in the environment or host to enable tracking over time. Recent studies have described genetic diversity of the Mycobacterium avium complex (MAC), of which MAP is a member, through pulsed-field gel electrophoresis, single sequence repeats, variable-number tandem repeats, genome rearrangements, single nucleotide polymorphisms and genomewide comparisons to identify insertions and deletions. Combinations of these methods can now provide discrimination sufficient for dependable strain tracking. These molecular epidemiology techniques are being applied to understand transmission of Johne's disease within dairy cattle herds as well as identify which strains predominate in wildlife.

Culture phenotypes and molecular characterization of Mycobacterium avium subsp. paratuberculosis isolates from small ruminants

In this study the suitability of different solid media was investigated for the isolation of Mycobacterium avium subsp. paratuberculosis (Map) in order to identify the optimum single or combination of media to permit the isolation of all strain types from small ruminants. A subset of these Map strains was then further characterized by molecular typing methods to assess the genetic diversity of Map strains in the study area (Northern Greece). Map strains were isolated from tissues and faeces of infected goats (n=52) and sheep (n=8) and were analysed for polymorphisms in IS1311 to classify the strain type as Type C or S. The study found that M7H11 supplemented with mycobactin j, OADC and new born calf serum (M7H11+Mj) is the best single choice of medium for the primary isolation of Map of both Type C and S from small ruminants. The combination of M7H11+Mj and Herrolds egg yolk medium supplemented with mycobactin j and sodium pyruvate allowed the detection of all Map isolates in this study. Nineteen Map isolates were characterised by pulsed-field gel electrophoresis and the isolates demonstrated significant genetic diversity. Twelve different SnaBI and 16 distinct SpeI profiles were detected of which 25 have not been described previously and are new profiles. The combination of both enzyme profiles gave 13 different multiplex profiles. Ten different multiplex profiles were detected in goats and three in sheep. One ovine isolate gave the same multiplex profile as a caprine isolate and two different profiles were found within a single goat herd.

Inter- and intra-subtype genotypic differences that differentiate Mycobacterium avium subspecies paratuberculosis strains

Background

Mycobacterium avium subspecies paratuberculosis (Map) is the aetiological agent of Johne’s disease or paratuberculosis and is included within the Mycobacterium avium complex (MAC). Map strains are of two major types often referred to as ‘Sheep’ or ‘S-type’ and ‘Cattle’ or ‘C-type’. With the advent of more discriminatory typing techniques it has been possible to further classify the S-type strains into two groups referred to as Type I and Type III. This study was undertaken to genotype a large panel of S-type small ruminant isolates from different hosts and geographical origins and to compare them with a large panel of well documented C-type isolates to assess the genetic diversity of these strain types. Methods used included Mycobacterial Interspersed Repetitive Units - Variable-Number Tandem Repeat analysis (MIRU-VNTR), analysis of Large Sequence Polymorphisms by PCR (LSP analysis), Single Nucleotide Polymorphism (SNP) analysis of gyr genes, Pulsed-Field Gel Electrophoresis (PFGE) and Restriction Fragment Length Polymorphism analysis coupled with hybridization to IS900 (IS900-RFLP) analysis.

Results

The presence of LSP^A4 and absence of LSP^A20 was confirmed in all 24 Map S-type strains analysed. SNPs within the gyr genes divided the S-type strains into types I and III. Twenty four PFGE multiplex profiles and eleven different IS900-RFLP profiles were identified among the S-type isolates, some of them not previously published. Both PFGE and IS900-RFLP segregated the S-type strains into types I and III and the results concurred with those of the gyr SNP analysis. Nine MIRU-VNTR genotypes were identified in these isolates. MIRU-VNTR analysis differentiated Map strains from other members of Mycobacterium avium Complex, and Map S-type from C-type but not type I from III. Pigmented Map isolates were found of type I or III.

Conclusion

This is the largest panel of S-type strains investigated to date. The S-type strains could be further divided into two subtypes, I and III by some of the typing techniques (IS900-RFLP, PFGE and SNP analysis of the gyr genes). MIRU-VNTR did not divide the strains into the subtypes I and III but did detect genetic differences between isolates within each of the subtypes. Pigmentation is not exclusively associated with type I strains.

A 16 kb naturally occurring genomic deletion including mce and PPE genes in Mycobacterium avium subspecies paratuberculosis isolates from goats with Johne's disease

In this study we characterise the genomic and transcriptomic variability of a natural deletion strain of Mycobacterium avium subspecies paratuberculosis (MAP) prevalent in Spanish Guadarrama goats. Using a pan-genome microarray including MAP and M. avium subspecies hominissuis 104 genomes (MAPAC) we demonstrate the genotype to be MAP Type II with a single deletion of 19 contiguous ORFs (16 kb) including a complete mammalian cell entry (mce7_1) operon and adjacent proline-glutamic acid (PE)/proline-proline-glutamic acid (PPE) genes. A deletion specific PCR test was developed and a subsequent screening identified four goat herds infected with the variant strain. Each was located in central Spain and showed epidemiological links suggestive of transmission between herds. A majority of animals infected with the variant manifested a paucibacillary form of the disease. Comparisons between virulent complete genome compliment strains isolated from multibacillary diseased goats and the MAP variant strain during entry into activated macrophages demonstrated an increased sensitivity in the variant to intracellular killing in human and ovine macrophages. As PPE and mce genes are associated with mycobacterial virulence and pathogenesis we investigated the interplay of these gene sets during cell entry using the MAPAC array. This showed significant differential transcriptome profiles compared to full genome complement MAP controls that included changes in other undeleted mce operons and PE/PPE genes, esx-like signalling operons and stress response/fatty acid metabolism pathways. This strain represents the first report of a MAP Type II genotype with significant natural genomic deletions which remains able to cause disease and is transmissible in goats.

Genome-wide analysis of the emerging infection with Mycobacterium avium subspecies paratuberculosis in the Arabian camels (Camelus dromedarius)

Mycobacterium avium subspecies paratuberculosis (M. ap) is the causative agent of paratuberculosis or Johne's disease (JD) in herbivores with potential involvement in cases of Crohn's disease in humans. JD is spread worldwide and is economically important for both beef and dairy industries. Generally, pathogenic ovine strains (M. ap-S) are mainly found in sheep while bovine strains (M. ap-C) infect other ruminants (e.g. cattle, goat, deer), as well as sheep. In an effort to characterize this emerging infection in dromedary/Arabian camels, we successfully cultured M. ap from several samples collected from infected camels suffering from chronic, intermittent diarrhea suggestive of JD. Gene-based typing of isolates indicated that all isolates belong to sheep lineage of strains of M. ap (M. ap-S), suggesting a putative transmission from infected sheep herds. Screening sheep and goat herds associated with camels identified the circulation of this type in sheep but not goats. The current genome-wide analysis recognizes these camel isolates as a sub-lineage of the sheep strain with a significant number of single nucleotide polymorphisms (SNPs) between sheep and camel isolates (∼1000 SNPs). Such polymorphism could represent geographical differences among isolates or host adaptation of M. ap during camel infection. To our knowledge, this is the first attempt to examine the genomic basis of this emerging infection in camels with implications on the evolution of this important pathogen. The sequenced genomes of M. ap isolates from camels will further assist our efforts to understand JD pathogenesis and the dynamic of disease transmission across animal species.

Exploring the zoonotic potential of Mycobacterium avium subspecies paratuberculosis through comparative genomics

A comparative genomics approach was utilised to compare the genomes of Mycobacterium avium subspecies paratuberculosis (MAP) isolated from early onset paediatric Crohn's disease (CD) patients as well as Johne's diseased animals. Draft genome sequences were produced for MAP isolates derived from four CD patients, one ulcerative colitis (UC) patient, and two non-inflammatory bowel disease (IBD) control individuals using Illumina sequencing, complemented by comparative genome hybridisation (CGH). MAP isolates derived from two bovine and one ovine host were also subjected to whole genome sequencing and CGH. All seven human derived MAP isolates were highly genetically similar and clustered together with one bovine type isolate following phylogenetic analysis. Three other sequenced isolates (including the reference bovine derived isolate K10) were genetically distinct. The human isolates contained two large tandem duplications, the organisations of which were confirmed by PCR. Designated vGI-17 and vGI-18 these duplications spanned 63 and 109 open reading frames, respectively. PCR screening of over 30 additional MAP isolates (3 human derived, 27 animal derived and one environmental isolate) confirmed that vGI-17 and vGI-18 are common across many isolates. Quantitative real-time PCR of vGI-17 demonstrated that the proportion of cells containing the vGI-17 duplication varied between 0.01 to 15% amongst isolates with human isolates containing a higher proportion of vGI-17 compared to most animal isolates. These findings suggest these duplications are transient genomic rearrangements. We hypothesise that the over-representation of vGI-17 in human derived MAP strains may enhance their ability to infect or persist within a human host by increasing genome redundancy and conferring crude regulation of protein expression across biologically important regions.

Multiplexed typing of Mycobacterium avium subsp. paratuberculosis types I, II, and III by Luminex xMAP suspension array

Differentiation among types I, II, and III is the primary step in typing Mycobacterium avium subsp. paratuberculosis. We propose an innovative approach based on detection of gyrase B (gyrB) gene polymorphisms by suspension array technology, with high discriminatory power and high-throughput potential.

Rapid identification and differentiation of Mycobacterium avium subspecies paratuberculosis types by use of real-time PCR and high-resolution melt analysis of the MAP1506 locus

High-resolution melt (HRM) analysis can identify sequence polymorphisms by comparing the melting curves of amplicons generated by real-time PCR amplification. We describe the application of this technique to identify Mycobacterium avium subspecies paratuberculosis types I, II, and III. The HRM approach was based on type-specific nucleotide sequences in MAP1506, a member of the PPE (proline-proline-glutamic acid) gene family.

Molecular characterization of Mycobacterium avium subspecies paratuberculosis Types II and III isolates by a combination of MIRU-VNTR loci

Mycobacterial interspersed repetitive units and variable number tandem repeats typing (MIRU-VNTR) is a useful technique that has been recently applied to characterize members of the Mycobacterium avium complex (MAC). The aim of this study was to examine the genetic variability among a collection of Spanish M. avium subspecies paratuberculosis (M. a. paratuberculosis) isolates with a combination of MIRU-VNTR loci. For this purpose we tested six MIRU-VNTR loci (MIRU-2, MIRU-3, VNTR-25, VNTR-32, VNTR-292 and VNTR-259) in 70 M. a. paratuberculosis isolates of Types II and III that were recovered from 22 Spanish localities during a nine-year period (1998-2007). The combination of five loci (MIRU-2, MIRU-3, VNTR-25, VNTR-32 and VNTR-259) enabled the differentiation of 12 allelic profiles, with a resulting Hunter and Gaston discriminatory index (HGDI) of 0.84. Moreover, we obtained MIRU-VNTR patterns that were unique for each of the M. a. paratuberculosis types analyzed (II and III); other patterns were host-related or restricted to geographic areas. Therefore, this MIRU-VNTR approach could be a useful sub-typing molecular tool in order to get a better sense of the epidemiology of Johne's disease.

Occurrence of Mycobacterium avium subspecies paratuberculosis across host species and European countries with evidence for transmission between wildlife and domestic ruminants

Background

Mycobacterium avium subspecies paratuberculosis (Map) causes an infectious chronic enteritis (paratuberculosis or Johne's disease) principally of ruminants. The epidemiology of Map is poorly understood, particularly with respect to the role of wildlife reservoirs and the controversial issue of zoonotic potential (Crohn's disease). Genotypic discrimination of Map isolates is pivotal to descriptive epidemiology and resolving these issues. This study was undertaken to determine the genetic diversity of Map, enhance our understanding of the host range and distribution and assess the potential for interspecies transmission.

Results

164 Map isolates from seven European countries representing 19 different host species were genotyped by standardized IS900 - restriction fragment length polymorphism (IS900-RFLP), pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphisms (AFLP) and mycobacterial interspersed repeat unit-variable number tandem repeat (MIRU-VNTR) analyses. Six PstI and 17 BstEII IS900-RFLP, 31 multiplex [SnaBI-SpeI] PFGE profiles and 23 MIRU-VNTR profiles were detected. AFLP gave insufficient discrimination of isolates for meaningful genetic analysis. Point estimates for Simpson's index of diversity calculated for the individual typing techniques were in the range of 0.636 to 0.664 but a combination of all three methods increased the discriminating power to 0.879, sufficient for investigating transmission dynamics. Two predominant strain types were detected across Europe with all three typing techniques. Evidence for interspecies transmission between wildlife and domestic ruminants on the same property was demonstrated in four cases, between wildlife species on the same property in two cases and between different species of domestic livestock on one property.

Conclusion

The results of this study showed that it is necessary to use multiple genotyping techniques targeting different sources of genetic variation to obtain the level of discrimination necessary to investigate transmission dynamics and trace the source of Map infections. Furthermore, the combination of genotyping techniques may depend on the geographical location of the population to be tested. Identical genotypes were obtained from Map isolated from different host species co-habiting on the same property strongly suggesting that interspecies transmission occurs. Interspecies transmission of Map between wildlife species and domestic livestock on the same property provides further evidence to support a role for wildlife reservoirs of infection.

Single nucleotide polymorphisms in the IS900 sequence of Mycobacterium avium subsp. paratuberculosis are strain type specific

Insertion sequence IS900 is used as a target for the identification of Mycobacterium avium subsp. paratuberculosis. Previous reports have revealed single nucleotide polymorphisms within IS900. This study, which analyzed the IS900 sequences of a panel of isolates representing M. avium subsp. paratuberculosis strain types I, II, and III, revealed conserved type-specific polymorphisms that could be utilized as a tool for diagnostic and epidemiological purposes.

Discovery of stable and variable differences in the Mycobacterium avium subsp. paratuberculosis type I, II, and III genomes by pan-genome microarray analysis

Mycobacterium avium subsp. paratuberculosis is an important animal pathogen widely disseminated in the environment that has also been associated with Crohn's disease in humans. Three M. avium subsp. paratuberculosis genomotypes are recognized, but genomic differences have not been fully described. To further investigate these potential differences, a 60-mer oligonucleotide microarray (designated the MAPAC array), based on the combined genomes of M. avium subsp. paratuberculosis (strain K-10) and Mycobacterium avium subsp. hominissuis (strain 104), was designed and validated. By use of a test panel of defined M. avium subsp. paratuberculosis strains, the MAPAC array was able to identify a set of large sequence polymorphisms (LSPs) diagnostic for each of the three major M. avium subsp. paratuberculosis types. M. avium subsp. paratuberculosis type II strains contained a smaller genomic complement than M. avium subsp. paratuberculosis type I and M. avium subsp. paratuberculosis type III genomotypes, which included a set of genomic regions also found in M. avium subsp. hominissuis 104. Specific PCRs for genes within LSPs that differentiated M. avium subsp. paratuberculosis types were devised and shown to accurately screen a panel (n = 78) of M. avium subsp. paratuberculosis strains. Analysis of insertion/deletion region INDEL12 showed deletion events causing a reduction in the complement of mycobacterial cell entry genes in M. avium subsp. paratuberculosis type II strains and significantly altering the coding of a major immunologic protein (MPT64) associated with persistence and granuloma formation. Analysis of MAPAC data also identified signal variations in several genomic regions, termed variable genomic islands (vGIs), suggestive of transient duplication/deletion events. vGIs contained significantly low GC% and were immediately flanked by insertion sequences, integrases, or short inverted repeat sequences. Quantitative PCR demonstrated that variation in vGI signals could be associated with colony growth rate and morphology.

Sequence polymorphisms in a surface PPE protein distinguish types I, II, and III of Mycobacterium avium subsp. paratuberculosis

In the last 2 decades, a variety of different molecular typing methods have been developed to differentiate strains of Mycobacterium avium subsp. paratuberculosis. The most successful techniques are based on insertion sequences, repetitive loci, comparative genomics, or single nucleotide polymorphisms. In the present study, we chose to examine whether a single M. avium subsp. paratuberculosis gene could serve as a means of differentiation of a variety of isolates. The MAP1506 gene locus encodes a member of the polymorphic PPE protein family that has putative roles relevant to M. avium subsp. paratuberculosis pathogenicity. The MAP1506 locus was sequenced from a collection of 58 M. avium subsp. paratuberculosis isolates from different sources, hosts, and typing profiles. Following sequence alignment and analysis, it was found that bovine (type II) strains of M. avium subsp. paratuberculosis consistently differed from ovine (type I) and intermediate (type III) strains in seven and eight nucleotides, respectively. Polymorphic regions of the MAP1506 locus were selected for analysis by denaturing gradient gel electrophoresis, allowing visual discrimination of the three subtypes of M. avium subsp. paratuberculosis isolates. This is the first report describing the use of PCR and denaturing gradient gel electrophoresis on a single gene as a method to distinguish types I, II, and III of M. avium subsp. paratuberculosis.