Genetic structure of Mycobacterium avium subsp. paratuberculosis population in cattle herds in Quebec as revealed by using a combination of multilocus genomic analyses.

Mycobacterium avium subsp. paratuberculosis is the etiological agent of paratuberculosis, a granulomatous enteritis affecting a wide range of domestic and wild ruminants worldwide. A variety of molecular typing tools are used to distinguish M. avium subsp. paratuberculosis strains, contributing to a better understanding of M. avium subsp. paratuberculosis epidemiology. In the present study, PCR-based typing methods, including mycobacterial interspersed repetitive units/variable-number tandem repeats (MIRU-VNTR) and small sequence repeats (SSR) in addition to IS1311 PCR-restriction enzyme analysis (PCR-REA), were used to investigate the genetic heterogeneity of 200 M. avium subsp. paratuberculosis strains from dairy herds located in the province of Quebec, Canada. The majority of strains were of the "cattle type," or type II, although 3 strains were of the "bison type." A total of 38 genotypes, including a novel one, were identified using a combination of 17 genetic markers, which generated a Simpson's index of genetic diversity of 0.876. Additional analyses revealed no differences in genetic diversity between environmental and individual strains. Of note, a spatial and spatiotemporal cluster was evidenced regarding the distribution of one of the most common genotypes. The population had an overall homogeneous genetic structure, although a few strains stemmed out of the consensus cluster, including the bison-type strains. The genetic structure of M. avium subsp. paratuberculosis populations within most herds suggested intraherd dissemination and microevolution, although evidence of interherd contamination was also revealed. The level of genetic diversity obtained by combining MIRU-VNTR and SSR markers shows a promising avenue for molecular epidemiology investigations of M. avium subsp. paratuberculosis transmission patterns.