Isolation of Mycobacterium avium subspecies paratuberculosis from Ugandan cattle and strain differentiation using optimised DNA typing techniques

Mycobacterium avium subspecies paratuberculosis (MAP) is the causative agent of paratuberculosis or Johne’s disease; a chronic intractable enteritis which affects many species of animals including cattle, goats, sheep and other domestic and wild ruminants [1]. Paratuberculosis poses a serious economic challenge where it occurs, especially in dairy cattle [2]. Some reports indicate that MAP may be involved in the causation of Crohn’s disease (CD), a human chronic enteropathy whose lesions resemble those of paratuberculosis [3]. It has also been reported that MAP from CD patients and animal species have similar genetic patterns [4]. The seroprevalence of paratuberculosis in Ugandan cattle has been estimated to be 8.8% [5]. However, the host range of MAP in Uganda is still unknown and so is the diversity of infecting strains. Furthermore, there are no reports on the diversity of MAP from any African country at the moment.

Molecular characterisation of MAP was initially based on Restriction Fragment Length Polymorphism using hybridization to MAP specific insertion sequence - IS900 (RFLP-IS900); and Pulsed Field Gel Electrophoresis (PFGE) [6, 7]. These methods are cumbersome, costly and also require large amounts of purified DNA which can only be obtained from cultures.

The first PCR based characterisation method based on restriction enzyme analysis of insertion sequence 1311 (IS1311 PCR-REA) and Heat shock protein 65 gene, used limited amounts of DNA but could not distinguish between most isolates [8, 9]. IS1311 PCR-REA is unable to distinguish between most strains except the simple classification of MAP as cattle (C), sheep (S) and bison (B) types [10, 11]. Further advances in molecular typing of mycobacteria introduced the use of mycobacterial interspersed repeat units (MIRU) [12, 13], variable number tandem repeats (VNTR) [13, 14] and short sequence repeats (SSR) [4, 15] that have enabled the discrimination between two isolates from the same animal [13]. An additional advantage of these systems is that they require minimal quantity of DNA. Amonsin et al. [16] identified 11 SSR loci, three of which have been used in several studies for genotyping [4, 15].

El-Sayed et al. [15] evaluated nine MIRU loci, three SSR loci and six VNTR loci, for characterisation of MAP using 34 isolates of German origin. They recommended the combined use of MIRU locus 2 and SSR loci 1, 2 and 8 (G1, G2 and GGT repeats), although MIRU locus 3 was also discriminatory. They found VNTR loci to be less discriminatory and not very suitable for MAP characterisation. However, other studies that used different sets of VNTR loci indicated that loci 25, 32 and 292 could have some useful discriminatory ability [4, 13, 14]. There is however conflicting information between Thibault et al. [14] and Castellanos et al. [13] as to whether VNTR 32 or VNTR 25 has a higher Hunter-Gaston index of diversity. Thibault et al. [14] studied 183 isolates from 10 different countries (USA, Venezuela, Argentina, France, Italy, Czech Republic, New Zealand, UK, Slovenia, and Netherlands), while Castellanos et al. [13] used 70 isolates from only Spain. The differences observed might have been due to geographical differences or genetic makeup of the cattle in the different study populations.

In this paper we report the existence of a wide range of MAP genotypes from cattle in Uganda as determined using MIRU-VNTR and SSR typing, in combination with IS1311 PCR-REA. This is the first report on characterisation of MAP from cattle in Uganda and from any African country. We also present some modifications to improve the current PCR procedures used in the amplification of some loci.

Twenty four isolates were obtained from culture of tissue and faecal materials. Twenty one of them were positive for IS900 PCR and IS1311 PCR yielding 229 bp and 608 bp respectively. None of the 21 isolates which were positive for IS900 was positive for IS901 PCR ruling out their classification as M. avium avium and M. avium silvaticum. Digestion of IS1311 PCR products with HinfI showed that seventeen isolates were of cattle type, two were of the bison type whereas two others had a novel restriction pattern and were classified as belonging to an unknown type designated here as X (Figure 1, Table 1). The cattle type has resulted in fragments of 323, 285, 218, 67 bases while the bison type has resulted in fragments of 323 and 285 bases. On the other hand type X has fragments of approximately 225, 175, 120, 90 base pairs estimated from the ladder.

Three MIRU 2 allele variants were found among the isolates (Table 1): the first type had seven repeats of approximately 300 bp and was found in one isolate; the second type had nine repeats of approximately 350 bp and was observed in 14 isolates and the third type carried 15 repeats of approximately 520 bp observed in five isolates (Figure 2). Two bands were observed in isolate MapUg3. Therefore, this isolate was considered as one of those isolates in which MIRU 2 amplification was not successful. Based on MIRU 3 profiles, two allele variants were found: one repeat (198 bp) in two isolates and five repeats (300 bp) in 18 isolates. Only 18 isolates were successfully amplified with VNTR 32 and all belonged to the same allele type with nine repeats (300 bp).

Regarding the SSR typing, SSR locus 1 amplification was successful in only 15 isolates with a 450 bp product. Sequencing of the PCR products resulted in two allele variants: 12 isolates with six G1 repeats and three isolates with seven G1 repeats. SSR locus 2 was successfully amplified from all the 21 isolates giving a 442 bp product. Sequencing of the PCR products revealed three allele variants: one with 10 G2 repeats in 14 isolates, a second with 11 G2 repeats in six isolates and a third one with 12 G2 repeat in one isolate. SSR locus 8 PCR resulted in approximately 355 bp products. With subsequent sequence analysis, this was the most informative of all loci for genotyping. Three allele profiles were identified bearing trinucleotide repeats: 2GGT in one isolate, 4GGT in 11 isolates and 3GGT with mutations, in nine isolates (GGTtgtGGT and GGTcgtcgt) (Table 1). The phylogenetic analysis of SSR locus 8 partial sequences was even more informative and separated the isolates into 8 genotypes based on their nucleotide polymorphisms (Figure 3). This separation of isolates was in general agreement to that revealed by the combination of MIRU, SSR and IS1311 PCR-REA typing results. However, the above mentioned combination of loci was able to differentiate into three genotypes the isolates MapUg2, MapUg6, MapUg8 and MapUg9, all of which were grouped together as a single genotype by the phylogenetic analysis of SSR locus 8 (Figure 3). Similarly, MapUg7 was also discriminated from isolates MapUg11-Ug20. As a result the combination of all analyses was more discriminative than the SSR locus 8 SNPs analysis alone.

In general, the combination of IS1311, MIRU 2, 3; VNTR 32; and SSR loci 1, 2 and 8, resolved the 21 MAP isolates into 11 distinct genotypes which were distributed in the five districts from which the samples were obtained. The details of the different isolates and their genotype profiles are shown in Table 1.

The best amplification reaction for MIRU, VNTR and SSR was obtained using Ethylene glycol at a concentration of 1.075 M (Figure 4). DMSO was the second best denaturant (at a concentration of 7.5%) followed by Propan-1, 2-diol (0.816 M), then Betaine, which was also better than a PCR reaction without any denaturant. The thermal profile was improved by increasing the cyclic denaturation temperature to 97 °C.

This study aimed at determining the strain diversity of MAP isolates responsible for paratuberculosis in Ugandan cattle and comparing them to strains reported from other countries. The results show that most of the MAP strains examined in the study (17 of 21) belonged to the cattle type of strains, while two isolates were of the bison type. A novel restriction pattern was revealed in two isolates (MapUg4 and MapUg2) indicating a type that has not been previously reported. This type was tentatively designated as X pending further characterisation. Cattle strains have been found in most areas where MAP has been detected in cattle. The bison type was first reported in the USA on cattle farms where bison infected with the same strain had been previously kept [10, 17, 18]. This strain has been reported to be widespread in India among buffaloes, cattle and goats [17, 19]. To our knowledge this is the first time the bison type is observed outside the USA and India. Since this is the first characterisation study of MAP in Uganda, it is not clear how widespread this type might be but it has the potential of spreading into more hosts than the cattle type [17, 20].

Based on the SSR locus 1, 15 isolates were divided into two allele groups having 6 and 7 G mononucleotide repeats respectively; while three allele profiles were revealed for SSR locus 2 (10 G, 11 G and 12 G) and four allele profiles for SSR locus 8: 4GGT, 2GGT and two polymorphic 3GGT repeats carrying mutations respectively, for all 21 isolates examined. Therefore, based on our results, the examination of the three loci is valuable in discriminating MAP isolates in this region, as has also been shown by El-Sayed et al. [15]. Although it has been reported that loci 1 and 8 were the most discriminative [20], in our study the three SSR loci complement each other and the highest level of discrimination is reached when all the three are used.

Harris et al. [21], found that the majority of the 211 isolates from the USA had more than 14 G1 and 9 G2 repeats for SSR loci 1 and 2 respectively and 5 GGT for SSR locus 8. On the other hand, El-Sayed et al. [15] found 7–14 G1 repeats, the majority being 7 G1 repeats, for SSR locus 1; 9–13 G2 repeats for SSR locus 2, the most common being 11 G among 34 German isolates. Regarding SSR locus 8, El-Sayed et al. [15] found that their German isolates had 4GGT and 5GGT repeats. In another study where SSR has been used, the most common profiles of SSR loci 1, 2 and 8 were >14 G-9 G-5GGT, >14-10-5GGT [21]. No study prior to this one has shown the occurrence of MAP strains with 6 G1 repeats in SSR locus 1, and 2GGT. Polymorphic 3GGT repeats has been observed in sheep strain [4] but the nature of the polymorphism was not defined.

According to Sevilla et al. [20] and Motiwala et al. [22], there is a relationship between the number of repeats in the three SSR loci of a particular strain and the host species. For instance, analysis of SSR locus 2 (G2 repeats) showed that 7 G were found in isolates from bison, impala, nyala, Thomson’s gazelles and goat, while 9 G were found in isolates from duikers, transcaspian urial and waterbucks. Elks were infected with isolates with 7 G and 13 G [19]. Sevilla et al. [20] observed that cattle isolates in Spain had 8, 9 and 11 G repeats in the same SSR locus, and that sheep and goats were mainly infected with isolates having 3GGT repeats on SSR locus 8. In contrast to the findings of El-Sayed et al. [15] who observed that all isolates with seven or greater number of G1 repeats had 5GGT repeats on SSR locus 8, this was not the case in the present study.

The significance of these SSR genotypes with regards to the potential differences in pathogenesis is not yet known but it suffices to say that the 7 G-4GGT profile seen in one isolate, is associated with MAP isolated from Crohn’s disease patients [4] and has also been found to be the most common type among cattle herds in Ohio, USA [23].

According to the present study, a combination of different loci comprising IS1311, MIRU (2, 3), VNTR 32 and SSR (loci 1, 2 and 8) discriminated 21 MAP isolates into 11 distinct strains. This shows very high strain diversity in Uganda. It implies that either the Ugandan isolates have evolved over a very long period of time or that they have all come from different regions in the recent past. Further studies to map out the distribution of these isolates and the historical origins of the cattle might clarify this question. It would also be of interest to determine if the different strains found in this study have any association with any particular breed or genotype of cattle and if the strains themselves differ in their pathogenecity. The most common genotype was number 10 with 10 isolates, with a 9-5-9-6 G-10 G-4GGT profile according to MIRU 2 and 3, VNTR 32 and SSR loci 1, 2 and 8 respectively (Table 1). Nine of these isolates came from abattoir specimens. Although we could not get reliable information on the districts of origin of the cattle from which these isolates were derived, we conjecture that they might be from western Uganda, since the same genotype was seen in a cow from Masindi district which is located in Western Uganda and also due to the fact that most of the cattle slaughtered in these two abattoirs come from western Uganda. The remaining genotype profiles had 1–3 isolates each and were evenly distributed in the districts of Wakiso, Mpigi and Luwero, where there are also several breeds of cattle and different husbandry practices.

Castellanos et al. [13] observed that type II isolates of MAP from Spain had five repeats at MIRU 3 locus, while type III had three repeats in contrast to some German type III isolates which had five repeats at MIRU 3 locus. In our study, 18 of the 20 isolates amplified, had five repeats on MIRU 3 locus indicating that they may be related to type II isolates. However, five of our isolates which had five repeats at MIRU 3 locus also had 15 repeats at MIRU 2. The 15 repeats at MIRU 2 locus has been associated with type I [13, 24]. As more strains continue to be isolated, the classification of MAP could become more complex than the simple type I, II and III system that was predominant in the last decade or new criteria will have to be devised to classify them into those types for epidemiological reasons. The isolates in this study differ from isolates reported in other studies especially with regard to their SSR profiles and should therefore be considered as novel strains.

Characterisation using SSR locus 8 SNP profile alone was only slightly less discriminatory than the combination of all three systems: MIRU-VNTR/SSR/IS1311 PCR-REA. Therefore phylogenetic analysis of SSR locus 8 could offer an opportunity for strain discrimination and follow up of the transmission patterns in the country, but it also failed to distinguish between several diverse isolates including the cattle, sheep and bison types found in GenBank (see Figure 3).

In our case MIRU 2 was found to be more discriminative than MIRU 3 and VNTR 32. This is in agreement with the findings of El-Sayed et al. [15] and Castellanos et al. [13]. It was also possible to amplify MIRU 2 from almost all the DNA templates that we had. Another positive aspect of the use of MIRU 2 is that it distinguishes MAP from M. intracellulare in which the locus is not amplified [25]. Finally VNTR 32 was the least discriminative and most difficult to amplify.

As stated in the materials and methods, in-silico analysis indicated high melting temperatures (>100 °C) for the GC-rich MIRU, VNTR and SSR templates. This was the possible reason for the failure of amplification we encountered in many cases using published primers and protocols [9, 12, 15, 19] (data not shown). In an effort to improve the efficiency of amplification, new specific primers were designed and several additives such as DMSO, Betaine, 1, 2, Propan-diol and Ethylene glycol that facilitate DNA denaturation [12, 26], were tested. The new primers were more effective than published primers both with and without the additives. With the new primers, we were able to obtain products from samples which had been too difficult to amplify despite the use of the additives. A cyclic denaturation step of 97 °C for 30 seconds was also adopted. The best results were obtained using Ethylene glycol. Previous studies on MIRUs and SSR have used DMSO and Betaine [12–14] but in our case those two denaturants were not very effective. Unfortunately despite the improvement of the PCR protocols we were not able to obtain amplicons for some of the isolates especially for SSR locus 1 and VNTR 32. Motiwala et al. [9] were also unable to amplify SSR locus 1 from some isolates.

There is great diversity among MAP isolates responsible for Paratuberculosis in Uganda. The greatest strain differentiation is obtained using a combination of different markers; however phylogenetic analysis of SSR locus 8 offers the single most discriminating tool for the differentiation of MAP isolates in Uganda. VNTR 32 is difficult to amplify and yet it offers no information on the diversity of MAP whereas, SSR locus 1 offers some information but failure of amplification is possible with some isolates. We conclude that the combined use of IS1311 PCR-REA, MIRU 2, 3, SSR 2, 8 and phylogenetic analysis of SSR locus 8 is more pragmatic in the characterisation of MAP isolates in Uganda but further studies may be required to validate this combination in other areas.