# Mycobacterium bovis in Burkina Faso: Epidemiologic and Genetic Links between Human and Cattle Isolates

Adama Sanou, Zekiba Tarnagda, Estelle Kanyala, Dezemon Zingué, Moumini Nouctara, Zakaria Ganamé, Adjima Combary, Hervé Hien, Mathurin Dembele, Antoinette Kabore, et al.

► To cite this version:

Adama Sanou, Zekiba Tarnagda, Estelle Kanyala, Dezemon Zingué, Moumini Nouctara, et al. Mycobacterium bovis in Burkina Faso: Epidemiologic and Genetic Links between Human and Cattle Isolates. PLoS Neglected Tropical Diseases, 2014, 8 (10), pp.e3142. 10.1371/journal.pntd.0003142. hal-02967869

HAL Id: hal-02967869

https://hal.science/hal-02967869

Submitted on 15 Oct 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# *Mycobacterium bovis* in Burkina Faso: Epidemiologic and Genetic Links between Human and Cattle Isolates

Adama Sanou¹,²*, Zekiba Tarnagda³,⁴, Estelle Kanyala⁴, Dezemon Zingué², Moumini Nouctara², Zakaria Ganamé², Adjima Combary⁵, Hervé Hien²,³, Mathurin Dembele⁵, Antoinette Kabore², Nicolas Meda²,⁴, Philippe Van de Perre⁶,⁷,⁸, Dorine Neveu⁷,⁸, Anne Laure Bañuls¹,⁹, Sylvain Godreuil⁶,⁷,⁸,⁹

1 MIVEGEC, UMR IRD 224-CNRS 5290-Montpellier Universités 1 et 2, Montpellier, France, 2 Centre Muraz, Bobo-Dioulasso, Burkina Faso, 3 Institut de Recherche en Sciences de la Santé (IRSS), Direction Régionale de l'Ouest (DRO), Bobo Dioulasso, Burkina Faso, 4 Université de Ouagadougou, Ouagadougou, Burkina Faso, 5 Programme National Tuberculose, Ouagadougou, Burkina Faso, 6 Centre Hospitalier Régional Universitaire (CHRU) de Montpellier, Département de Bactériologie-Virologie, Montpellier, France, 7 Université Montpellier 1, Montpellier, France, 8 INSERM U 1058, Infection by HIV and by Agents with Mucocutaneous Tropism: From Pathogenesis to Prevention, Montpellier, France

<figure><img src="image_2.png" /></figure>

## Abstract

**Background:** In sub-Saharan Africa, bovine tuberculosis (bTB) is a potential hazard for animals and humans health. The goal of this study was to improve our understanding of bTB epidemiology in Burkina Faso and especially *Mycobacterium bovis* transmission within and between the bovine and human populations.

**Methodology/principal findings:** Twenty six *M. bovis* strains were isolated from 101 cattle carcasses with suspected bTB lesions during routine meat inspections at the Bobo Dioulasso and Ouagadougou slaughterhouses. In addition, 7 *M. bovis* strains were isolated from 576 patients with pulmonary tuberculosis. Spoligotyping, RDAf1 deletion and MIRU-VNTR typing were used for strains genotyping. The isolation of *M. bovis* strains was confirmed by spoligotyping and 12 spoligotype signatures were detected. Together, the spoligotyping and MIRU-VNTR data allowed grouping the 33 *M. bovis* isolates in seven clusters including isolates exclusively from cattle (5) or humans (1) or from both (1). Moreover, these data (genetic analyses and phenetic tree) showed that the *M. bovis* isolates belonged to the African 1 (Af1) clonal complex (81.8%) and the putative African 5 (Af5) clonal complex (18.2%), in agreement with the results of RDAf1 deletion typing.

**Conclusions/Significance:** This is the first detailed molecular characterization of *M. bovis* strains from humans and cattle in Burkina Faso. The distribution of the two Af1 and putative Af5 clonal complexes is comparable to what has been reported in neighbouring countries. Furthermore, the strain genetic profiles suggest that *M. bovis* circulates across the borders and that the Burkina Faso strains originate from different countries, but have a country-specific evolution. The genetic characterization suggests that, currently, *M. bovis* transmission occurs mainly between cattle, occasionally between cattle and humans and potentially between humans. This study emphasizes the bTB risk in cattle but also in humans and the difficulty to set up proper disease control strategies in Burkina Faso.

**Citation:** Sanou A, Tarnagda Z, Kanyala E, Zingué D, Nouctara M, et al. (2014) *Mycobacterium bovis* in Burkina Faso: Epidemiologic and Genetic Links between Human and Cattle Isolates. PLoS Negl Trop Dis 8(10): e3142. doi:10.1371/journal.pntd.0003142

**Editor:** Joseph M. Vinetz, University of California, San Diego, United States of America

**Received:** March 20, 2014; **Accepted:** July 23, 2014; **Published:** October 2, 2014

**Copyright:** © 2014 Sanou et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability:** The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was funded by the French National Agency for Research on AIDS and viral hepatitis (ANRS), the Global Fund and the Ministry of Animal Resources and Fishery of Burkina Faso. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Introduction

*Mycobacterium bovis* is the causative agent of bovine tuberculosis (bTB) in a broad spectrum of hosts, such as cattle, goats, sheep and wild animals to which it can be transmitted through the oral or respiratory route [1]. Humans also can acquire *M. bovis* generally through the aerogenous route when in close contact with infected animals, but also by consuming unpasteurized dairy products from infected animals and through the skin when handling infected carcasses [2,3,4].

In sub-Saharan Africa, bTB is a serious problem for livestock production but also a health risk for humans as most human

populations live in close contact with domestic animals in which the disease is highly prevalent and imperfectly controlled [5]. Therefore, bTB has a deleterious economic burden, although this has not been quantified in Africa as yet [5,6]. In Burkina Faso, little is known about bTB epidemiology and national strategies for disease control are almost non-existent [7]. Tuberculin testing of livestock is not routinely performed and bTB screening is limited to visually checking the meat in abattoirs. Despite the high prevalence of bTB in cattle and the presence of *M. bovis* in 26.5% of unpasteurized milk samples, the zoonotic transmission of bTB is also poorly known [8]. In Burkina Faso, cattle breeding relies mainly on extensive transhumance and is the prerogative of few

<footnote>* Email: sanouadama2003@yahoo.fr</footnote>

<footnote>¶ These authors contributed equally to this work.</footnote>
## Author Summary

Bovine tuberculosis is an infectious disease caused by *Mycobacterium bovis* in livestock and wild animals. Humans can acquire this germ by aerogenous route when in close contact with infected animals, or by consuming unpasteurized dairy products from infected animals and also through the skin when handling infected carcasses. For the present study in Burkina Faso, *M. bovis* strains were collected from slaughtered animals during routine veterinary inspection at the slaughterhouses of Bobo Dioulasso and Ouagadougou and also from patients with suspected pulmonary tuberculosis. The isolates were genetically characterized using three techniques: spoligotyping, MIRU-VNTR and RDAf1 deletion analysis. Our results highlight two aspects of *M. bovis* epidemiology that are crucial for disease control: i) *M. bovis* circulates between Burkina Faso and its neighbouring countries and ii) *M. bovis* is transmitted mainly between cattle, but also between cattle and humans, and potentially between humans. This study stresses the need to develop an efficient strategy to control *M. bovis* transmission, but also the difficulty to implement control measures because of the complex epidemiology of bovine tuberculosis in Burkina Faso.

ethnic groups [7]. These populations live in close and permanent contact with their livestock and consume raw and unpasteurized dairy products that could favour *M. bovis* transmission [8]. Moreover, in slaughterhouses where cattle are slaughtered, butchers wear minimal protective clothing and handle infected offals with bare hands [8]. These close contacts are an important source of zoonotic transmission [7,8,9]. Although culture and species identification of the *M. tuberculosis* complex are not routinely performed in Burkina Faso, previous studies have suggested that *M. bovis* is present in 0.4% to 1.4% of isolates from patients with pulmonary tuberculosis [10,11].

The main goal of this work was thus to genetically characterize the *M. bovis* population in Burkina Faso in order to improve our understanding of bTB epidemiology and the circulation of *M. bovis* within and between the bovine and human populations in this country. We used spoligotyping and mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) analysis, because a combination of them is a powerful tool for the genetic characterization of *M. bovis* [12,13]. We also explored the presence or absence of a specific chromosomal region called RDAf1, which is a region of difference in *M. bovis* [14].

## Materials and Methods

### Human and cattle samples

The tissues showing macroscopic lesions compatible with bTB were collected from slaughtered cattle carcasses during the post-mortem inspection at the slaughterhouses of Ouagadougou and Bobo Dioulasso (the two largest cities of Burkina Faso) between May and October 2011. The cattle slaughtered in these two slaughterhouses are mainly originated from neighbouring villages of Ouagadougou and Bobo Dioulasso, as well as the main areas of cattle production throughout the country and the cattle production system is mainly pastoralism. Samples were collected and transported in sterile containers at 4°C to the Mycobacteria Laboratory of the Muraz Centre (in Bobo Dioulasso) for analysis.

Smear positive sputum specimens of patients with suspected pulmonary TB were collected in the framework of two studies in Burkina Faso: i) a nationwide survey on anti-tuberculosis drug resistance between 2009 and 2011 (funded by the Global Fund); ii) and a regional study in the Hauts Bassins area between 2011 and 2013 (funded by the French National Agency for Research on AIDS and Viral Hepatitis, ANRS 1224 project, "Impact of HIV/ *Mycobacterium tuberculosis* co-infection on the dynamics of tuberculosis transmission in Burkina Faso").

### Isolation and identification of mycobacteria

Bovine tissue samples were processed for mycobacteria isolation following the standard procedures described by the World Organization for Animal Health [15]. Briefly, tissues with tuberculous lesions were dissected into pieces using sterile scissors and forceps, and were then crushed using sterile sea-sand, mortar and pestle. The homogenate of each sample was recovered into 50 ml sterile tube with 10 ml of sterile distilled water, and the obtained solution was homogenized on a vortex mixer for few minutes. After 5 minutes of settling, two milliliters of supernatant were decontaminated with 10 ml of NaOH at 4% according to the protocol described by Petroff's method. Patients' sputum specimens were treated according to the Petroff's method too. The obtained suspensions were inoculated in four Lowenstein-Jensen (LJ) slants, two of which were supplemented with 0.2% of sodium pyruvate. Isolates were identified as mycobacteria and as *M. bovis* species by Ziehl-Neelsen (ZN) staining for Acid Fast Bacilli (AFB) and a conventional biochemical method, previously described by Ledru et al. [11].

### Genotyping

DNA extraction from mycobacteria isolates and high-throughput spoligotyping on Luminex 100 (Luminex Corp., TX) were performed as previously described [16,17]. The obtained data were compared with those of the international databases SpoDB4.0 [18] and http://www.Mbovis.org [19]. The new spoligotypes were submitted to the http://www.Mbovis.org database, and new SB numbers (spoligotype codes) were assigned accordingly.

Human and cattle isolates were also genotyped by PCR amplification of 26 MIRU-VNTR loci: ETR A, B, C, D, E; QUB-11a, 11b, 26, 4156, 3232; MIRU 2, 10, 16, 20, 23, 24, 26, 27, 39, 40 and Mtub 04, 21, 29, 30, 34, 39 [20,21]. We used multiplex PCR and capillary electrophoresis-based sequencers (ABI 3730-XL), as previously described [22]. PCR fragment sizing and assignment of the different MIRU-VNTR alleles were done using Genemapper, version 4.0 (PE Applied Biosystems). The results for each of the 26 loci were combined into 26-digit allelic profiles [22].

Finally, the presence of the RDAf1 deletion was determined using a multiplex PCR method with a set of three primers followed by agarose gel electrophoresis according to Müller et al. [14].

### Genetic diversity and population structure analyses

Several diversity indices, including the genotypic diversity (Gd = the number of different genotypes divided by the total number of samples using the combination of MIRU-VNTR and Spoligotyping data), the allelic diversity per locus and the mean genetic diversity (Hs) were calculated. The population structure was explored by estimating the Fst (index of genetic differentiation between samples) value (0 = no differentiation and 1 = fixation of alternative alleles). The allelic diversity, the Hs and the Fst were calculated using F-STAT, version 2.9.3 with the 26 MIRU-VNTR loci [23].
## Phenetic tree and statistical analyses

Genetic relationships among isolates were built with the UPGMA (unweighted pair group method with arithmetic average) and NJ (Neighbour Joining) clustering methods using the MIRU-VNTR and spoligotyping data. The Phylip and Populations packages [24] were used for tree elaboration based on the Nei's distance, and Treedyn for tree visualization and annotation [25].

## Ethical considerations

The recruitment of human patients and the collection of bovine samples were done according to protocols approved by the Ethics Committee for Health Research in Burkina Faso (2007-031; June 28, 2009 and 2010-049; 7 July 2010) and by the Ministry of Animal Resources and Fishery. All patients and cattle owners provided written informed consent. The bovine study was conducted according to guidelines recommended by the Government of Burkina Faso (KITI n° AN VII 114 FP-AGRI-EL portant réglementation de la santé publique vétérinaire au Burkina Faso).

## Results

### *M. bovis* strains in human and cattle samples

Among the sputum samples collected between 2009 and 2013, the presence of *M. bovis* was biochemically confirmed in 5/269 (1.85%) samples from the nationwide survey and in 2/307 (0.65%) samples from the regional study respectively. Patients originated from cities that were quite distant one from each other (Table 1). Out of the 6 patients (6/7) with available HIV serology, only one (1/6) was HIV-positive.

Among the 1499 cattle carcasses inspected between May and October 2011, suspicious TB lesions were detected in 101 (6.74%) and 48/101 (47.5%) had a mycobacterial culture positive for AFB. Of the 48 strains isolated, 26 were biochemically identified as *M. bovis*. The remaining 22 samples were either nontuberculous mycobacteria (2/22), or *M. africanum* (3/22), or *M. tuberculosis* (2/22), or were contaminated (9/22) or had an insufficient growth (6/22).

## Genetic characterization

Spoligotyping confirmed species identification for all 33 *M. bovis* isolates. Among the 12 spoligotype signatures obtained (Table 2), only five (SB1398, SB0300, SB0857, SB0944, SB1439)

were already described in the http://www.Mbovis.org database. New codes (SB2282, SB2283, SB2284, SB2285, SB2286, SB2287 and SB2288) have been assigned to the 7 spoligotypes absent in the database. Based on their spoligotype signature, 25 isolates (75.8%) were clustered in 4 groups. The largest cluster included 17 strains bearing the SB0944 spoligotype profile; two clusters had only three strains each (SB0300 and SB2286 spoligotypes) and the last cluster included two strains with the SB1398 spoligotype. Each of the other eight spoligotype signatures (SB0857, SB1439, SB2282, SB2283, SB2284, SB2285, SB2287 and SB2288) was found in a single *M. bovis* isolate (Fig. 1 and Table 2).

By MIRU-VNTR typing, only one (human sample h7, spoligotype = SB2284, pattern P4) of the 33 isolates could not be genotyped. MIRU-VNTR typing produced 24 distinct patterns (Fig. 1 and Table 2). Fifteen isolates (46.9%) were grouped in seven clusters and the other 17 isolates (53.1%) had a specific pattern each. Moreover analysis of the 26 MIRU-VNTR loci, showed that nine loci were monomorphic in the 32 *M. bovis* strains analysed (see Table 3), indicating the presence of moderate polymorphism. The five most discriminatory loci were ETR A, ETR B, QUB-11a, QUB-26 and MIRU 26. The number of alleles for the MIRU-VNTR loci ranged from 1 to 6, with a mean of 2.5 (Table 3). The combination of spoligotyping and MIRU-VNTR typing revealed that fourteen isolates (42.4%) were grouped in seven clusters and the other 18 isolates (57.6%) had a unique pattern. The seven clusters (14 strains) included isolates from cattle (c4-c5, P3; c17-c13, P9; c20-c21, P18; c22-c23, P22; c1-c11, P25), humans (h5-h6, P7) and from humans and cattle (c19-h4, P20) (see Fig. 1 and Table 2). The MIRU-VNTR pattern P25 was split up in two spoligotypes and these results were confirmed by retesting using the two genotyping methods (Table 2).

The RDAfl deletion was detected in 27/33 isolates (81.8%). For these samples, as expected, the spoligotyping signatures revealed the absence of spacer 30.

## Phenetic analysis

We obtained comparable trees with UPGMA and NJ methods (data not shown). Only the UPGMA is presented here in order to facilitate the confrontation of tree with the spoligotyping data (Fig. 1). From the tree, we could distinguish two groups (Fig. 1). Group I had six strains (18.2%) that were characterized by the absence of spacers 4 and 5, the presence of spacer 30 in the spoligotype signatures and of the RDAfl region. This group of

<table><caption>**Table 1.** Socio-demographic informations about hosts and *M. bovis* isolates.</caption><thead><tr><td>Strain ID</td><td>Host</td><td>Geographical location</td><td>HIV status<sup>a</sup></td><td>Year of isolation</td><td>Number</td></tr></thead><tbody><tr><td>c1-c7</td><td>cattle</td><td>Bobo-Dioulasso</td><td></td><td>2011</td><td>7<sup>b</sup></td></tr><tr><td>c8-c26</td><td>cattle</td><td>Ouagadougou</td><td></td><td>2011</td><td>19<sup>b</sup></td></tr><tr><td>h1</td><td>human</td><td>Solenzo</td><td>HIV-1+</td><td>2009</td><td>1<sup>c</sup></td></tr><tr><td>h2</td><td>human</td><td>Ouagadougou</td><td>unknown</td><td>2010</td><td>1<sup>c</sup></td></tr><tr><td>h3</td><td>human</td><td>Ouagadougou</td><td>HIV-</td><td>2010</td><td>1<sup>c</sup></td></tr><tr><td>h4</td><td>human</td><td>Koupéla</td><td>HIV-</td><td>2011</td><td>1<sup>c</sup></td></tr><tr><td>h5</td><td>human</td><td>Bobo-Dioulasso</td><td>HIV-</td><td>2011</td><td>1<sup>c</sup></td></tr><tr><td>h6</td><td>human</td><td>Bobo-Dioulasso</td><td>HIV-</td><td>2011</td><td>1<sup>d</sup></td></tr><tr><td>h7</td><td>human</td><td>Bobo-Dioulasso</td><td>HIV-</td><td>2013</td><td>1<sup>d</sup></td></tr></tbody></table>

<footnote><sup>a</sup>specific to human hosts,</footnote>

<footnote><sup>b</sup>bovine study,</footnote>

<footnote><sup>c</sup>nationwide survey,</footnote>

<footnote><sup>d</sup>regional study.</footnote>

<footnote>doi:10.1371/journal.pntd.0003142.t001</footnote>
<table><caption>Table 2. Cont.</caption><thead><tr><th>SB number</th><th>Spoligotype signature<sup>1</sup></th><th>MIRU-VNTR (patterns, codes and number of strains)<sup>2</sup></th><th>Strains</th><th>Clone [no (%)]<sup>3</sup></th></tr></thead><tbody><tr><td></td><td></td><td>3 5 5 3 11 3 1 3 1 2 3 2 4 2 4<br>3 2 2 3 5 4 3 2 (P22, n=2)</td><td>c22, c23</td><td></td></tr><tr><td></td><td></td><td>3 6 5 3 11 3 1 3 1 2 2 3 2 4 2 4<br>3 2 2 3 5 4 3 2 (P23, n=1)</td><td>c6</td><td></td></tr><tr><td></td><td></td><td>4 5 5 3 10 3 1 3 2 2 3 2 4 2 4 3<br>2 2 3 5 3 3 2 (P24, n=1)</td><td>c26</td><td></td></tr><tr><td></td><td></td><td>4 5 5 3 10 3 1 2 0 2 2 3 2 4 2 4<br>3 2 2 3 5 4 3 2 (P25, n=2)</td><td>c1, c11</td><td></td></tr><tr><td>SB1439</td><td></td><td>4 5 5 3 10 3 1 2 0 2 2 3 2 4 2 4<br>3 2 2 3 5 4 3 2 (P25, n=1)</td><td>c10</td><td>Af1 [1 (3%)]</td></tr><tr><td>SB2288</td><td></td><td>4 5 5 3 11 3 1 4 1 2 2 3 2 4 2 4<br>3 2 2 3 5 4 3 2 (P17, n=1)</td><td>h3</td><td>Af1 [1 (3%)]</td></tr><tr><td>Total</td><td>12 Spoligotype signatures</td><td>25 MIRU-VNTR patterns</td><td>33 (100%)</td><td>Af5 [6 (18.2%)]<br>Af1 [27 (81.8%)]</td></tr></tbody></table>

<footnote n="1">¹■ presence of spacer; ☐, absence of spacer.</footnote>

<footnote n="2">²MIRU-VNTR loci: ETR A, ETR B, ETR C, ETR D, ETR E, QUB-11a, QUB-11b, QUB-3232, QUB-26, QUB-4156, MIRU 2, MIRU 10, MIRU 16, MIRU 20, MIRU 23, MIRU 24, MIRU 26, MIRU 27, MIRU 39, MIRU 40, M tub 04, M tub 21, M tub 29, M tub 30, M tub 34, M tub 39. NA = Not Available.</footnote>

<footnote n="3">³Af1 = African 1 clonal complex; Af5 = putative African 5 clonal complex.</footnote>

<footnote>doi:10.1371/journal.pntd.0003142.t002</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1. UPGMA tree based on the MIRU-VNTR (26 loci) and spoligotyping data.** 1, SB number = name of spoligotype based on http://www.bovis.org database nomenclature; 2, RDAf1 = Genomic deletion specific to Af1 clonal complex; 3, The MIRU-VNTR patterns are detailed in table 2. doi:10.1371/journal.pntd.0003142.g001</figcaption></figure>

strains appeared similar to the one described by Müller et al. [26], provisionally called African 5 (Af5) clonal complex. Even if these strains revealed specific spoligotype signatures as described above, no genomic deletion allowed to characterize this group as a well defined clonal complex [14,26,27,28]. In this study, this group of strains is thus called “putative Af5 clonal complex”. Group II included 27 strains that were assigned to the African 1 (Af1) clonal complex based on two criteria: (i) the absence of spacer 30 in the spoligotype; and (ii) the presence of the specific RDAf1 deletion in the genome (Fig. 1) [14].

## Genetic diversity and population structure of *M. bovis* in Burkina Faso

The mean genetic diversity ($H_s$) and the genotypic diversity ($G_d$) were respectively 0.187 and 0.79. As expected, the genetic differentiation between Af1 and the putative Af5 strains was high and significant [group II (n = 27) versus group I (n = 5); $F_{st} = 0.35$; $p<0.05$]. In addition to the RDAf1 deletion and spoligotype signatures, three MIRU-VNTR loci (QUB-3232, QUB-26 and MIRU 26) allowed assigning the isolates to the Af1 or the putative Af5 clonal complex because the number of repetitions (n) is ≥5 for putative Af5 and n≤1 for Af1 with QUB-3232 locus, n≥5 for putative Af5 and n≤4 for Af1 with locus QUB-26, n = 7 or 1 for putative Af5, and 1<n<7 for Af1 with MIRU 26. Nevertheless, these differences should be confirmed on a larger sample with isolates from different regions.

## Discussion

We present here the first detailed molecular characterization of *M. bovis* strains from humans and cattle in Burkina Faso. The 6.8% prevalence of bTB in cattle recorded in the present study on the basis of observable lesions in carcasses is lower than in previously published works. A 19% prevalence was reported for the Bobo Dioulasso slaughterhouse, fifteen years ago [8].

Nevertheless, our finding is in agreement with results in other Sahelian countries [6,29,30,31]. Despite this still high prevalence of bTB in cattle in our study, the prevalence of *M. bovis* in human TB was low (0.6–1.85%), as that reported in other studies in this country [7,8,10,11]. However, prevalences could have been underestimated because only patients with pulmonary TB were included in the present and previous studies. Classically, *M. bovis* is mainly responsible for extra-pulmonary TB [8]. In any case, the presence of glycerol in LJ medium may have negative impact on the yield of *M. bovis* growth. The use of LJ medium with 0.5% of sodium pyruvate and without glycerol could allow the optimization of this mycobacterial species isolation in Burkina Faso.

In our study, the combination of spoligotyping and RDAf1 deletion analysis showed that all *M. bovis* strains belonged either to the Af1 clonal complex, also found in other countries of West-Central Africa, such as Mali, Cameroon, Nigeria, Chad and Niger [14,32], or to the putative Af5 clonal complex, also previously described in Mali [26]. In Burkina Faso, the putative Af5 strains are geographically located in the Western Region (Bobo Dioulasso and Solenzo), an area bordering Mali. The presence of common or related genotypes between Burkina Faso and Mali can be explained by the transhumance activities between these countries and by the transit of Malian livestock on their way to the South, for instance to Ghana and Nigeria. Spoligotype signatures belonging to the putative Af5 clonal complex have also been reported in Europe [14,33]. The putative Af5 clonal complex could have been introduced in these African regions from Europe, possibly via North Africa [26,33,34]. In the Af1 clonal complex, the SB0944 spoligotype signature is defined as the most recent common ancestor (progenitor) and is the most frequent pattern within this group (40% in Chad, 46.1% in Nigeria and 62.7% in Cameroon). It was also the most abundant (52%) in our study [14]. The spread of the Af1 clonal complex over this large area of West-Central Africa could be explained by the long distance transhumance for livestock production in the Sahel, mainly practised by
<table><caption>**Table 3.** Allelic diversity of the 26 MIRU-VNTR loci in *M. bovis* isolates from humans and livestock in Burkina Faso.*</caption><thead><tr><th rowspan="2">Locus</th><th colspan="3">Number of alleles</th><th colspan="3">Allelic diversity</th></tr><tr><th>Global (n=32)</th><th>Af1 (n=27)</th><th>Af5 (n=5)</th><th>Global (n=32)</th><th>Af1 (n=27)</th><th>Af5 (n=5)</th></tr></thead><tbody><tr><td>ETR A</td><td>4</td><td>3</td><td>2</td><td>0.66</td><td>0.6</td><td>0.4</td></tr><tr><td>ETR B</td><td>4</td><td>4</td><td>2</td><td>0.54</td><td>0.44</td><td>0.4</td></tr><tr><td>ETR C</td><td>3</td><td>3</td><td>1</td><td>0.12</td><td>0.15</td><td>0</td></tr><tr><td>ETR D</td><td>2</td><td>1</td><td>2</td><td>0.06</td><td>0</td><td>0.4</td></tr><tr><td>ETR E</td><td>2</td><td>2</td><td>1</td><td>0.12</td><td>0.14</td><td>0</td></tr><tr><td>QUB-11a</td><td>5</td><td>4</td><td>2</td><td>0.60</td><td>0.63</td><td>0.4</td></tr><tr><td>QUB-11b</td><td>2</td><td>2</td><td>1</td><td>0.31</td><td>0.36</td><td>0</td></tr><tr><td>QUB-3232</td><td>5</td><td>3</td><td>4</td><td>0.34</td><td>0.15</td><td>0.9</td></tr><tr><td>QUB-26</td><td>6</td><td>4</td><td>3</td><td>0.76</td><td>0.70</td><td>0.7</td></tr><tr><td>QUB-4156</td><td>2</td><td>2</td><td>1</td><td>0.20</td><td>0.22</td><td>0</td></tr><tr><td>MIRU 2</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>MIRU 10</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>MIRU16</td><td>2</td><td>2</td><td>1</td><td>0.18</td><td>0.2</td><td>0</td></tr><tr><td>MIRU 20</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>MIRU 23</td><td>2</td><td>1</td><td>2</td><td>0.12</td><td>0</td><td>0.6</td></tr><tr><td>MIRU 24</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>MIRU 26</td><td>6</td><td>5</td><td>2</td><td>0.71</td><td>0.62</td><td>0.4</td></tr><tr><td>MIRU 27</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>MIRU 39</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>MIRU 40</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Mtub 04</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Mtub 21</td><td>4</td><td>4</td><td>2</td><td>0.24</td><td>0.21</td><td>0.4</td></tr><tr><td>Mtub 29</td><td>2</td><td>2</td><td>1</td><td>0.12</td><td>0.14</td><td>0</td></tr><tr><td>Mtub 30</td><td>2</td><td>2</td><td>2</td><td>0.27</td><td>0.07</td><td>0.4</td></tr><tr><td>Mtub 34</td><td>3</td><td>3</td><td>1</td><td>0.12</td><td>0.14</td><td>0</td></tr><tr><td>Mtub 39</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Mean</td><td>2.5</td><td>2.15</td><td>1.5</td><td>0.21</td><td>0.18</td><td>0.19</td></tr></tbody></table>

*Excluding one strain of the putative African 5 clonal complex that hasn't MIRU-VNTR data.
Af1 = African 1 clonal complex, Af5 = putative African 5 clonal complex.
doi:10.1371/journal.pntd.0003142.t003

the Fulani pastoralists [26]. This points out the difficulty to develop an efficient strategy to stop bTB transmission.

As expected, MIRU-VNTR typing revealed more polymorphisms than spoligotyping with a high genotypic diversity, but a low genetic diversity because the majority of the analysed loci presented a low mean allelic diversity. Despite a significant genetic differentiation, the low number of isolates analyzed in this study does not allow any assumption about the chronological emergence of these two groups of strains. Nevertheless, to explain the predominance of the Af1 clonal complex in many countries, Müller and al. have suggested that Af1 might have a selective advantage compared to the putative Af5 [14].

When we compared the *M. bovis* population of Burkina Faso with the populations from Mali, Chad, Nigeria and Cameroon described in Müller et al. [14], we observed a significant genetic differentiation (data not shown). Studies conducted in Chad and Nigeria showed inter-country variation in terms of discriminatory power of MIRU-VNTR loci. These studies and the present one had eight markers in common and only ETR A and ETR B had very high discriminatory power, while MIRU 2 had a very low or no discriminatory power [12,35]. On the basis of the polymorphic

loci, each country presents a specific pool of genotypes, although they also share several genotypes with the neighbouring countries. As suggested by Müller et al., these different genetic patterns could be explained by specific evolutionary processes (such as genetic drift and/or selection pressure) depending on the ecosystem set after the emergence and spread of Af1 and putative Af5 progenitors in each country [14].

From an epidemiological point of view, despite the absence of reliable and accurate conventional data the finding that specific spoligotype signatures/MIRU-VNTR patterns were shared only by bovine isolates or by bovine and human isolates suggests a recent transmission within the cattle population and between cattle and humans. Moreover, two patients were infected by Af1 strains with the same genotype (spoligotype plus MIRU-VNTR results). The sputum samples from these two patients were processed at the same medical centre (TB Diagnosis and Treatment Centre of Dafra, Bobo Dioulasso) and at the Mycobacteria Laboratory of the Muraz Centre in 2011, but not at the same time, thus excluding intra-laboratory contamination. Nevertheless, the epidemiological link could not be definitively established for this cluster. Different routes of contamination could be responsible for the infection: i)
contamination by a common animal or food-borne source; ii) human to human transmission. Indeed inter-human transmission cannot be excluded since intra-familial and community based transmissions of human TB cases due to Afl strains of *M. bovis* have been already suspected [36,37].

In conclusion, our study shows that two groups of *M. bovis* circulate in Burkina Faso; a major group belonging to the Afl clonal complex and a minor group belonging to the putative Afl5 clonal complex. Furthermore, the comparison with data from other African regions indicates an inter-country transmission associated with a country-specific evolution. Finally, the clusters suggest current transmission that occurs mainly within cattle populations, less frequently between cattle and humans and possibly between humans. This study points out the difficulty to develop an efficient national control strategy of bTB in Burkina Faso.

## References

1. Neill SD, Pollock JM, Bryson DB, Hanna J (1994) Pathogenesis of Mycobacterium bovis infection in cattle. Vet Microbiol 40: 41-52.

2. Kaneene JB, Miller R, de Kantor IN, Thoen CO (2010) Tuberculosis in wild animals. Int J Tuberc Lung Dis 14: 1508-1512.

3. Belchior I, Seabra B, Duarte R (2011) Primary inoculation skin tuberculosis by accidental needle stick. BMJ Case Rep 2011: pii: bcr1120103496. doi: 10.1136/bcr.11.2010.3496.

4. Gavier-Widen D, Chambers MA, Palmer N, Newell DG, Hewinson RG (2001) Pathology of natural Mycobacterium bovis infection in European badgers (Meles meles) and its relationship with bacterial excretion. Vet Rec 148: 299-304.

5. Ayele WY, Neill SD, Zinstag J, Weiss MG, Pavlik I (2004) Bovine tuberculosis: an old disease but a new threat to Africa. Int J Tuberc Lung Dis 8: 924-937.

6. Diguimbaye-Djaibe C, Hilty M, Ngandolo R, Mahamat HH, Pfyffer GE, et al. (2006) Mycobacterium bovis isolates from tuberculous lesions in Chadian zebu carcasses. Emerg Infect Dis 12: 769-771.

7. Rey JL, Villon A, Saliou P, Gidel R (1986) [Tuberculosis infection in a cattle-breeding region in Sahelian Africa]. Ann Soc Belg Med Trop 66: 235-243.

8. Vekemans M, Cartoux M, Diagbouga S, Dembele M, Kone B, et al. (1999). Potential source of human exposure to Mycobacterium bovis in Burkina Faso, in the context of the HIV epidemic. Clin Microbiol Infect 5: 617-621.

9. Gidel R, Albert JP, Lefevre M, Menard M, Retif M (1969) [Mycobacteria of animal origin isolated by the Muraz Center from 1965 to 1968: techniques of isolation and identification; results]. Rev Elev Med Vet Pays Trop 22: 495-508.

10. Gomgnimbou MK, Refregier G, Diagbouga SP, Adama S, Kabore A, et al. (2012) Spoligotyping of Mycobacterium africanum, Burkina Faso. Emerg Infect Dis 18: 117-119.

11. Ledru S, Cauchoux B, Yameogo M, Zoubga A, Lamande-Chiron J, et al. (1996) Impact of short-course therapy on tuberculosis drug resistance in South-West Burkina Faso. Tuber Lung Dis 77: 429-436.

12. Hilty M, Diguimbaye C, Schelling E, Baggi F, Tanner M, et al. (2005) Evaluation of the discriminatory power of variable number tandem repeat (VNTR) typing of Mycobacterium bovis strains. Vet Microbiol 109: 217-222.

13. Smith NH, Gordon SV, de la Rua-Domenech R, Clifton-Hadley RS, Hewinson RG (2006) Bottlenecks and broomsticks: the molecular evolution of Mycobacterium bovis. Nat Rev Microbiol 4: 670-681.

14. Muller B, Hilty M, Berg S, Garcia-Pelayo MC, Dale J, et al. (2009) African 1, an epidemiologically important clonal complex of Mycobacterium bovis dominant in Mali, Nigeria, Cameroon, and Chad. J Bacteriol 191: 1951-1960.

15. OIE (2008) Manual of diagnostic tests and vaccines for terrestrial animals (mammals, birds and bees). 1: 598p.

16. van Soolingen D, de Haas PE, Hermans PW, van Embden JD (1994) DNA fingerprinting of Mycobacterium tuberculosis. Methods Enzymol 235: 196-205.

17. Zhang J, Abadia E, Refregier G, Tafaj S, Boschiroli ML, et al. (2010) Mycobacterium tuberculosis complex CRISPR genotyping: improving efficiency, throughput and discriminative power of 'spoligotyping' with new spacers and a microbead-based hybridization assay. J Med Microbiol 59: 285-294.

18. Brudey K, Driscoll JR, Rigouts L, Proddinger WM, Gori A, et al. (2006) Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol 6: 23.

19. Smith NH, Upton P (2012) Naming spoligotype patterns for the RD9-deleted lineage of the Mycobacterium tuberculosis complex; www.Mbovis.org. Infect Genet Evol 12: 873-876.

## Acknowledgments

We would like to thank INSERM, IRD, CNRS and the Ministry of Animal Resources and Fishery of Burkina Faso for technical support. We thank Elisabetta Andermarcher for assistance in preparing and editing the manuscript.

## Author Contributions

Conceived and designed the experiments: ZT EK AC MD ZD HH NM PVdP DN ALB SG. Performed the experiments: AS EK MN ZG AK. Analyzed the data: AS ALB DN SG. Contributed reagents/materials/ analysis tools: NM ALB DZ SG. Contributed to the writing of the manuscript: AS ALB SG.

20. Frothingham R, Meeker-O'Connell WA (1998) Genetic diversity in the Mycobacterium tuberculosis complex based on variable numbers of tandem DNA repeats. Microbiology 144 (Pt 5): 1189-1196.

21. Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, et al. (2000) Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome. Mol Microbiol 36: 762-771.

22. Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, et al. (2006) Proposal for Standardization of Optimized Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeat Typing of Mycobacterium tuberculosis V. J Clin Microbiol 44: 4498-4510.

23. Goudet J, Buchi L (2006) The effects of dominance, regular inbreeding and sampling design on QST, an estimator of population differentiation for quantitative traits. Genetics 172: 1337-1347.

24. Felsenstein J (1985) Phylogenies and the Comparative Method. The American Naturalist 125, No 1: 1-15.

25. Chevenet F, Brun C, Banuls AL, Jacq B, Christen R (2006) TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics 7: 439.

26. Muller B, Steiner B, Bonfoh B, Fane A, Smith NH, et al. (2008) Molecular characterization of Mycobacterium bovis isolated from cattle slaughtered at the Bamako abattoir in Mali. BMC Vet Res 4: 26.

27. Smith NH, Berg S, Dale J, Allen A, Rodriguez S, et al. (2011) European 1: a globally important clonal complex of Mycobacterium bovis. Infect Genet Evol 11: 1340-1351.

28. Berg S, Garcia-Pelayo MC, Muller B, Hailu E, Asimwe B, et al. (2011) African 2, a clonal complex of Mycobacterium bovis epidemiologically important in East Africa. J Bacteriol 193: 670-678.

29. Asseged B, Woldesenbet Z, Yimer E, Lemma E (2004) Evaluation of abattoir inspection for the diagnosis of Mycobacterium bovis infection in cattle at Addis Ababa abattoir. Trop Anim Health Prod 36: 537-546.

30. Mekibeb A, Fulasa TT, Firdessa R, Hailu E (2013) Prevalence study on bovine tuberculosis and molecular characterization of its causative agents in cattle slaughtered at Addis Ababa municipal abattoir, Central Ethiopia. Trop Anim Health Prod 45: 763-769.

31. Tekul A, Asseged B, Yimer E, Gebeyehu M, Woldesenbet Z (2004) Tuberculous lesions not detected by routine abattoir inspection: the experience of the Hossana municipal abattoir, southern Ethiopia. Rev Sci Tech 23: 957-964.

32. Boukary AR, Thys E, Rigouts L, Matthys F, Berkvens D, et al. (2012) Risk factors associated with bovine tuberculosis and molecular characterization of Mycobacterium bovis strains in urban settings in Niger. Transbound Emerg Dis 59: 490-502.

33. Haddad N, Ostyn A, Karoui C, Masselot M, Thorel MF, et al. (2001) Spoligotype diversity of Mycobacterium bovis strains isolated in France from 1979 to 2000. J Clin Microbiol 39: 3623-3632.

34. Sahraoui N, Muller B, Guetarni D, Boulahbal F, Yala D, et al. (2009) Molecular characterization of Mycobacterium bovis strains isolated from cattle slaughtered at two abattoirs in Algeria. BMC Vet Res 5: 4.

35. Jenkins AO, Cadmus SI, Venter EH, Pourcel C, Hauk Y, et al. (2011) Molecular epidemiology of human and animal tuberculosis in Ibadan, Southwestern Nigeria. Vet Microbiol 151: 139-147.

36. Evans JT, Smith EG, Banerjee A, Smith RM, Dale J, et al. (2007) Cluster of human tuberculosis caused by Mycobacterium bovis: evidence for person-to-person transmission in the UK. Lancet 369: 1270-1276.

37. Godreuil S, Jeziorski E, Banuls AL, Fraisse T, Van de Perre P, et al. (2010) Intrafamilial cluster of pulmonary tuberculosis due to Mycobacterium bovis of the African 1 clonal complex. J Clin Microbiol 48: 4680-4683.