# Transmission network of deer-borne Mycobacterium bovis infection revealed by a WGS approach.

Lorraine Michelet, Cyril Condé, Maxime Branger, Thierry Cochard, Franck

Biet, Maria Laura Boschiroli

► To cite this version:

Lorraine Michelet, Cyril Condé, Maxime Branger, Thierry Cochard, Franck Biet, et al. Transmission network of deer-borne Mycobacterium bovis infection revealed by a WGS approach. Microorganisms, 2019, 7 (12), 11 p. 10.3390/microorganisms7120687. [https://hal.inrae.fr/hal-02622745](https://hal.inrae.fr/hal-02622745)

HAL Id: hal-02622745

https://hal.inrae.fr/hal-02622745

Submitted on 26 May 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
<figure><img src="image_2.png" /></figure>

Article

# Transmission Network of Deer-Borne *Mycobacterium bovis* Infection Revealed by a WGS Approach

Lorraine Michelet<fnref n="1" />,*, Cyril Conde<fnref n="2" />, Maxime Branger<fnref n="2" />,*, Thierry Cochard<fnref n="2" />, Franck Biet<fnref n="2" /> and Maria Laura Boschiroli<fnref n="1" />

Received: 28 October 2019; Accepted: 9 December 2019; Published: 12 December 2019

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

check for
updates

**Abstract:** Bovine tuberculosis (TB) is a zoonotic disease, mainly caused by *Mycobacterium bovis*. France was declared officially TB free in 2001, however, the disease persists in livestock and wildlife. Among wild animals, deer are particularly susceptible to bovine TB. Here, a whole genome sequence (WGS) analysis was performed on strains with the same genetic profile—spoligotype SB0121, Multiple Loci VNTR Analysis (MLVA) 6 4 5 3 11 2 5 7—isolated from different types of outbreaks, including from deer or cattle herds, or zoological or hunting parks where the presence of infected deer was a common trait in most of them. The results of the phylogeny based on the SNP calling shows that two sub-clusters co-exist in France, one related to deer bred to be raised as livestock, and the other to hunting parks and zoos. The persistence over almost 30 years of sporadic cases due to strains belonging to these clusters highlights the deficiency in the surveillance of captive wildlife and the need for better monitoring of animals, especially before movement between parks or herds.

**Keywords:** bovine tuberculosis; *Mycobacterium bovis*; deer herd; animal park; WGS; SNP calling

## 1. Introduction

Bovine tuberculosis (bTB), mainly caused by *Mycobacterium bovis*, is an important re-emergent zoonotic disease in Europe [1]. Although France has been officially bTB free (OTF) since 2001, the persistence of the disease in livestock and its occurrence in some areas of wildlife are matters of great concern [2]. This ancient disease is no longer considered just a cattle-borne problem, but a concern for multihost communities that include wildlife species such as wild boar (*Sus scrofa*), red deer (*Cervus elaphus*), and badgers (*Meles meles*) [3].

In Europe, red and fallow deer (*Dama dama*) are frequently reported to be infected by *M. bovis* or *M. caprae* [4]. In the wild population, red deer are considered to play an important role in the dissemination of the disease as an important number of infected animals disclose generalized patterns of clinical tuberculosis [5]. Farmed deer are also very susceptible to *M. bovis*, as a high incidence of animals with generalized infection is commonly observed in deer herds [6]. Gross lesions are generally observed in lymph nodes of the head, thorax and abdomen [4]; most frequently in the retropharyngeal lymph nodes for red deer, and the thoracic region for fallow deer [7].

Interaction among different wildlife species, and between wildlife and human beings, may now be more frequent due to the increasing amount of captive wildlife in zoos or other types of animal parks [8]. The zoonotic transmission of *Mycobacterium tuberculosis* complex (MTBC) pathogens has

<footnote n="1">¹ Paris-Est University, National Reference Laboratory for Tuberculosis, Animal Health Laboratory, French Agency for Food, Environmental and Occupational Health and Safety (Anses), 94701 Maisons-Alfort CEDEX, France; lorraine.michelet@anses.fr</footnote>

<footnote n="2">² Infectiologie et Santé Publique (ISP), Institut National de Technologie Agronomique (INRA), Université de Tours, UMR 1282, 37380 Nouzilly, France; cyril.conde@inra.fr (C.C.); maxime.branger@inra.fr (M.B.); thierry.cochard@inra.fr (T.C.); franck.biet@inra.fr (F.B.)</footnote>

<footnote>* Correspondence: maria-laura.boschiroli@anses.fr; Tel.: +33149771321</footnote>
been reported worldwide: e.g., an *M. bovis* infection in zookeepers exposed to a white rhinoceros [9], an *M. pinnipedii* infection in animal handlers linked to a sea lion colony [10], and an *M. tuberculosis* infection of staff members exposed to elephants [11,12]. Infrastructure where the public has potential contact with captive wildlife is a serious emerging concern for zoonotic transmission, especially for immunocompromised individuals [8]. The same concern arose after the discovery in an Iranian zoo of an *M. bovis* infection in free-roaming fallow deer that had been in close contact with the public [13].

MTBC members are excellent microorganisms for conducting molecular epidemiology studies on, by virtue of their low mutation rate nature. In the case of bTB, two genotyping methods are widely used—spoligotyping [14] and Multiple Loci VNTR Analysis (MLVA) typing [15]—both of which establish differentiation between *M. bovis* strains, taking into account relative genetic stability. The combination of these two techniques allows for a very fine differentiation of strains. In France, 497 different genotypes could be identified between 1978 and 2013, in a collection of *M. bovis* strains isolated from cattle and wild animals [2,16]. The vast diversity of the molecular profiles obtained by combining these two techniques makes it possible, by comparing the profiles of isolated strains of cattle and wild animals, to determine the origin of the infection of a very large number of outbreaks and to highlight possible inter-species transmission [17]. However, in areas where the disease is still highly prevalent, unique dominant genetic profiles-per-zone are shared by almost all isolates, which makes the reconstitution of a transmission chain impossible [17,18]. Molecular techniques with finer resolutions are necessary in such contexts. Analyses of mycobacterial whole genome sequences (WGS), which detect genomic changes at a very small scale, hold promise as a strain discrimination alternative [19]. This technique has particularly been used to better understand the human-to-human transmission of different mycobacteria, such as *M. leprae* [20] or *M. tuberculosis* [21]. WGS data have also been used to reconstruct transmission scenarios of *M. bovis* infection [22,23].

The aim of this study was to use WGS analysis in order to investigate the transmission network of *M. bovis* with a particular genotype, mainly identified in infected red deer belonging to French farms, zoos, or hunting or other type of parks.

## 2. Material and Methods

### 2.1. Ethical Statement

In France, bTB is a notifiable disease with an eradication program in cattle and surveillance in free range wildlife. The official methods for the diagnosis of this disease are culture, PCR and histopathology. Therefore, all the samples included in this study are issued from animals analyzed within an official context. No purposeful killing of animals was performed for this study. All samplings were in complete agreement with national and European regulations. No ethical approval was necessary.

### 2.2. Bacterial Strains

This study included 24 field *M. bovis* isolates belonging to the French bTB National Reference Laboratory (NRL, Anses) collection (Table 1). The source population was animals whose samples were analyzed because (i) they presented non-negative skin tests and/or γ-interferon tests results and were culled for diagnostic purposes; (ii) they were culled as a result of the total or partial slaughter of their herd of origin after confirmation of the herd's infection, or (iii) they presented macroscopic bTB-like lesions at routine abattoir inspection. Bacterial culture was performed following the protocol established by the French NRL (NF U 47-104) for the isolation of *M. bovis*. The identity of the isolated *Mycobacterium tuberculosis* complex colonies was confirmed by DNA amplification as described by Hénault et al. [24], and *M. bovis* was confirmed by spoligotyping and VNTR typing as described below. For livestock breakdowns, a herd-based criterion was used for strain choice: at least 1 strain/outbreak was selected. Spoligotyping was performed as described by Zhang et al. [14], using TB-SPOL kits purchased from Beamedex® (Beamedex® SAS, Orsay, France) on BioPLex200/Luminex 200® as described by Hauer et al. [2]. Spoligotypes have been named according to an agreed international convention
(www.mbovis.org). A Multi-Locus Variable numbers tandem repeats Analysis (MLVA) profile identification [25] was performed by Genoscreen (Lille, France), using PCR amplification targeting genetic loci including mycobacterial MIRU-VNTR. Analysis was based on 8 loci, ETR A (VNTR2165), ETR B (VNTR2461), ETR C (VNTR577), ETR D (MIRU4 or VNTR580), QUB 11a (VNTR2163a), QUB 11b (VNTR2163b), QUB 26 (VNTR4052) and QUB 3232 (VNTR3232), chosen in the framework of a European consortium based on their degree of polymorphism and their ability to discriminate local strains [26].
<table><caption>**Table 1.** Description of the 24 sequenced *M. bovis* strains.</caption><thead><tr><th></th><th></th><th>ID</th><th>Species</th><th>Year of Isolation</th><th>Context</th><th>Spoligotype</th><th>VNTR Profile</th><th>Accession Number</th></tr></thead><tbody><tr><td>1</td><td>1979-02442</td><td>S79-020</td><td>Wild boar</td><td>1979</td><td>Park</td><td>GB54/SB0121</td><td>6 3 5 3 11 1 5 7</td><td>ERS3860435</td></tr><tr><td>2</td><td>1990-08827</td><td>F90-035</td><td>Red deer</td><td>1990</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860436</td></tr><tr><td>3</td><td>1990-09372</td><td>S90-020</td><td>Cattle</td><td>1990</td><td>French herd</td><td>F70/SB0295</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860437</td></tr><tr><td>4</td><td>1993-04379</td><td>F93-039</td><td>Cattle</td><td>1993</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860438</td></tr><tr><td>5</td><td>1993-06259</td><td>S93-007</td><td>Bison</td><td>1993</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860439</td></tr><tr><td>6</td><td>1993-11575</td><td>F93-076</td><td>Red deer</td><td>1993</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860440</td></tr><tr><td>7</td><td>1993-11583</td><td>S93-013</td><td>Monkey</td><td>1993</td><td>Park</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860441</td></tr><tr><td>8</td><td>1995-00583</td><td>S95-007</td><td>Cattle</td><td>1995</td><td>Foreign-born animal</td><td>F070/SB0295</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860442</td></tr><tr><td>9</td><td>1996-07003</td><td>S96-046</td><td>Cattle</td><td>1996</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860443</td></tr><tr><td>10</td><td>1996-08822</td><td>S96-050</td><td>Cattle</td><td>1996</td><td>Foreign-born animal</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860444</td></tr><tr><td>11</td><td>1997-00766</td><td>S97-004</td><td>Cattle</td><td>1997</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860445</td></tr><tr><td>12</td><td>1998-12646</td><td>S98-099</td><td>Bison</td><td>1998</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860446</td></tr><tr><td>13</td><td>00-06572</td><td>S00-206</td><td>Cattle</td><td>2000</td><td>Foreign-born animal</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860447</td></tr><tr><td>14</td><td>07-03151</td><td>S07-524</td><td>Wild boar</td><td>2007</td><td>Park</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860448</td></tr><tr><td>15</td><td>08-00202</td><td>S08-044</td><td>Red deer</td><td>2008</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860449</td></tr><tr><td>16</td><td>08-02744</td><td>S08-481</td><td>Cattle</td><td>2008</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860450</td></tr><tr><td>17</td><td>D11-02906</td><td>S11-496</td><td>Cattle</td><td>2011</td><td>Foreign-born animal</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860451</td></tr><tr><td>18</td><td>D12-01540</td><td>S12-164</td><td>Wild boar</td><td>2012</td><td>Park</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860452</td></tr><tr><td>19</td><td>D12-02012</td><td>S12-295</td><td>Cattle</td><td>2012</td><td>Foreign-born animal</td><td>F070/SB0295</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860453</td></tr><tr><td>20</td><td>D14-00543</td><td>S14-130</td><td>Cattle</td><td>2014</td><td>Foreign-born animal</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860454</td></tr><tr><td>21</td><td>D14-02802</td><td>S14-680</td><td>Cattle</td><td>2014</td><td>French herd</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860455</td></tr><tr><td>22</td><td>D17-02544</td><td>S17-321</td><td>Fallow deer</td><td>2017</td><td>Park</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860456</td></tr><tr><td>23</td><td>D17-03607</td><td>S17-426</td><td>Wapiti</td><td>2017</td><td>Park</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860457</td></tr><tr><td>24</td><td>D18-00301</td><td>S18-037</td><td>Bison</td><td>2018</td><td>Park</td><td>GB54/SB0121</td><td>6 4 5 3 11 2 5 7</td><td>ERS3860458</td></tr></tbody></table>
## 2.3. Whole Genome Sequencing and Variant Calling

Whole genome paired end (2 × 250 bp) sequencing was performed on Rapid Run HiSeq 2500 (Illumina) using NextEra XT libraries (Illumina, San Diego, CA, United States), and indexed following the manufacturer's recommendations (Illumina, San Diego, CA, United States) by Genoscreen (Lille, France). The theoretical coverage was expected to be higher than 50×. Raw sequences were aligned to the reference genome AF2122/97 (Genbank accession NC_002945.4) using BioNumerics v7.6 (Applied Maths). A wgSNP analysis was performed with the strict SNP filtering template, i.e., removing positions with at least one unreliable base, ambiguous base, or gap; removing nondiscriminatory positions; with a minimum of 10× coverage; with a minimum distance of 12 nucleotides between SNPs. The evolutionary history of this dataset was inferred on MEGA7 software [27], using the Maximum Likelihood method based on the General Time Reversible (GTR) model (model selected based on jmodeltest v2.1.10 [28]). A bootstrap analysis was performed with 500 replicates.

WGS data were submitted to the EBI/EMBL. The accession Numbers are listed in Table 1.

## 3. Results

Among the 86 strains available at the NRL with the genotype of interest, i.e., a spoligotype profile belonging to the European 2 complex (either GB54/SB0121 or F70/SB0295) and a unique VNTR profile (6 4 5 3 11 2 5 7), 23 strains were selected from different herds or other origins, to represent all potential outbreaks (Table 1). They were isolated between 1979 and 2018 from 12 cattle (*Bos taurus*), 3 bison (*Bison bonasus*), 3 wild boar (*Sus scrofa*), 3 red deer (*Cervus elaphus*), 1 fallow deer (*Dama dama*), 1 wapiti (*Cervus canadensis*) and 1 nonhuman primate (species not specified). The strain isolated in 1979 from a wild boar was selected as an outgroup, as it shared the same spoligotype but a slightly different VNTR profile (6 3 5 3 11 1 5 7).

These isolates were analyzed by WGS. In total, we identified 770 SNPs compared to the *M. bovis* AF2122/97 reference genome (supplementary Table S1). The phylogenetic analysis highlighted different clusters, all enrooted by the wild boar strain isolated in 1979 (Figures 1 and 2). A distance matrix of SNPs was generated using MEGA (Table 2) and the number of SNPs between each strain was determined. A first cluster supported by a bootstrap value of 100% is defined by 47 SNPs. This cluster is sub-divided into two sub-clusters: the first one, which is phylogenetically well supported (cluster 1.1, 98%), containing strains isolated from the cattle or deer of livestock herds; the second one, which is less robust (cluster 1.2, 67%), contains strains isolated from zoo animals or belonging to hunting parks (Figure 2). A low number of SNPs defined these sub-clusters—four and one respectively. Cluster 2 groups all the strains with spoligotype F70/SB0295, isolated from one French cow and from foreign-born animals, and is defined by 24 SNPs. Cluster 3, defined by 51 SNPs, contains strains isolated from foreign-born animals with spoligotype GB54/SB0121.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Molecular Phylogenetic tree of the 24 strains using the Maximum Likelihood method. The evolutionary history was inferred on MEGA by using the Maximum Likelihood method based on the General Time Reversible model on a final dataset of 770 SNP positions for 24 strains [1]. The number above each branch corresponds to the number of SNPs. Bootstrap values (%) based on 500 replications are indicated in italics after a slash if greater than 75%. A violet square represents strains isolated from French herds. A red square represents strains isolated from foreign-born animals. A green square represents strains isolated from park animals.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Zoomed-in view of the phylogenetic tree of cluster 1. The number above each branch corresponds to the number of SNPs. Bootstrap values (%) based on 500 replications are indicated in italics after a slash if greater than 75%. A violet square represents strains isolated from French herds. A green square represents strains isolated from park animals.</figcaption></figure>

The 16 French strains gathered in cluster 1 seem to be closer (between 0 and 46 SNPs) than the strains of cluster 2 and 3. Indeed, some of them have very few genetic differences. Two strains isolated from a bison and a wapiti from the same animal park are identical. The WGS analysis allows us to confirm a link between two strains isolated from one wild boar in 2007, another in 2012 and a strain isolated from a fallow deer in 2017. The total number of SNPs for cluster 1 is 107, 50 for cluster 1.1, and 57 for cluster 1.2. On the other hand, the number of differences in the remaining clusters is higher with fewer strains, i.e., 163 SNPs for cluster 2 and 199 SNPs for cluster 3.
<table><caption>**Table 2.** Estimates of Evolutionary Divergence between Sequences. The number of base differences between sequences are shown.</caption><thead><tr><th></th><th>F90-035</th><th>F93-039</th><th>S93-007</th><th>F93-076</th><th>S96-046</th><th>S97-004</th><th>S98-099</th><th>S08-044</th><th>S08-481</th><th>S14-680</th><th>S93-013</th><th>S07-524</th><th>S12-164</th><th>S17-321</th><th>S17-426</th><th>S18-037</th><th>S90-020</th><th>S95-007</th><th>S12-295</th><th>S96-050</th><th>S00-206</th><th>S11-496</th><th>S14-130</th></tr></thead><tbody><tr><td>Cluster 1.1</td><td>F90-035</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>F93-039</td><td>4</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>S93-007</td><td>2</td><td>6</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>F93-076</td><td>10</td><td>10</td><td>12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>S96-046</td><td>5</td><td>3</td><td>7</td><td>11</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>S97-004</td><td>6</td><td>4</td><td>8</td><td>12</td><td>3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>S98-099</td><td>8</td><td>8</td><td>10</td><td>14</td><td>9</td><td>10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>S08-044</td><td>7</td><td>11</td><td>9</td><td>17</td><td>12</td><td>13</td><td>15</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>S08-481</td><td>11</td><td>15</td><td>13</td><td>21</td><td>16</td><td>17</td><td>19</td><td>12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.1</td><td>S14-680</td><td>11</td><td>11</td><td>13</td><td>13</td><td>12</td><td>13</td><td>15</td><td>18</td><td>22</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.2</td><td>S93-013</td><td>23</td><td>23</td><td>25</td><td>29</td><td>24</td><td>25</td><td>17</td><td>30</td><td>34</td><td>30</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.2</td><td>S07-524</td><td>27</td><td>27</td><td>29</td><td>33</td><td>28</td><td>29</td><td>21</td><td>34</td><td>38</td><td>34</td><td>26</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.2</td><td>S12-164</td><td>28</td><td>28</td><td>30</td><td>34</td><td>29</td><td>30</td><td>22</td><td>35</td><td>39</td><td>35</td><td>27</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.2</td><td>S17-321</td><td>32</td><td>32</td><td>34</td><td>38</td><td>33</td><td>34</td><td>26</td><td>39</td><td>43</td><td>39</td><td>31</td><td>9</td><td>10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.2</td><td>S17-426</td><td>35</td><td>35</td><td>37</td><td>41</td><td>36</td><td>37</td><td>29</td><td>42</td><td>46</td><td>42</td><td>34</td><td>38</td><td>39</td><td>43</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 1.2</td><td>S18-037</td><td>35</td><td>35</td><td>37</td><td>41</td><td>36</td><td>37</td><td>29</td><td>42</td><td>46</td><td>42</td><td>34</td><td>38</td><td>39</td><td>43</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 2</td><td>S90-020</td><td>134</td><td>134</td><td>136</td><td>140</td><td>135</td><td>136</td><td>128</td><td>141</td><td>145</td><td>141</td><td>135</td><td>139</td><td>140</td><td>144</td><td>147</td><td>147</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 2</td><td>S95-007</td><td>139</td><td>139</td><td>141</td><td>145</td><td>140</td><td>141</td><td>133</td><td>146</td><td>150</td><td>146</td><td>140</td><td>144</td><td>145</td><td>149</td><td>152</td><td>152</td><td>87</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 2</td><td>S12-295</td><td>145</td><td>145</td><td>147</td><td>151</td><td>146</td><td>147</td><td>139</td><td>152</td><td>156</td><td>152</td><td>146</td><td>150</td><td>151</td><td>155</td><td>158</td><td>158</td><td>93</td><td>98</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 3</td><td>S96-050</td><td>151</td><td>151</td><td>153</td><td>157</td><td>152</td><td>153</td><td>145</td><td>158</td><td>162</td><td>158</td><td>152</td><td>156</td><td>157</td><td>161</td><td>164</td><td>164</td><td>147</td><td>152</td><td>158</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 3</td><td>S00-206</td><td>159</td><td>159</td><td>161</td><td>165</td><td>160</td><td>161</td><td>153</td><td>166</td><td>170</td><td>166</td><td>160</td><td>164</td><td>165</td><td>169</td><td>172</td><td>172</td><td>155</td><td>160</td><td>166</td><td>70</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 3</td><td>S11-496</td><td>164</td><td>164</td><td>166</td><td>170</td><td>165</td><td>166</td><td>158</td><td>171</td><td>175</td><td>171</td><td>165</td><td>169</td><td>170</td><td>174</td><td>177</td><td>177</td><td>160</td><td>165</td><td>171</td><td>65</td><td>83</td><td></td><td></td><td></td><td></td></tr><tr><td>Cluster 3</td><td>S14-130</td><td>159</td><td>159</td><td>161</td><td>165</td><td>160</td><td>161</td><td>153</td><td>166</td><td>170</td><td>166</td><td>160</td><td>164</td><td>165</td><td>169</td><td>172</td><td>172</td><td>155</td><td>160</td><td>166</td><td>68</td><td>78</td><td>81</td><td></td><td></td><td></td></tr><tr><td>Outgroup</td><td>S79-020</td><td>301</td><td>301</td><td>303</td><td>307</td><td>302</td><td>303</td><td>295</td><td>308</td><td>312</td><td>308</td><td>302</td><td>306</td><td>307</td><td>311</td><td>314</td><td>314</td><td>319</td><td>324</td><td>330</td><td>336</td><td>344</td><td>349</td><td>344</td><td></td><td></td></tr></tbody></table>
## 4. Discussion and Conclusions

This investigation, based on genomic data, allowed us to reveal the circulation of closely related *M. bovis* strains that seem to be confined to deer farms and animal parks. Although these types of strains have been regularly found since the 1990s, no geographical link seems to exist between most herd cases. WGS analysis showed a main cluster in France, divided into two very close sub-clusters (with just five SNP differences)—cluster 1.1 related to cattle and deer herds reared as livestock, and cluster 1.2 related to recreational or hunting animal parks. This analysis confirms the infection by the same strain of a bison and a wapiti from the same park. It also shows a link between the strain isolated recently in a fallow deer and two wild boar strains in 2007 and 2012 [29]. An epidemiological link exists between these three cases: the wild boars belonged to the same hunting park and the fallow deer, found infected but in another animal park, was originally issued from the same hunting park as that of the wild boar cases. Altogether, the common point in all these French cases was the presence of deer in the parks/herds or close to cattle herds. This genotype has been circulating in France for the last 30 years, however such strains have been isolated sporadically which clearly shows that there is a lack of appropriate surveillance in animal parks. Given that primates are not maintenance hosts of *M. bovis*, the case described in the monkey—acting as a sentinel of the infection—could be the result of the under-detection of another wild ungulate case in the park. The low number of screening tests, such as the skin test, conducted on deer is the result of several factors, such as the difficulty to perform the test on these wild animals, the potential cross-reaction with environmental mycobacteria, or the inherent technical variability of the test [30].

Our data highlight the limit of classical genotyping techniques that can be overcome by the use of whole genome SNP analysis. Strains that shared the same genotype (spoligotype-VNTR) from French and foreign-born animals were genetically different and very distant.

The use of WGS revealed an unsuspected transmission network between hunting parks and zoos that needs to be further explored. These results reinforce the suspicion of an infectious deer-breeding network that could be at the origin of the bacillus dissemination. Indeed, this transmission could have occurred through inter-breeding animal exchanges. Our results raise awareness about the infectious risk of wild animals in captivity and the consequences of trading non-appropriately-tested animals. P22 ELISA, a currently available, highly specific and sensitive new test, could contribute to a more accurate diagnosis of TB in this species [30].

This study raises the problem of the zoonotic transmission of the disease, especially for cases in animal parks where close contact between free-roaming animals and the public is allowed. The fact that deer often develop generalized TB, and thus that they probably excrete the bacillus at significant rates, increases the risk to humans. This study therefore shows the need for the strengthened monitoring of these animals with regular ante-mortem tests (as for livestock) before allowing animal movement between parks.

**Supplementary Materials:** Supplementary materials can be found at http://www.mdpi.com/2076-2607/7/12/687/s1.

**Author Contributions:** Conceptualization, M.L.B., F.B. and L.M.; Software, C.C., M.B., T.C., L.M.; Formal Analysis, C.C., M.B., T.C., L.M.; Investigation, M.L.B. and L.M.; Writing—Original Draft Preparation, M.L.B. and L.M.; Writing—Review & Editing, C.C., M.B., T.C., F.B.

**Funding:** This research received no external funding.

**Acknowledgments:** We thank Victoria Boschiroli for her useful comments on this article.

**Conflicts of Interest:** The authors declare no conflict of interest.
References

1. Muller, B.; Durr, S.; Alonso, S.; Hattendorf, J.; Laisse, C.J.; Parsons, S.D.; van Helden, P.D.; Zinsstag, J. Zoonotic *Mycobacterium bovis*-induced tuberculosis in humans. *Emerg. Infect. Dis.* **2013**, *19*, 899–908. [CrossRef] [PubMed]

2. Hauer, A.; De Cruz, K.; Cochard, T.; Godreuil, S.; Karoui, C.; Henault, S.; Bulach, T.; Banuls, A.L.; Biet, F.; Boschiroli, M.L. Genetic evolution of *Mycobacterium bovis* causing tuberculosis in livestock and wildlife in France since 1978. *PLoS ONE* **2015**, *10*, e0117103. [CrossRef] [PubMed]

3. Gortázar, C.; Delahay, R.J.; McDonald, R.A.; Boadella, M.; Wilson, G.J.; Gavier-Widen, D.; Acevedo, P. The status of tuberculosis in European wild mammals. *Mammal Rev.* **2012**, *42*, 193–206. [CrossRef]

4. Griffin, J.F.; Mackintosh, C.G. Tuberculosis in deer: Perceptions, problems and progress. *Vet. J.* **2000**, *160*, 202–219. [CrossRef] [PubMed]

5. Zanella, G.; Duvauchelle, A.; Hars, J.; Moutou, F.; Boschiroli, M.L.; Durand, B. Patterns of lesions of bovine tuberculosis in wild red deer and wild boar. *Vet. Rec.* **2008**, *163*, 43–47. [CrossRef] [PubMed]

6. Mackintosh, C.G.; de Lisle, G.W.; Collins, D.M.; Griffin, J.F. Mycobacterial diseases of deer. *N. Z. Vet. J.* **2004**, *52*, 163–174. [CrossRef]

7. Martin-Hernando, M.P.; Torres, M.J.; Aznar, J.; Negro, J.J.; Gandia, A.; Gortazar, C. Distribution of lesions in red and fallow deer naturally infected with *Mycobacterium bovis*. *J. Comp. Pathol.* **2010**, *142*, 43–50. [CrossRef]

8. Miller, M.; Olea-Popelka, F. One Health in the shrinking world: Experiences with tuberculosis at the human-livestock-wildlife interface. *Comp. Immunol. Microbiol. Infect. Dis.* **2013**, *36*, 263–268. [CrossRef]

9. Dalovisio, J.R.; Stetter, M.; Mikota-Wells, S. Rhinoceros' rhinorrhea: Cause of an outbreak of infection due to airborne *Mycobacterium bovis* in zookeepers. *Clin. Infect. Dis.* **1992**, *15*, 598–600. [CrossRef]

10. Kiers, A.; Klarenbeek, A.; Mendelts, B.; Van Soolingen, D.; Koeter, G. Transmission of *Mycobacterium pinnipedii* to humans in a zoo with marine mammals. *Int. J. Tuberc. Lung Dis.* **2008**, *12*, 1469–1473.

11. Montali, R.J.; Mikota, S.K.; Cheng, L.I. *Mycobacterium tuberculosis* in zoo and wildlife species. *Rev. Sci. Tech.* **2001**, *20*, 291–303. [CrossRef] [PubMed]

12. Murphree, R.; Warkentin, J.V.; Dunn, J.R.; Schaffner, W.; Jones, T.F. Elephant-to-human transmission of tuberculosis, 2009. *Emerg. Infect. Dis.* **2011**, *17*, 366–371. [CrossRef] [PubMed]

13. Gharib Mombeni, E.; Mosavari, N.; Moradi Gravand, M.; Amir Rezai, A.; Keshavarz, R.; Tadayon, K.; Bakhshi, R.; Behmanesh, R. First Report of *Mycobacterium bovis* Isolation from a European Fallow Deer (*Dama Dama*) in Iran. *Iran. J. Public Health* **2016**, *45*, 814–816. [PubMed]

14. Zhang, J.; Abadia, E.; Refregier, G.; Tafaj, S.; Boschiroli, M.L.; Guillard, B.; Andremont, A.; Ruimy, R.; Sola, C. *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.* **2010**, *59*, 285–294. [CrossRef] [PubMed]

15. Skuce, R.A.; McDowell, S.W.; Mallon, T.R.; Luke, B.; Breadon, E.L.; Lagan, P.L.; McCormick, C.M.; McBride, S.H.; Pollock, J.M. Discrimination of isolates of *Mycobacterium bovis* in Northern Ireland on the basis of variable numbers of tandem repeats (VNTRs). *Vet. Rec.* **2005**, *157*, 501–504. [CrossRef]

16. Hauer, A.; Michelet, L.; De Cruz, K.; Cochard, T.; Branger, M.; Karoui, C.; Henault, S.; Biet, F.; Boschiroli, M.L. MIRU-VNTR allelic variability depends on *Mycobacterium bovis* clonal group identity. *Infect. Genet. Evol.* **2016**, *45*, 165–169. [CrossRef]

17. Boschiroli, M.L.; Michelet, L.; Hauer, A.; De Cruz, K.; Courcoul, A.; Hénault, S.; Palisson, A.; Karoui, C.; Biet, F.; Zanella, G. Tuberculose bovine en France: Cartographie des souches de *Mycobacterium bovis* entre 2000–2013. *Bull. Epid. St. Anim. Alim.* **2015**, *70*, 2–8.

18. Cavalerie, L.; Courcoul, A.; Boschiroli, M.L.; Réveillaud, E.; Gay, P. Bovine Tuberculosis in France in 2014: A stable situation. *Bull. Epid. St. Anim. Alim.* **2015**, *71*, 4–11.

19. Hauer, A.; Michelet, L.; Cochard, T.; Branger, M.; Nunez, J.; Boschiroli, M.L.; Biet, F. Accurate Phylogenetic Relationships Among *Mycobacterium bovis* Strains Circulating in France Based on Whole Genome Sequencing and Single Nucleotide Polymorphism Analysis. *Front. Microbiol.* **2019**, *10*, 955. [CrossRef]

20. Turankar, R.P.; Lavania, M.; Chaitanya, V.S.; Sengupta, U.; Darlong, J.; Darlong, F.; Siva Sai, K.S.; Jadhav, R.S. Single nucleotide polymorphism-based molecular typing of *M. leprae* from multicase families of leprosy
patients and their surroundings to understand the transmission of leprosy. *Clin. Microbiol. Infect.* **2014**, *20*, O142–O149. [CrossRef]

21. Roetzer, A.; Diel, R.; Kohl, T.A.; Ruckert, C.; Nubel, U.; Blom, J.; Wirth, T.; Jaenicke, S.; Schuback, S.; Rusch-Gerdes, S.; et al. Whole genome sequencing versus traditional genotyping for investigation of a *Mycobacterium tuberculosis* outbreak: A longitudinal molecular epidemiological study. *PLoS Med.* **2013**, *10*, e1001387. [CrossRef] [PubMed]

22. Biek, R.; O'Hare, A.; Wright, D.; Mallon, T.; McCormick, C.; Orton, R.J.; McDowell, S.; Trewby, H.; Skuce, R.A.; Kao, R.R. Whole genome sequencing reveals local transmission patterns of *Mycobacterium bovis* in sympatric cattle and badger populations. *PLoS Pathog.* **2012**, *8*, e1003008. [CrossRef] [PubMed]

23. Trewby, H.; Wright, D.; Breadon, E.L.; Lycett, S.J.; Mallon, T.R.; McCormick, C.; Johnson, P.; Orton, R.J.; Allen, A.R.; Galbraith, J.; et al. Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis. *Epidemics* **2016**, *14*, 26–35. [CrossRef] [PubMed]

24. Henault, S.; Karoui, C.; Boschiroli, M.L. A PCR-based method for tuberculosis detection in wildlife. *Dev. Biol.* **2006**, *126*, 23–132, discussion 325–126.

25. Roring, S.; Scott, A.; Brittain, D.; Walker, I.; Hewinson, G.; Neill, S.; Skuce, R. Development of variable-number tandem repeat typing of *Mycobacterium bovis*: Comparison of results with those obtained by using existing exact tandem repeats and spoligotyping. *J. Clin. Microbiol.* **2002**, *40*, 2126–2133. [CrossRef]

26. Supply, P.; Allix, C.; Lesjean, S.; Cardoso-Oelemann, M.; Rusch-Gerdes, S.; Willery, E.; Savine, E.; de Haas, P.; van Deutekom, H.; Roring, S.; et al. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of *Mycobacterium tuberculosis*. *J. Clin. Microbiol.* **2006**, *44*, 4498–4510. [CrossRef]

27. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870–1874. [CrossRef]

28. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. *Nat. Methods* **2012**, *9*, 772. [CrossRef]

29. Richomme, C.; Rivière, J.; Hars, J.; Boschiroli, M.L.; Gueneau, E.; Fediaevsky, A.; Dufour, H. Tuberculose bovine: Infection de sangliers dans un parc de chasse. *Bull. Epid. St. Anim. Alim.* **2013**, *56*, 14–16.

30. Thomas, J.; Infantes-Lorenzo, J.A.; Moreno, I.; Romero, B.; Garrido, J.M.; Juste, R.; Dominguez, M.; Dominguez, L.; Gortazar, C.; Risalde, M.A. A new test to detect antibodies against *Mycobacterium tuberculosis* complex in red deer serum. *Vet. J.* **2019**, *244*, 98–103. [CrossRef]

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

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).