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

Mycobacterium bovis in Burkina Faso: Epidemiologic and Genetic Links between Human and Cattle Isolates Adama Sanou 1,2 *, Zekiba Tarnagda 3,4, Estelle Kanyala 4, Dezemon Zingué 2, Moumini Nouctara 2, Zakaria Ganamé 2, Adjima Combary 5, Hervé Hien 2,3, Mathurin Dembele 5, Antoinette Kabore 2, Nicolas Meda 2,4, Philippe Van de Perre 6,7,8, Dorine Neveu 7,8, Anne Laure Bañuls 1., Sylvain Godreuil 6,7,8. 1 MIVEGEC, UMRIRD 224-CNRS 5290-Montpellier Universités 1 et 2, Montpellier, France, 2 CentreMuraz, Bobo-Dioulasso,Burkina Faso,3 Institut de Recherche en Sciences de la Santé (IRSS), DirectionRégionale de l Ouest (DRO),Bobo Dioulasso, Burkina Faso,4 Université de Ouagadougou, Ouagadougou, Burkina Faso, 5 Programme NationalTuberculose, 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 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. * Email: sanouadama2003@yahoo.fr. These authors contributed equally to this work. 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 PLOS Neglected Tropical Diseases www.plosntds.org 1 October 2014 Volume 8 Issue 10 e3142

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 veterinarian 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 postmortem 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 4uC 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 SpolDB4.0[18] andhttp://www.mbovis.org [19].Thenew 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 (H s ) were calculated. The population structure was explored by estimating the F st (index of genetic differentiation between samples) value (0 = no differentiation and 1 = fixation of alternative alleles). The allelic diversity, the H s and the F st were calculated using F-STAT, version 2.9.3 with the 26 MIRU-VNTR loci [23]. PLOS Neglected Tropical Diseases www.plosntds.org 2 October 2014 Volume 8 Issue 10 e3142

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 nu 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 singlem. 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 RDAf1 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 RDAf1 region. This group of Table 1. Socio-demographic informations about hosts and M. bovis isolates. Strain ID Host Geographical location HIV status a Year of isolation Number c1 c7 cattle Bobo-Dioulasso 2011 7 b c8 c26 cattle Ouagadougou 2011 19 b h1 human Solenzo HIV-1+ 2009 1 c h2 human Ouagadougou unknown 2010 1 c h3 human Ouagadougou HIV- 2010 1 c h4 human Koupéla HIV- 2011 1 c h5 human Bobo-Dioulasso HIV- 2011 1 c h6 human Bobo-Dioulasso HIV- 2011 1 d h7 human Bobo-Dioulasso HIV- 2013 1 d a specific to human hosts, b bovine study, c nationwide survey, d regional study. doi:10.1371/journal.pntd.0003142.t001 PLOS Neglected Tropical Diseases www.plosntds.org 3 October 2014 Volume 8 Issue 10 e3142

Table 2. Spoligotypes, MIRU-VNTR patterns and clonal complex identification of the M. bovis strains isolated in Burkina Faso. MIRU-VNTR (patterns, codes and SB number Spoligotype signature 1 number of strains) 2 Strains Clone [no (%)] 3 SB2282 &%%%%&&&%&&&&&&%&&&%&&& &&&&&&&&&&&%%&&%%%%% 4253311356122324273 22225432(P1,n=1) c7 Af5 [1 (3%)] SB2283 &%%%%&&&%&&&&&&%&&&&&&& &&&&&&&&&&&%&&&%%%%% 7252312311712232621 322235332(P2,n=1) h1 Af5 [1 (3%)] SB1398 &&%%%&&&%&&&&&&%&&&&&&& &&&&&&&&&&&%%&&%%%%% 7253311310512232421 322235332(P3,n=2) c4, c5 Af5 [2 (6%)] SB2284 &%%%%&&&%&&&&&%%&&&&&&& &&&&&&&&&&&%%&&%%%%% NA (P4, n = 1) h7 Af5 [1 (3%)] SB2285 &%%%%&&&%&&&&&&%&&&&&&& &&&&&&&&&&&%%&&%%%%% 7253311355122326213 22235332(P5,n=1) c2 Af5 [1 (3%)] SB0857 &&%&&&&&%&&&&&&%&&&&&&& &&&&&&%%%%%%%&&%%%%% 4553311312122324253 22235432(P6,n=1) c15 Af1 [1 (3%)] SB0300 &&%&&%&&%&&&&&&%&&&&&&& 545329214122324263 22235432(P7,n=2) h5, h6 Af1 [2 (6%)] 545339214122324263 22235422(P8,n=1) c24 Af1 [1 (3%)] SB2286 &&%&&%&&%&&&&%&%&&&&&&& 555339214122224253 22235432(P9,n=2) c13, c17 Af1 [2 (6%)] 555331021412222425 322235432(P10,n=1) c12 Af1 [1 (3%)] SB2287 &&%&&%%&%&&&&&&%&&&&&&& &&&&&&%&&&&%%&&%%%%% 455331031312232424 322235432(P11,n=1) c25 Af1 [1 (3%)] SB0944 &&%&&&&&%&&&&&&%&&&&&&& 345331130312232424 322215432(P12,n=1) h2 Af1 [17 (52%)] 365338313122324233 22235432(P13,n=1) c3 455331131212232423 322245412(P14,n=1) c9 455331131212232425 322225432(P15,n=1) c16 453331131222324253 22235432(P16,n=1) c8 455331131312232424 322234432(P18,n=2) c21, c20 455331131212232424 322235432(P19,n=1) c18 455331131312232424 322235432(P20,n=2) h4, c19 452331131312232424 322235432(P21,n=1) c14 PLOS Neglected Tropical Diseases www.plosntds.org 4 October 2014 Volume 8 Issue 10 e3142

Table 2. Cont. MIRU-VNTR (patterns, codes and SB number Spoligotype signature 1 number of strains) 2 Strains Clone [no (%)] 3 355331131312232424 322235432(P22,n=2) c22, c23 365331131312232424 322235432(P23,n=1) c6 455331031322324243 22235332(P24,n=1) c26 455331031202232424 322235432(P25,n=2) c1, c11 SB1439 &&%&&&&&%&&&&&%%&&&&&&& 455331031202232424 322235432(P25,n=1) c10 Af1 [1 (3%)] SB2288 &&%&&&&&%&&&&&&%&&%&&&& 455331131412232424 322235432(P17,n=1) h3 Af1 [1 (3%)] Total 12 Spoligotype signatures 25 MIRU-VNTR patterns 33 (100%) Af5 [6 (18.2%)] Af1 [27 (81.8%)] 1 &, presence of spacer; %, absence of spacer. 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, Mtub 04, Mtub 21, Mtub 29, Mtub 30, Mtub 34, Mtub 39. NA = Not Available. 3 Af1 = African 1 clonal complex, Af5 = putative African 5 clonal complex. doi:10.1371/journal.pntd.0003142.t002 PLOS Neglected Tropical Diseases www.plosntds.org 5 October 2014 Volume 8 Issue 10 e3142

Figure 1. UPGMA tree based on the MIRU-VNTR(26 loci) and spoligotyping data. 1, SB number = name of spoligotype based on http://www. M.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 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 (Gd) 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 PLOS Neglected Tropical Diseases www.plosntds.org 6 October 2014 Volume 8 Issue 10 e3142

Table 3. Allelic diversity of the 26 MIRU-VNTR loci in M. bovis isolates from humans and livestock in Burkina Faso.* Locus Number of alleles Allelic diversity Global (n = 32) Af1 (n = 27) Af5 (n = 5) Global (n = 32) Af1 (n = 27) Af5 (n = 5) ETR A 4 3 2 0.66 0.6 0.4 ETR B 4 4 2 0.54 0.44 0.4 ETR C 3 3 1 0.12 0.15 0 ETR D 2 1 2 0.06 0 0.4 ETR E 2 2 1 0.12 0.14 0 QUB-11a 5 4 2 0.60 0.63 0.4 QUB-11b 2 2 1 0.31 0.36 0 QUB-3232 5 3 4 0.34 0.15 0.9 QUB-26 6 4 3 0.76 0.70 0.7 QUB-4156 2 2 1 0.20 0.22 0 MIRU 2 1 1 1 0 0 0 MIRU 10 1 1 1 0 0 0 MIRU16 2 2 1 0.18 0.2 0 MIRU 20 1 1 1 0 0 0 MIRU 23 2 1 2 0.12 0 0.6 MIRU 24 1 1 1 0 0 0 MIRU 26 6 5 2 0.71 0.62 0.4 MIRU 27 1 1 1 0 0 0 MIRU 39 1 1 1 0 0 0 MIRU 40 1 1 1 0 0 0 Mtub 04 1 1 1 0 0 0 Mtub 21 4 4 2 0.24 0.21 0.4 Mtub 29 2 2 1 0.12 0.14 0 Mtub 30 2 2 2 0.27 0.07 0.4 Mtub 34 3 3 1 0.12 0.14 0 Mtub 39 1 1 1 0 0 0 Mean 2.5 2.15 1.5 0.21 0.18 0.19 *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) PLOS Neglected Tropical Diseases www.plosntds.org 7 October 2014 Volume 8 Issue 10 e3142

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 Af1 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 Af1 clonal complex and a minor group belonging to the putative Af5 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, Zinsstag 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 cattlebreeding 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: technics 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, Cauchoix 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, Prodinger 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,{. J Clin Microbiol 44: 4498 4510. 23. Goudet J, Buchi L (2006) The effects of dominance, regular inbreeding and sampling design on Q(ST), 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 characterisation 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, Asiimwe 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. Teklul 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-toperson 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. PLOS Neglected Tropical Diseases www.plosntds.org 8 October 2014 Volume 8 Issue 10 e3142