Stability of genotyping target sequences of Mycobacterium avium subsp. paratuberculosis upon cultivation on different media, in vitro- and in vivo passage, and natural infection

Genomic epidemiology of Mycobacterium avium subsp. paratuberculosis isolates from Canadian dairy herds provides evidence for multiple infection events

Mycobacterium avium subsp. paratuberculosis (MAP) is the pathogen responsible for paratuberculosis or Johne's Disease (JD) in ruminants, which is responsible for substantial economic losses worldwide. MAP transmission primarily occurs through the fecal-oral route, and the introduction of an MAP infected animal into a herd is an important transmission route. In the current study, we characterized MAP isolates from 67 cows identified in 20 herds from the provinces of Quebec and Ontario, Canada. Whole genome sequencing (WGS) was performed and an average genome coverage (relative to K-10) of ∼14.9 fold was achieved. The total number of SNPs present in each isolate varied from 51 to 132 and differed significantly between herds. Isolates with the highest genetic variability were generally present in herds from Quebec. The isolates were broadly separated into two main clades and this distinction was not influenced by the province from which they originated. Analysis of 8 MIRU-VNTR loci and 11 SSR loci was performed on the 67 isolates from the 20 dairy herds and publicly available references, notably major genetic lineages and six isolates from the province of Newfoundland and Labrador. All 67 field isolates were phylogenetically classified as Type II (C-type) and according to MIRU-VNTR, the predominant type was INMV 2 (76.1%) among four distinct patterns. Multilocus SSR typing identified 49 distinct INMV SSR patterns. The discriminatory index of the multilocus SSR typing was 0.9846, which was much higher than MIRU-VNTR typing (0.3740). Although multilocus SSR analysis provides good discriminatory power, the resolution was not informative enough to determine inter-herd transmission. In select cases, SNP-based analysis was the only approach able to document disease transmission between herds, further validated by animal movement data. The presence of SNPs in several virulence genes, notably for PE, PPE, mce and mmpL, is expected to explain differential antigenic or pathogenetic host responses. SNP-based studies will provide insight into how MAP genetic variation may impact host-pathogen interactions. Our study highlights the informative power of WGS which is now recommended for epidemiological studies and to document mixed genotypes infections.

Phylogenetic and Multiple-Locus Variable number tandem repeat analysis of Mycobacterium avium subsp. paratuberculosis isolates from Argentina

Paratuberculosis is a worldwide chronic enteric disease of ruminants, caused by Mycobacterium avium subsp. paratuberculosis (MAP). While MAP has been widely investigated all around the world, little is known about the different strains that circulate in each country. This study describes the genetic diversity of MAP isolates from different bovine and deer herds from Argentina, analyzed by Multiple-Locus Variable number tandem repeat Analysis (MLVA), as well as the phylogenetic relatedness between geographically distant isolates through Whole Genome Sequencing (WGS) and core-genome analysis. A total of 90 MAP isolates were analyzed. The results showed seven different MLVA genotypes, with almost 75% of them belonging to pattern INMV 1, described in all the herds studied. WGS results suggested the presence of a common INMV 1 strain circulating throughout the country. Our results allow confirming the coexistence of different strains in time and space and the mixed infections identified in some animals. These observations suggest the absence of animal monitoring prior to introduction to the herds and the need for a control program in the country. This study represents the first to report WGS of MAP strains in Argentina.

Transmission patterns of a Mycobacterium avium subsp. paratuberculosis clone within a single heard investigated by Whole Genome Sequencing

Mycobacterium avium subsp. paratuberculosis (MAP) is characterized by a low genomic rate of mutation. Current subtyping tools, such as Mini-Micro-satellite analyses, do to have not sufficient discriminatory power to disclose MAP's evolution on small spatial and temporal scales. The aim of the study was to investigate the population structure of MAP inside a single dairy herd using whole genome sequencing (WGS) approaches. For this purpose, the genomes of 43 field isolates, recovered from the faeces of 36 cows of the same dairy herd from 2012 to 2016, were sequenced by WGS. The isolates' genomes showed a low number (43) of polymorphic sites (SNPs), confirming the clonal origin of the herd infection. However, despite the limited genomic diversity found in WGS, the phylogenetic analysis was discriminatory enough to detect the presence of different genomic clades and sub-clades inside the herd population. In addition, the phylodynamic reconstruction showed the existence of an ancestor clade from which the other clades and sub-clades originated. Moreover, by reconstructing the putative within-herd transmission networks using WGS data, we demonstrated that: (i) in a herd where MAP is endemic, multiple isolates recovered from a single animal and differing from each other by few (three/four) SNPs can originate from different transmission or passive shedding events and not from intra-host evolution; and (ii) variability of minisatellites coupled with a few microsatellites does not represent reliable tracers of within-herd infection chains. Our findings show that WGS, coupled with relevant epidemiological information, represents a valuable tool to work out fine epidemiological and micro-evolutionary relationships such as those at herd-level scale.

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.

Genetic diversity of bovine Mycobacterium avium subsp. paratuberculosis discriminated by IS1311 PCR-REA, MIRU-VNTR, and MLSSR genotyping

The aim of this study was to describe the genetic diversity of Mycobacterium avium subsp. paratuberculosis (MAP) obtained from individual cows in Korea. Twelve MAP-positive fecal DNA samples and 19 MAP isolates were obtained from 10 cattle herds located in 5 provinces in Korea. In addition, 5 MAP isolates obtained from the Czech Republic and Slovakia and 3 isolates from Australia were genotyped for comparison with the domestic isolates. The most prevalent strains in Korea were of the “bison-type” genotype (23 of 31 fecal DNA/isolates) and were distributed nationwide. The remaining MAP isolates (8) and all of the foreign isolates were identified as “cattle-type”. The bison-type strains which were discriminated only as INMV 68 in variable-number tandem repeats of mycobacterial interspersed repetitive units (MIRU-VNTR) typing. Multilocus short sequence repeat (MLSSR) typing differentiated the bison-type strains into 3 different subtypes. The cattle-type strains were divided into 3 subtypes by MIRU-VNTR and 8 subtypes by MLSSR. The allelic diversities in the MIRU-VNTR and MLSSR results were calculated as 0.567 and 0.866, respectively. These results suggest that MIRU-VNTR typing cannot provide a sufficient description of the epidemiological situation of MAP. Therefore, an alternative method, such as MLSSR, is needed for typing of MAP strains to elucidate the molecular epidemiology of MAP infections. Overall, this study is the first epidemiological survey report in Korea using both MIRU-VNTR and MLSSR typing methods, and it has provided basic data necessary to elucidate the characteristics of MAP infections in Korea.

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.

Mycobacterium aviumsubsp. paratuberculosis isolated from wild red deer (Cervus elaphus) in Northern Italy

Paratuberculosis (or Johne's disease) is an infectious disease which affects mainly ruminants and it is caused by Mycobacterium avium subsp. paratuberculosis (MAP). During a culling program (years 2011-2015) aimed at controlling the red deer (Cervus elaphus) population in Stelvio National Park (Italian Alps), where paratuberculosis was already described in this species, 382 tissue samples from the Lombardy Region and 102 fecal specimens from the Autonomous Province of Bolzano were analyzed by PCR. Of these, 77 samples (20.16%) from the Lombardy area and 19 specimens (18.63%) from the Bolzano area resulted PCR positive. The cultural test was carried out on PCR positive samples (n = 96), enabling the isolation of 19 MAP field strains which were genotyped using MIRU-VNTR typing and Short Sequence repeats (SSRs). Our results suggest that all isolates share an identical VNTR profile corresponding to the INMV1 genotype. The only variation was on the locus SSR2, but the utility of this last locus has already been questioned because of its instability. Overall, these data suggest a common clonal origin and host adaptation during the diffusion of paratuberculosis in this population. Finally, this profile is the same as that which has already been described in the cattle population in Northern Italy, suggesting a possible inter-species disease transmission pattern from wildlife to domestic ruminants and vice versa.

Knowledge gaps that hamper prevention and control of Mycobacterium avium subspecies paratuberculosis infection

In the last decades, many regional and country-wide control programmes for Johne's disease (JD) were developed due to associated economic losses, or because of a possible association with Crohn's disease. These control programmes were often not successful, partly because management protocols were not followed, including the introduction of infected replacement cattle, because tests to identify infected animals were unreliable, and uptake by farmers was not high enough because of a perceived low return on investment. In the absence of a cure or effective commercial vaccines, control of JD is currently primarily based on herd management strategies to avoid infection of cattle and restrict within-farm and farm-to-farm transmission. Although JD control programmes have been implemented in most developed countries, lessons learned from JD prevention and control programmes are underreported. Also, JD control programmes are typically evaluated in a limited number of herds and the duration of the study is less than 5 year, making it difficult to adequately assess the efficacy of control programmes. In this manuscript, we identify the most important gaps in knowledge hampering JD prevention and control programmes, including vaccination and diagnostics. Secondly, we discuss directions that research should take to address those knowledge gaps.

Evaluation of associations between genotypes of Mycobacterium avium subsp. paratuberculsis and presence of intestinal lesions characteristic of paratuberculosis

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of paratuberculosis affecting ruminants worldwide. Depending on the MAP-Type (MAP-C or MAP-S, cattle or sheep type), strains differ in virulence and host preference. There is not yet any strong evidence indicating that individual field strains of the same MAP-subgroup exhibit differences in virulence. The aim of this study was to evaluate a potential association between the genotype of individual field strains belonging to the MAP-C group and the presence of macroscopic intestinal lesions characteristic of paratuberculosis in the infected animals. 88 MAP-C isolates were sampled from clinically healthy cows at slaughter. Cows were grouped as A (n=46) with, and B (n=42) without macroscopic intestinal lesions. Sampled cows from both the A and B groups came from different farms and had a similar age distribution. MAP isolates were characterized by MIRU-VNTR and IS900-RFLP analysis. Resulting genotypes were examined for an association with the presence of macroscopic intestinal lesions characteristic of paratuberculosis. MAP isolates from groups A and B exhibited similar strain diversity: 20 and 18 combined genotypes, altogether 32 genotypes. Six of these genotypes were detected in both groups. Although no association was found between individual combined genotypes and presence of macroscopic intestinal lesions, IS900-RFLP-(BstEII)-Type-C1 (the most common type worldwide) was found more often in group A (p<0.01). The data give only weak indication for the existence of differences in virulence among MAP-cattle type isolates. Differences in the development and severity of lesions may rather depend on unknown host factors or inoculation dose. Virulence properties of IS900-RFLP-(BstEII)-Type-C1 isolates should be examined in more detail.

Examination of Mycobacterium avium subspecies paratuberculosis mixed genotype infections in dairy animals using a whole genome sequencing approach

Many pathogenic mycobacteria are known to cause severe disease in humans and animals. M. avium subspecies paratuberculosis (Map) is the causative agent of Johne’s disease—a chronic wasting disease affecting ruminants such as cattle and sheep, responsible for significant economic losses in the dairy and beef industries. Due to the lack of treatment options or effective vaccines, mitigating losses can be difficult. In addition, the early stages of Map infection may occur in asymptomatic hosts that continue to shed viable bacteria in their faeces, leading to the infection of other healthy animals. Using multi-locus short sequence repeat (ML-SSR) analysis we previously reported that individual Johne’s positive dairy cattle from farms across the island of Newfoundland were infected by Map with multiple SSR-types simultaneously. The occurrence of multiple mixed genotype infections has the potential to change pathogen and disease dynamics as well as reduce the efficacy of treatments and vaccines. Therefore, we conducted whole genome sequencing (WGS) and single nucleotide polymorphism (SNP) analysis on a subset of these isolates for a more in-depth examination. We also implemented a PCR assay using two discriminatory SNPs and demonstrated the incidence of a mixed infection by three genotypically diverse Map isolates in a single animal. In addition, results show that WGS and SNP analysis can provide a better understanding of the relationship between Map isolates from individual and different animals. In the future such studies on the occurrence of mixed genotype infections could potentially lead to the identification of variable pathogenicity of different genotypes and allow for better tracking of Map isolates for epidemiological studies.

New polymorphisms within the variable number tandem repeat (VNTR) 7 locus of Mycobacterium avium subsp. paratuberculosis

Variable number tandem repeat (VNTR) is a frequently employed typing method of Mycobacterium avium paratuberculosis (MAP) isolates. Based on whole genome sequencing in a previous study, allelic diversity at some VNTR loci seems to over- or under-estimate the actual phylogenetic variance among isolates. Interestingly, two closely related isolates on one farm showed polymorphism at the VNTR 7 locus, raising concerns about the misleading role that it might play in genotyping. We aimed to investigate the underlying basis of VNTR 7-polymorphism by analyzing sequence data for published genomes and field isolates of MAP and other M. avium complex (MAC) members. In contrast to MAP strains from cattle, strains from sheep displayed an "imperfect" repeat within VNTR 7, which was identical to respective allele types in other MAC genomes. Subspecies- and strain-specific single nucleotide polymorphisms (SNPs) and two novel (16 and 56 bp) repeats were detected. Given the combination of the three existing repeats, there are at least five different patterns for VNTR 7. The present findings highlight a higher polymorphism and probable instability of VNTR 7 locus that needs to be considered and challenged in future studies. Until then, sequencing of this locus in future studies is important to correctly assign the underlying allele types.(1).

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.

Typing of Mycobacterium avium subspecies paratuberculosis isolates from Newfoundland using fragment analysis

Short Sequence Repeat (SSR) typing of Mycobacterium avium subspecies paratuberculosis (Map) isolates is one of the most commonly used method for genotyping this pathogen. Currently used techniques have challenges in analyzing mononucleotide repeats >15 bp, which include some of the Map SSRs. Fragment analysis is a relatively simple technique, which can accurately measure the size of DNA fragments and can be used to calculate the repeat length of the target SSR loci. In the present study, fragment analysis was used to analyze 4 Map SSR loci known to provide sufficient discriminatory power to determine the relationship between Map isolates. Eighty-five Map isolates from 18 animals from the island of Newfoundland were successfully genotyped using fragment analysis. To the best of our knowledge, this is the first report on Map SSR diversity from Newfoundland dairy farms. Previously unreported Map SSR-types or combinations were also identified during the course of the described work. In addition, multiple Map SSR-types were isolated from a single animal in many cases, which is not a common finding.

Ecology and genomic features of infection with Mycobacterium avium subspecies paratuberculosis in Egypt

Mycobacterium avium subspecies paratuberculosis (M. paratuberculosis) is the causative agent of paratuberculosis, or Johne's disease, in cattle, with potential involvement in cases of Crohn's disease in humans. Johne's disease is found worldwide and is economically important for both beef and dairy industries. In an effort to characterize this important infection in Egypt, we analysed the ecological and genomic features of recent isolates of M. paratuberculosis. In this report, we examined 26 Holstein dairy herds distributed throughout Egypt, from 2010 to 2013. Using PCR analysis of faecal samples, we estimated a mean herd-level prevalence of 65.4 %, with animal-level infection that reached a mean of 13.6 % among animals suffering from diarrhoea. Whole genome sequencing of field isolates identified numerous single nucleotide polymorphisms among field isolates relative to the standard M. paratuberculosis K10 genome. Interestingly, the virulence of M. paratuberculosis isolates from Egypt revealed diverse virulence phenotypes in the murine model of paratuberculosis, with significant differences in tissue colonization, particularly during the chronic stage of infection. Overall, our analysis confirmed that Johne's disease is a newly identified problem in Egypt and indicated that M. paratuberculosis has potentially diverse genotypes that impact its virulence. Further ecological mapping and genomic analysis of M. paratuberculosis will enhance our understanding of the transmission and evolutionary dynamics of this pathogen under natural field conditions.

Variable-number tandem repeats genotyping used to aid and inform management strategies for a bovine Johne’s disease incursion in tropical and subtropical Australia

The application of variable-number tandem repeats (VNTR) genotyping of Mycobacterium avium subsp. paratuberculosis isolates to assist in investigating incidents of bovine Johne’s disease in a low-prevalence region of Australia is described in the current study. Isolates from a response to detection of bovine Johne’s disease in Queensland were compared with strains from national and international sources. The tandem application of mycobacterial interspersed repetitive unit (MIRU) and multilocus short sequence repeats (MLSSR) genotyping identified 2 strains, 1 that infected cattle on multiple properties with trace-forward histories from a common infected property, and 1 genotypically different strain recovered from a single property. The former strain showed an identical genotype to an isolate from India. Neither strain showed a genotypic link to regions of Australia with a higher prevalence of the disease. Genotyping has indicated incursions from 2 independent sources. This intelligence has informed investigations into potential routes of entry and the soundness of ongoing control measures, and supported strategy and policy decisions regarding management of Mycobacterium avium subsp. paratuberculosis incursions for Queensland.

Facts, myths and hypotheses on the zoonotic nature of Mycobacterium avium subspecies paratuberculosis

Mycobacterium avium subspecies paratuberculosis (MAP) is the causative agent of paratuberculosis (Johne's disease [JD]), a chronic granulomatous enteritis in ruminants. JD is one of the most widespread bacterial diseases of domestic animals with significant economic impact. The histopathological picture of JD resembles that of Crohn's disease (CD), a human chronic inflammatory bowel disease of still unresolved aetiology. An aetiological relevance of MAP for CD has been proposed. This and the ambiguity of other published epidemiological findings raise the question whether MAP represents a zoonotic agent. In this review, we will discuss evidence that MAP has zoonotic capacity.

Improved short-sequence-repeat genotyping of Mycobacterium avium subsp. paratuberculosis by using matrix-assisted laser desorption ionization-time of flight mass spectrometry

Accurate sequence analysis of mononucleotide repeat regions is difficult, complicating the use of short sequence repeats (SSRs) as a tool for bacterial strain discrimination. Although multiple SSR loci in the genome of Mycobacterium avium subsp. paratuberculosis allow genotyping of M. avium subsp. paratuberculosis isolates with high discriminatory power, further characterization of the most discriminatory loci is limited due to inherent difficulties in sequencing mononucleotide repeats. Here, a method was evaluated using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) as an alternative to Sanger sequencing to further differentiate the dominant mycobacterial interspersed repetitive-unit (MIRU)-variable-number tandem-repeat (VNTR) M. avium subsp. paratuberculosis type (n = 37) in Canadian dairy herds by targeting a highly discriminatory mononucleotide SSR locus. First, PCR-amplified DNA was digested with two restriction enzymes to yield a sufficiently small fragment containing the SSR locus. Second, MALDI-TOF MS was performed to identify the mass, and thus repeat length, of the target. Sufficiently intense, discriminating spectra were obtained to determine repeat lengths up to 15, an improvement over the limit of 11 using traditional sequencing techniques. Comparison to synthetic oligonucleotides and Sanger sequencing results confirmed a valid and reproducible assay that increased discrimination of the dominant M. avium subsp. paratuberculosis MIRU-VNTR type. Thus, MALDI-TOF MS was a reliable, fast, and automatable technique to accurately resolve M. avium subsp. paratuberculosis genotypes based on SSRs.