# Molecular Epidemiology of Rabies in Wild Canidae in Tunisia

Zied Bouslama, Habib Kharmachi, Nourhene Basdouri, Jihen Ben Salem, Samia Ben Maiez, Mariem Handous, Mohamed Saadi, Abdeljalil Ghram, Imed Turki

► To cite this version:

Zied Bouslama, Habib Kharmachi, Nourhene Basdouri, Jihen Ben Salem, Samia Ben Maiez, et al. Molecular Epidemiology of Rabies in Wild Canidae in Tunisia. Viruses, 2021, 13 (12), pp.2473. 10.3390/v13122473. pasteur-03545143

HAL Id: pasteur-03545143

https://riip.hal.science/pasteur-03545143

Submitted on 27 Jan 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
Article

# Molecular Epidemiology of Rabies in Wild Canidae in Tunisia

Zied Bouslama<sup>1,2,*</sup>, Habib Kharmachi<sup>1</sup>, Nourhene Basdouri<sup>1</sup>, Jihen Ben Salem<sup>1</sup>, Samia Ben Maiez<sup>1</sup>, Mariem Handous<sup>1 ID</sup>, Mohamed Saadi<sup>1</sup>, Abdeljalil Ghram<sup>3</sup> and Imed Turki<sup>4</sup>

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

check for
updates

**Citation:** Bouslama, Z.; Kharmachi, H.; Basdouri, N.; Ben Salem, J.; Ben Maiez, S.; Handous, M.; Saadi, M.; Ghram, A.; Turki, I. Molecular Epidemiology of Rabies in Wild Canidae in Tunisia. *Viruses* **2021**, *13*, 2473. https://doi.org/10.3390/v13122473

Academic Editors: Paola De Benedictis, Wanda Markotter and Stefania Leopardi

Received: 30 September 2021

Accepted: 2 November 2021

Published: 10 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 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 (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Rabies is a viral zoonosis that is transmissible to humans via domestic and wild animals. There are two epidemiological cycles for rabies, the urban and the sylvatic cycles. In an attempt to study the epidemiological role of wild canidae in rabies transmission, the present study aimed to analyze the genetic characteristics of virus isolates and confirm prior suggestions that rabies is maintained through a dog reservoir in Tunisia. Virus strains isolated from wild canidae were subject to viral sequencing, and Bayesian phylogenetic analysis was performed using Beast2 software. Essentially, the virus strains isolated from wild canidae belonged to the Africa-1 clade, which clearly diverges from fox-related strains. Our study also demonstrated that genetic characteristics of the virus isolates were not as distinct as could be expected if a wild reservoir had already existed. On the contrary, the geographic landscape is responsible for the genetic diversity of the virus. The landscape itself could have also acted as a natural barrier to the spread of the virus.

**Keywords:** rabies; wild canidae; Tunisia; phylogeny

## 1. Introduction

The rabies virus belongs to the *Lyssavirus* genus. It is a negative-sense single-stranded RNA virus with a genome size of approximately 12 kb. It contains five genes, namely the nucleoprotein (N), the phosphoprotein (P), the matrix protein (M), the glycoprotein (G) and the large protein or polymerase (L), and a non-coding region [1,2]. Among these, the nucleoprotein (N) is the most conserved gene, and the phosphoprotein (P) and glycoprotein (G) are the most variable ones [3]. To date, 17 different viral species have been identified [4]. Rabies is a viral infectious zoonotic disease endemic in many regions of the world, particularly Asia, Africa, and South America. In Tunisia, rabies is enzootic and endemic. Rabies virus (RABV) can be associated with a variety of mammalian hosts that maintain independent epidemiological cycles within geographical locations, resulting in the differentiation of the virus into several lineages [5,6]. RABV strains that circulate in dogs (*Canis lupus familiaris*) are responsible for more than 99% of all human cases worldwide [7,8]. However, different wild mammals have been identified as RABV reservoirs, such as bat-eared foxes, black-backed jackals, and yellow mongooses in Southern Africa and red foxes in Europe and the Middle East. In different parts of the world, these species maintain sylvatic transmission cycles that can be independent from the primary urban cycle in which the dog is the main reservoir [9,10].

Rabies viruses isolated from foxes in Europe belong to the same cosmopolitan lineage associated with dogs. However, it has unique genetic characteristics and thus forms a separate clade [11]. In North Africa, and particularly in Tunisia, the molecular characterization

<footnote n="1">¹ Laboratory for Rabies Diagnostics, Institute Pasteur of Tunis, Belvedere, Tunis 1002, Tunisia; h.kharmachi@yahoo.fr (H.K.); nourhene.basdouri@yahoo.com (N.B.); jihen.bensalem@pasteur.tn (J.B.S.); samia.benmaiez@pasteur.tn (S.B.M.); mariem handous@pasteur.tn (M.H.); mohamed.saadi@pasteur.tn (M.S.)</footnote>

<footnote n="2">² Faculty of Sciences, Université Tunis El Manar, Tunis 2092, Tunisia</footnote>

<footnote n="3">³ Laboratory of Epidemiology and Veterinary Microbiology, LR 16 IPT 03, Institut Pasteur de Tunis, Université Tunis El Manar, Tunis 2092, Tunisia; abdeljalil.ghram@pasteur.tn</footnote>

<footnote n="4">⁴ Service des Maladies Contagieuses, Ecole Nationale de Médecine Vétérinaire-Sidi Thabet, Université Manouba, Sidi Thabet 2020, Tunisia; turkiimed58@yahoo.fr</footnote>

<footnote>* Correspondence: ziedbouslama80@gmail.com</footnote>
of rabies virus isolated from dogs [12], domestic animals [13] and humans [14] allowed for the sole identification of canine rabies clade Africa 1-a, belonging to the cosmopolitan lineage, showing that dogs are the main reservoir species in the area [12,15]. However, no study has yet been conducted on rabies viruses isolated from wild animals in Tunisia to investigate the existence of a sylvatic circle. The laboratory for rabies diagnostics at Institut Pasteur of Tunis (IPT) is the national reference laboratory and the only one authorized for rabies diagnostics in Tunisia. Yearly, an increasing number of samples are tested, exceeding 1000 samples per year since 2014. According to the national commission of rabies control, samples submitted to the laboratory are mostly from dogs identified as suspected cases, while a maximum of three samples per year are collected from wildlife [16]. A suspected case is identified as an animal that died after presenting neurological symptoms or after having bitten a human. In turn, among animal rabies cases reported yearly in Tunisia, dogs represent 60%, while wild animals only represent 0.3 % [15]. For example, among 383 rabies cases reported in Tunisia in 2018, 236 cases were identified in dogs and 3 were human cases. It is known that a certain number of cases can be recorded in wild canids even when the country adheres to mass dog vaccination. Indeed, rabies cases are broadly recorded in wild canids among North Africa, including in Tunisia, even if these host species are poorly studied for their potential role in rabies transmission. On the contrary, studies performed in Middle Eastern countries such as Iraq, Jordan, and Syria confirmed the coexistence of urban and sylvatic rabies cycles [17]. The objective of this study was the molecular characterization of rabies viruses isolated in Tunisia from wild canids, specifically foxes and jackals, in order to compare them with variants that were previously reported from dogs in the area. In addition, we also analyzed and discussed the temporal evolution of variants in comparison with North African isolates.

## 2. Materials and Methods

### 2.1. Origin of Virus Isolates

Animals suspected to be dead from rabies infection were collected by the veterinary services and sent to the IPT laboratory to undergo necropsy, brain extraction and rabies testing using a FAT test, as recommended by the OIE and the World Health Organization (WHO) [18,19]. All samples that tested positive for rabies were stored in the rabies laboratory biobank. In total, we selected 11 brain samples collected from wild animals, namely five foxes and six jackals, and 14 brain samples isolated from domestic dogs in the same location and during the same period of time (same delegation, same month), as shown in Table 1. These samples were characterized and integrated into the phylogenetic analysis together with additional sequences downloaded from GenBank.

<table><caption>**Table 1.** List of virus isolates selected for molecular characterization. Information on location and sampling date, corresponding to the animal's death, were extracted from the laboratory's database.</caption><thead><tr><th>Isolate</th><th>Animal Species</th><th>Collection Date</th><th>Location (Governorate)</th><th>Location (Delegation)</th><th>Latitude</th><th>Longitude</th><th>Nucleoprotein Sequence Accession Number</th><th>Phosphoprotein Sequence Accession Number</th></tr></thead><tbody><tr><td>15841</td><td>Fox</td><td>10 April 2017</td><td>Nabeul</td><td>Nabeul</td><td>36.4911995</td><td>10.6626997</td><td>OK275685</td><td>-</td></tr><tr><td>15397</td><td>Fox</td><td>1 November 2016</td><td>Nabeul</td><td>El Mida</td><td>36.773047</td><td>10.7656</td><td>OK275686</td><td>OK275684</td></tr><tr><td>11284</td><td>Fox</td><td>17 November 2012</td><td>Ben Arous</td><td>Mhamdia</td><td>36.647136</td><td>10.066214</td><td>OK275688</td><td>-</td></tr><tr><td>9484</td><td>Fox</td><td>10 October 2008</td><td>Medenine</td><td>Sidi Makhlouf</td><td>33.499538</td><td>10.473707</td><td>OK275694</td><td>OK275673</td></tr><tr><td>9393</td><td>Fox</td><td>22 July 2008</td><td>Gabes</td><td>Mareth</td><td>33.663899</td><td>10.348131</td><td>OK275692</td><td>OK275670</td></tr><tr><td>9413</td><td>Jackal</td><td>10 August 2008</td><td>Medenine</td><td>Beni Khdache</td><td>33.436471</td><td>10.206275</td><td>OK275691</td><td>OK275671</td></tr><tr><td>13618</td><td>Jackal</td><td>16 February 2015</td><td>Kasserine</td><td>Sbeitla</td><td>35.241052</td><td>9.131288</td><td>OK275687</td><td>OK275681</td></tr><tr><td>13119</td><td>Jackal</td><td>17 September 2014</td><td>Siliana</td><td>Kesra</td><td>35.8359985</td><td>9.4758101</td><td>OK275693</td><td>OK275679</td></tr><tr><td>10582</td><td>Jackal</td><td>5 June 2011</td><td>Kasserine</td><td>Thala</td><td>35.557643</td><td>8.680733</td><td>OK275689</td><td>OK275675</td></tr><tr><td>9438</td><td>Jackal</td><td>29 August 2008</td><td>Medenine</td><td>Beni Khdache</td><td>33.436471</td><td>10.206275</td><td>OK275690</td><td>OK275672</td></tr><tr><td>8112</td><td>Jackal</td><td>23 February 2006</td><td>El Kef</td><td>Dahmani</td><td>36.021820</td><td>8.906841</td><td>OK275695</td><td>OK275668</td></tr><tr><td>12596</td><td>Dog</td><td>3 November 2014</td><td>Siliana</td><td>Kesra</td><td>35.8359985</td><td>9.4758101</td><td>-</td><td>OK275678</td></tr></tbody></table>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Isolate</th><th>Animal Species</th><th>Collection Date</th><th>Location (Governorate)</th><th>Location (Delegation)</th><th>Latitude</th><th>Longitude</th><th>Nucleoprotein Sequence Accession Number</th><th>Phosphoprotein Sequence Accession Number</th></tr></thead><tbody><tr><td>14109</td><td>Dog</td><td>24 July 2015</td><td>Kasserine</td><td>Sbeitla</td><td>35.2153015</td><td>9.0153503</td><td>-</td><td>OK275682</td></tr><tr><td>10910</td><td>Dog</td><td>2 May 2012</td><td>Ben Arous</td><td>Fouchana</td><td>36.6683998</td><td>10.2124996</td><td>-</td><td>OK275677</td></tr><tr><td>15341</td><td>Dog</td><td>24 October 2016</td><td>Nabeul</td><td>Menzel Temime</td><td>36.796965</td><td>10.853034</td><td>-</td><td>OK275683</td></tr><tr><td>13239</td><td>Dog</td><td>24 October 2014</td><td>Siliana</td><td>Rouhia</td><td>35.7809982</td><td>9.0813599</td><td>-</td><td>OK275680</td></tr><tr><td>10584</td><td>Dog</td><td>2 June 2011</td><td>Kasserine</td><td>El Hammar</td><td>35.3737984</td><td>9.0671101</td><td>-</td><td>OK275676</td></tr><tr><td>9118</td><td>Dog</td><td>17 Mars 2008</td><td>Medenine</td><td>Ben Guerdane</td><td>33.2186012</td><td>10.7330999</td><td>-</td><td>OK275669</td></tr><tr><td>8064</td><td>Dog</td><td>1 November 2006</td><td>Zaghouan</td><td>Nadhour</td><td>36.0915985</td><td>10.0277004</td><td>-</td><td>OK275667</td></tr><tr><td>9767</td><td>Dog</td><td>2 April 2009</td><td>Medenine</td><td>Medenine Nord</td><td>33.3704987</td><td>10.4882002</td><td>-</td><td>OK275674</td></tr><tr><td>10792</td><td>Dog</td><td>31 January 2011</td><td>Jendouba</td><td>Ouled Mliz</td><td>36.4488983</td><td>8.5616598</td><td>-</td><td>OK275665</td></tr><tr><td>11747</td><td>Dog</td><td>11 June 2013</td><td>Jendouba</td><td>Sloul</td><td>36.7193985</td><td>8.6416702</td><td>-</td><td>OK275664</td></tr><tr><td>12900</td><td>Dog</td><td>13 June 2014</td><td>Jendouba</td><td>Ouechtata</td><td>36.5028000</td><td>8.2497797</td><td>-</td><td>OK275663</td></tr><tr><td>13961</td><td>Dog</td><td>5 June 2015</td><td>El Kef</td><td>Eddir</td><td>36.2057991</td><td>8.7846699</td><td>-</td><td>OK275666</td></tr></tbody></table>

## 2.2. RNA Extraction and Amplification

RNA was extracted using TRIzol reagent (invitrogen®) as described in previous studies [20]. Brain samples were homogenized in 1 mL of TRIzol reagent and 0.2 mL of chloroform. After centrifugation (11,000 rpm during 15 min), the supernatant was mixed with 0.5 mL of isopropanol. After a second centrifugation (11,000 rpm during 10 min), the collected supernatant was discarded, and the pellet was washed using 1ml of ethanol. The tubes were then dried and eluted in 50 µL of RNase free water. In the case of wild canids, we characterized the partial Nucleoprotein gene with a hemi-nested PCR protocol using the SuperScript III One step RT-PCR kit (invitrogen) for the first step and the platinum Taq (Invitrogen) in the second step, following the manufacturer's instructions. A partial sequence of the phosphoprotein gene was characterized for both wild canids and domestic dogs with in one step using an Affinity Script One-Step RT-PCR kit (Agilent) as recommended by the manufacturer. The primers used for both analyses are presented in Table 2.

<table><caption>**Table 2. List of primers.**</caption><thead><tr><th>Primer Name</th><th>Primer Sequence</th><th>Sense</th><th>Target Gene</th><th>Reference</th></tr></thead><tbody><tr><td>Rab N1</td><td>AACACCTCTACAATGGATGCCGACAA</td><td>Forward</td><td>Nucleoprotein</td><td>Nadin Davis S, 1998</td></tr><tr><td>Rab N5</td><td>GGATTGAC (AG)<br>AAGATCTTGCTCAT</td><td>Reverse</td><td>Nucleoprotein</td><td>Nadin Davis S, 1998</td></tr><tr><td>Rab Nfor</td><td>TTGT(AG) GA (TC)<br>CAATATGAGTACAA</td><td>Forward</td><td>Nucleoprotein</td><td>Nadin Davis S, 1998</td></tr><tr><td>Rab Nrev</td><td>CCGGCTCAAACATTCTTCTTA</td><td>Reverse</td><td>Nucleoprotein</td><td>Nadin Davis S, 1998</td></tr><tr><td>P for</td><td>GAACCATCCCAAAYATG</td><td>Forward</td><td>Phosphoprotein</td><td>Wang L. et al., 2013 [21]</td></tr><tr><td>P rev</td><td>CTATCTTGCGCAGAAARTTCAT</td><td>Reverse</td><td>Phosphoprotein</td><td>Wang L. et al., 2013</td></tr></tbody></table>

## 2.3. DNA Sequencing

PCR products were purified with ExoSAP and sequenced using the BigDye Terminator v3.1 cycle sequencing kit (Applied Biosystems, Waltham, MA, USA). Each amplicon was sequenced using one forward and one reverse primer to determine consensus sequences reaching up to 751 bp (N gene) and 1017 bp (P gene), using an ABI PRISM 3500 Genetic Analyzer (Applied Biosystems). Original sequences were trimmed and edited using BioEdit software. Then, they were aligned along with a set of sequences downloaded from GenBank using MAFFT online [22]. All sequences were deposited in GenBank (accession numbers OK275663-OK275695).
## 2.4. Maximum Likelihood Phylogeny

In order to determine taxonomy and clade membership for each virus isolate, we performed maximum likelihood (ML) phylogenetic analysis on the nucleoprotein dataset [23]. We first determined the best-fit substitution model by Smart Model Selection (SMS) [24], and certain GTR model + Gamma 4 categories were selected as the best parameters based on the lower Akaike Information Criterion (AIC). We then implemented ML phylogeny in PhyML with 1000 bootstrap replicates, using PhyML online software available at "http://www.atgc-montpellier.fr/phyml/, accessed on 26 January 2021". The resulting tree was edited and annotated on FigTree 1.7 [25].

## 2.5. Bayesian Phylogeny

To determine the time to common ancestry among RABV clades circulating in Tunisia, we conducted a Bayesian MCMC phylogenetic analysis on the phosphoprotein dataset [26]. Initial substitution rates were estimated for use as tree priors by using Tempest [27]. A substitution model was averaged by using the bModelTest Package in Beast2 [28]. A Bayesian phylogenetic tree was constructed in Beast2 using a TN93 model with a discrete Gamma distribution among variations and a relaxed molecular clock model [29]. A Bayesian skyline plot model was specified as a tree prior. The analysis was set to run for 30 million iterations. ESS values were examined using Tracer software. A Maximum Clade Credibility (MCC) tree was constructed after discarding the first 10% burn in using Tree Annotator V2.6.2. Rates of nucleotide substitution and the time to most recent common ancestor (tMRCA) were estimated. Time tree was edited and visualized using iTOL online software available at: "https://itol.embl.de/, accessed on 18 March 2021". The spatial distribution of variants was edited using the ggmap package v3.0.0 [30].

## 2.6. Reservoir Species Analysis

A pairwise genetic distance matrix was calculated using the 'Ape' package on R software and then visualized using the 'ggplot2' package on R software. The dataset included P sequences obtained from wild canids and dogs sharing similar sampling dates and locations. To identify the most likely animal species reservoir for each Clade, a second analysis was carried out on Beast2 using the animal host species as a discrete trait, as specified in a previous study [16]. The analysis was set to run for 30 million MCMC iterations. An MCC tree with traits was set and annotated using TreeAnnotatorV2.6.2, then visualized and edited in Fig-tree 1.7.

# 3. Results

## 3.1. Virus Taxonomy and Clade Membership

After trimming and alignment, the obtained N gene sequences were up to 481 bp in length. Original N sequences were obtained for six jackals and five foxes. A maximum likelihood inferred tree was divided into six different lineages, identified as Asian, Cosmopolitan, Africa-2, Africa-3, Africa-4, and Arctic-related lineages (Figure 1), all of which are of the RABV virus species. Cosmopolitan lineage was divided into four clades, which were Africa-1a, Africa-1b, Africa-1c, and European Fox-related strains (Figure 1). All Tunisian viruses sequenced in the present study clustered with the clade Africa-1a within the Cosmopolitan lineage, together with the canine-related African rabies strains.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Phylogenetic analysis of rabies virus isolated from wild canidae in Tunisia. Cosmopolitan lineage (highlighted in purple) comprises Africa-1a (light blue), Africa-1b (light green), Africa-1c (orange), and fox-related strains (grey).</figcaption></figure>

## 3.2. *Comparison of Wild Canidae Virus Isolates with Tunisian Canine Variant Strains*

We compared viruses associated with foxes and jackals in Tunisia with variants previously identified in Tunisian dogs, which were downloaded from GenBank. The resulting ML tree is displayed in Figure 2, showing that variants found in Tunisian wild canids cluster together with those identified in dogs. The phylogenetic tree was further divided in two sub-clades that were independent from the host species but related to the area of sampling, identified as northwest (NW) and northeast-center-south (NCS). All viruses isolated from wild canids clustered with high support within the northeast-center-south (NCS) sub-clade, except for the isolate 8112, associated with a jackal found in 2006 in the El Kef governorate (northwest of Tunisia). Another interesting result was that isolate 10582, associated with a jackal that originated from the Kasserine governorate (center west) in 2011, was mostly related with the isolate DogTN/Ks98 (isolated from Kasserine in 1998). Although they belonged to the NCS sub-clade, they did slightly deviate from other NCS isolates. This deviation was supported by a node with high bootstrap values (100%). The old Tunisian isolates, downloaded from GenBank (Tunisia 1986), also belonged to the NW sub-clade.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Phylogenetic analysis of rabies virus isolates in Tunisia.</figcaption></figure>

### 3.3. Reservoir Species Analysis

To better understand virus transmission between dogs and wild canids, we performed a pairwise genetic distance analysis followed by a Bayesian analysis using a dataset composed of sequences of the phosphoprotein gene isolated from wild canids and dogs that shared the same sampling times and locations. Original P sequences were obtained for six jackals, three foxes and fourteen dogs. A pairwise distance matrix, calculated and visualized on R software, showed 100% similarity (0% dissimilarity) between viral sequences associated with animals that shared the same sampling location, irrespective of the animal species (Figure 3). Similarly, the MCC tree built using the host species as a discrete trait (Figure 4) confirmed that variants sharing similar sampling dates and locations cluster together independently of the animal species. The origin (ancestry) of most of the clusters defined within the tree was represented by a canine isolate.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Heatmap representing the pairwise genetic distance between viral sequences associated with wild and domestic animals.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Species tree with traits representing the animal species probably associated with the ancestral node.</figcaption></figure>

3.4. *Time Scale of Rabies Virus Evolution in Tunisia*

In order to estimate the time scale of divergence of the Tunisian rabies variants, we
performed a Bayesian analysis using the phosphoprotein dataset, including sequences from
the study, plus two sequences from Algeria, and eight additional sequences from Tunisia
(Table 3) that were downloaded from GenBank. Molecular clock analyses implemented
in Tempest confirmed a positive correlation between time and substitution events and
estimated a slope rate of 5.4 × 10⁻⁴ substitution/site/year that was thus set as a prior for
the mean clock rate. The MCC tree resulting from evolutionary analyses is displayed in
Figure 5. The overall tree structure confirms that Tunisian isolates are divided into two
phylogenetic groups, NW and NCS, that diverged a few decades ago; the most recent
common ancestor (tMRCA) between these sub-clades was estimated to occur in the year
1957 (Figure 5). The NCS variants are further divided in two clusters, NCS-1 and NCS-2.
RABV variants belonging to the cluster NCS-1 are located in the northeast, the center-
east and the southeast of the country, while cluster NCS-2 includes viruses found in the
northwest and the center-west (Figure 6). Interestingly, the two Algerian variants included
in this study clustered within the cluster NCS-2. The time of the most recent common
ancestor (tMRCA) between NCS's clusters was estimated to occur in the year 1989 (95%
HPD: 1980–2000), around thirty years after the divergence between the two sub-clades. The
virus isolate number 8112, associated with a Jackal in El Kef in 2006, was highly divergent
from all other sequences of the P gene alignment, although it clustered with NW variant
based on phylogenetic analyses performed on the N gene. Its most recent common ancestor
was estimated to occur in the year 1892.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Time scale divergence of the Tunisian rabies phylogenetic variants.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Spatial distribution of branches of the NCS variants.</figcaption></figure>
<table><caption>**Table 3.** List of sequences downloaded from GenBank.</caption><thead><tr><th>Isolate</th><th>GenBank Accession Number</th><th>Host Species</th><th>Location (Country)</th><th>Reference</th></tr></thead><tbody><tr><td>ALG/08-153</td><td>GU798549.1</td><td>Dog</td><td>Algeria</td><td>(Talbi et al., 2010)</td></tr><tr><td>ALG/08-200</td><td>GU798548.1</td><td>Human</td><td>Algeria</td><td>(Talbi et al., 2010)</td></tr><tr><td>8676TUN/1986</td><td>GU798390.1</td><td>Human</td><td>Tunisia</td><td>(Talbi et al., 2010)</td></tr><tr><td>8727TUN/1986</td><td>GU798391.1</td><td>Human</td><td>Tunisia</td><td>(Talbi et al., 2010)</td></tr><tr><td>86075TUN/1986</td><td>GU798389.1</td><td>Human</td><td>Tunisia</td><td>(Talbi et al., 2010)</td></tr><tr><td>86130TUN/1986</td><td>GU798392.1</td><td>Human</td><td>Tunisia</td><td>(Talbi et al., 2010)</td></tr><tr><td>86131TUN/1986</td><td>GU798393.1</td><td>Human</td><td>Tunisia</td><td>(Talbi et al., 2010)</td></tr><tr><td>86127TUN/1986</td><td>GU798395.1</td><td>Human</td><td>Tunisia</td><td>(Talbi et al., 2010)</td></tr><tr><td>86128TUN/1986</td><td>GU798394.1</td><td>Human</td><td>Tunisia</td><td>(Talbi et al., 2010)</td></tr><tr><td>Ariana2</td><td>MK981888.1</td><td>Dog</td><td>Tunisia</td><td>[31] (Bonnaud et al., 2019)</td></tr><tr><td>DogTN/Bj01</td><td>EU643525.1</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Bz98</td><td>EU643527.1</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Bz99</td><td>EU643528.1</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Jd01</td><td>EU643534.1</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Jd96</td><td>EU643535.1</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Jd98</td><td>EU643536.1</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Ar98</td><td>EU643555</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Gb96</td><td>EU643532</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Gb98</td><td>EU643533</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Jd99</td><td>EU643537</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Kf96</td><td>EU643540</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Kr02</td><td>EU643538</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Kr03</td><td>EU643553</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Kr98</td><td>EU643539</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Ks00</td><td>EU643541</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Ks98</td><td>EU643543</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Md95</td><td>EU643545</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Md99</td><td>EU643546</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Nb00</td><td>EU643548</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Ta02</td><td>EU643550</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Tn03</td><td>EU643554</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Tz96</td><td>EU643551</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr><tr><td>DogTN/Zg00</td><td>EU643552</td><td>Dog</td><td>Tunisia</td><td>(Amouri et al., 2011)</td></tr></tbody></table>

4. Discussion

Our study provided a genetic characterization of rabies viruses isolated from wild canids in Tunisia and broadly in North Africa, using a stepwise approach, including viral isolation, sequencing, and phylodynamics reconstruction. The study highlighted interesting spatial and temporal dynamics of rabies virus evolution that, in Tunisia, are independent of the host species. Phylogenetic analysis of the nucleoprotein sequences
confirmed that all Tunisian isolates originated from wild canids belong to the Africa-1a clade of the Cosmopolitan lineage, commonly associated with dogs in several African localizations, especially in the north and the east [3,12,32]. These isolates do not form a distinct cluster within the clade but are mostly related to dog sequences. The inclusion of animal species as a discrete trait in the phylogenetic analysis suggests that all clusters defined in Tunisian canines and wild canids have a rabies virus from domestic dogs as an ancestor. These results differ from what is known for the clade associated with red foxes (*Vulpes vulpes*) in Europe that, despite belonging to the same cosmopolitan lineage, show unique characteristics likely resulting from the adaptation to this wild host [9] and clearly diverge from the root of all North African variants [33]. Each viral variant has its specificity for a host reservoir encrypted in its genetic settings [11,31]. Our results also suggest that there is not any independent spillover of rabies virus in wild canids in Tunisia. This is not the case for the neighboring Middle Eastern region (Iran, Irak, Syria, Jordan) where both urban and wild rabies cycles exist in the same country [16]. When aligning the Tunisian wild canids sequences with those of the Tunisian dog strains downloaded from GenBank, the ML phylogeny showed that all isolates that have similar sampling dates and locations strongly cluster together independently of the animal species. Specifically, the current genetic diversity of the Tunisian isolates is represented by two major and geographically distinct phylogenetic sub-clades, previously defined by Amouri et al. (1986) as north-east-center-south (NCS) and northwest (NW) [12]. Our analyses showed that very old Tunisian isolates (Tunisia 1986) and two Algerian sequences (ALG/08-153 and ALG/08-200) clustered together with the Tunisian isolates. In addition, the distribution of Tunisian isolates (Figure 6) showed two clusters within the NCS sub-clade that are spatially distinct. Variants from the newly identified cluster, NCS2, are exclusively located in the western parts of the country, precisely the northwest and the center-west. Interestingly, this cluster also included variants isolated from the northeast of Algeria (close to the northwest of Tunisia), suggesting a possible epidemiological link between these neighboring countries. Previous studies that analyzed rabies strains in Algeria, Morocco, and Tunisia suggested the grouping of viruses according to their country of origin [17]. On the other hand, we confirmed the genetic distinction of two sub-clades and further genetic clusters that strictly associated with different geographical areas and are evolving independently. Our data are consistent with results from a preliminary investigation undertaken on dog rabies cases in Tunisia (2011–2016) and support distinct epidemiological cycles of dog rabies in the northwest and northeast [34]. Viral variants associated with wild canids follow the same pattern seen for dogs, suggesting that the geography rather than the animal species is responsible for the phylogenetic structure. In addition, wild canids located in the governorate of Kasserine (NCS) can be infected with variants from both clusters NCS1 and NCS2, according to their city of origin. In particular, we detected variants from cluster NCS-1 in Sbeitla and El Hammar and from NCS-2 in Thala. Thala is geographically separated from Sbeitla and El Hammar by the presence of Chaambi and Semmama Mountains (Figure 6), supporting the fact that geographical barriers have a critical role in the evolution of rabies virus. Multiple studies have identified the role of natural landscape features such as mountain ranges and waterways as natural barriers to virus spread [35–37]. Specifically, elevation has been proven to act as a barrier to the dissemination of the rabies virus [35], given that human settlements, and therefore dog density, are less common at high elevations [38]. The Tunisian landscape is characterized by the presence of the dorsal ridge, the Tunisian part of the Atlas Mountains (including Jebal Chaambi and Jebal Semmama) that cross along the country from the center-west to the northeast. On both sides of the dorsal ridge, on the plains of the northwest and the center-east, scattered farms are widely spread. Domestic animal and dog presence is also influenced by such landscapes [39]. Thus, the dorsal ridge may act as a natural barrier, responsible for the geographical dispersion of both branches of the NCS variants. In addition to the latter hypothesis, if a spillover existed in wildlife, the viruses isolated from wild animals living on both sides of the same mountain would have had close genetic
characteristics. On the contrary, mountain chains seem to act as a natural barriers to
virus spread.

In the present study, we performed molecular clock analysis using sequences of
phosphoprotein but not nucleoprotein, because R squared analysis only showed positive
correlation on the P dataset. In addition, the higher genetic variability of the phospho-
protein gene makes it a better choice for phylodynamics reconstruction [3]. In fact, the
estimated substitution rate using Tempest and further confirmed by Tracer was 5 × 10⁻⁴
substitutions per site per year. This rate did not differ widely with previous studies of
lyssavirus evolution [3,40]. The relaxed molecular clock analysis suggested that the diver-
sification of the two sub-clades of Tunisian rabies variants (NCS/NW) occurred during the
second half of the twentieth century at the latest. Genetic clusters within the NCS sub-clade
diverged around 30 years ago; 95% HPD of the tMRCA estimates for the two nodes not
overlapping. Additionally, the NCS-2 cluster is genetically related to both Algerian rabies
isolates, probably introduced to Algeria before 1990. These age estimates are consistent
with previous analyses of RABV in North Africa, where the most recent common ancestor
for all North African RABV was estimated to have existed during the period 1878–1945 [17].
Thus, the diversification of the two Tunisian sub-clades (NCS/NW) is more recent than
that of all the North African clades. Our main limitation was the inability to access recent
North African sequences in GenBank. This might have affected our accuracy in tMRCA
estimation. In addition, because of the unavailability of recent Algerian sequences from
the northeast (near the Tunisian border), it was not possible to reliably distinguish whether
the virus isolates were the result of introduction of the Tunisian variants into Algeria.

**5. Conclusions**

Throughout this work, we demonstrated that rabies viruses isolated from dogs and
wild canids are not genetically distinct but rather evolve according to the geographical
area. Viruses from both sides of the dorsal ridge of the Atlas Mountains have evolved
independently, most likely through an urban cycle maintained by dogs. It can therefore be
recommended to the veterinary services that the fight against rabies in Tunisia remains
dependent on fighting dog-mediated rabies. Vaccination campaigns focusing on dogs
should be maintained and undertaken in coordination with laboratory-based surveillance,
dog population management, vaccination awareness and post-exposure prophylaxis.

**Author Contributions:** Conceptualization, Z.B. and H.K.; data curation, Z.B.; formal analysis, Z.B.; methodology, Z.B. and H.K.; resources (Biobank), N.B., J.B.S., S.B.M., M.H., and M.S.; software, Z.B.; supervision, H.K. and I.T.; validation, H.K. and I.T.; writing—original draft, Z.B.; writing—review and editing, H.K., A.G., and I.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Tunisian Ministry of higher education and scientific research under the grant number 'LR16-IPT03' and by Institute Pasteur of Tunis through the project "Phylogeography of rabies in Tunisia" under the grant number 'PCI-2017-06'.

**Institutional Review Board Statement:** The study protocol was approved by the Ethics Committee (Comité d'Ethique Bio-Médicale) of the Institute Pasteur de Tunis. File number 2018/13/lboRage approved on 8 June 2018.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** The authors would like to thank Jaber A. Belkhiria (One Health Institute, School of Veterinary Medicine, University of California Davis, Davis, CA, USA) for his time and efforts in editing the manuscript.

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

1. Bourhy, H.; Kissi, B.; Tordo, N. Molecular Diversity of the Lyssavirus Genus. *Virology* **1993**, *194*, 70–81. [CrossRef]

2. Tordo, N.; Poch, O.; Ermine, A.; Keith, G. Primary Structure of Leader RNA and Nucleoprotein Genes of the Rabies Genome: Segmented Homology with VSV. *Nucleic Acids Res.* **1986**, *14*, 2671–2683. [CrossRef]

3. Troupin, C.; Dacheux, L.; Tanguy, M.; Sabeta, C.; Blanc, H.; Bouchier, C.; Vignuzzi, M.; Duchene, S.; Holmes, E.C.; Bourhy, H. Large-Scale Phylogenomic Analysis Reveals the Complex Evolutionary History of Rabies Virus in Multiple Carnivore Hosts. *PLoS Pathog.* **2016**, *12*, e1006041. [CrossRef] [PubMed]

4. International Committee on Taxonomy of Viruses (ICTV). 2021. Available online: https://talk.ictvonline.org/taxonomy/ (accessed on 21 October 2021).

5. Mollentze, N.; Biek, R.; Streicker, D.G. The Role of Viral Evolution in Rabies Host Shifts and Emergence. *Curr. Opin. Virol.* **2014**, *8*, 68–72. [CrossRef]

6. Rupprecht, C.E.; Hanlon, C.A.; Hemachudha, T. Rabies Re-Examined. *Lancet Infect. Dis.* **2002**, *2*, 327–343. [CrossRef]

7. Bourhy, H.; Reynes, J.-M.; Dunham, E.J.; Dacheux, L.; Larrous, F.; Huong, V.T.Q.; Xu, G.; Yan, J.; Miranda, M.E.G.; Holmes, E.C. The Origin and Phylogeography of Dog Rabies Virus. *J. Gen. Virol.* **2008**, *89*, 2673–2681. [CrossRef] [PubMed]

8. Fooks, A.R.; Banyard, A.C.; Horton, D.L.; Johnson, N.; McElhinney, L.M.; Jackson, A.C. Current Status of Rabies and Prospects for Elimination. *Lancet* **2014**, *384*, 1389–1399. [CrossRef]

9. Blancou, J. Ecology and Epidemiology of Fox Rabies. *Rev. Infect. Dis.* **1988**, *10*, S606–S609. [CrossRef]

10. Oliveira, R.N.; Freire, C.C.; Iamarino, A.; Zanotto, P.M.; Pessoa, R.; Sanabani, S.S.; de Souza, S.P.; Castilho, J.G.; Batista, H.B.C.R.; Carnieli, P., Jr.; et al. Rabies Virus Diversification in Aerial and Terrestrial Mammals. *Genet. Mol. Biol.* **2020**, *43*, e20190370. [CrossRef] [PubMed]

11. Blancou, J.; Aubert, M.F. Transmission of rabies virus: Importance of the species barrier. *Bull. Acad. Natl. Med.* **1997**, *181*, 301–311; discussion 311–312. [PubMed]

12. Amouri, I.K.; Kharmachi, H.; Djebbi, A.; Saadi, M.; Hogga, N.; Zakour, L.B.; Ghram, A. Molecular Characterization of Rabies Virus Isolated from Dogs in Tunisia: Evidence of Two Phylogenetic Variants. *Virus Res.* **2011**, *158*, 246–250. [CrossRef] [PubMed]

13. Hogga, N. *Etude Moléculaire des Isolats de Virus Rabique Chez les Animaux Domestiques*; Mémoire de Master de Microbiologie; Faculté des Sciences de Tunis: Tunis, Tunisia, 2007.

14. Lamine, H. *Caractérisation Moléculaire des Isolats du Virus Rabique à Partir des Cas de Rage Humaine en Tunisie*; Mémoire de fin d'études; Institut de Biotechnologie de Monastir: Monastir, Tunisia, 2007.

15. Talbi, C.; Lemey, P.; Suchard, M.A.; Abdelatif, E.; Elharrak, M.; Jalal, N.; Faouzi, A.; Echevarría, J.E.; Morón, S.V.; Rambaut, A.; et al. Phylodynamics and Human-Mediated Dispersal of a Zoonotic Virus. *PLoS Pathog.* **2010**, *6*, e1001166. [CrossRef]

16. Available online: http://www.rage.tn/upload/1516305067.pdf (accessed on 30 September 2021).

17. Horton, D.L.; McElhinney, L.M.; Freuling, C.M.; Marston, D.A.; Banyard, A.C.; Goharrriz, H.; Wise, E.; Breed, A.C.; Saturday, G.; Kolodziejek, J.; et al. Complex Epidemiology of a Zoonotic Disease in a Culturally Diverse Region: Phylogeography of Rabies Virus in the Middle East. *PLoS Negl. Trop. Dis.* **2015**, *9*, e0003569. [CrossRef] [PubMed]

18. Available online: https://www.oie.int/fr/ce-que-nous-faisons/normes/codes-et-manuels/acces-en-ligne-au-manuel-terrestre/ (accessed on 26 October 2021).

19. World Health Organization. WHO Expert Consultation on Rabies: Third Report; WHO Technical Report Series; World Health Organization: Geneva, Switzerland, 2018; ISBN 978-92-4-121021-8.

20. Nadin-Davis, S.A. Polymerase Chain Reaction Protocols for Rabies Virus Discrimination. *J. Virol. Methods* **1998**, *75*, 1–8. [CrossRef]

21. Wang, L.; Wu, H.; Tao, X.; Li, H.; Rayner, S.; Liang, G.; Tang, Q. Genetic and Evolutionary Characterization of RABVs from China Using the Phosphoprotein Gene. *Virol. J.* **2013**, *10*, 14. [CrossRef] [PubMed]

22. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. *Mol. Biol. Evol.* **2013**, *30*, 772–780. [CrossRef]

23. Guindon, S.; Dufayard, J.-F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. *Syst. Biol.* **2010**, *59*, 307–321. [CrossRef] [PubMed]

24. Lefort, V.; Longueville, J.-E.; Gascuel, O. SMS: Smart Model Selection in PhyML. *Mol. Biol. Evol.* **2017**, *34*, 2422–2424. [CrossRef]

25. FigTree. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 12 January 2021).

26. Bouckaert, R.; Vaughan, T.G.; Barido-Sottani, J.; Duchêne, S.; Fourment, M.; Gavryushkina, A.; Heled, J.; Jones, G.; Kühnert, D.; Maio, N.D.; et al. BEAST 2.5: An Advanced Software Platform for Bayesian Evolutionary Analysis. *PLoS Comput. Biol.* **2019**, *15*, e1006650. [CrossRef]

27. Rambaut, A.; Lam, T.T.; Max Carvalho, L.; Pybus, O.G. Exploring the Temporal Structure of Heterochronous Sequences Using TempEst (Formerly Path-O-Gen). *Virus Evol.* **2016**, *2*, vew007. [CrossRef] [PubMed]

28. Bouckaert, R.R.; Drummond, A.J. BModelTest: Bayesian Phylogenetic Site Model Averaging and Model Comparison. *BMC Evol. Biol.* **2017**, *17*, 42. [CrossRef] [PubMed]

29. Drummond, A.J.; Ho, S.Y.W.; Phillips, M.J.; Rambaut, A. Relaxed Phylogenetics and Dating with Confidence. *PLoS Biol.* **2006**, *4*, e88. [CrossRef]

30. Kahle, D.; Wickham, H. Ggmap: Spatial Visualization with Ggplot2. *R J.* **2013**, *5*, 144–161. [CrossRef]
31. Bonnaud, E.M.; Troupin, C.; Dacheux, L.; Holmes, E.C.; Monchatre-Leroy, E.; Tanguy, M.; Bouchier, C.; Cliquet, F.; Barrat, J.; Bourhy, H. Comparison of Intra- and Inter-Host Genetic Diversity in Rabies Virus during Experimental Cross-Species Transmission. *PLoS Pathog.* **2019**, *15*, e1007799. [CrossRef]

32. Omodo, M.; Ar Gouilh, M.; Mwiine, F.N.; Okurut, A.R.A.; Nantima, N.; Namatovu, A.; Nakanjako, M.F.; Isingoma, E.; Arinaitwe, E.; Esau, M.; et al. Rabies in Uganda: Rabies Knowledge, Attitude and Practice and Molecular Characterization of Circulating Virus Strains. *BMC Infect. Dis.* **2020**, *20*, 200. [CrossRef]

33. Bourhy, H.; Kissi, B.; Audry, L.; Smreczak, M.; Sadkowska-Todys, M.; Kulonen, K.; Tordo, N.; Zmudzinski, J.F.; Holmes, E.C. Ecology and Evolution of Rabies Virus in Europe. *J. Gen. Virol.* **1999**, *80 Pt 10*, 2545–2557. [CrossRef] [PubMed]

34. Bouslama, Z.; Belkhiria, J.A.; Turki, I.; Kharmachi, H. Spatio-Temporal Evolution of Canine Rabies in Tunisia, 2011–2016. *Prev. Vet. Med.* **2020**, *185*, 105195. [CrossRef]

35. Brunker, K.; Lemey, P.; Marston, D.A.; Fooks, A.R.; Lugelo, A.; Ngeleja, C.; Hampson, K.; Biek, R. Landscape Attributes Governing Local Transmission of an Endemic Zoonosis: Rabies Virus in Domestic Dogs. *Mol. Ecol.* **2018**, *27*, 773–788. [CrossRef]

36. Saito, M.; Oshitani, H.; Orbina, J.R.C.; Tohma, K.; de Guzman, A.S.; Kamigaki, T.; Demetria, C.S.; Manalo, D.L.; Noguchi, A.; Inoue, S.; et al. Genetic Diversity and Geographic Distribution of Genetically Distinct Rabies Viruses in the Philippines. *PLoS Negl. Trop. Dis.* **2013**, *7*, e2144. [CrossRef]

37. Tohma, K.; Saito, M.; Kamigaki, T.; Tuason, L.T.; Demetria, C.S.; Orbina, J.R.C.; Manalo, D.L.; Miranda, M.E.; Noguchi, A.; Inoue, S.; et al. Phylogeographic Analysis of Rabies Viruses in the Philippines. *Infect. Genet. Evol.* **2014**, *23*, 86–94. [CrossRef]

38. Cohen, J.E.; Small, C. Hypsographic Demography: The Distribution of Human Population by Altitude. *PNAS* **1998**, *95*, 14009–14014. [CrossRef] [PubMed]

39. Maurer, G. Agriculture in the Rif and Tell Mountains of North Africa. *Mt. Res. Dev.* **1992**, *12*, 337–347. [CrossRef]

40. Holmes, E.C.; Woelk, C.H.; Kassis, R.; Bourhy, H. Genetic Constraints and the Adaptive Evolution of Rabies Virus in Nature. *Virology* **2002**, *292*, 247–257. [CrossRef] [PubMed]