Genetic Diversity among Mycobacterium bovis Isolates: a Preliminary Study of Strains from Animal and Human Sources

Tuberculosis is caused by Mycobacterium tuberculosis, a slow-growing acid-fast bacterium which shares microbiological and biochemical characteristics with M. bovis. Together with other closely related organisms (M. microti, M. africanum, andM. canetti), these bacteria are referred to as the M. tuberculosis complex. M. bovis has the broadest host range of any member of the complex, causing disease in a wide range of mammals, including humans (21).

Members of the complex are traditionally characterized on the basis of a set of tests of microbiological and biochemical phenotypes, but these are time-consuming and consequently not always performed. Differentiation and drug susceptibility testing of M. tuberculosis complex members solely by conventional means may be complicated by the slow growth of the organisms in culture and their propensity to change growth characteristics (7). Species identification may sometimes be difficult with multiple-drug-resistant organisms; a nosocomial outbreak of multiple-drug-resistant M. tuberculosis was originally attributed incorrectly to M. bovis on the basis of phenotype (3, 14), which emphasises the benefits of genotyping methods, and indeed, a growing list of these is available.

Polymorphisms that have been used to discriminate between M. bovis and M. tuberculosis include single base changes in the pyrazinamidase gene (pncA) and in the oxyRpseudogene. M. bovis strains customarily contain guanine at position 169 in the pncA gene, and M. tuberculosis isolates contain cytosine at the same position. (23). This polymorphism renders the majority of M. bovis isolates resistant to treatment with pyrazinamide. TheoxyR pseudogene in M. bovis is characterized by adenine at position 285, and M. tuberculosis is characterized at the same position by guanine (24).

In the classification scheme proposed for the division of the members of the tuberculosis complex into three major groups defined by single nucleotide polymorphisms in the catalase (katG463)- and gyrase A (gyrA95)-encoding genes (25), M. bovis (and M. microti, M. africanum, and some M. tuberculosis isolates) displays polymorphisms characteristic of group I, having the genotypekatG463(Leu) gyrA95 (Thr). M. bovisand group I M. tuberculosis isolates also show greater differences in the oxyR pseudogene, which has been advanced as evidence for the ancestral nature of group I organisms (25).

M. bovis and M. bovis BCG vaccine strains have undergone a series of genomic deletion events compared to M. tuberculosis reference strains (9). M. microti and M. africanum appear to fall between theM. bovis and BCG strains and the reference strains in their deletion patterns (9). Associated with one of these deletions, namely, RD5, is the mtp40 element, a 3-kb fragment which encodes a section of a phospholipase C enzyme (9). This fragment is missing from the majority ofM. bovis isolates and provides the basis for a PCR to distinguish M. tuberculosis from M. bovis(17).

Differences in the direct repeat (DR) locus, which are exploited by the DNA fingerprinting technique of spoligotyping, have also been used to discriminate between members of the M. tuberculosis complex as well as between host-specific strains of M. bovis(5, 15, 26). M. bovis characteristically lacks spacers 39 to 43 in the spoligotype system. Compared to other members of the M. tuberculosis complex, M. microti has a much smaller number of variable spacers in the DR region (K. Kremer, D. van Soolingen, J. van Embden, S. Hughes, J. Inwald, and G. Hewinson, Letter, J. Clin. Microbiol.36:2793–2794, 1998).

Reports of studies applying biomolecular assays have occasionally described strains with unusual genotypic profiles between those ofM. tuberculosis and M. bovis. One variant isolated from seals and other pinnipeds (22, 26) displayspncA and oxyR polymorphisms identical to those ofM. tuberculosis but shares other biochemical and genetic properties with M. bovis. Caprine strains also appear to occupy an intermediate genotypic position, with pncA169 (C), typical of M. tuberculosis, and oxyR285 (A), characteristic of M. bovis (1, 6, 13). These strains also possess multiple copies of the repetitive element IS6110, a characteristic generally associated with M. tuberculosis (12, 27). These subtypes of M. bovis carry a zoonotic risk of disease for humans (13, 27), which adds clinical importance to their pyrazinamide-susceptible phenotype (20).

The present study, using a panel of diverse M. bovisisolates of human and animal origins, combines a number of these approaches and was undertaken to assess molecular assays as diagnostic tools as well as the extent of genotypic diversity in this taxon.

M. bovis BCG Pasteur and the H37Rv isolate of M. tuberculosis were used as the reference strains. A total of sevenM. tuberculosis complex strains, including some from exotic species, were provided by the Veterinary Laboratory Agency (VLA), New Haw, Surrey United Kingdom. The strains were chosen to typify infections from geographic hot spot regions of bovine tuberculosis in the United Kingdom and from unrelated outbreaks outside these areas (Table 1). Thus, bovine and badger isolates were from cases in Cornwall, United Kingdom, the fur seal isolate was from an outbreak at a zoo in southern Britain, and the goat and llama isolates were from mid Wales, United Kingdom.

Six isolates of M. bovis strains that cause human disease were obtained from the Public Health Laboratory Service at Dulwich, London, United Kingdom (Table 2). Species identification by standard biochemical assays was performed at each source laboratory (4).

DNA was extracted by a standard lysozyme-proteinase K-phenol chloroform procedure, followed by ethanol precipitation (10). To maximize the integrity of the genomic DNA, cultures were heat killed by boiling for 5 min at 100°C (2).

Genotyping of the various M. bovis isolates was performed by using PCR assays directed towards a series of loci with defined specificity patterns. PCR was performed as described previously, with a number of measures taken to avoid contamination (18). All oligonucleotide primers were obtained from Sigma-Genosys Ltd., Cambridge, United Kingdom. Table 3 lists the sequences of these primers and the parameters of the various methods used. Gel electrophoretic analysis of PCR products and automated DNA sequencing were performed as previously reported (18).

The mycobacterial genomic regions amplified were the RNA polymerase beta gene (rpoB), a hypervariable region of the 16S rRNA gene (16), the pyrazinamidase gene (pncA), and the oxyR pseudogene (oxyR), which contains species-defining polymorphisms (6, 23, 24). Analysis ofkatG463 and gyrA95 for grouping was performed according to the Sreevatsan model (25). The isolates were examined for the presence or absence of the mtp40 element in the plcA gene. Additionally, primers spanning deletion region 7 (RD7) (9) were designed to amplify a PCR product from strains which have undergone this event. For consistency, we have adopted the numerical classification of deletion regions proposed by Gordon and colleagues (9). Spoligotyping (15) was used for further strain identification.

All clinical and veterinary isolates showed microbiological and biochemical characteristics typical of M. bovis, including resistance to pyrazinamide, susceptibility to thiophene-2-carboxylic acid hydrazide (TCH), and preferential growth on media containing pyruvate (Tables 1 and 2).

The results from screening a panel of M. bovis isolates of human and animal origin by using a range of genotypic criteria emphasize the complex phylogeny of such strains. In all, five different genotypes were distinguished among the 13 isolates (Table4).

As expected, all of the isolates had the katG463 (CTG)gyrA95 (ACC) genotype characteristic of a group 1 organism (Table 5) (25). The M. tuberculosis H37Rv reference strain showed the predicted group 3 genotype. Furthermore, all the isolates as well as the reference strain (M. bovis BCG) generated the expected 211-bp product characteristic of the RD7 deletion event (9). As anticipated, no product was obtained from the M. tuberculosis laboratory strain H37Rv.

Diverse patterns were seen in the pncA and oxyRloci. Five of the veterinary isolates had the pncA169(G)-oxyR285 (A) genotype characteristic of M. bovis (Table 5). Typifying these, and consistent with their derivation from a single outbreak, the bovine and badger 1 isolates were identical according to all the criteria used. The strain is AF 2122/97, which was recently sequenced at the Sanger Centre and is now awaiting full description and commentary.

The typical M. bovis group included the feline isolate, which is a finding in contrast to those of previous reports of the acquisition of M. microti infections by cats through the hunting and consumption of small rodents (11).

Two isolates (from badger 2 and a fur seal) were found to have theM. tuberculosis pncA169 (C)-oxyR285 (G) genotype.

Three different patterns were observed among the animal isolates with genotypes between those of M. tuberculosis and M. bovis. Unexpectedly, three of the veterinary strains generated positive PCR products for mtp40, which is generally missing from M. bovis. These strains included the two isolates with the M. tuberculosis pncA169(C)-oxyR285 (G) genotype as well as the caprine isolate which typifies M. bovis at these loci. Unlike the strains isolated from Spanish goats (1, 20), our goat isolate demonstrated the M. bovis-resistant allelepncA169 (G). This strain, originally isolated from a goat in mid Wales in 1996, shared a spoligotype identical with that recovered from a captive llama in Gwent, South Wales, United Kingdom, in 1999. The two strains differed in their mtp40 status (Table 5). This finding illustrates the fact that spoligotyping is limited in its ability to differentiate between epidemiologically related strains; identical spoligotypes do not exclude differences at other loci (19).

As found previously, the fur seal isolate retained the RD5-associatedmtp40 element and showed polymorphisms at thepncA and oxyR loci characteristic of M. tuberculosis (22, 27). In contrast to the RD5 region, which was retained, the RD7 region had been deleted from all the isolates analyzed, including that from the fur seal, which suggests that this deletion event may have occurred early in the evolution of the M. tuberculosis complex. Undoubtedly, knowledge of deletion region subsets in M. bovis strains will play an increasingly important part in understanding the phenotype, pathogenesis, and evolution of individual isolates. The seal strain showed a typical M. bovis spoligotype in that it lacked spacers 39 to 43 in the DR region (Fig.1). With minor differences, the pattern was very similar to the two examples previously recorded for seals (27).

• Open in new tab
• Download powerpoint

Fig. 1.

Spoligotyping hybridization patterns obtained with animal and human isolates of M. bovis. The M. tuberculosis strain H37Rv and M. bovis BCG patterns are shown for reference. Filled and open circles indicate the presence and absence of spacers, respectively.

The spoligotype for badger 2 was characteristic of M. microti. The pattern is distinctive, with hybridization to only 3 of the possible 43 spacers of the DR region (spacers 4, 37, and 38) (Fig. 1). This pattern was previously found in a cat and a cow from the West Country, United Kingdom, is consistent with the transfer of disease between these species, and has been reported as part of a study of isolates from the United Kingdom and Holland (Kremer et al., letter). A second, similar pattern found in voles shows hybridization only to spacers 37 and 38 (Kremer et al., letter).

Five of the six human isolates and four of the animal isolates had classic M. bovis profiles in that they possessed the pyrazinamide-resistant pncA allele and lacked the RD5 region and the 3′ spoligotype spacers. The spoligotype patterns of these human isolates did not match any patterns of the animal isolates or any patterns in the M. bovis database kept at the VLA. Considering the ages of the patients (Table 2), it is likely that disease was the result of the reactivation of an infection acquired earlier in life.

In marked contrast, the isolate from patient 1 (Table 2) showed a pattern immediately suggestive of M. tuberculosis because of the retention of 3′-end spacers 40 to 43 (Fig. 1). This isolate was cultured from a patient of Indian origin who had been living in England for over 40 years and who presented with pulmonary tuberculosis. In spite of having the pncA169 (C) genotype, this strain displayed a pyrazinamidase-resistant phenotype in microbiological assays and deletion of the RD5 and RD7 regions, consistent withM. bovis. Spoligotyping data compiled on disease-producing strains of M. bovis in the United Kingdom have revealed very occasionally the presence of spacers in the 3′ spoligotype region; to date, a total of three other isolates with two or more spacers have been recorded (VLA, unpublished observations). Strains possessing spacers 4 and 5 of these 3′-terminal spacers have also been reported as appearing in isolates from patients born in Asia (5). Patient isolate 24 as described by this group (5) resembled patient isolate 1 in the present study in being negative for the mtp40 fragment.

Finally, the rpoB and 16S rRNA sequences were screened for differences as described in an earlier study of various mycobacterial species (8). Among the atypical mycobacterial species examined, intraspecies diversity in these genes was highest in M. smegmatis and M. gordonae (8). The study recorded a lack of genetic diversity among 63 M. tuberculosis isolates but did not include M. bovis or other members of the complex. Despite the diverse origins of theM. bovis isolates in the present study, no differences were found in either 738 bp of the rpoB gene or 208 bp of a hypervariable region of the 16S rRNA region, with all isolates displaying sequences identical to those of the reference strains. These findings therefore extend the observations made by Gingeras et al. (8) to members of the M. bovis lineage and are consistent with the lack of variation at these loci in the M. tuberculosis complex.

It was anticipated that the selection of strains from different animal sources might help uncover associations between genotype and host specificity in addition to exploring the scope of genotypic diversity. Our results are broadly consistent with those of previous reports which found that strains with a genotype between those of M. tuberculosis and M. bovis are often implicated in infections of fur seals and goats, but more studies at the molecular level are needed to confirm this. Investigations with multiple genetic markers can confound expectations: the cat in our study was not infected with M. microti, as might be expected, whereas one of the badgers was; the caprine subtype was not the cause of disease in the goat. Thus, while the current range of genotypic markers provides a useful resource for studying evolution within the M. tuberculosis complex, these basic observations suggest that they are probably not directly linked to the genetic determinants underlying host specificity.

• Copyright © 2001 American Society for Microbiology