# Deciphering the role of host species for two Mycobacterium bovis genotypes from the European 3 clonal complex circulation within a cattle-badger-wild boar multihost system

Laetitia Canini, Gabriela Modenesi, Aurélie Courcoul, Maria-laura Boschiroli, Benoit Durand, Lorraine Michelet

► To cite this version:

Laetitia Canini, Gabriela Modenesi, Aurélie Courcoul, Maria-laura Boschiroli, Benoit Durand, et al. Deciphering the role of host species for two Mycobacterium bovis genotypes from the European 3 clonal complex circulation within a cattle-badger-wild boar multihost system. MicrobiologyOpen, 2023, 12 (1), pp.e1331. 10.1002/mbo.3.1331. anses-03909482

HAL Id: anses-03909482

https://anses.hal.science/anses-03909482

Submitted on 21 Dec 2022

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
ORIGINAL ARTICLE

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

WILEY

# Deciphering the role of host species for two *Mycobacterium bovis* genotypes from the European 3 clonal complex circulation within a cattle-badger-wild boar multihost system

Laetitia Canini<fnref n="1" />

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

Gabriela Modenesi<fnref n="1" />

Aurélie Courcoul<fnref n="1" />

Maria-Laura Boschiroli<fnref n="2" />

Benoit Durand<fnref n="1" />

Lorraine Michelet<fnref n="2" />

### Correspondence

Laetitia Canini, Paris-Est University, Epidemiology Unit, Laboratory for Animal Health, Anses, Maisons-Alfort, France.
Email: laetitia.canini@anses.fr

### Present address

Gabriela Modenesi, Santé publique France, French National Public Health Agency, Regional Unit (Cire), Ile-de-France, France.
Aurélie Courcoul, Oniris, Nantes 44300, France.

### Funding information

Ministère de l'Agriculture et de l'Alimentation, Grant/Award Number: RFSA 170131

## Abstract

Bovine tuberculosis is a common disease affecting cattle and wildlife worldwide. *Mycobacterium bovis* circulation in wildlife decreases the efficacy of surveillance and control programs in cattle. Strains of the European 3 clonal complex are the most frequent in France. The aim of our work was hence to investigate the role played by cattle and wildlife species in the circulation of two *M. bovis* European 3 strains circulation. WGS of *M. bovis* strains collected between 2010 and 2017 in two distinct areas (Nouvelle-Aquitaine region, NAq, and Côte-d'Or département, CdO), from badgers, wild boars, and cattle were used in an evolutionary model to infer the transition between the three species. We computed host species transition and persistence between two consecutive nodes and the average number of transitions per tree. In total, 144 and 218 samples were collected respectively in CdO and NAq. In CdO, three between-species transition rates stood out: from cattle to badgers, from badgers to wild boars, and from wild boars to cattle. In NAq an additional fourth transition rate was identified: from badgers to cattle. However, host transition remained a rare event. Our results suggest that wild boars could be an intermediary host between badgers and cattle in the circulation of the studied strains in CdO and NAq. Our results also highlight the differences between these two areas, suggesting that the transition pattern does not only depend on the host species and other ecological, landscape and anthropic factors are important.

## KEYWORDS

bovine tuberculosis, multi-host system, *Mycobacterium bovis*, phylodynamics

# 1 | INTRODUCTION

Bovine tuberculosis (bTB) is a disease affecting cattle and wildlife worldwide (Bovine tuberculosis (2021)). *Mycobacterium bovis* can infect a large variety of wildlife hosts (Fitzgerald & Kaneene, 2013),

which differ from country to country. *M. bovis* was detected in European badgers (*Meles meles*) in the UK (Rivière et al., 2014), Ireland, and continental Europe, in wild boars (*Sus scrofa*) in continental Europe (Rivière et al., 2014), in red foxes (*Vulpes vulpes*) in France (Michelet et al., 2018), in cervids, more specifically red deer

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2022 The Authors. *MicrobiologyOpen* published by John Wiley & Sons Ltd.

<footnote n="1">¹Epidemiology Unit, Laboratory for Animal Health, Anses, Paris-Est University, Maisons-Alfort, France</footnote>

<footnote n="2">²Tuberculosis National Reference Laboratory, Bacterial Zoonosis unit, Laboratory for Animal Health, Anses, Paris-Est University, Maisons-Alfort, France</footnote>
(Cervus elaphus) and roe deer (Capreolus capreolus) in continental Europe (Rivière et al., 2014), Sika deer (Cervus nippon) recently identified in Ireland (Kelly et al., 2021) and white-tailed deer (Odocoileus virginianus) and elk (Cervus canadensis) in North America or brush-tailed possums (Trichosurus vulpecula) in New Zealand (Fitzgerald & Kaneene, 2013). The circulation of M. bovis in wildlife hampers control programs when implemented.

In France, a control program in line with the European Union (EU) Directive 64/432/EEC has been implemented starting in 1954 to eradicate bTB in cattle farms. This program led to a rapid decrease in herd incidence (Michelet et al., 2020), resulting in a disease-free status in 2001 when herd prevalence was below 0.1% for six consecutive years (Decision 2001/26/EC). The bTB-free status alleviates the control measures for export and therefore ensures French cattle farming competitiveness.

There are two surveillance programs in France: one for cattle and one for wildlife. The cattle surveillance program relies on three components. First, at slaughterhouses at the national level, all carcasses are systematically inspected, with incisions of specific tissues (lungs, retropharyngeal, tracheobronchial, and mediastinal lymph nodes). Samples from suspect lesions are sent for bTB confirmation by polymerase chain reaction (PCR) or bacteriology to certified laboratories. Second, depending on the epidemiological situation in the *département* (intermediate administrative divisions of France), periodic systematic screening of all animals over 6 weeks of age is performed at a regularity ranging from yearly to none. Finally, based on the epidemiological investigation of bTB-infected farms, targeted screening is performed on animals before they depart from at-risk farms (Delavenne et al., 2020a).

The wildlife surveillance program created in 2011, called Sylvatub, focuses on red deer, roe deer, wild boars, and badgers (Rivière et al., 2014). It is based on event-based surveillance (e.g. carcass inspection of hunted wild boars and cervids), enhanced event-based surveillance (e.g., carcass inspection of badgers found dead on the roadside), and programmed surveillance (e.g. badger trapping or hunted wild boar direct diagnosis). The implementation of these modalities depends on the level of surveillance defined at the *département* level, which in turn depends on their epidemiological situation. Animals are necropsied and appropriate samples are taken for bTB diagnosis by PCR and culture and molecular typing when confirmed positive (Réveillaud et al., 2018).

In France, the epidemiological situation is heterogeneous. Few areas concentrate most national outbreaks, such as Côte-d'Or *département* in the East of France or the Nouvelle-Aquitaine region in the South-West. As a result, in these two areas, surveillance of cattle was biennial until 2018 and is annual since 2018 in the Dordogne *département*; and for wildlife, the program includes event-based, enhanced event-based, and programmed surveillance (Delavenne et al., 2020b).

In addition, molecular typing revealed that M. bovis strains circulating in these two regions are specific per area. For instance, SB0120 VNTR profile 5 5 4 3 11 4 5 6 (SB0120-NAq) is mainly found in an area of Nouvelle-Aquitaine overlapping Dordogne, Haute-Vienne, Charente and Charente-Maritime *départements*, while

SB0120 VNTR profile 5 5 4 3 11 4 5 6 (SB0120-CdO) is mainly found in Burgundy, especially Côte-d'Or *département* (Michelet et al., 2020). To differentiate strains of these dominant genotypes within these regions at a finer scale, a more discriminating method is required such as whole genome sequencing (WGS), which has shown higher resolution (Crispell et al., 2017; Price-Carter et al., 2018).

Despite being phylogenetically close genotypes belonging to the European 3 (Eu3) clonal complex (also described as lineage La1.2 (Zwyer et al., 2021)) affecting the same host species (Hauer et al., 2015, 2019), the epidemiologic situations were contrasted in these two areas. While the number of cattle outbreaks has been steadily decreasing in CdO, the number of detected cases has been increasing in NAq (Delavenne et al., 2020a). In wildlife, a similar trend was observed in badgers, with the apparent prevalence decreasing from 8.1% to 4.2% in Côte-d'Or between 2013 and 2014 (as identified by culture) and 2016–2017 (as identified by PCR) while it increased in Nouvelle-Aquitaine from 2.7% to 5.3% (Réveillaud et al., 2018). However, the apparent prevalence of wild boars decreased in both areas (from 3.1% to 2.2% in Côte-d'Or and from 4.1% to 2.7% in Nouvelle-Aquitaine during the same time intervals) (Réveillaud et al., 2018).

The drivers of these epidemic dynamics are still unclear. In particular, the role played by each species in these two multi-host systems remains to be determined. Previous studies on bTB have focused on interactions between two species: primarily cattle and badgers (Biek et al., 2012; Bouchez-Zacria et al., 2018; Crispell et al., 2019; Rossi et al., 2020; Trewby et al., 2014) but also cattle and possums (Crispell et al., 2017), cattle and elk (Salvador et al., 2019), cattle and cervids (Crispell et al., 2020), or wild boars and cervids (Zanella et al., 2008). However, species such as badgers and wild boars have different life traits: while badgers are sedentary and have a life expectancy of about 14 years, wild boars can travel long distances and are often hunted before they are 4 or 5 years old (Byrne et al., 2014; Podgórski et al., 2013). In addition, the amount and the timing of M. bovis shedding by the different species, hence their infectiousness, may differ. Thus, the roles played by different species in the multi-host system could be different.

The aim of our work was therefore to investigate the role played by each species in the circulation of two Eu3 M. bovis lineages using genomic data. To do so, we defined two study areas in CdO and NAq, in which samples collected from cattle, badgers, and wild boars have been analyzed by WGS. We then used these data to model the evolutionary history of the pathogen and to infer host species of ancestors. This allowed us to analyze the transitions between species in a three-species multi-host system.

## 2 | MATERIAL AND METHODS

### 2.1 | bTB detection

All bTB detections which have been declared to the Directorate General on Food Safety (Direction Générale de l'Alimentation - DGA)
between 2010 and 2017 included herd identification for cattle and date of bTB confirmation. The National Reference Laboratory (NRL) (Anses, Maisons-Alfort) verified the data set for consistency with the samples that are analyzed at the lab for bTB confirmation.

Suspect cattle had been identified either by skin tests (in the cervical region, using single intradermal tuberculin test (SITT) or single intradermal comparative tuberculin test (SICTT)) provided by different surveillance protocols (periodic screening, epidemiological investigations, pre-movement of cattle (Delavenne et al., 2020a)) or following the detection of lesions at the slaughterhouse. For each detection, the presence of *M. bovis* was confirmed by PCR and/or bacterial culture (Delavenne et al., 2020a). Isolated strains are genotyped by spoligotyping and MIRU-VNTR. Only isolates belonging to the SB0120 spoligotype were included, more specifically isolates of genotype SB0120-CdO in CdO and isolates of genotype SB0120-NAq in NAq.

Wildlife animals shot during hunting or found dead were subject to a necropsy and samples were collected to detect mycobacteria by PCR and culture. All strains isolated from badgers and wild boars identified during the wildlife surveillance program in the study area were included in our study. We decided to exclude samples collected from red deer, roe deer, and red foxes since their limited number (n = 2, 1, and 5, respectively) would have altered parameter inference (Réveillaud, 2013).

## 2.2 | Study areas selection

In CdO and NAq, study areas were defined according to the following criteria: (1) the municipality with the most isolates was included; (2) the final study area was well delimited (i.e., it was mostly surrounded by municipalities without detected cases); (3) the final study area was compact (i.e., with a limited number of municipalities without cases detected within it). The municipalities with detected cases were municipalities where infected wildlife was identified or municipalities with pastures belonging to farms with outbreaks. The pastures were identified with the 2013–2018 French graphic parcel register (Relevé

Parcellaire Graphique, RPG) provided by the French Ministry of Agriculture. We call below "outbreak" the official declaration of one or several bTB-infected animals in a given farm. A farm could hence have several outbreaks (in case of breakdown after the farm has recovered the bTB-free status). We limited the number of samples to be analyzed per outbreak to three. Indeed, in NAq, the number of available isolates per outbreak varied between one and 26 and was <3 in most outbreaks (84%). As two or three different VNTR profiles were identified in some outbreaks (2.5%), bounding the number of isolates per outbreak to three allowed limiting the number of isolates to sequences, while allowing for the detection of different sequences in the same outbreak.

CdO was the pilot area and the municipalities were chosen manually according to the detection of infected wildlife or localization of pastures belonging to infected farms. In CdO, the selected area consisted of 38 municipalities, covering 499 km², in which 144 isolates (only one *M. bovis* strain per individual animal) were selected between 2009 and 2014. The number of samples to be analyzed was limited to three per outbreak. The host species were cattle (n = 77, from 74 outbreaks), badgers (n = 52) and wild boars (n = 15) (Figure 1, left panel) (Table A1).

For NAq, we decided to formalize the municipality selection process with a procedure that reflects the decisions taken to select the municipalities in CdO. First, for each municipality, the number of samples available was defined as the total of available samples collected from the outbreak and the wildlife. We then defined the "starting zone" as the municipality with the most available samples and the neighboring municipalities. We iteratively added one municipality randomly selected from the list of neighboring municipalities with at least one sample. The added municipality was thus adjacent to the previously defined area. Iterations stopped once a predefined maximal number of samples to include (n_max) was reached. We ran this iterative process 1000 times. Among the 1000 areas thus obtained, we selected the best-delineated area, minimizing the ratio of the total length of borders with noninfected municipalities to the perimeter of the selected area.

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** Temporal distribution of sequences for each species (blue for badgers, green for cattle, and red for wild boars) included in the study areas of Côte-d'Or (left) and Nouvelle-Aquitaine (right).</figcaption></figure>
The resulting selected area in NAq consisted of 76 municipalities, covering 1540 km², in which $n_{max}$ = 219 samples collected between 2010 and 2017 were selected. The host species were cattle (n = 161, from 95 outbreaks), badgers (n = 41) and wild boars (n = 17) (Figure 1, right panel). In CdO, the maximum number of samples was collected in 2014 (n = 51), the minimal number in 2009 (n = 11) with no clear trend, whereas in NAq, the number of selected samples varied between 8 in 2010 and 69 in 2017 with a trend toward an increase with time (Figure 1).

## 2.3 | Sequencing and alignment

Thermolysates of selected isolates were sequenced by Illumina sequencing (paired-end 2\*250 bp) at Genoscreen (Pasteur Institute, Lille) for CdO and Illumina sequencing (paired-end 2\*150 bp) at the Paris Brain Institute (ICM) for NAq. Sequencing quality was controlled using FASTQC with an acceptability Phred score threshold of 30. Sequence alignment and Single-nucleotide polymorphism (SNP) calling were computed at the NRL using the Mb3601 reference strain (Branger et al., 2020) on Bionumerics software, version 7.6 (AppliedMath, Belgium). Identified SNPs were selected according to strict criteria of the wgSNP module: (1) they had to be present on at least five reads in both forward and reverse direction, (2) 12 base pairs had to separate them, (3) they were not present in repetitive regions of the genome, and (4) ambiguous SNPs (at least one unreliable (N) base, ambiguous (non ATCG) base or gap) were not included. SNPs were then used to reconstruct a maximum parsimony tree on Bionumerics to identify genetic outliers.

## 2.4 | Bayesian phylogenetic modeling

We used BEAST2 (Bayesian evolutionary analysis by sampling trees) 2.6.4 to model bTB evolutionary (Bouckaert et al., 2019). The sequences for each region were analyzed separately.

## 2.5 | The structured coalescent model

The differences in surveillance protocols between cattle and wildlife induce sampling biases. Indeed, wildlife cannot be exhaustively monitored, while all slaughtered cattle are tested for bTB. To take into account this sampling bias, we used the approximation of the structured coalescent as implemented in the Mascot (Marginal Approximation of the Structured COalescent) 2.1.2 package (Müller et al., 2018). Indeed structured coalescent population models, contrary to migration models, are less susceptible to sampling bias (Müller et al., 2017). We assumed constant between-species transition rates in times and used the Bayesian stochastic search variable selection (BSSVS) procedure to select only transition rates that explain the transition of *M. bovis* between the different species (Lemey et al., 2009). This procedure was designed to limit the number of transition rates to infer only those adequately explaining the

diffusion between the subpopulations (here the host species) (Lemey et al., 2009). The estimated transition rates were backward in time, however, to avoid confusion we present the transition rates as the forward transition rates thereafter.

## 2.6 | Substitution and molecular clock model selection

To select the best-fitting model, the marginal likelihood (ML) was computed using the nested sampling algorithm as implemented in the NS 1.1.0 package (Russel et al., 2019). Models were then compared two by two by computing the Bayes factor (BF) as $BF = \log(ML2)-\log(ML1)$, where ML1 and ML2 are the ML of models 1 and 2, respectively. The level of support was considered overwhelming when $|BF| > 150$, strong if $20 < |BF| \le 150$, positive if $3 < |BF| \le 20$ and hardly worth mentioning if $1 \le |BF| \le 3$ (Kass & Raftery, 1995). The model with the largest ML was favored. We tested the substitution models (JC69, HKY, and GTR) and three molecular clock models (strict, uncorrelated exponential, and uncorrelated lognormal relaxed molecular clock).

## 2.7 | Parameters inference

To infer the parameters, we set the Markov chain Monte Carlo (MCMC) length to 50 million, the burn-in to 10%, and sampled every 5000 iterations. Four replicates were performed. For each replicate, we inspected each parameter trace, looking for a stationary distribution, as well as the effective sample size (ESS). We considered ESS > 200 would guarantee that samples were independent and only selected models for each of the inferred parameters. We combined the four replicates with LogCombiner v2.6.6 with a lower sampling frequency to analyze thereafter 10,000 parameters. All trees were plotted using the ggtree package in R (Yu et al., 2017).

## 2.8 | Host species transition

To study host species transitions, we resampled 1001 trees. We recorded for each node of each tree the predicted host species, its probability, and its height. We considered that the host species was known if its probability was >0.9 and unknown otherwise. We then computed the number of host species transitions between two consecutive nodes as well as the number of host species persistence (i.e., when the descendant host species is the same as the parental host species). If the host species for the parent and/or the descendant node was unknown, then the transition was considered unknown as well. We computed the host species transition and persistence at the lineage level. It is noteworthy that two consecutive nodes do not necessarily represent two distinct hosts; indeed these two nodes could represent strains hosted by the same individual or conversely strains that could have been transmitted to one or several other individuals. We computed the average number of transitions
per tree for each parental node time as the sum of each transition or persistence divided by the total number of trees.

## 3 | RESULTS

From the 144 SB0120-CdO strains selected in CdO and from the 219 SB0120-NAq strains selected in NAq, 123 and 290 SNPs were identified relative to the Mb3601 reference strain, respectively. When multiple isolates were available for a given outbreak, two or three unique SNPs sequences were detected in 39.7% of CdO outbreaks and 51.7% of NAq outbreaks.

### 3.1 | Bayesian phylogenetic estimates

According to the parameters, trace inspection showing stationary distribution, and ESS < 200, all models fitted the data well. The best-

fitted model, was an HKY substitution model, with an uncorrelated lognormal molecular clock for both study areas (Table A2). The corresponding inferred maximum clade credibility trees are shown in Figure 2.

According to the best-fitting model for CdO, the most recent common ancestor (MRCA) of the studied isolates circulated in 1997 (with a 95% credibility interval between 1989 and 2002), as suggested by the median tree height (17.5 years, 95% highest posterior density (HPD): [12.6, 25.1]). It is however unclear whether the MRCA was hosted by wildlife or cattle (median host probability: 0.49 (95% HPD = [0, 0.97]) for wild boar, 0.35 (95% HPD = [0, 0.97]) for cattle, and 0.05 (95% HPD = [0, 0.33]) for badger). In NAq, the best-fitting model predicted that the MRCA of the studied isolates circulated in 1991 (with a 95% credibility interval between 1975 and 1999), according to the median tree height (25.4 years, 95% HPD: [18.1, 42.0]). While the HPD intervals are overlapping, the predicted MRCA would have been hosted in wildlife (host probability: 0.55 (95% HPD = [0.02; 0.94]) for wild boar, vs. 0.40 (95% HPD = [0.01; 0.94]).

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** Phylogenetic trees for Côte-d'Or (left) and Nouvelle-Aquitaine (right). The tips are shown by triangles whose colors represent the species. The pies represent the median probability that the internal nodes were hosted by the three different species (badger in blue, cattle in green and wild boar in red). Pies were drawn for nodes with posterior probabilities >0.9.</figcaption></figure>
0.92)) for badger and 0.02 (95% HPD = [0; 0.11]) for cattle). The mean molecular clock was significantly (p < 0.001, Wilcoxon test) smaller in CdO with a median estimate of 0.42 substitution/genome/year (95% HPD: [0.31, 0.54]) versus 0.57 substitution/genome/year (95% HPD: [0.44, 0.71]) for NAq.

For both study areas, we inferred that the effective population size was the largest for badgers (CdO: median of 9.4, 95% HPD: [4.9, 14.6]; NAq: median of 21.2, 95% HPD: [13.4, 29.6]) before cattle (CdO: median of 3.8, 95% HPD: [1.6, 6.6]; NAq: median of 5.0, 95% HPD: [1.8, 9.6]) and wild boars (CdO: median of 2.7, 95% HPD: [0.3, 7.0]; Naq: median of 2.0, 95% HPD: [0.3, 5.3]). This means that two lineages hosted by badgers are expected to coalesce more slowly than two lineages hosted by cattle or wild boar because the coalescent rate is the inverse of the effective population size.

## 3.2 | Host species transition

Table 1 shows the proportion of transition rates selected (i.e., different from zero) by the BSSVS procedure as well as the inferred transition rates median and HPDs when these rates were selected. We defined the selection threshold as the proportion of transition rates selected by the BSSVS procedure. For example, a selection threshold of 0.5 for a given transition rate means that this transition rate was selected by the BSSVS procedure in half of the outputs. Considering an arbitrary selection threshold of 0.80, three transition rates were selected in CdO, namely from cattle to badgers, from badgers to wild boars, and from wild boars to cattle, and four in NAq which were the same as in CdO with the additional transition rate from badgers to cattle.

In both study areas, the inferred transition rates showed large variations. It is therefore difficult to conclude the frequency of each transition event. Globally, these results suggest that *M. bovis* migrated from cattle to badgers, from badgers to wild boars, and from wild boars to cattle. In addition, in NAq a transition back from badgers to cattle was also predicted.

The proportion of unknown events, when the host species is considered unknown for a probability <0.9, is more important in NAq

(0.89) than in CdO (0.36). For the known events, the majority is species persistence for both study areas with 0.90 in CdO (with 89.2% of persistence events being cattle persistence, 8.4% being badger persistence and 2.4% being wild boar persistence) and 0.88 in NAq (with 50.7% of persistence events being cattle persistence, 19.3% being badger persistence and 30.0% being wild boar persistence). Consistently with the inferred transition rate, the most frequent between-species transition was cattle-to-badger which represented 95.2% of all between-species transition events in CdO, and wild boar-to-cattle which represented 83.4% of all between-species transition events in NAq. These results suggest that between-species transitions remain rare events. Figure 3 shows the evolution with time of the average number of lineages and the proportion of transitions and persistence events per tree for both study areas. Moreover, it is worth mentioning that persistence events changed over time. Indeed, for both study areas, cattle persistence was most frequent and while badger persistence is identified onward starting as early as 1979 in NAq, it is gaining importance starting in 2008 only in CdO. On the opposite, cattle persistence is identified first in 1990 in CdO while only in 2009 in NAq.

## 4 | DISCUSSION

In this study, we studied in two distinct study areas the transition of two different *M. bovis* strains belonging to the Eu3 clonal group among three host species, namely cattle, badgers, and wild boars. For this purpose, we used a structured coalescent model to infer transition rates between these three subpopulations. We showed that while the transition events remain rare events, our model predictions suggest that wild boars may be intermediaries for the transmission of *M. bovis* from badgers to cattle in both CdO and NAq.

We selected the same Bayesian evolutionary models for both study areas. More specifically, we selected an HKY substitution model and a lognormal molecular clock. The HKY substitution model has been previously used to model *M. bovis* phylogeny in cattle and brushtail possums in New Zealand (Crispell et al., 2017) and elks,

<table><caption>TABLE 1 Median transition rates for both study areas, Côte-d'Or (CdO) and Dordogne/Haute-Vienne (NAq)</caption><thead><tr><th rowspan="2">Transition rate</th><th colspan="3">CdO</th><th colspan="3">NAq</th></tr><tr><th>Selected</th><th>Median</th><th>HPD</th><th>Selected</th><th>Median</th><th>HPD</th></tr></thead><tbody><tr><td>From cattle to badgers</td><td>0.99</td><td>0.27</td><td>0.08-0.54</td><td>0.80</td><td>0.17</td><td>0.03-0.36</td></tr><tr><td>From wild boars to badgers</td><td>0.58</td><td>0.09</td><td>0.01-0.47</td><td>0.41</td><td>0.11</td><td>0.01-0.30</td></tr><tr><td>From badgers to cattle</td><td>0.63</td><td>0.09</td><td>0.01-0.40</td><td>0.88</td><td>0.56</td><td>0.11-1.32</td></tr><tr><td>From wild boars to cattle</td><td>0.94</td><td>0.24</td><td>0.03-0.84</td><td>0.99</td><td>1.03</td><td>0.31-2.57</td></tr><tr><td>From badgers to wild boars</td><td>0.97</td><td>0.24</td><td>0.03-0.83</td><td>0.90</td><td>0.42</td><td>0.08-1.16</td></tr><tr><td>From cattle to wild boars</td><td>0.60</td><td>0.13</td><td>0.01-0.62</td><td>0.42</td><td>0.25</td><td>0.01-1.07</td></tr></tbody></table>

Note: "Selected" represents the proportion transition rate selected by the BSSVS procedure (i.e., non-zero). The median and HPD were computed for the selected transition rates.

Abbreviation: HPD, 95% highest probability distribution.
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** The proportion of transitions per tree and per year for CdO (left) and NAq (right). Computed from 1001 sampled posterior trees. Ca stands for Cattle, Ba for Badger, and Wb for Wild boar. Unknown transitions are considered when the posterior probability of the parent or descendant node is <0.9. The black line represents the average number of lineages per tree (right axis).</figcaption></figure>

white-tailed deer, and cattle in Michigan, USA (Salvador et al., 2019). We estimated two significantly different substitution rates for both regions. Both estimates were higher than previously estimated in Northern Ireland from cattle and badgers (Biek et al., 2012; Trewby et al., 2016) and in Michigan from cattle and elk (Salvador et al., 2019). However, the substitution rates estimated from cattle and possums in New Zealand are in the same order of magnitude as the one we estimated from NAq (Crispell et al., 2017) and the one estimated from the South-West of France from cattle and badgers (Duault et al., 2022) is in the same order of magnitude as the one we estimated for CdO. While, as noted by Duault et al., the differences between the studies could result from the different *M. bovis* lineage or the sampled species (Duault et al., 2022), this does not explain the significantly higher substitution rates inferred in NAq than in CdO. Indeed, for both study areas, the same host species were sampled during the same period. However, even if the spoligotype was the same in both study areas, the VNTR profiles differed. Hauer et al. showed that while SB0120 spoligotype was found all over France, specific VNTR profiles that spread locally were identified (such as 5 5 4 3 11 4 5 6 in CdO and 5 3 5 3 9 4 5 6 in NAq) (Hauer et al., 2015). While belonging to the Eu3 clonal group and sharing the same spoligotype, these strains were identified on different phylogenetic branches (Hauer et al., 2019) suggesting that VNTR profile and substitution rates could represent specific lineage characteristics.

To infer the internal nodes for host species, we used a structured coalescent method, which limits the impact of sampling bias (Müller et al., 2017). Indeed, contrary to the migration method (Lemey et al., 2009), structured coalescent methods do not assume that the migration process and tree-generating process are independent. According to the selected transition rates between badgers and cattle, the relationship between both species was clearly different between the two study areas: transition rate from cattle to badgers was more frequently predicted than from badgers to cattle in CdO, contrary to NAq where both transition rates were predicted in more than 80% of the outputs. This could relate to the implementation of bTB biosecurity measures toward wildlife that were more stringent in CdO than in NAq. Recent studies including samples collected from two species, namely badger and cattle, are in favor of the transition from badger to cattle (Crispell et al., 2019; Duault et al., 2022; van Tonder et al., 2021), however in one of these studies, while the overall badger-to-cattle transition rate was higher than the cattle-to-badger transition rate, a more refined analysis on transmission cluster levels revealed that for 4/12 clusters, cattle-to-badger transition rates were higher than badger-to-cattle transition ones (van Tonder et al., 2021). Moreover, an epidemiological study reconstructing the contact network between badger setts and cattle farms concluded the intermediary role of badgers in *M. bovis* transmission in the South-West of France (Bouchez-Zacria et al., 2018). In addition,
according to the French graphic parcel register (Relevé parcellaire Graphique, https://www.geoportail.gouv.fr), the landscape in NAq is more fragmented with numerous interfaces between pastures and wooded areas, than in CdO. Fragmented landscapes were shown to be associated with lower adult badger population densities, which was corroborated by their lower densities in NAq than in CdO (Jacquier et al., 2021). However, the increased interfaces between pastures and wooded areas could increase the contact rates between badgers and cattle in NAq (Bouchez-Zacria et al., 2017) and explain the transition back and forth between these two species in this area compared to CdO. Also, bTB apparent prevalence in badgers might have increased between 2014 and 2016–2017 in NAq, while it has decreased in CdO during the same time interval (Réveillaud et al., 2018). The lack of evidence for the badger-to-cattle transition could also result from a lower statistical power due to the different sampling schemes in CdO and NAq.

It is noteworthy that we studied for the first time a multi-host system including wild boars in addition to badgers and cattle. We showed that among the selected transition rate in both NAq and CdO, the transition rate from wild boar to cattle was identified, suggesting a relatively frequent transition from wild boar to cattle and the intermediary role played by wild boars between badgers and cattle. This highlight the important role played by wild boars in the spread of *M. bovis* to other species. This could be related to the large distances traveled by wild boars (Podgórski et al., 2013) compared to badgers (Byrne et al., 2014). In addition, cattle movement, even if less frequent, could also enhance *M. bovis* spread between distant farms (Palisson et al., 2016). Furthermore, in both study areas host species persistence was identified as the main event while between-species transition represented less than 12% of known events, highlighting that between-species transition remains a rare event.

To summarize, we showed that while wild boars played the role of intermediary host between badgers and cattle, the role of badgers differed between both regions: in CdO badgers were intermediaries from cattle to wild boars, whereas in NAq, badgers transmitted to both cattle and wild boars. Several factors could explain these differences. Some of them are inherent to our study design, such as the lineage that can lead to different substitution rates, the temporal depth, or the number of strains collected, while other relates to the landscape or the population density of wildlife species.

Our work has several limitations. First, we excluded red foxes, red deer, and roe deer samples from our study because of the limited number of available samples, which prevented us from evaluating their role in this multi-host system. BTB-infected foxes (*Vulpes vulpes*) were detected in Dordogne in 2015 (Michelet et al., 2018). A further study showed that in Dordogne, Landes, and Charente, bTB prevalence in foxes ranged between 5% and 10%, similar to that observed in badgers and wild boars (Richomme et al., 2020). However, due to the lack of information concerning red foxes, we cannot conclude about their role. Second, the sampling procedure varies between cattle and wildlife, with nearly exhaustive testing of cattle while samples collected from wildlife depend on events such as hunting or road kills leading to an underestimation of *M. bovis*

infection prevalence in wildlife. In addition, wildlife carcasses are subject to contamination and deterioration lowering the culture sensitivity and thus poorer statistical representativeness in our sample compared to cattle (Rivière et al., 2015). While inference with structured coalescent models is less altered by sampling bias, unsampled demes (such as the red foxes or cervids in our study) could reshape our results. Thirdly, our approach did not allow the inclusion of spatialized information that could describe the localization of sampled isolates or the different sizes of the home range for badgers and wild boars, or the change of pastures for cattle. Finally, a large number of transition events was labeled as unknown. This is particularly true for NAq where nearly 90% of transition could not be labeled. The results concerning host species persistence and between-species transition should therefore be considered with caution.

## 5 | CONCLUSIONS

In conclusion, using a Bayesian evolutionary model, we inferred transition rates between cattle, badgers, and wild boars. Although this approach does not allow us to quantify within-species transmission, our result shed light on the wild boar role, which appears to act as an intermediary between badgers and cattle in the circulation of two different *Eu3 M. bovis* in two distinct study areas.

## AUTHOR CONTRIBUTIONS

**Laetitia Canini:** Data curation (equal); formal analysis (equal); methodology (equal); validation (equal); visualization (equal); writing – original draft (equal). **Gabriela Modenesi:** Data curation (supporting); formal analysis (supporting); visualization (supporting); writing – original draft (supporting). **Aurélie Courcoul:** Conceptualization (equal); funding acquisition (equal); writing – original draft (equal).

**Maria-Laura Boschiroli:** Conceptualization (equal); project administration (equal); resources (equal); supervision (equal); writing – original draft (equal). **Benoit Durand:** Conceptualization (equal); funding acquisition (equal); supervision (equal); writing – original draft (equal). **Lorraine Michelet:** Conceptualization (equal); formal analysis (equal); funding acquisition (equal); investigation (equal); methodology (equal); project administration (equal); writing – original draft (equal).

## ACKNOWLEDGMENTS

The French Ministry of Agriculture financed the sampling and the analyses in the framework of the RFSA call on TB projects (Anses-DGAI credit agreement no 170131).

## CONFLICT OF INTEREST

None declared.

## DATA AVAILABILITY STATEMENT

All data are provided in full in the results section of this paper apart from all WGS data which are available in NCBI GenBank under
BioProject PRJEB46102 for NAq: https://www.ncbi.nlm.nih.gov/bioproject/PRJEB46102 and PRJEB46417 for CdO: https://www.ncbi.nlm.nih.gov/bioproject/PRJEB46417. The individual isolates can be accessed under the following BioSample accession numbers: SAMEA8955321 - SAMEA8955540 for NAq and SAMEA8987071 - SAMEA8987214 for CdO.

## ETHICS STATEMENT

None required.

## ORCID

Laetitia Canini http://orcid.org/0000-0001-7643-1219

## REFERENCES

Biek, R., O'Hare, A., Wright, D., Mallon, T., McCormick, C., Orton, R. J., McDowell, S., Trewby, H., Skuce, R. A., & Kao, R. R. (2012). Whole genome sequencing reveals local transmission patterns of *Mycobacterium bovis* in sympatric cattle and badger populations. *PLoS Pathogens*, 8(11), e1003008.

Bouchez-Zacria, M., Courcoul, A., & Durand, B. (2018). The distribution of bovine tuberculosis in cattle farms is linked to cattle trade and badger-mediated contact networks in South-Western France, 2007–2015. *Front Vet Sci [Internet]. Disponible sur https://www.frontiersin.org/articles/10.3389/fvets.2018.00173/full*

Bouchez-Zacria, M., Courcoul, A., Jabert, P., Richomme, C., & Durand, B. (2017). Environmental determinants of the *Mycobacterium bovis* concomitant infection in cattle and badgers in France. *European Journal of Wildlife Research*, 63(5), 74.

Bouckaert, R., Vaughan, T. G., Barido-Sottani, J., Duchêne, S., Fourment, M., Gavryushkina, A., Heled, J., Jones, G., Kühnert, D., De Maio, N., Matschiner, M., Mendes, F. K., Müller, N. F., Ogilvie, H. A., du Plessis, L., Popinga, A., Rambaut, A., Rasmussen, D., Siveroni, I., ... Drummond, A. J. (2019). BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. *PLoS Computational Biology*, 15(4), e1006650.

Bovine tuberculosis: global distribution and implementation of prevention and control measures according to WAHIS data [Internet]. OIE Bulletin. [cité 27 mai 2021]. Disponible sur https://oiebulletin.com/?panorama=3-01-tb-wahis-en

Branger, M., Loux, V., Cochard, T., Boschiroli, M. L., Biet, F., & Michelet, L. (2020). The complete genome sequence of mycobacterium bovis Mb3601, a SB0120 spoligotype strain representative of a new clonal group. *Infection, Genetics and Evolution*, 82, 104309.

Byrne, A. W., Quinn, J. L., O'Keeffe, J. J., Green, S., Paddy Sleeman, D., Wayne Martin, S., & Davenport, J. (2014). Large-scale movements in European badgers: Has the tail of the movement kernel been underestimated? *Journal of Animal Ecology*, 83(4), 991–1001.

Crispell, J., Benton, C. H., Balaz, D., De Maio, N., Ahkmetova, A., Allen, A., et al. (2019). Combining genomics and epidemiology to analyse bidirectional transmission of *Mycobacterium bovis* in a multi-host system. *eLife*, 8, e45833.

Crispell, J., Cassidy, S., Kenny, K., McGrath, G., Warde, S., Cameron, H., et al. (2020). *Mycobacterium bovis* genomics reveals transmission of infection between cattle and deer in Ireland. *Microb Genomics [Internet]*, 6(8). https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7641417/

Crispell, J., Zadoks, R. N., Harris, S. R., Paterson, B., Collins, D. M., de-Lisle, G. W., Livingstone, P., Neill, M. A., Biek, R., Lycett, S. J., Kao, R. R., & Price-Carter, M. (2017). Using whole genome sequencing to investigate transmission in a multi-host system: bovine tuberculosis in New Zealand. *BMC Genomics*, 18(1), 180.

Delavenne, C., Desvaux, S., Boschiroli, M. L., Carles, S., Chaigneau, P., Dufour, B., Durand, B., Gache, K., Garapin, F., Girard, S., Jabert, P.,

Keck, N., Réveillaud, É., Rivière, J., Dupuy, C., & Chevalier, F. (2020). Surveillance de la tuberculose due à *Mycobacterium bovis* en France métropolitaine en 2018: Résultats et indicateurs de fonctionnement. *Bull Epidémiologique*, 92(1), 1–9.

Delavenne, C., Pandolfi, F., Girard, S., Réveillaud, É., Jabert, P., Boschiroli, M. L., Dommergues, L., Garapin, F., Keck, N., Martin, F., Moussu, M., Philizot, S., Rivière, J., Tourette, I., Calavas, D., Dupuy, C., Dufour, B., & Chevalier, F. (2020). Tuberculose bovine: bilan et évolution de la situation épidémiologique entre 2015 et 2017 en France métropolitaine. *Bull Epidémiologique*, 91(12), 1–9.

Duault, H., Michelet, L., Boschiroli, M. L., Durand, B., & Canini, L. (2022). A Bayesian evolutionary model towards understanding wildlife contribution to F4-family *Mycobacterium bovis* transmission in the South-West of France. *Veterinary Research*, 53(1), 28.

Fitzgerald, S. D., & Kaneene, J. B. (2013). Wildlife reservoirs of bovine tuberculosis worldwide: Hosts, pathology, surveillance, and control. *Veterinary Pathology*, 50(3), 488–499.

Hauer, A., De Cruz, K., Cochard, T., Godreuil, S., Karoui, C., Henault, S., Bulach, T., Bañuls, A. L., Biet, F., & Boschiroli, M. L. (2015). Genetic evolution of *Mycobacterium bovis* causing tuberculosis in livestock and wildlife in France since 1978. *PLoS One*, 10(2), e0117103.

Hauer, A., Michelet, L., Cochard, T., Branger, M., Nunez, J., Boschiroli, M. L., & Biet, F. (2019). Accurate phylogenetic relationships among mycobacterium bovis strains circulating in France based on whole genome sequencing and single nucleotide polymorphism analysis. *Frontiers in Microbiology*, 10, 955.

Jacquier, M., Vandel, J. M., Léger, F., Duhayer, J., Pardonnet, S., Say, L., Devillard, S., & Ruette, S. (2021). Breaking down population density into different components to better understand its spatial variation. *BMC Ecology and Evolution*, 21(1), 82.

Kass, R. E., & Raftery, A. E. (1995). Bayes factors. *Journal of the American Statistical Association*, 90(430), 773–795.

Kelly, D. J., Mullen, E., & Good, M. (2021). Bovine tuberculosis: The emergence of a new wildlife maintenance host in Ireland. *Frontiers in Veterinary Science*, 8, 632525. https://www.frontiersin.org/article/10.3389/fvets.2021.632525

Lemey, P., Rambaut, A., Drummond, A. J., & Suchard, M. A. (2009). Bayesian phylogeography finds its roots. *PLoS Computational Biology*, 5(9), e1000520.

Michelet, L., De Cruz, K., Hénault, S., Tambosco, J., Richomme, C., Réveillaud, É., Gares, H., Moyen, J. L., & Boschiroli, M. L. (2018). *Mycobacterium bovis* infection of Red Fox, France. *Emerging Infectious Diseases*, 24(6), 1150–1153.

Michelet, L., Durand, B., & Boschiroli, M. L. (2020). Tuberculose bovine: Bilan génotypique de *M. bovis* à l'origine des foyers bovins entre 2015 et 2017 en France métropolitaine. *Bull Épidémiologique*, 91(13), 1–4.

Müller, N. F., Rasmussen, D., & Stadler, T. (2018). MASCOT: Parameter and state inference under the marginal structured coalescent approximation. *Bioinformatics*, 34(22), 3843–3848.

Müller, N. F., Rasmussen, D. A., & Stadler, T. (2017, November). The structured coalescent and its approximations. *Molecular Biology and Evolution*, 34(11), 2970–2981.

Palisson, A., Courcoul, A., & Durand, B. (2016). Role of cattle movements in bovine tuberculosis spread in France between 2005 and 2014. *PLoS One*, 11(3), e0152578.

Podgórski, T., Baś, G., Jędrzejewska, B., Sönnichsen, L., Śnieżko, S., Jędrzejewski, W., & Okarma, H. (2013). Spatiotemporal behavioral plasticity of wild boar (*Sus scrofa*) under contrasting conditions of human pressure: Primeval forest and metropolitan area. *Journal of Mammalogy*, 94(1), 109–119.

Price-Carter, M., Brauning, R., de Lisle, G. W., Livingstone, P., Neill, M., Sinclair, J., Paterson, B., Atkinson, G., Knowles, G., Crews, K., Crispell, J., Kao, R., Robbe-Austerman, S., Stuber, T., Parkhill, J., Wood, J., Harris, S., & Collins, D. M. (2018). Whole genome sequencing for determining the source of *Mycobacterium bovis*.
infections in livestock herds and wildlife in New Zealand. *Frontiers in Veterinary Science*, **5**, 272.

Réveillaud, É. (2013). Sylvatub un dispositif de surveillance de la tuberculose bovine dans la faune sauvage. *Nouv Prat Vét*, **6**(25), 8.

Réveillaud, É., Desvaux, S., Boschiroli, M. L., Hars, J., Faure, É., Fediaevsky, A., Cavalerie, L., Chevalier, F., Jabert, P., Poliak, S., Tourette, I., Hendrikx, P., & Richomme, C. (2018). Infection of wildlife by *Mycobacterium bovis* in France assessment through a national surveillance system, sylvatub. *Frontiers in Veterinary Science*, **5**, 262.

Richomme, C., Réveillaud, E., Moyen, J. L., Sabatier, P., De Cruz, K., Michelet, L., & Boschiroli, M. L. (2020). Mycobacterium bovis infection in Red foxes in four animal tuberculosis endemic areas in France. *Microorganisms*, **8**(7), 1070.

Rivière, J., Carabin, K., Le Strat, Y., Hendrikx, P., & Dufour, B. (2014). Bovine tuberculosis surveillance in cattle and free-ranging wildlife in EU Member states in 2013: A survey-based review. *Veterinary Microbiology*, **173**(3), 323-331.

Rivière, J., Le Strat, Y., Dufour, B., & Hendrikx, P. (2015). Sensitivity of bovine tuberculosis surveillance in wildlife in France: A scenario tree approach. *PLoS One*, **10**(10), e0141884.

Rossi, G., Crispell, J., Balaz, D., Lycett, S. J., Benton, C. H., Delahay, R. J., & Kao, R. R. (2020). Identifying likely transmissions in *Mycobacterium bovis* infected populations of cattle and badgers using the Kolmogorov forward equations. *Scientific Reports*, **10**(1), 21980.

Russel, P. M., Brewer, B. J., Klaere, S., & Bouckaert, R. R. (2019). Model selection and parameter inference in phylogenetics using nested sampling. *Systematic Biology*, **68**(2), 219-233.

Salvador, L. C. M., O'Brien, D. J., Cosgrove, M. K., Stuber, T. P., Schooley, A. M., Crispell, J., Church, S. V., Gröhn, Y. T., Robbe-Austerman, S., & Kao, R. R. (2019). Disease management at the wildlife-livestock interface: Using whole-genome sequencing to study the role of elk in *Mycobacterium bovis* transmission in Michigan, USA. *Molecular Ecology*, **28**(9), 2192-2205.

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., Herzyk, P., Skuce, R. A., Biek, R., & Kao, R. R. (2016). Use of bacterial

whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis. *Epidemics*, **14**, 26-35.

Trewby, I. D., Young, R., McDonald, R. A., Wilson, G. J., Davison, J., Walker, N., Robertson, A., Doncaster, C. P., & Delahay, R. J. (2014). Impacts of removing badgers on localised counts of Hedgehogs. *PLoS One*, **9**(4), e95477.

van Tonder, A. J., Thornton, M. J., Conlan, A. J. K., Jolley, K. A., Goolding, L., Mitchell, A. P., Dale, J., Palkopoulou, E., Hogarth, P. J., Hewinson, R. G., Wood, J. L. N., & Parkhill, J. (2021). Inferring *Mycobacterium bovis* transmission between cattle and badgers using isolates from the randomised badger culling trial. *PLoS Pathogens*, **17**(11), e1010075.

Yu, G., Smith, D. K., Zhu, H., Guan, Y., & Lam, T. T. Y. (2017). Ggtree: An r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. *Methods in Ecology and Evolution*, **8**(1), 28-36.

Zanella, G., Durand, B., Hars, J., Moutou, F., Garin-Bastuji, B., Duvauchelle, A., Fermé, M., Karoui, C., & Boschiroli, M. L. (2008). *Mycobacterium bovis* in wildlife in France. *Journal of Wildlife Diseases*, **44**(1), 99-108.

Zwyer, M., Çavusoglu, C., Ghielmetti, G., Pacciarini, M. L., Scaltriti, E., Soolingen, D. V., et al. (2021). A new nomenclature for the livestock-associated *Mycobacterium tuberculosis* complex based on phylogenomics. *Open Res Eur*, **1**(100), 25.

**How to cite this article:** Canini, L., Modenesi, G., Courcoul, A., Boschiroli, M.-L., Durand, B., & Michelet, L. (2022). Deciphering the role of host species for two *Mycobacterium bovis* genotypes from the European 3 clonal complex circulation within a cattle-badger-wild boar multihost system. *MicrobiologyOpen*, **12**, e1331.

https://doi.org/10.1002/mbo3.1331
# APPENDIX

<table><caption>TABLE A1 Accession numbers</caption><thead><tr><td>DA</td><td>Sample</td><td>Collection date</td><td>Host scientific name</td><td>Biosample</td><td>Accession</td></tr></thead><tbody><tr><td>10-00086-00</td><td>10Z000099</td><td>2009-11-09</td><td>Bos taurus</td><td>SAMEA8955321</td><td>ERR6198117</td></tr><tr><td>10-00566-00</td><td>10Z002121</td><td>2010-02-12</td><td>Bos taurus</td><td>SAMEA8955322</td><td>ERR6198204</td></tr><tr><td>10-01026-00</td><td>10Z005043</td><td>2010-03-26</td><td>Bos taurus</td><td>SAMEA8955323</td><td>ERR6198207</td></tr><tr><td>D-10-02190</td><td>10Z008544</td><td>2010-09-18</td><td>Sus scrofa</td><td>SAMEA8955324</td><td>ERR6198208</td></tr><tr><td>D-10-02191</td><td>10Z008545</td><td>2010-09-10</td><td>Bos taurus</td><td>SAMEA8955325</td><td>ERR6198209</td></tr><tr><td>D-10-02192</td><td>10Z008546</td><td>2010-08-15</td><td>Sus scrofa</td><td>SAMEA8955326</td><td>ERR6198210</td></tr><tr><td>D-10-02315</td><td>10Z009059</td><td>2010-09-12</td><td>Sus scrofa</td><td>SAMEA8955327</td><td>ERR6198216</td></tr><tr><td>D-10-02451</td><td>10Z009733</td><td>2010-10-18</td><td>Meles meles</td><td>SAMEA8955328</td><td>ERR6198228</td></tr><tr><td>D-11-00916</td><td>11Z001983</td><td>2011-01-04</td><td>Bos taurus</td><td>SAMEA8955329</td><td>ERR6198247</td></tr><tr><td>D-11-01150</td><td>11Z002149</td><td>2011-02-08</td><td>Bos taurus</td><td>SAMEA8955330</td><td>ERR6198291</td></tr><tr><td>D-11-01158</td><td>11Z002157</td><td>2011-02-27</td><td>Meles meles</td><td>SAMEA8955331</td><td>ERR6198317</td></tr><tr><td>D-11-01374</td><td>11Z002528</td><td>2011-02-22</td><td>Bos taurus</td><td>SAMEA8955332</td><td>ERR6198349</td></tr><tr><td>D-11-01375</td><td>11Z002529</td><td>2011-03-01</td><td>Bos taurus</td><td>SAMEA8955333</td><td>ERR6198365</td></tr><tr><td>D-11-01905</td><td>11Z003570</td><td>2011-04-12</td><td>Bos taurus</td><td>SAMEA8955334</td><td>ERR6198374</td></tr><tr><td>D-11-01992</td><td>11Z003734</td><td>2011-05-03</td><td>Bos taurus</td><td>SAMEA8955335</td><td>ERR6198375</td></tr><tr><td>D-11-02082</td><td>11Z003894</td><td>2011-05-05</td><td>Bos taurus</td><td>SAMEA8955336</td><td>ERR6198376</td></tr><tr><td>D-11-03278</td><td>11Z006170</td><td>2011-10-10</td><td>Bos taurus</td><td>SAMEA8955337</td><td>ERR6198377</td></tr><tr><td>D-12-00020</td><td>12Z000029</td><td>2011-10-10</td><td>Meles meles</td><td>SAMEA8955338</td><td>ERR6198378</td></tr><tr><td>D-12-01240</td><td>12Z002503</td><td>2012-01-11</td><td>Meles meles</td><td>SAMEA8955339</td><td>ERR6198379</td></tr><tr><td>D-12-01581</td><td>12Z003372</td><td>2012-01-17</td><td>Bos taurus</td><td>SAMEA8955340</td><td>ERR6198389</td></tr><tr><td>D-12-01582</td><td>12Z003373</td><td>2013-01-22</td><td>Meles meles</td><td>SAMEA8955341</td><td>ERR6198390</td></tr><tr><td>D-12-01589</td><td>12Z003430</td><td>2012-03-12</td><td>Sus scrofa</td><td>SAMEA8955342</td><td>ERR6198391</td></tr><tr><td>D-12-01597</td><td>12Z003438</td><td>2012-01-27</td><td>Bos taurus</td><td>SAMEA8955343</td><td>ERR6198393</td></tr><tr><td>D-12-01601</td><td>12Z003442</td><td>2012-01-24</td><td>Bos taurus</td><td>SAMEA8955344</td><td>ERR6198394</td></tr><tr><td>D-12-01697</td><td>12Z003519</td><td>2012-02-14</td><td>Bos taurus</td><td>SAMEA8955345</td><td>ERR6198395</td></tr><tr><td>D-12-01854</td><td>12Z004182</td><td>2012-02-13</td><td>Meles meles</td><td>SAMEA8955346</td><td>ERR6198396</td></tr><tr><td>D-12-01929</td><td>12Z004424</td><td>2012-02-14</td><td>Bos taurus</td><td>SAMEA8955347</td><td>ERR6198398</td></tr><tr><td>D-12-02001</td><td>12Z004494</td><td>2012-03-16</td><td>Bos taurus</td><td>SAMEA8955348</td><td>ERR6198399</td></tr><tr><td>D-12-02007</td><td>12Z004505</td><td>2012-02-09</td><td>Bos taurus</td><td>SAMEA8955349</td><td>ERR6198400</td></tr><tr><td>D-12-02120</td><td>12Z004741</td><td>2012-03-23</td><td>Bos taurus</td><td>SAMEA8955350</td><td>ERR6198401</td></tr><tr><td>D-12-02121</td><td>12Z004742</td><td>2012-02-28</td><td>Bos taurus</td><td>SAMEA8955351</td><td>ERR6198442</td></tr><tr><td>D-12-02129</td><td>12Z004749</td><td>2012-03-13</td><td>Bos taurus</td><td>SAMEA8955352</td><td>ERR6198443</td></tr><tr><td>D-12-02131</td><td>12Z004751</td><td>2012-03-13</td><td>Bos taurus</td><td>SAMEA8955353</td><td>ERR6198444</td></tr><tr><td>D-12-02240</td><td>12Z004928</td><td>2012-03-16</td><td>Bos taurus</td><td>SAMEA8955354</td><td>ERR6198445</td></tr><tr><td>D-12-02242</td><td>12Z004930</td><td>2012-03-27</td><td>Bos taurus</td><td>SAMEA8955355</td><td>ERR6198447</td></tr><tr><td>D-12-02292</td><td>12Z004995</td><td>2012-03-07</td><td>Meles meles</td><td>SAMEA8955356</td><td>ERR6198448</td></tr><tr><td>D-12-02293</td><td>12Z004996</td><td>2012-03-22</td><td>Bos taurus</td><td>SAMEA8955357</td><td>ERR6198449</td></tr></tbody></table>

(Continues)
<table><caption>TABLE A1 (Continued)</caption><thead><tr><td>DA</td><td>Sample</td><td>Collection date</td><td>Host scientific name</td><td>Biosample</td><td>Accession</td></tr></thead><tbody><tr><td>D-12-02296</td><td>12Z004999</td><td>2012-04-02</td><td>Meles meles</td><td>SAMEA8955358</td><td>ERR6199027</td></tr><tr><td>D-12-02372</td><td>12Z005163</td><td>2012-04-06</td><td>Meles meles</td><td>SAMEA8955359</td><td>ERR6201794</td></tr><tr><td>D-12-02472</td><td>12Z005339</td><td>2012-04-05</td><td>Meles meles</td><td>SAMEA8955360</td><td>ERR6201795</td></tr><tr><td>D-12-02749</td><td>12Z005660</td><td>2012-05-05</td><td>Meles meles</td><td>SAMEA8955361</td><td>ERR6201796</td></tr><tr><td>D-12-02971</td><td>12Z006182</td><td>2012-05-22</td><td>Bos taurus</td><td>SAMEA8955362</td><td>ERR6201797</td></tr><tr><td>D-12-03180</td><td>12Z006994</td><td>2012-05-09</td><td>Bos taurus</td><td>SAMEA8955363</td><td>ERR6201798</td></tr><tr><td>D-12-03307</td><td>12Z007111</td><td>2013-05-24</td><td>Meles meles</td><td>SAMEA8955364</td><td>ERR6201800</td></tr><tr><td>D-12-03422</td><td>12Z007219</td><td>2012-07-23</td><td>Bos taurus</td><td>SAMEA8955365</td><td>ERR6201801</td></tr><tr><td>D-12-03423</td><td>12Z007220</td><td>2012-07-20</td><td>Bos taurus</td><td>SAMEA8955366</td><td>ERR6201802</td></tr><tr><td>D-13-00134</td><td>13Z000269</td><td>2012-11-16</td><td>Meles meles</td><td>SAMEA8955367</td><td>ERR6201804</td></tr><tr><td>D-13-00505</td><td>13Z001275</td><td>2013-01-03</td><td>Meles meles</td><td>SAMEA8955368</td><td>ERR6201806</td></tr><tr><td>D-13-00830</td><td>13Z001729</td><td>2013-01-25</td><td>Bos taurus</td><td>SAMEA8955369</td><td>ERR6201807</td></tr><tr><td>D-13-00831</td><td>13Z001730</td><td>2013-01-25</td><td>Bos taurus</td><td>SAMEA8955370</td><td>ERR6201808</td></tr><tr><td>D-13-00956</td><td>13Z001898</td><td>2013-02-01</td><td>Bos taurus</td><td>SAMEA8955371</td><td>ERR6201809</td></tr><tr><td>D-13-00962</td><td>13Z001903</td><td>2013-02-15</td><td>Bos taurus</td><td>SAMEA8955372</td><td>ERR6201810</td></tr><tr><td>D-13-01061</td><td>13Z001990</td><td>2013-01-02</td><td>Bos taurus</td><td>SAMEA8955373</td><td>ERR6201811</td></tr><tr><td>D-13-01238</td><td>13Z002208</td><td>2013-02-20</td><td>Bos taurus</td><td>SAMEA8955374</td><td>ERR6201812</td></tr><tr><td>D-13-01245</td><td>13Z002215</td><td>2013-02-15</td><td>Bos taurus</td><td>SAMEA8955375</td><td>ERR6201813</td></tr><tr><td>D-13-01303</td><td>13Z002277</td><td>2013-03-08</td><td>Bos taurus</td><td>SAMEA8955376</td><td>ERR6201814</td></tr><tr><td>D-13-01556</td><td>13Z002699</td><td>2013-03-29</td><td>Bos taurus</td><td>SAMEA8955377</td><td>ERR6201815</td></tr><tr><td>D-13-01562</td><td>13Z002705</td><td>2013-04-12</td><td>Bos taurus</td><td>SAMEA8955378</td><td>ERR6201817</td></tr><tr><td>D-13-01565</td><td>13Z002708</td><td>2013-04-22</td><td>Bos taurus</td><td>SAMEA8955379</td><td>ERR6201818</td></tr><tr><td>D-13-01566</td><td>13Z002709</td><td>2013-04-22</td><td>Bos taurus</td><td>SAMEA8955380</td><td>ERR6201819</td></tr><tr><td>D-13-02048</td><td>13Z003478</td><td>2013-06-05</td><td>Bos taurus</td><td>SAMEA8955381</td><td>ERR6201820</td></tr><tr><td>D-13-02100</td><td>13Z003643</td><td>2013-04-23</td><td>Bos taurus</td><td>SAMEA8955382</td><td>ERR6201821</td></tr><tr><td>D-13-02101</td><td>13Z003644</td><td>2013-06-07</td><td>Bos taurus</td><td>SAMEA8955383</td><td>ERR6201822</td></tr><tr><td>D-13-02225</td><td>13Z003769</td><td>2013-03-22</td><td>Bos taurus</td><td>SAMEA8955384</td><td>ERR6201823</td></tr><tr><td>D-13-02226</td><td>13Z003770</td><td>2013-06-13</td><td>Bos taurus</td><td>SAMEA8955385</td><td>ERR6201825</td></tr><tr><td>D-13-02368</td><td>13Z004012</td><td>2013-06-13</td><td>Bos taurus</td><td>SAMEA8955386</td><td>ERR6201826</td></tr><tr><td>D-13-02401</td><td>13Z004051</td><td>2013-06-13</td><td>Bos taurus</td><td>SAMEA8955387</td><td>ERR6201828</td></tr><tr><td>D-13-02436</td><td>13Z004129</td><td>2013-04-23</td><td>Bos taurus</td><td>SAMEA8955388</td><td>ERR6201829</td></tr><tr><td>D-13-02620</td><td>13Z004667</td><td>2013-06-13</td><td>Bos taurus</td><td>SAMEA8955389</td><td>ERR6201830</td></tr><tr><td>D-13-02704</td><td>13Z004785</td><td>2013-08-06</td><td>Bos taurus</td><td>SAMEA8955390</td><td>ERR6201831</td></tr><tr><td>D-13-03134</td><td>13Z005631</td><td>2013-06-13</td><td>Bos taurus</td><td>SAMEA8955391</td><td>ERR6201832</td></tr><tr><td>D-13-03137</td><td>13Z005634</td><td>2013-08-22</td><td>Sus scrofa</td><td>SAMEA8955392</td><td>ERR6201833</td></tr><tr><td>D-13-03140</td><td>13Z005637</td><td>2013-09-15</td><td>Sus scrofa</td><td>SAMEA8955393</td><td>ERR6201835</td></tr><tr><td>D-13-03227</td><td>13Z005787</td><td>2013-09-03</td><td>Bos taurus</td><td>SAMEA8955394</td><td>ERR6201836</td></tr><tr><td>D-13-03897</td><td>13Z010139</td><td>2013-10-29</td><td>Bos taurus</td><td>SAMEA8955395</td><td>ERR6201837</td></tr><tr><td>D-13-03898</td><td>13Z010140</td><td>2013-11-15</td><td>Bos taurus</td><td>SAMEA8955396</td><td>ERR6201838</td></tr><tr><td>D-13-03899</td><td>13Z010141</td><td>2013-10-13</td><td>Sus scrofa</td><td>SAMEA8955397</td><td>ERR6201839</td></tr></tbody></table>
<table><caption>TABLE A1 (Continued)</caption><thead><tr><td>DA</td><td>Sample</td><td>Collection date</td><td>Host scientific name</td><td>Biosample</td><td>Accession</td></tr></thead><tbody><tr><td>D-13-03900</td><td>13Z010142</td><td>2013-11-10</td><td>Sus scrofa</td><td>SAMEA8955398</td><td>ERR6201840</td></tr><tr><td>D-14-00007</td><td>14Z000006</td><td>2013-11-03</td><td>Sus scrofa</td><td>SAMEA8955399</td><td>ERR6201841</td></tr><tr><td>D-14-00013</td><td>14Z000010</td><td>2013-11-20</td><td>Bos taurus</td><td>SAMEA8955400</td><td>ERR6201842</td></tr><tr><td>D-14-00016</td><td>14Z000013</td><td>2013-11-19</td><td>Bos taurus</td><td>SAMEA8955401</td><td>ERR6201843</td></tr><tr><td>D-14-00019</td><td>14Z000015</td><td>2013-11-19</td><td>Bos taurus</td><td>SAMEA8955402</td><td>ERR6201844</td></tr><tr><td>D-14-00020</td><td>14Z000016</td><td>2013-11-19</td><td>Bos taurus</td><td>SAMEA8955403</td><td>ERR6201845</td></tr><tr><td>D-14-01072</td><td>14Z003302</td><td>2013-12-22</td><td>Meles meles</td><td>SAMEA8955404</td><td>ERR6201847</td></tr><tr><td>D-14-01074</td><td>14Z003304</td><td>2014-01-13</td><td>Meles meles</td><td>SAMEA8955405</td><td>ERR6201848</td></tr><tr><td>D-14-01623</td><td>14Z004188</td><td>2014-03-28</td><td>Bos taurus</td><td>SAMEA8955406</td><td>ERR6201849</td></tr><tr><td>D-14-01757</td><td>14Z004539</td><td>2014-03-13</td><td>Bos taurus</td><td>SAMEA8955407</td><td>ERR6201850</td></tr><tr><td>D-14-01758</td><td>14Z004540</td><td>2014-03-12</td><td>Bos taurus</td><td>SAMEA8955408</td><td>ERR6201851</td></tr><tr><td>D-14-02086</td><td>14Z005338</td><td>2014-04-28</td><td>Bos taurus</td><td>SAMEA8955409</td><td>ERR6201852</td></tr><tr><td>D-14-02099</td><td>14Z005351</td><td>2014-04-23</td><td>Bos taurus</td><td>SAMEA8955410</td><td>ERR6201853</td></tr><tr><td>D-14-02112</td><td>14Z005364</td><td>2014-04-23</td><td>Bos taurus</td><td>SAMEA8955411</td><td>ERR6201854</td></tr><tr><td>D-14-02116</td><td>14Z005365</td><td>2014-04-09</td><td>Bos taurus</td><td>SAMEA8955412</td><td>ERR6201855</td></tr><tr><td>D-14-02259</td><td>14Z005576</td><td>2014-05-06</td><td>Meles meles</td><td>SAMEA8955413</td><td>ERR6201857</td></tr><tr><td>D-14-02574</td><td>14Z006089</td><td>2014-05-25</td><td>Meles meles</td><td>SAMEA8955414</td><td>ERR6201858</td></tr><tr><td>D-14-02811</td><td>14Z006436</td><td>2014-06-23</td><td>Bos taurus</td><td>SAMEA8955415</td><td>ERR6201860</td></tr><tr><td>D-14-02871</td><td>14Z006581</td><td>2014-07-24</td><td>Bos taurus</td><td>SAMEA8955416</td><td>ERR6201861</td></tr><tr><td>D-14-03006</td><td>14Z006798</td><td>2014-06-05</td><td>Meles meles</td><td>SAMEA8955417</td><td>ERR6201862</td></tr><tr><td>D-14-03118</td><td>14Z007047</td><td>2014-06-05</td><td>Meles meles</td><td>SAMEA8955419</td><td>ERR6201863</td></tr><tr><td>D-14-03380</td><td>14Z007947</td><td>2014-09-10</td><td>Bos taurus</td><td>SAMEA8955420</td><td>ERR6201864</td></tr><tr><td>D-14-03463</td><td>14Z008062</td><td>2014-09-10</td><td>Bos taurus</td><td>SAMEA8955421</td><td>ERR6201865</td></tr><tr><td>D-14-03644</td><td>14Z008376</td><td>2014-09-17</td><td>Bos taurus</td><td>SAMEA8955422</td><td>ERR6201866</td></tr><tr><td>D-14-03697</td><td>14Z008498</td><td>2014-08-24</td><td>Sus scrofa</td><td>SAMEA8955423</td><td>ERR6201868</td></tr><tr><td>D-14-04057</td><td>14Z009345</td><td>2014-09-21</td><td>Sus scrofa</td><td>SAMEA8955424</td><td>ERR6201869</td></tr><tr><td>D-15-00525</td><td>15Z001806</td><td>2015-01-14</td><td>Bos taurus</td><td>SAMEA8955425</td><td>ERR6201870</td></tr><tr><td>D-15-00613</td><td>15Z002145</td><td>2015-01-23</td><td>Bos taurus</td><td>SAMEA8955426</td><td>ERR6201871</td></tr><tr><td>D-15-00828</td><td>15Z002493</td><td>2015-02-13</td><td>Bos taurus</td><td>SAMEA8955427</td><td>ERR6201872</td></tr><tr><td>D-15-00970</td><td>15Z002980</td><td>2015-02-17</td><td>Bos taurus</td><td>SAMEA8955428</td><td>ERR6201873</td></tr><tr><td>D-15-01174</td><td>15Z004066</td><td>2015-03-06</td><td>Bos taurus</td><td>SAMEA8955429</td><td>ERR6201874</td></tr><tr><td>D-15-01176</td><td>15Z004068</td><td>2015-03-20</td><td>Bos taurus</td><td>SAMEA8955430</td><td>ERR6201875</td></tr><tr><td>D-15-01301</td><td>15Z004334</td><td>2015-01-25</td><td>Meles meles</td><td>SAMEA8955431</td><td>ERR6201877</td></tr><tr><td>D-15-01302</td><td>15Z004335</td><td>2015-02-27</td><td>Bos taurus</td><td>SAMEA8955432</td><td>ERR6201878</td></tr><tr><td>D-15-01311</td><td>15Z004344</td><td>2015-03-30</td><td>Meles meles</td><td>SAMEA8955433</td><td>ERR6201879</td></tr><tr><td>D-15-01312</td><td>15Z004345</td><td>2015-04-15</td><td>Bos taurus</td><td>SAMEA8955434</td><td>ERR6201881</td></tr><tr><td>D-15-01346</td><td>15Z004377</td><td>2015-03-12</td><td>Meles meles</td><td>SAMEA8955435</td><td>ERR6201882</td></tr><tr><td>D-15-01460</td><td>15Z004939</td><td>2015-03-13</td><td>Meles meles</td><td>SAMEA8955436</td><td>ERR6201883</td></tr><tr><td>D-15-01557</td><td>15Z005359</td><td>2015-03-22</td><td>Meles meles</td><td>SAMEA8955437</td><td>ERR6201884</td></tr><tr><td>D-15-01665</td><td>15Z005833</td><td>2015-04-17</td><td>Bos taurus</td><td>SAMEA8955438</td><td>ERR6201885</td></tr></tbody></table>

(Continues)
<table><caption>TABLE A1 (Continued)</caption><thead><tr><td>DA</td><td>Sample</td><td>Collection date</td><td>Host scientific name</td><td>Biosample</td><td>Accession</td></tr></thead><tbody><tr><td>D-15-01667</td><td>15Z005835</td><td>2013-05-21</td><td>Meles meles</td><td>SAMEA8955439</td><td>ERR6201886</td></tr><tr><td>D-15-01759</td><td>15Z006243</td><td>2015-05-18</td><td>Bos taurus</td><td>SAMEA8955440</td><td>ERR6201887</td></tr><tr><td>D-15-01919</td><td>15Z006492</td><td>2015-04-08</td><td>Meles meles</td><td>SAMEA8955441</td><td>ERR6201888</td></tr><tr><td>D-15-01922</td><td>15Z006495</td><td>2015-05-18</td><td>Bos taurus</td><td>SAMEA8955442</td><td>ERR6201889</td></tr><tr><td>D-15-01960</td><td>15Z006541</td><td>2015-05-12</td><td>Meles meles</td><td>SAMEA8955443</td><td>ERR6201890</td></tr><tr><td>D-15-02070</td><td>15Z006760</td><td>2015-06-17</td><td>Meles meles</td><td>SAMEA8955444</td><td>ERR6201891</td></tr><tr><td>D-15-03139</td><td>15Z010210</td><td>2015-10-14</td><td>Bos taurus</td><td>SAMEA8955445</td><td>ERR6201893</td></tr><tr><td>D-15-03409</td><td>15Z010962</td><td>2015-09-13</td><td>Sus scrofa</td><td>SAMEA8955446</td><td>ERR6201894</td></tr><tr><td>D-16-00172</td><td>16Z000595</td><td>2015-12-15</td><td>Bos taurus</td><td>SAMEA8955447</td><td>ERR6201895</td></tr><tr><td>D-16-00301</td><td>16Z000964</td><td>2015-12-13</td><td>Sus scrofa</td><td>SAMEA8955448</td><td>ERR6201896</td></tr><tr><td>D-16-00302</td><td>16Z000965</td><td>2015-12-26</td><td>Sus scrofa</td><td>SAMEA8955449</td><td>ERR6201897</td></tr><tr><td>D-16-00361</td><td>16Z001097</td><td>2015-12-15</td><td>Bos taurus</td><td>SAMEA8955450</td><td>ERR6201898</td></tr><tr><td>D-16-01021</td><td>16Z002887</td><td>2016-02-15</td><td>Bos taurus</td><td>SAMEA8955451</td><td>ERR6201899</td></tr><tr><td>D-16-01467</td><td>16Z004107</td><td>2016-03-11</td><td>Bos taurus</td><td>SAMEA8955452</td><td>ERR6201901</td></tr><tr><td>D-16-01566</td><td>16Z004271</td><td>2016-03-25</td><td>Bos taurus</td><td>SAMEA8955453</td><td>ERR6201902</td></tr><tr><td>D-16-01567</td><td>16Z004272</td><td>2016-03-20</td><td>Meles meles</td><td>SAMEA8955454</td><td>ERR6201903</td></tr><tr><td>D-16-01569</td><td>16Z004274</td><td>2016-03-10</td><td>Bos taurus</td><td>SAMEA8955455</td><td>ERR6201906</td></tr><tr><td>D-16-01664</td><td>16Z004560</td><td>2016-03-22</td><td>Bos taurus</td><td>SAMEA8955456</td><td>ERR6201907</td></tr><tr><td>D-16-01669</td><td>16Z004566</td><td>2016-03-08</td><td>Bos taurus</td><td>SAMEA8955457</td><td>ERR6201908</td></tr><tr><td>D-16-01672</td><td>16Z004569</td><td>2016-02-17</td><td>Meles meles</td><td>SAMEA8955458</td><td>ERR6201909</td></tr><tr><td>D-16-01723</td><td>16Z004699</td><td>2016-03-08</td><td>Meles meles</td><td>SAMEA8955459</td><td>ERR6201910</td></tr><tr><td>D-16-01847</td><td>16Z005204</td><td>2016-04-12</td><td>Bos taurus</td><td>SAMEA8955460</td><td>ERR6201915</td></tr><tr><td>D-16-01848</td><td>16Z005206</td><td>2016-04-05</td><td>Bos taurus</td><td>SAMEA8955461</td><td>ERR6201937</td></tr><tr><td>D-16-01849</td><td>16Z005207</td><td>2016-04-06</td><td>Bos taurus</td><td>SAMEA8955462</td><td>ERR6201945</td></tr><tr><td>D-16-01851</td><td>16Z005209</td><td>2016-03-25</td><td>Meles meles</td><td>SAMEA8955463</td><td>ERR6201960</td></tr><tr><td>D-16-01852</td><td>16Z005210</td><td>2016-04-11</td><td>Bos taurus</td><td>SAMEA8955464</td><td>ERR6201972</td></tr><tr><td>D-16-02216</td><td>16Z006257</td><td>2017-03-27</td><td>Meles meles</td><td>SAMEA8955465</td><td>ERR6201978</td></tr><tr><td>D-16-02405</td><td>16Z006707</td><td>2016-05-27</td><td>Bos taurus</td><td>SAMEA8955466</td><td>ERR6201986</td></tr><tr><td>D-16-02601</td><td>16Z007121</td><td>2016-05-10</td><td>Meles meles</td><td>SAMEA8955467</td><td>ERR6201998</td></tr><tr><td>D-16-02833</td><td>16Z007476</td><td>2016-06-20</td><td>Bos taurus</td><td>SAMEA8955468</td><td>ERR6202004</td></tr><tr><td>D-16-02834</td><td>16Z007477</td><td>2016-06-24</td><td>Bos taurus</td><td>SAMEA8955469</td><td>ERR6202040</td></tr><tr><td>D-16-02835</td><td>16Z007478</td><td>2016-06-27</td><td>Bos taurus</td><td>SAMEA8955470</td><td>ERR6202401</td></tr><tr><td>D-16-02836</td><td>16Z007479</td><td>2016-07-01</td><td>Bos taurus</td><td>SAMEA8955471</td><td>ERR6203392</td></tr><tr><td>D-16-02924</td><td>16Z007681</td><td>2016-06-20</td><td>Bos taurus</td><td>SAMEA8955472</td><td>ERR6204677</td></tr><tr><td>D-16-03225</td><td>16Z008197</td><td>2016-07-21</td><td>Bos taurus</td><td>SAMEA8955473</td><td>ERR6206178</td></tr><tr><td>D-16-03320</td><td>16Z008443</td><td>2016-07-17</td><td>Meles meles</td><td>SAMEA8955474</td><td>ERR6208726</td></tr><tr><td>D-17-00049</td><td>17Z000058</td><td>2016-12-01</td><td>Bos taurus</td><td>SAMEA8955475</td><td>ERR6208976</td></tr><tr><td>D-17-00050</td><td>17Z000059</td><td>2016-12-01</td><td>Bos taurus</td><td>SAMEA8955476</td><td>ERR6209044</td></tr><tr><td>D-17-00051</td><td>17Z000060</td><td>2016-12-01</td><td>Bos taurus</td><td>SAMEA8955477</td><td>ERR6209102</td></tr><tr><td>D-17-00055</td><td>17Z000064</td><td>2016-11-01</td><td>Sus scrofa</td><td>SAMEA8955478</td><td>ERR6209192</td></tr></tbody></table>
<table><caption>TABLE A1 (Continued)</caption><thead><tr><td>DA</td><td>Sample</td><td>Collection date</td><td>Host scientific name</td><td>Biosample</td><td>Accession</td></tr></thead><tbody><tr><td>D-17-00269</td><td>17Z000349</td><td>2016-12-06</td><td>Bos taurus</td><td>SAMEA8955479</td><td>ERR6209278</td></tr><tr><td>D-17-00270</td><td>17Z000350</td><td>2016-12-06</td><td>Bos taurus</td><td>SAMEA8955480</td><td>ERR6209301</td></tr><tr><td>D-17-00270</td><td>17Z000351</td><td>2016-12-06</td><td>Bos taurus</td><td>SAMEA8955481</td><td>ERR6209309</td></tr><tr><td>D-17-00540</td><td>17Z000898</td><td>2016-12-01</td><td>Bos taurus</td><td>SAMEA8955482</td><td>ERR6209323</td></tr><tr><td>D-17-00563</td><td>17Z000914</td><td>2017-01-09</td><td>Bos taurus</td><td>SAMEA8955483</td><td>ERR6209335</td></tr><tr><td>D-17-00641</td><td>17Z001046</td><td>2017-01-20</td><td>Bos taurus</td><td>SAMEA8955484</td><td>ERR6209343</td></tr><tr><td>D-17-00710</td><td>17Z001177</td><td>2017-01-09</td><td>Bos taurus</td><td>SAMEA8955485</td><td>ERR6209344</td></tr><tr><td>D-17-00849</td><td>17Z001576</td><td>2017-02-03</td><td>Bos taurus</td><td>SAMEA8955486</td><td>ERR6209345</td></tr><tr><td>D-17-00850</td><td>17Z001577</td><td>2017-02-03</td><td>Bos taurus</td><td>SAMEA8955487</td><td>ERR6209346</td></tr><tr><td>D-17-00851</td><td>17Z001578</td><td>2017-02-02</td><td>Bos taurus</td><td>SAMEA8955488</td><td>ERR6209347</td></tr><tr><td>D-17-00854</td><td>17Z001581</td><td>2017-02-02</td><td>Bos taurus</td><td>SAMEA8955489</td><td>ERR6209348</td></tr><tr><td>D-17-01066</td><td>17Z002146</td><td>2017-02-07</td><td>Bos taurus</td><td>SAMEA8955490</td><td>ERR6209349</td></tr><tr><td>D-17-01071</td><td>17Z002151</td><td>2017-02-10</td><td>Bos taurus</td><td>SAMEA8955491</td><td>ERR6209350</td></tr><tr><td>D-17-01146</td><td>17Z002194</td><td>2017-02-15</td><td>Bos taurus</td><td>SAMEA8955492</td><td>ERR6209351</td></tr><tr><td>D-17-01148</td><td>17Z002196</td><td>2017-02-09</td><td>Bos taurus</td><td>SAMEA8955493</td><td>ERR6209352</td></tr><tr><td>D-17-01149</td><td>17Z002197</td><td>2017-02-10</td><td>Bos taurus</td><td>SAMEA8955494</td><td>ERR6209353</td></tr><tr><td>D-17-01150</td><td>17Z002198</td><td>2017-02-07</td><td>Bos taurus</td><td>SAMEA8955495</td><td>ERR6209422</td></tr><tr><td>D-17-01286</td><td>17Z002553</td><td>2017-01-18</td><td>Bos taurus</td><td>SAMEA8955496</td><td>ERR6209423</td></tr><tr><td>D-17-01306</td><td>17Z002558</td><td>2017-02-27</td><td>Bos taurus</td><td>SAMEA8955497</td><td>ERR6209424</td></tr><tr><td>D-17-01412</td><td>17Z002740</td><td>2017-02-07</td><td>Bos taurus</td><td>SAMEA8955498</td><td>ERR6209443</td></tr><tr><td>D-17-01416</td><td>17Z002744</td><td>2017-03-08</td><td>Bos taurus</td><td>SAMEA8955499</td><td>ERR6209444</td></tr><tr><td>D-17-01517</td><td>17Z002865</td><td>2017-03-23</td><td>Bos taurus</td><td>SAMEA8955500</td><td>ERR6209445</td></tr><tr><td>D-17-01610</td><td>17Z003109</td><td>2017-03-01</td><td>Bos taurus</td><td>SAMEA8955501</td><td>ERR6209491</td></tr><tr><td>D-17-01611</td><td>17Z003110</td><td>2017-02-10</td><td>Bos taurus</td><td>SAMEA8955502</td><td>ERR6209492</td></tr><tr><td>D-17-01695</td><td>17Z003235</td><td>2017-04-05</td><td>Bos taurus</td><td>SAMEA8955503</td><td>ERR6209493</td></tr><tr><td>D-17-01696</td><td>17Z003236</td><td>2017-03-28</td><td>Bos taurus</td><td>SAMEA8955504</td><td>ERR6209524</td></tr><tr><td>D-17-01698</td><td>17Z003238</td><td>2017-03-21</td><td>Bos taurus</td><td>SAMEA8955505</td><td>ERR6209526</td></tr><tr><td>D-17-01705</td><td>17Z003245</td><td>2017-03-27</td><td>Bos taurus</td><td>SAMEA8955506</td><td>ERR6209527</td></tr><tr><td>D-17-01929</td><td>17Z003568</td><td>2017-03-30</td><td>Bos taurus</td><td>SAMEA8955507</td><td>ERR6209542</td></tr><tr><td>D-17-01931</td><td>17Z003570</td><td>2017-04-13</td><td>Bos taurus</td><td>SAMEA8955508</td><td>ERR6209573</td></tr><tr><td>D-17-02078</td><td>17Z004077</td><td>2017-04-21</td><td>Bos taurus</td><td>SAMEA8955509</td><td>ERR6209576</td></tr><tr><td>D-17-02080</td><td>17Z004080</td><td>2017-04-20</td><td>Bos taurus</td><td>SAMEA8955510</td><td>ERR6210127</td></tr><tr><td>D-17-02089</td><td>17Z004088</td><td>2017-04-20</td><td>Bos taurus</td><td>SAMEA8955511</td><td>ERR6210128</td></tr><tr><td>D-17-02090</td><td>17Z004128</td><td>2017-04-19</td><td>Bos taurus</td><td>SAMEA8955512</td><td>ERR6210129</td></tr><tr><td>D-17-02088</td><td>17Z004275</td><td>2017-04-20</td><td>Bos taurus</td><td>SAMEA8955513</td><td>ERR6210130</td></tr><tr><td>D-17-02260</td><td>17Z004423</td><td>2017-03-29</td><td>Bos taurus</td><td>SAMEA8955514</td><td>ERR6210147</td></tr><tr><td>D-17-02262</td><td>17Z004425</td><td>2017-04-19</td><td>Bos taurus</td><td>SAMEA8955515</td><td>ERR6210148</td></tr><tr><td>D-17-02264</td><td>17Z004428</td><td>2017-05-02</td><td>Bos taurus</td><td>SAMEA8955516</td><td>ERR6210149</td></tr><tr><td>D-17-02335</td><td>17Z004525</td><td>2017-04-24</td><td>Bos taurus</td><td>SAMEA8955517</td><td>ERR6210150</td></tr><tr><td>D-17-02438</td><td>17Z005047</td><td>2017-04-26</td><td>Bos taurus</td><td>SAMEA8955518</td><td>ERR6210158</td></tr></tbody></table>

(Continues)
<table><caption>TABLE A1 (Continued)</caption><thead><tr><td>DA</td><td>Sample</td><td>Collection date</td><td>Host scientific name</td><td>Biosample</td><td>Accession</td></tr></thead><tbody><tr><td>D-17-02580</td><td>17Z005388</td><td>2017-04-19</td><td>Bos taurus</td><td>SAMEA8955519</td><td>ERR6210161</td></tr><tr><td>D-17-02586</td><td>17Z005394</td><td>2017-05-09</td><td>Bos taurus</td><td>SAMEA8955520</td><td>ERR6210169</td></tr><tr><td>D-17-02587</td><td>17Z005395</td><td>2017-05-18</td><td>Bos taurus</td><td>SAMEA8955521</td><td>ERR6210175</td></tr><tr><td>D-17-02698</td><td>17Z005495</td><td>2017-05-22</td><td>Bos taurus</td><td>SAMEA8955522</td><td>ERR6210188</td></tr><tr><td>D-17-02699</td><td>17Z005496</td><td>2017-05-22</td><td>Bos taurus</td><td>SAMEA8955523</td><td>ERR6210195</td></tr><tr><td>D-17-02775</td><td>17Z005671</td><td>2017-05-30</td><td>Bos taurus</td><td>SAMEA8955524</td><td>ERR6210201</td></tr><tr><td>D-17-02902</td><td>17Z005849</td><td>2017-06-23</td><td>Bos taurus</td><td>SAMEA8955525</td><td>ERR6210207</td></tr><tr><td>D-17-03220</td><td>17Z006451</td><td>2017-03-13</td><td>Meles meles</td><td>SAMEA8955526</td><td>ERR6210212</td></tr><tr><td>D-17-03217</td><td>17Z006454</td><td>2017-06-12</td><td>Bos taurus</td><td>SAMEA8955527</td><td>ERR6210216</td></tr><tr><td>D-17-03251</td><td>17Z006498</td><td>2017-07-17</td><td>Bos taurus</td><td>SAMEA8955528</td><td>ERR6210223</td></tr><tr><td>D-17-03430</td><td>17Z006770</td><td>2017-07-17</td><td>Bos taurus</td><td>SAMEA8955529</td><td>ERR6210230</td></tr><tr><td>D-17-03431</td><td>17Z006771</td><td>2017-07-15</td><td>Meles meles</td><td>SAMEA8955530</td><td>ERR6210236</td></tr><tr><td>D-17-03656</td><td>17Z007283</td><td>2017-08-06</td><td>Bos taurus</td><td>SAMEA8955531</td><td>ERR6210252</td></tr><tr><td>D-17-03910</td><td>17Z007631</td><td>2017-08-24</td><td>Bos taurus</td><td>SAMEA8955532</td><td>ERR6210258</td></tr><tr><td>D-17-03912</td><td>17Z007633</td><td>2017-08-10</td><td>Bos taurus</td><td>SAMEA8955533</td><td>ERR6210259</td></tr><tr><td>D-17-03913</td><td>17Z007634</td><td>2017-08-10</td><td>Bos taurus</td><td>SAMEA8955534</td><td>ERR6210266</td></tr><tr><td>D-17-04408</td><td>17Z008463</td><td>2017-08-24</td><td>Bos taurus</td><td>SAMEA8955535</td><td>ERR6210267</td></tr><tr><td>D-17-04412</td><td>17Z008467</td><td>2017-07-09</td><td>Meles meles</td><td>SAMEA8955536</td><td>ERR6210268</td></tr><tr><td>D-17-04421</td><td>17Z008476</td><td>2017-08-14</td><td>Meles meles</td><td>SAMEA8955537</td><td>ERR6210269</td></tr><tr><td>D-17-04423</td><td>17Z008478</td><td>2017-04-27</td><td>Meles meles</td><td>SAMEA8955538</td><td>ERR6210270</td></tr><tr><td>D-17-04936</td><td>17Z010330</td><td>2017-10-01</td><td>Sus scrofa</td><td>SAMEA8955539</td><td>ERR6210271</td></tr><tr><td>D-17-05074</td><td>17Z010552</td><td>2017-10-29</td><td>Sus scrofa</td><td>SAMEA8955540</td><td>ERR6210272</td></tr></tbody></table>
<table><caption>TABLE A2 Model selection</caption><thead><tr><th>Model</th><th>Substitution model</th><th>Molecular clock model</th><th>Marginal likelihood</th><th>SD</th><th>BF</th><th>2√SD<sub>2</sub><sup>2</sup> + SD<sub>1</sub><sup>2</sup></th><th>Reference model</th></tr></thead><tbody><tr><td colspan="8"><strong>CdO</strong></td></tr><tr><td>M1CdO</td><td>JC69</td><td>Strict</td><td>-24554</td><td>0.09</td><td>-</td><td>-</td><td>-</td></tr><tr><td>M2CdO</td><td>HKY</td><td>Strict</td><td>-24225</td><td>1.00</td><td>329</td><td>2.00</td><td>M1CdO</td></tr><tr><td>M3CdO</td><td>GTR</td><td>Strict</td><td>Does not converge</td><td></td><td></td><td></td><td></td></tr><tr><td>M4CdO</td><td>HKY</td><td>Uncorrelated lognormal</td><td>-24222</td><td>0.84</td><td>3</td><td>2.61</td><td>M2CdO</td></tr><tr><td>M5CdO</td><td>HKY</td><td>Uncorrelated exponential</td><td>-24223</td><td>0.91</td><td>-1</td><td>2.48</td><td>M4CdO</td></tr><tr><td colspan="8"><strong>NAq</strong></td></tr><tr><td>M1NAq</td><td>JC69</td><td>Strict</td><td>-87642</td><td>0.18</td><td>-</td><td>-</td><td>-</td></tr><tr><td>M2NAq</td><td>HKY</td><td>Strict</td><td>-3400</td><td>9.41</td><td>84242</td><td>18.82</td><td>M1NAq</td></tr><tr><td>M3NAq</td><td>GTR</td><td>Strict</td><td>Does not converge</td><td></td><td></td><td></td><td></td></tr><tr><td>M4NAq</td><td>HKY</td><td>Uncorrelated lognormal</td><td>-3350</td><td>8.87</td><td>50</td><td>25.86</td><td>M2NAq</td></tr><tr><td>M5NAq</td><td>HKY</td><td>Uncorrelated exponential</td><td>-87616</td><td>1.04</td><td>-84266</td><td>17.9</td><td>M4NAq</td></tr></tbody></table>

The marginal likelihood was computed with the Nested Sampling algorithm implemented in the BEAST NS package, with 10 particles. For all models, an unstructured coalescent population model was used. BF stands for Bayes Factor and was computed as ML1-ML2 if ML1 is the marginal likelihood of model 1 and ML2 marginal likelihood of model 2. To distinguish between both MLs, the difference should be greater than $2\sqrt{SD_2^2 - SD_1^2}$, where SD1 stands for the standard deviation of model 1 and SD2 standard deviation of model 2.

Abbreviations: BF, Bayes factor; ML1, marginal likelihood of Model 1; ML2, marginal likelihood of Model 2; SD1, standard deviation of Model 1; SD2, SD2 standard deviation of Model 2.