Relative frequency of 4 major strain types of Mycobacterium avium ssp. paratuberculosis in Canadian dairy herds using a novel single nucleotide polymorphism-based polymerase chain reaction

Johne's disease is a worldwide concern, as it causes huge economic losses. The etiological agent, Mycobacterium avium ssp. paratuberculosis (MAP), has limited genetic diversity, impeding efforts to understand transmission and distribution of strain types. Whole-genome sequencing was previously performed on a representative set of MAP isolates from Canadian dairy herds and 9 divergent clades were identified. Four clades were of particular interest, as they were either MAP types rarely reported in North America, or they represented a substantial proportion of isolates recovered from dairy farms in Canada. One clade included type I/III isolates, whereas the remaining clades included type II isolates. Variant sites in the MAP genome are often separated by thousands of base pairs, limiting use of single nucleotide polymorphism (SNP)-based genotyping on a single genomic region. Therefore, a SNP-PCR assay was developed to facilitate interrogation of 5 SNP in 2 distant regions of the genome, linking them together in a single PCR reaction for subsequent Sanger sequencing. This high-throughput assay enabled discrimination of 602 MAP isolates from 264 herds (from all 10 provinces). More than 1 isolate was cultured from 133 herds, 14 of which included multiple subtypes. A previously identified dominant type included 87% of isolates, whereas the Bison type was more widespread than previously reported. The latter type and isolates from a second clade of interest were overrepresented in Québec and Saskatchewan, respectively. In conclusion, the distribution and relative frequency of MAP subtypes within Canadian dairy herds were assessed using a novel SNP-based typing assay. These findings will contribute to understanding the clinical relevance and transmission dynamics of MAP in this population and elsewhere.