# Serological Evidence of Lyssaviruses among Bats on Southwestern Indian Ocean Islands

Julien Mélade, Stewart Mcculloch, Beza Ramasindrazana, Erwan Lagadec, Magali Turpin, Hervé Pascalis, Steven M. Goodman, Wanda Markotter, Koussay Dellagi

► To cite this version:

Julien Mélade, Stewart Mcculloch, Beza Ramasindrazana, Erwan Lagadec, Magali Turpin, et al. Serological Evidence of Lyssaviruses among Bats on Southwestern Indian Ocean Islands. PLoS ONE, 2016, 11 (8), pp.e0160553. 10.1371/journal.pone.0160553. hal-01452894

HAL Id: hal-01452894

https://univ-reunion.hal.science/hal-01452894v1

Submitted on 28 Jun 2018

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 CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Serological Evidence of Lyssaviruses among Bats on Southwestern Indian Ocean Islands

Julien Mélade¹,²,³, Stewart McCulloch⁴, Beza Ramasindrazana¹,²,³,⁵,⁶, Erwan Lagadec¹,²,³, Magali Turpin¹,²,³, Hervé Pascalis¹,²,³, Steven M. Goodman⁵,⁷, Wanda Markotter⁴*, Koussay Dellagi¹,²,³*

¹ Centre de Recherche et de Veille sur les Maladies Emergentes dans l'Océan Indien (CRVOI), Plateforme de Recherche CYROI, Sainte Clotilde, La Réunion, France, ² Université de La Réunion, UMR PIMIT «Processus Infectieux en Milieu Insulaire Tropical», INSERM U1187, CNRS 9192, IRD 249. Plateforme de Recherche CYROI, Saint Denis, La Réunion, France, ³ Institut de Recherche pour le Développement (IRD), Sainte-Clotilde, La Réunion, France, ⁴ Center for Viral Zoonoses, Department of Medical Virology, School of Medicine, Faculty of Health Sciences, University of Pretoria, Pretoria, South Africa, ⁵ Association Vahatra, Antananarivo, Madagascar, ⁶ Institut Pasteur de Madagascar, Antananarivo, Madagascar, ⁷ Field Museum of Natural History, Chicago, Illinois, United States of America

* koussay.dellagi@gmail.com (KD); wanda.markotter@up.ac.za (WM)

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

OPEN ACCESS

**Citation:** Mélade J, McCulloch S, Ramasindrazana B, Lagadec E, Turpin M, Pascalis H, et al. (2016) Serological Evidence of Lyssaviruses among Bats on Southwestern Indian Ocean Islands. PLoS ONE 11(8): e0160553. doi:10.1371/journal.pone.0160553

**Editor:** Pierre Roques, CEA, FRANCE

**Received:** March 31, 2016

**Accepted:** July 21, 2016

**Published:** August 8, 2016

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

**Data Availability Statement:** Data are available via Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.gp4h2.

**Funding:** This work was supported by funding from CRVOI and from the European Regional Development Funds FEDER POCT "Pathogènes associés à la Faune Sauvage Océan Indien" project (N° 31189), and FEDER-POCT ParamyxOI project (N° 33857), as well as the National Research Foundation, South Africa (UID 91496, and 85786). Sampling on Mayotte was supported by ECOSAN/ CNRS-INEE "Batman". B. Ramasindrazana's postdoctoral fellowship at CRVOI was funded by "RUN-Emerge, European Commission FP7 Regpot

## Abstract

We provide serological evidence of lyssavirus circulation among bats on southwestern Indian Ocean (SWIO) islands. A total of 572 bats belonging to 22 species were collected on Anjouan, Mayotte, La Réunion, Mauritius, Mahé and Madagascar and screened by the Rapid Fluorescent Focus Inhibition Test for the presence of neutralising antibodies against the two main rabies related lyssaviruses circulating on the African continent: *Duvenhage lyssavirus* (DUVV) and *Lagos bat lyssavirus* (LBV), representing phylogroups I and II, respectively. A total of 97 and 42 sera were able to neutralise DUVV and LBV, respectively. No serum neutralised both DUVV and LBV but most DUVV-seropositive bats (n = 32/220) also neutralised *European bat lyssavirus 1* (EBLV-1) but not *Rabies lyssavirus* (RABV), the prototypic lyssavirus of phylogroup I. These results highlight that lyssaviruses belonging to phylogroups I and II circulate in regional bat populations and that the putative phylogroup I lyssavirus is antigenically closer to DUVV and EBLV-1 than to RABV. Variation between bat species, roost sites and bioclimatic regions were observed. All brain samples tested by RT-PCR specific for lyssavirus RNA were negative.

## Introduction

Lyssaviruses (order: *Mononegavirales*, family: *Rhabdoviridae*) are RNA viruses with single-stranded, negative-sense genomes approximately 12 kb in length [1]. Fifteen distinct viral species, classified into three phylogroups, have been identified to date within the *Lyssavirus* genus [2, 3]. These are the prototypic *Rabies lyssavirus* (RABV) and the 14 genetically related species referred to as rabies-related viruses that cluster into different phylogroups; *Duvenhage lyssavirus* (DUVV), *European bat lyssavirus* type 1 (EBLV-1) and type 2 (EBLV-2), *Australian bat*
Capacity program", the "Fonds de Coopération régionale, Prefecture de La Réunion", and The Field Museum of Natural History, Chicago, through the Dr. Ralph and Marian Falk Medical Research Trust.

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

*lyssavirus* (ABLV), *Irkut lyssavirus* (IRKV), *Aravan lyssavirus* (ARAV), *Khujand lyssavirus* (KHUV) and *Bokeloh bat lyssavirus* (BBLV) constitute phylogroup I; *Lagos bat lyssavirus* (LBV), *Mokola lyssavirus* (MOKV) and *Shimoni bat lyssavirus* (SHIV) form phylogroup II; *Ikoma lyssavirus* (IKOV) and *West Caucasian bat lyssavirus* (WCBV) make up phylogroup III, whilst the newly described *Lleida bat lyssavirus* (LLEBV) potentially constituting phylogroup IV [3, 4]. With the exceptions of MOKV and IKOV, all lyssaviruses have been isolated from bats [5, 6].

RABV, the type species of the genus, is the principal causative agent of human rabies, a fatal acute encephalitis claiming an estimated 59,000 lives annually, primarily affecting Africa and Asia [7]. The natural cyclic maintenance of RABV worldwide is principally via dogs and wild canids [8], while bats are seen as a reservoir only in the New World, where the virus is known to infect insectivorous and hematophagous species [9, 10]. Nowadays, the majority of cases of human exposure to rabies in the New World is related to contact with rabid bats [10, 11]. In contrast, among the 14 rabies-related viruses, only EBLV-1 and 2, ABLV, DUVV, IRKV and MOKV have caused sporadic human deaths [12].

Terrestrial rabies is endemic on Madagascar where no native canids occur and dogs are both maintenance reservoirs and vectors. The other islands of the SWIO are deemed terrestrial rabies free and no case of RABV infection being recorded in non-volant mammals [13]. With respect to rabies related viruses, only one study conducted between 2005 and 2010 on Madagascar, recorded antibodies to LBV and EBLV-1, being positive in frugivorous bats (*Eidolon dupreanum* and *Pteropus rufus*) and negative in the six tested insectivorous bat species (n = 92) [14]. No data has been reported on rabies related lyssavirus infections in bats from other SWIO islands [13].

We report the results of a serosurvey conducted on bats from SWIO islands based on the screening of neutralising antibodies against two rabies related viruses, namely LBV and DUVV, known to circulate among African bat populations [15, 16]. Our data show evidence that at least two lyssaviruses circulate among bats in the insular ecosystems of the SWIO region. This study is the first comprehensive analysis of bat exposure to lyssaviruses in this region and aims to strengthen available information.

# Materials and Methods

## Specimen collection

From March 2010 to March 2015, 653 bats representing 22 species from six families were collected from the islands of Anjouan and Mayotte (Comoros Archipelago), La Réunion and Mauritius (Mascarenes), Mahé (Seychelles Archipelago) and Madagascar (Fig 1). Both insectivorous and frugivorous bat species were captured using harp traps, hand nets and mist nets. Individual bats were identified by morphological characteristics [17] and in most cases by comparison to voucher specimens in The Field Museum of Natural History (FMNH) and The Université d'Antananarivo, Département de Biologie Animale (UADBA).

Animals were manipulated in accordance with the guidelines for the handling of wild mammals [18], and all field protocols were designed to strictly follow the terms of the research permits issued by national authorities in the different countries this study took place (S1 Text). Details on ethical clearance are provided (S2 Text). This study benefited from the sampling efforts conducted in the context of an ongoing long-term international project to catalogue the terrestrial vertebrate fauna of Madagascar [19]. The tissue samples collected during these inventories were optimized by different research programs and already reported on in different publications [20–25]. Bats trapped on Mayotte were released at the capture site after sampling. A total of 38 *Pteropus seychellensis* from Mahé and 12 *Pteropus rufus* from Madagascar, all
<figure><img src="image_4.png" /><figcaption>**Fig 1. Geographic distribution of bats collected on the southwestern Indian Ocean islands of Anjouan (Union of Comoros), Mayotte (France), Madagascar, La Réunion (France), Mahé (Seychelles Archipelago) and Mauritius.** n: number of bats sampled.</figcaption></figure>

doi:10.1371/journal.pone.0160553.g001

destined for the bush meat trade, were obtained alive from local hunters, individual bats were processed within five minutes after death and tissue samples collected as for other islands (see below). Bat specimens sampled on the other islands were preserved as voucher specimens and curated in the FMNH and UADBA. Details concerning the sex and age (based on wing bone ossification and dental eruption patterns) [26] for each captured bat were recorded and are presented in database deposited on Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.gp4h2. The names of capture sites were abbreviated (e.g. "Mangajou" as MGJ) and reported in S1 Table. An open-source Geographic Information System software, QGIS [27], was used to generate Fig 1.

Blood was collected either by veinipuncture of the brachial vein for bats that were released at the site of capture or by cardiac puncture for bats that were dispatched (in some cases, needles and syringes were rinsed with heparin or EDTA to prevent blood coagulation during the procedure). Sera were recovered by centrifugation at 3100 g for 5 min. Brain samples were collected from vouchered individuals. Tissue samples and sera were immediately stored in liquid nitrogen in the field and then transferred to -80°C freezers upon arrival at the laboratory.

## Detection of lyssavirus neutralising antibodies

Detection and titration of lyssavirus neutralising antibodies was performed using the miniaturised Rapid Fluorescent Focus Inhibition Test (RFFIT) [15]. The challenge viruses were the
two most represented African bat associated lyssaviruses: DUVV (isolate DUVVSA06) and LBV (isolate LBVAFR1999) belonging to the phylogroups I and II, respectively. In some experiments, in order to delineate the antigenic relatedness of the putative locally circulating lyssavirus with other members of the *Lyssavirus* genus, sera were tested against two other representatives of phylogroup I: RABV (CVS-11) (n = 528) and EBLV-1 (isolate RV20) (n = 220).

Sera were heated at 56°C for 30 min to inactivate complement and tested at 1:5, 1:25, 1:125 and 1:625 dilutions. Each neutralisation assay was performed in triplicate and was read and validated by two different investigators. Serum samples producing significant viral neutralisation (i.e. > 50% reduction of fluorescent foci) at a dilution of 1:25 or higher, were considered positive for all tested viruses. Control wells were prepared with three dilutions of the challenge virus at 0.5FFD<sub>50</sub>, 5FFD<sub>50</sub> and 50FFD<sub>50</sub>, and a well containing media and cell culture without virus was used as negative control. Mouse neuroblastoma cells at 2.10<sup>6</sup> cells/ml were added to each well. Cultures were incubated in a humidified incubator at 37°C in the presence of 5% CO<sub>2</sub> for 20 to 24 hours depending on the virus used. Slides were then acetone fixed and stained with 10 µl of polyclonal fluorescein isothiocyanate (FITC) conjugated immunoglobulin (Onderstepoort Veterinary Institute, South Africa) at a dilution of 1:300 with Evans Blue (0.5% in phosphate buffered saline (PBS)) for 30 min at 37°C. Slides were examined for fluorescent foci using an epifluorescence microscope. Results were validated if i) the negative control chamber (i.e. no virus) showed a confluent layer of cells and no fluorescent focus, ii) fluorescence foci were observed in the 0.5FFD<sub>50</sub>, 5FFD<sub>50</sub> and 50FFD<sub>50</sub> chambers (i.e. virus multiplication).

Of the 653 serum samples, 81 were excluded from analysis either because they displayed a cell-adhesion inhibitory-type cytotoxicity (n = 45) or because the volume of sera was insufficient to do the assay in triplicate (n = 36). *In finé*, the RFFIT could be performed for 572 sera. Results for each individual bat are presented in database deposited on Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.gp4h2.

## Lyssaviruses RNA Detection

Although the Direct Fluorescent Antibody test (DFA) is the gold standard method to detect rabies antigen in brain tissue of rabid animals, we preferred to use the ultra-sensitive real time reverse transcription PCR, which can detect as low as 10 RNA copies, since all animals collected in our study were apparently “healthy” and we anticipated very low viral antigen expression in brain tissue [28]. Lyssavirus RNA detection was performed using a real time reverse transcription PCR (QuantiTect Probe RT-PCR Kit, QIAGEN, Cat. number 204443) with specific primers designed to target the conserved region of the nucleoprotein gene of lyssaviruses, including RABV, DUVV, EBLV-1 and EBLV-2 from phylogroup I, LBV and MOKV from phylogroup II [28].

## Statistical analysis

Statistical analyses were performed using Pearson chi-square (χ<sup>2</sup>) or Fisher's exact tests in R software for calculation procedures (95% confidence intervals with a continuity correction). Probability values < 0.05 were considered statistically significant. All statistical computations were conducted using R version 3.0.0 [29].

## Results

A total of 572 sera from 421 insectivorous bats and 151 frugivorous bats belonging to six families and 22 species, were tested by RFFIT for neutralising antibodies against two African bat lyssaviruses namely DUVV and LBV representing phylogroups I and II, respectively. Sera from 42 individual bats out of 572 (7.3%) were capable of neutralising LBV and sera from 97 out of
540 (17.9%) were capable of neutralising DUVV. The cut off value for positivity was at 1:25 serum dilution but most bats seropositive for DUVV reacted at a dilution of 1:125 or higher (85.5% and 81.8% of sera from insectivorous and frugivorous bats, respectively (S1 Fig)). The same held true for bats seropositive for LBV, with 66.5% and 60.0% of sera from insectivorous and frugivorous bats, respectively, reacting at a dilution of 1:125 or higher (S1 Fig). Table 1 shows the levels of reactivity for the two challenge viruses segregated by island, bat family and species.

In order to delineate the antigenic relatedness of the putative lyssavirus that circulates in the SWIO region with other viruses of the genus, sera were tested against two other representatives of phylogroup I: RABV (CVS11) (n = 528) and EBLV-1 (isolate RV20) (n = 220). Of the 97 sera neutralising DUVV, 32 cross-neutralised EBLV-1 and not one cross-neutralised RABV. Of note, 13 sera out of 220 only neutralised EBLV-1 and not DUVV.

Bats seropositive to either of the two test viruses belonged to the six different bat families represented in our collections. Among the 22 bat species sampled from the different islands, four tested seropositive for DUVV only, three for LBV only and nine for DUVV and LBV (in different individuals). Six species were negative for all tested viruses. The seropositivity rates for DUVV and LBV did not differ significantly between islands, or between insectivorous and frugivorous species, or between males and females or between sub-adults and adults (p > 0.05). However, statistically significant differences exist between species and are described below.

Among the 14 insectivorous bat species studied on Madagascar, the highest seropositivity rate for DUVV was detected in *Miniopterus griveaudi* (34.5%; n = 10/29) (p < 0.001) and the highest seropositivity rate for LBV in *Myotis goudoti* (36.4%; n = 4/11) (p = 0.005). Of the three Malagasy frugivorous bat species, all endemic to the island, two (*Pteropus rufus* and *Rousettus madagascariensis*) tested seropositive for DUVV (45.5% n = 5/11; and 12.1%, n = 4/33; respectively), with the highest rates in *Pteropus rufus* (p = 0.01). Only *Rousettus madagascariensis* tested seropositive for LBV (20.0%; n = 7/35). The third frugivorous bat, *Eidolon dupreanum*, tested seronegative (n = 0/9) for LBV and DUVV.

Significant differences in seropositivity rates were recorded between the four bioclimatic zones and three altitudinal levels of Madagascar delineated herein, but only for DUVV. The highest seropositivity rates for DUVV were recorded in the sub-humid region (27.3%; n = 15/55) (p = 0.003) and in the highest elevation zone (29.5%; n = 13/44, p = 0.02). No significant difference was observed for seroreactivity for LBV among bioclimatic zones and altitudes. Significant differences were observed among localities (with rates above 50.0% of seropositivity for DUVV recorded in ANJOZ2 and BELO) (p < 0.001) but not between sites containing multiple species as compared to sites with single species nor between bat populations collected in caves, synanthropic roost sites or forests. Interestingly, bats belonging to the same species collected at different localities, displayed significantly different seropositivity rates. For example, *Miniopterus griveaudi* and *Rousettus madagascariensis* collected at multi-species sites (i.e. ANJHB) had significantly higher seropositivity rates for DUVV (75.0%) and LBV (40.0%), respectively, whereas significantly lower rates at monospecific sites (p < 0.05): in ANJHK1, 9.1% for DUVV in *Miniopterus griveaudi* (n = 1/11) and at ANJHB, 0.0% for LBV in *Rousettus madagascariensis* (n = 0/11) (p < 0.001). At the site level, specific patterns of lyssavirus neutralisation were recorded. At 21 positive sites on Madagascar, the neutralisation pattern appeared homogeneous among the seropositive bats at the specific site. In 11 localities bats seroneutralised DUVV only, at three localities bats cross-neutralised DUVV and EBLV-1 and at three localities bats seroneutralised LBV only. The four other positive localities combine animals, which expressed any of these three profiles of lyssavirus neutralisation (S2 Fig).
<table><caption>Table 1. Seroreactivity of bats from SWIO islands to lyssaviruses as assessed by RFFIT.</caption><thead><tr><th rowspan="2"></th><th colspan="4">Neutralisation of only one virus</th><th>Neutralisation of two viruses</th></tr><tr><th>LBV</th><th>DUVV</th><th>EBLV-1</th><th>CVS-11</th><th>DUVV+EBLV-1</th></tr></thead><tbody><tr><td><strong>Anjouan</strong></td><td><strong>2/21</strong></td><td><strong>2/8</strong></td><td><strong>0/0</strong></td><td><strong>0/8</strong></td><td><strong>0/0</strong></td></tr><tr><td>Miniopteridae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Miniopterus griveaudi</td><td>0/2</td><td>0/1</td><td>0/0</td><td>0/8</td><td>0/0</td></tr><tr><td>Molossidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Chaerephon pusillus</td><td>2/19</td><td>2/7</td><td>0/0</td><td>0/0</td><td>0/0</td></tr><tr><td><strong>Madagascar</strong></td><td><strong>23/301</strong></td><td><strong>54/286</strong></td><td><strong>1/170</strong></td><td><strong>0/277</strong></td><td><strong>5/70</strong></td></tr><tr><td>Hipposideridae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Hipposideros commersoni</td><td>0/6</td><td>1/6</td><td>1/4</td><td>0/6</td><td>1/4</td></tr><tr><td>Miniopteridae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Miniopterus cf. ambohitrensis</td><td>0/17</td><td>4/15</td><td>0/4</td><td>0/15</td><td>0/4</td></tr><tr><td>Miniopterus gleni</td><td>0/4</td><td>0/4</td><td>0/0</td><td>0/4</td><td>0/0</td></tr><tr><td>Miniopterus griveaudi</td><td>2/33</td><td>10/29</td><td>0/6</td><td>0/23</td><td>0/6</td></tr><tr><td>Miniopterus mahafaliensis</td><td>0/11</td><td>0/11</td><td>0/0</td><td>0/9</td><td>0/0</td></tr><tr><td>Miniopterus sororculus</td><td>0/3</td><td>0/2</td><td>0/1</td><td>0/23</td><td>0/1</td></tr><tr><td>Molossidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Chaerephon atsinanana</td><td>0/16</td><td>0/15</td><td>0/0</td><td>0/16</td><td>0/0</td></tr><tr><td>Chaerephon leucogaster</td><td>0/28</td><td>9/28</td><td>0/6</td><td>0/26</td><td>1/6</td></tr><tr><td>Mops leucostigma</td><td>2/20</td><td>0/19</td><td>0/4</td><td>0/20</td><td>0/4</td></tr><tr><td>Mops midas</td><td>0/7</td><td>0/7</td><td>0/0</td><td>0/7</td><td>0/0</td></tr><tr><td>Mormopterus jugularis</td><td>3/58</td><td>17/58</td><td>0/8</td><td>0/58</td><td>0/8</td></tr><tr><td>Otomops madagascariensis</td><td>3/20</td><td>4/18</td><td>0/11</td><td>0/20</td><td>0/11</td></tr><tr><td>Pteropodidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Eidolon dupreanum</td><td>0/9</td><td>0/9</td><td>0/1</td><td>0/9</td><td>0/1</td></tr><tr><td>Pteropus rufus</td><td>0/12</td><td>5/11</td><td>0/5</td><td>0/10</td><td>3/5</td></tr><tr><td>Rousettus madagascariensis</td><td>7/35</td><td>4/33</td><td>0/16</td><td>0/33</td><td>0/16</td></tr><tr><td>Rhinonycteridae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Triaenops menamena</td><td>2/11</td><td>0/11</td><td>0/2</td><td>0/7</td><td>0/2</td></tr><tr><td>Vespertilionidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Myotis goudoti</td><td>4/11</td><td>0/10</td><td>0/2</td><td>0/11</td><td>0/2</td></tr><tr><td><strong>Mahé (Seychelles)</strong></td><td><strong>4/40</strong></td><td><strong>6/40</strong></td><td><strong>1/17</strong></td><td><strong>0/40</strong></td><td><strong>1/17</strong></td></tr><tr><td>Pteropodidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Pteropus seychellensis</td><td>4/40</td><td>6/40</td><td>1/17</td><td>0/40</td><td>1/17</td></tr><tr><td><strong>Mauritius</strong></td><td><strong>7/67</strong></td><td><strong>19/67</strong></td><td><strong>2/29</strong></td><td><strong>0/62</strong></td><td><strong>0/29</strong></td></tr><tr><td>Molossidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mormopterus acetabulosus</td><td>5/31</td><td>6/31</td><td>1/8</td><td>0/29</td><td>0/8</td></tr><tr><td>Pteropodidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Pteropus niger</td><td>2/36</td><td>13/36</td><td>1/21</td><td>0/33</td><td>0/21</td></tr><tr><td><strong>Mayotte</strong></td><td><strong>2/22</strong></td><td><strong>2/22</strong></td><td><strong>0/22</strong></td><td><strong>0/22</strong></td><td><strong>0/22</strong></td></tr><tr><td>Molossidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Chaerephon pusillus</td><td>1/3</td><td>1/3</td><td>0/3</td><td>0/3</td><td>0/3</td></tr><tr><td>Pteropodidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Pteropus seychellensis</td><td>2/19</td><td>1/19</td><td>0/19</td><td>0/19</td><td>0/19</td></tr><tr><td><strong>La Réunion</strong></td><td><strong>3/121</strong></td><td><strong>14/117</strong></td><td><strong>9/82</strong></td><td><strong>0/119</strong></td><td><strong>26/82</strong></td></tr><tr><td>Molossidae</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mormopterus francoismoutoui</td><td>3/121</td><td>14/117</td><td>9/82</td><td>0/119</td><td>26/82</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1.** (Continued)</caption><thead><tr><th rowspan="2"></th><th colspan="4">Neutralisation of only one virus</th><th>Neutralisation of two viruses</th></tr><tr><th>LBV</th><th>DUVV</th><th>EBLV-1</th><th>CVS-11</th><th>DUVV+EBLV-1</th></tr></thead><tbody><tr><th>Grand total</th><td>42/572</td><td>97/540</td><td>13/220</td><td>0/528</td><td>32/220</td></tr></tbody></table>

Columns 2–5 show sera that neutralise LBV only, DUVV only, EBLV-1 only and CVS-11 only. Column six represents sera that cross neutralise both DUVV and EBLV-1 among the 220 sera tested against the two viruses. Fractions indicate at the numerator, the number of bats neutralising lyssaviruses and at the denominator the number of tested samples.

doi:10.1371/journal.pone.0160553.t001

The reactivity profile of lyssavirus antibodies differed on the various regional islands (S3–S7 Figs) and, in most cases, positive sites combined different animals that neutralised either DUVV, EBLV-1, or cross neutralised DUVV and EBLV-1 and animals that neutralised LBV only. With the exception of Mayotte (S4 Fig), most sites on La Réunion (S5 Fig), Anjouan (S3 Fig), Mahé (S7 Fig) and Mauritius (S6 Fig), combined different animals that neutralised either phylogroup I viruses or LBV only.

In addition to the serology tests described above, brain samples from 550 of these individuals were screened for the presence of lyssavirus nucleic acid, all of which yielded negative results.

## Discussion

Islands of the southwestern Indian Ocean (SWIO) region, comprise one of the world's most biodiverse hotspots [30], hosting a rich and unique chiropteran fauna. For example, on Madagascar, of the 45 recognized bat species, 36 are endemic to the island [31, 32]; the neighbouring SWIO islands have notably lower measures of bat species diversity [19]. Several studies have reported on seroprevalence to lyssaviruses in bats from Africa and Asia [16, 33–39], but until now limited information was available for SWIO islands, located between these two continents. As compared to continents, island ecosystems have specificities and peculiarities that may influence the epidemiology of infection in humans and animals [40]. Indeed, our study reveals serological evidence of wide spread lyssaviruses exposure among frugivorous and insectivorous bats on these islands (including those deemed rabies free). At a regional scale, our study considerably expands the geographic coverage and bat species investigated on six different islands, when compared to the single study previously published on two Malagasy fruit bats, *Eidolon dupreanum* and *Pteropus rufus* [14].

Our data reveal that on SWIO islands, antibodies to lyssaviruses are present at high rates, both in frugivorous and insectivorous bats as already reported from Asia [37, 38]. The rates of seropositivity reported in our study varied from 10.0% to 44.4% in bat species for which the sampling was sufficient (n > 20 individuals). The high positive rates reported herein have also been found in other studies outside the region of insectivorous bat species, such as 60.8% (n = 45/74) for *Taphozous* sp. and 50.0% (n = 63/126) for *Hipposideros larvatus* in northern Vietnam [37], 59.0% (n = 62/105) for *Scotophilus kuhlii* and 54.1% for *Philetor brachypterus* in the Philippines [39]. Similarly, high rates have been reported in frugivorous bat species, such as 51.9% (n = 14/27) for *Pteropus hypomelanus* in the Philippines [39].

Although the high detection rates of lyssavirus infection in some bat species may indicate that they are particularly susceptible to virus transmission, prevalence data should be interpreted with prudence; they represent only a momentary picture of the immunity to *Lyssavirus* in the study sample at the specific moment they were collected and in some cases based on a limited number of individuals. Seroprevalence levels may change over time on a given island and across bat species, which may depend on age and sex, antibody kinetics, individual B cell
memory responses, co-infections, seasonal dynamics of infection, etc. These aspects clearly indicate caution to not over interpret the quantitative data, specifically in the absence of comparative analysis of longitudinal studies that are complicated to conduct in wild animal populations. Our results revealed significant differences related to bat species, roost sites and bioclimatic regions on Madagascar, and further suggest the role of abiotic or biotic factors (e.g. climate, altitude, sex, etc.) in lyssavirus circulation and dynamics as previously reported for Malagasy bats and their hosted paramyxoviruses [25].

Our data highlight that 12 insectivorous bat species and three out of the four frugivorous bat species from the SWIO region were exposed to lyssaviruses related to or similar to DUVV or LBV. Moreover, we recorded antibodies reacting to DUVV or LBV in nine bat species indicating an unrestricted circulation of these two viral species. Based on the sero-neutralisation profiles, at least one circulating *Lyssavirus* appears referable to phylogroup I, since most of the positive sera were able to neutralise DUVV or EBLV-1. In fact, at the regional level, the highest rates of seroreactivity were recorded against DUVV (17.9%), followed by a dual reactivity with DUVV and EBLV (14.5%), and a lower rate with LBV (7.3%). This reactivity profile is conserved across the different island populations tested in this study.

These results also support the view that a phylogroup II lyssavirus, closely related to, if not LBV itself, is circulating in the region. This is in agreement with the study of Reynes and colleagues [14], who reported that 24% of sera from Malagasy *Eidolon dupreanum* neutralised only LBV. In our series, the highest seropositivity rates against LBV across the different island and bat species were recorded on Mayotte and Mauritius (13.6% and 10.4%, respectively) and in populations of *Myotis goudoti* and *Rousettus madagascariensis* on Madagascar (36.4% and 20.0%, respectively).

The use of EBLV-1 and Rabies lyssavirus (CVS-11) as additional challenge viruses in the RFFIT, allowed delineation of cross reactivity profiles (Table 1). The observation that bat sera neutralizing DUVV occasionally cross neutralise EBLV-1 was expected, as the two viruses belong to the same phylogroup I. However, it strikingly contrasts with the absence of any cross neutralisation with CVS-11 (RABV), which also belongs to the same phylogroup [4, 33]. Cross neutralisation between lyssaviruses, within either phylogroup I or II has been previously reported [4, 41]. Experimental studies in animals also supported cross protection induced by RABV vaccines to other lyssaviruses of phylogroup I [42–45]. The lyssavirus glycoprotein carries most of the antigenic epitopes and regions targeted by neutralising antibodies [46]. Antibodies against conserved G protein epitopes shared by two lyssaviruses may allow cross neutralisation of these two lyssaviruses [4].

A brief comment is needed to explain why bat sera neutralizing DUVV/EBLV1 do not cross neutralize CVS11, even though they all belong to the same phylogroup. Cross-antigenicity is best detected by hyperimmune antisera generated after separate active immunization of mice or rabbits for each virus [4]. It is not required that a natural infection (such as in bats naturally infected by *Lyssavirus*) should constantly elicit antibodies specific for cross-antigenic epitopes, even if the latter do exist. This will depend on their intrinsic antigenicity and the strength of the immune antibody response after natural infection. Moreover, even if such antibodies were generated at some point after natural infection, they may have subsequently disappeared, due to their own kinetics compared to that of non-cross neutralizing antibodies. Antibodies to *Lyssavirus* detected in bat sera, persisting long after remote infection, most likely represent only a fraction of the repertoire of specific memory B cells.

Thus, the profile of cross neutralisation recorded in our study may indicate closer antigenic relativity of the putative lyssavirus(es) circulating in the SWIO with DUVV and EBLV-1 than with RABV. This might be particularly true in an insectivorous bat species from La Réunion, *Mormopterus francoismoutoui*, for which high rates of reactivity for EBLV-1 only and cross reactivity with DUVV and EBLV-1 were detected.
Madagascar and islands making up the Comoros archipelago are closer to the African continent than the more outlying islands of the Mascarenes and those sampled in the Seychelles archipelago. Among the lyssaviruses belonging to phylogroups I and II, only DUVV, SHIV and LBV, have been isolated to date from African bats. Hence, DUVV and LBV might be the putative lyssaviruses actually circulating on SWIO islands. Alternatively, the latter might correspond to closely related and unidentified new lyssavirus species. However, this will require successful isolation of these viruses and further molecular and genetic characterisation for their detailed diagnosis.

Although the main aim of this study was to detect serologic traces of exposure to *Lyssavirus* amongst bats of the SWIO islands [16], we also tested available brain samples for the presence of lyssavirus nucleic acids. We did not detect lyssavirus RNA in the brain of any insectivorous or frugivorous bat. This is not an unexpected result as no abnormal morbidity or mortality among sampled bats was observed during our study. This result indicates that the potential risk for humans or other animals to be exposed to an actively infected bat is low; however, the serological evidence indicates that rabies related lyssaviruses do circulate in bat populations in the SWIO. In the Old World, very few cases of human rabies due to rabies-related bat lyssaviruses are known: to date three fatal cases have been reported for DUVV [47], whilst no human case due to LBV has been confirmed. However, cases of rabies in Africa and Asia are known to be under estimated and when reported, the causative lyssavirus is seldom characterized. In the SWIO region, an increase awareness of the risk of rabies transmission is recommended after contact with bats, specifically via biting or scratching or saliva contact on open wounds or mucosal membranes. Any exposure should be immediately followed by appropriate post exposure prophylactic actions as recommended by the WHO guidelines [48].

**Supporting Information**

**S1 Fig. Distribution of neutralising antibody titers to DUVV and LBV in bats from southwestern Indian Ocean islands as assessed by the Rapid Fluorescent Focus Inhibition Test.**

Results represent the proportion of sera from insectivorous and frugivorous bat species titering at 1/25 (the cut off value of positivity), 1/125, and 1/625 among the sera tested positive for each challenge virus.

(TIF)

**S2 Fig. Sites recorded with samples seropositive to lyssaviruses and lyssaviruses neutralisation profiles on Madagascar.** Abbreviations next to illustrated islands indicate names of capture sites (e.g. ANKPK for "Ankapoka") reported in S1 Table. Coloured site names and in the chart correspond to the lyssaviruses antibodies detected in bats: sera neutralising DUVV only are in blue; those neutralising LBV only are in red, those cross-neutralising DUVV + EBLV-1 are in yellow. Sites which combine animals which sera neutralised DUVV only, EBLV-1 only, cross-neutralised DUVV and EBLV-1 and animals neutralised LBV only are in brown. The six red squares indicate the provincial capitals.

(TIF)

**S3 Fig. Sites recorded with samples seropositive to lyssaviruses and lyssaviruses neutralisation profiles on Anjouan.** Abbreviations next to illustrated islands indicate names of capture sites (e.g. COA for "College d'Ouani") reported in S1 Table. Coloured site names and in the chart correspond to the lyssaviruses antibodies detected in bats: sera neutralising DUVV only are in blue and sites which combine animals which sera neutralised DUVV only and LBV only are in brown.

(TIF)
**S4 Fig.** Sites recorded with samples seropositive to lyssaviruses and lyssaviruses neutralisation profiles on Mayotte. Abbreviations next to illustrated islands indicate names of capture sites (e.g. MGJ for "Mangajou") reported in S1 Table. Coloured site names and in the chart correspond to the lyssaviruses antibodies detected in bats: sera neutralising DUVV only are in blue and those neutralising LBV only are in red.
(TIF)

**S5 Fig.** Sites recorded with samples seropositive to lyssaviruses and lyssaviruses neutralisation profiles on La Réunion. Abbreviations next to illustrated islands indicate names of capture sites (e.g. GTB for "Grotte de Trois Bassin") reported in S1 Table. Coloured site names and in the chart correspond to the lyssaviruses antibodies detected in bats: sera neutralising DUVV only are in blue and sites which combine animals which sera cross-neutralised DUVV + EBLV-1 or LBV only are in brown.
(TIF)

**S6 Fig.** Sites recorded with samples seropositive to lyssaviruses and lyssaviruses neutralisation profiles on Mauritius. Abbreviations next to illustrated islands indicate names of capture sites (e.g. CSC for "Cascavelle") reported in S1 Table. Coloured site names and in the chart correspond to the lyssaviruses antibodies detected in bats: sera neutralising DUVV only are in blue; those cross-neutralising DUVV+EBLV-1 are in yellow. Sites which combine animals which sera neutralised DUVV only or EBLV-1 only and animals neutralised LBV only are in brown.
(TIF)

**S7 Fig.** Sites recorded with samples seropositive to lyssaviruses and lyssaviruses neutralisation profiles on Mahé. Abbreviations next to illustrated islands indicate names of capture sites (e.g. FRL for "FairyLand") reported in S1 Table. Coloured site names and in the chart correspond to the lyssaviruses antibodies detected in bats: sera neutralising. Sites which combine animals which sera neutralised DUVV only, EBLV-1 only, cross-neutralised DUVV+EBLV-1 and animals neutralised LBV only are in brown.
(TIF)

**S1 Table.** Information from sites on Anjouan, Madagascar, Mahé, Mauritius, Mayotte and La Réunion, indicated as: the island name, the site name (and abbreviation), GPS coordinates, type of habitat and number of species recorded at each site (one or more than one).
(DOC)

**S1 Text.** Authorisations for sampling of bats on different southwestern Indian Ocean islands. List of authorisations required for bat sampling from Madagascar, La Réunion, Mayotte, Anjouan, Mahé and Mauritius authorities.
(DOC)

**S2 Text.** Ethical clearance information. Ethical procedure applied for bat sampling on Madagascar, La Réunion, Mayotte, Anjouan, Mahé and Mauritius.
(DOCX)

## Acknowledgments

For kindly providing research permits, we are grateful on Madagascar to Département de Biologie Animale, Université d'Antananarivo; Direction du Système des Aires Protégées et Direction Générale de l'Environnement et des Forêts and Madagascar National Parks; on La Réunion to Direction de l'Environnement, de l'Aménagement et du Logement de La Réunion;
on Mauritius to National Park and Conservation; on Mayotte to Direction de l'Environnement,
de l'Aménagement et du Logement de Mayotte; on Anjouan to Centre National de Documen-
tation et de Recherche Scientifique (CNDRS) and Université des Comores; and on Mahé to
Direction of Wildlife, Trade and Conservation Section and the Ministry of Environment and
Energy of Republic of Seychelles. We thank the Institut Pasteur de Madagascar for their assis-
tance in storing samples and exportation. On Mayotte, we are grateful to Association Natura-
listes Environnement et Patrimoine de Mayotte for their assistance in locating bat roost sites
and to COOPADEM for providing laboratory space. We are grateful to Dr. Rajendraprasad
Sookhareea for his technical assistance in sampling bats on Mauritius. We thank the Secretary
of State for Environment, Food and Rural Affairs, London through the Animal and Plant
Health Agency of New Haw, Addlestone, Surrey, representing the European Virus Archive for
the supply of EBLV-1.

Author Contributions

**Conceptualization:** SMG WM KD.

**Data curation:** JM.

**Formal analysis:** JM WM KD.

**Funding acquisition:** WM KD.

**Investigation:** JM.

**Methodology:** SMG WM KD.

**Project administration:** WM KD.

**Resources:** SM BR EL MT SMG WM KD.

**Software:** JM.

**Supervision:** WM KD.

**Validation:** JM SM.

**Visualization:** JM.

**Writing - original draft:** JM SM HP SMG WM KD.

**Writing - review & editing:** JM SM BR EL HP SMG WM KD.

References

1. Bourhy H, Kissi B, Tordo N. Molecular diversity of the *Lyssavirus* genus. Virology. 1993 May; 194 (1):70-81. PMID: 8386891

2. International Committee on Taxonomy of Viruses (ICTV). Virus Taxonomy: 2014 Release.

3. Banyard AC, Evans JS, Luo TR, Fooks AR. Lyssaviruses and bats: Emergence and zoonotic threat. Viruses. 2014 Aug; 6(8):2974-90. doi: 10.3390/v6082974 PMID: 25093425

4. Badrane H, Bahloul C, Perrin P, Tordo N. Evidence of two Lyssavirus phylogroups with distinct patho- genicity and immunogenicity. J Virol. 2001 Apr; 75(7):3268-76. PMID: 11238853

5. Kgaladi J, Wright N, Coertse J, Markotter W, Marston D, Fooks AR, et al. Diversity and epidemiology of *Mokola* virus. PLoS Negl Trop Dis. 2013 Oct; 7(10):e2511. doi: 10.1371/journal.pntd.0002511 PMID: 24205423

6. Marston DA, Horton DL, Ngeleja C, Hampson K, McElhinney LM, Banyard AC, et al. *Ikoma* *Lyssavirus*, highly divergent novel *Lyssavirus* in African civet. Emerg Infect Dis. 2012 Apr; 18(4):664-7. doi: 10.3201/eid1804.111553 PMID: 22469151
7. Hampson K, Coudeville L, Lembo T, Sambo M, Kieffer A, Attlan M, et al. Estimating the global burden of endemic canine rabies. PLoS Negl Trop Dis. 2015 May; 9(5):e0003786. doi: 10.1371/journal.pntd.0003786 PMID: 25961848

8. Knobel DL, Cleaveland S, Coleman PG, Fèvre EM, Meltzer MI, Miranda ME, et al. Re-evaluating the burden of rabies in Africa and Asia. Bull WHO 2005 May; 83(5):360–8. PMID: 15976877

9. Escobar LE, Peterson AT, Favi M, Yung V, Medina-Vogel G. Bat-borne rabies in Latin America. Rev Inst Med Trop Sao Paulo. 2015 Jan-Feb; 57(1):63–72. doi: 10.1590/S0036-46652015000100009 PMID: 25651328

10. Webber BJ, Ayers KJ, Winterton BS, Yun HC, Cropper TL, Foster Jr, et al. Assessment of rabies exposure risk in a group of U.S. Air Force basic trainees—Texas, January 2014. MMWR Morb Mortal Wkly Rep. 2014 Aug; 63(34):749–52. PMID: 25166925

11. Johnson N, Aréchiga-Ceballos N, Aguilar-Setien A. Vampire bat rabies: ecology, epidemiology and control. Viruses. 2014 Apr; 6(5):1911–28. doi: 10.3390/v6051911 PMID: 24784570

12. WHO. Consultation on Rabies. Second report. Tech. Rep. Ser. 2013:1–139.

13. Andriamandimby SF, Héraud JM, Ramiandrasoa R, Ratsitorahina M, Rasambainarivo JH, Dacheux L, et al. Surveillance and control of rabies in La Réunion, Mayotte, and Madagascar. Vet Res. 2013 Sep; 44:77. doi: 10.1186/1297-9716-44-77 PMID: 24016204

14. Reynes JM, Andriamandimby SF, Razafitrimo GM, Razainirina J, Jeanmaire EM, Bourhy H, et al. Laboratory surveillance of rabies in humans, domestic animals, and bats in Madagascar from 2005 to 2010. Adv Prev Med. 2011; 2011:727821. doi: 10.4061/2011/727821 PMID: 21991442

15. Kuzmin IV, Niezgoda M, Franka R, Agwanda B, Markotter W, Beagley JC, et al. *Lagos Bat Virus* in Kenya. J Clin Microbiol. 2008 Apr; 46(4):1451–61. Epub 2008 Feb 27. doi: 10.1128/JCM.00016-08 PMID: 18305130

16. Markotter W, Monadjem A, Nel LH. Antibodies against *Duvenhage Virus* in insectivorous bats in Swaziland. J Wildl Dis. 2013 Oct; 49(4):1000–3. doi: 10.7589/2012-10-257 PMID: 24502729

17. Racey PA. Reproductive assessment in bats. In: Kunz TH, editor. Ecological and behavioral methods for the study of bats. Washington, DC: Smithsonian Institution Press. 1988; pp. 31–45.

18. Gannon WL, Sikes RS, American Society of Mammalogists. Guidelines of the American Society of Mammalogists for the use of wild mammals in research. J. Mammal. 2007 June; 88(3):809–23.

19. Goodman S.M. and Raherilalao MJ. Atlas of selected land vertebrates of Madagascar [Atlas d'une sélection de vertébrés terrestres de Madagascar] Association Vahatra, Antananarivo. 2013; 290 pp.

20. Wilkinson DA, Temmam S, Lebarbenchon C, Lagadec E, Chotte J, Guillebaud J, et al. Identification of novel paramyxoviruses in insectivorous bats of the Southwest Indian Ocean. Virus Res. 2012 Dec; 170(1–2):159–63. doi: 10.1016/j.virusres.2012.08.022 PMID: 22982204

21. Wilkinson DA, Mélade J, Dietrich M, Ramasindrazana B, Soarimalala V, Lagadec E, et al. Highly diverse morbillivirus-related paramyxoviruses in wild fauna of the southwestern Indian Ocean Islands: evidence of exchange between introduced and endemic small mammals. J Virol. 2014 Aug; 88 (15):8268–77. doi: 10.1128/JVI.01211-14 PMID: 24829336

22. Dietrich M, Wilkinson DA, Soarimalala V, Goodman SM, Dellagi K, Tortosa P. Diversification of an emerging pathogen in a biodiversity hotspot: *Leptospira* in endemic small mammals of Madagascar. Mol Ecol. 2014 Jun; 23(11):2783–96. doi: 10.1111/mec.12777 PMID: 24784171

23. Ramasindrazana B, Dellagi K, Lagadec E, Randrianarivelojosia M, Goodman SM, Tortosa P. Diversity, host specialization, and geographic structure of filarial nematodes infecting Malagasy bats. PLoS One. 2016 Jan; 11(1):e0145709. doi: 10.1371/journal.pone.0145709 PMID: 26751792

24. Gomard Y, Dietrich M, Wieseke N, Ramasindrazana B, Lagadec E, Goodman SM, et al. Malagasy bats shelter a considerable genetic diversity of pathogenic *Leptospira* suggesting notable host-specificity patterns. FEMS Microbiol Ecol. 2016 Apr; 92(4): pii: fiw037. 10.1093/femsec/fiw037.

25. Mélade J, Wieseke N, Ramasindrazana B, Flores O, Lagadec E, Gomard Y, et al. An eco-epidemiological study of Morbilli-related paramyxovirus infection in Madagascar bats reveals host-switching as the dominant macro-evolutionary mechanism. Sci. Rep. 2016 Apr 12; 6:23752. doi: 10.1038/srep23752 PMID: 27068130

26. Anthony ELP. Age determination in bats. In: Kunz TH, editor. Ecological and behavioral methods for the study of bats. Washington, DC: Smithsonian Institution Press. 1988; pp. 47–58.

27. Quantum GIS Development Team. Quantum GIS Geographic Information System. Open Source Geospatial Foundation Project. 2015. http://qgis.osgeo.org.

28. Coertse J, Weyer J, Nel L.H, Markotter W. Improved PCR methods for the detection of African rabies and rabies-related *Lyssaviruses*. J. Clin. Microbiol. 2010 Nov; 48(11):3949–55. doi: 10.1128/JCM.01256-10 PMID: 20810772
29. R. Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2013. ISBN 3-900051-07-0. Available: http://www.R-project.org/.

30. Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GA, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000 Feb; 403(6772):853-8. PMID: 10706275

31. Goodman S.M. The bats of Madagascar. [Les chauves-souris de Madagascar] Association Vahatra, Antananarivo. 2011; p. 129.

32. Goodman SM, Ramasindrazana B, Naughton K, Appleton B. Description of a new species of the *Mini-opterus aelleni* group (Chiroptera: Miniopteridae) from upland areas of central and northern Madagascar. Zootaxa. 2015 Mar; 3936(4):538-58. doi: 10.11646/zootaxa.3936.4.4 PMID: 25947452

33. Kuzmin IV, Turmelle AS, Agwanda B, Markotter W, Niezgoda M, Breiman RF, et al. Commerson's leaf-nosed bat (*Hipposideros commersoni*) is the likely reservoir of Shimoni bat virus. Vector Borne Zoonotic Dis. 2011 Nov; 11(11):1465-70. doi: 10.1089/vbz.2011.0663 PMID: 21867415

34. Kuzmin IV, Niezgoda M, Franka R, Agwanda B, Markotter W, Beagley JC, et al. Possible emergence of *West Caucasian bat virus* in Africa. Emerg Infect Dis. 2008 Dec; 14(12):1887-9. doi: 10.3201/eid1412.080750 PMID: 19046512

35. Hayman DT, Fooks AR, Horton D, Suu-Ire R, Breed AC, Cunningham AA, et al. Antibodies against *Lagos Bat Virus* in Megachiroptera from West Africa. Emerg Infect Dis. 2008 Jun; 14(6):926-8. doi: 10.3201/eid1406.071421 PMID: 18507903

36. Kemp GE, Causey OR, Moore DL, Odelola A, Fabiyi A. Mokola virus. Further studies on IbAn 27377, a new rabies-related etiologic agent of zoonosis in Nigeria. Am J Trop Med Hyg. 1972 May; 21(3):356-9. PMID: 5025622

37. Nguyen AT, Nguyen TT, Noguchi A, Nguyen DV, Ngo GC, Thong VD, et al. Bat lyssaviruses, Northern Vietnam. Emerg Infect Dis. 2014 Jan; 20(1):161-3. doi: 10.3201/eid2001.130813 PMID: 24377728

38. Reynes JM, Molia S, Audry L, Hout S, Ngin S, Walston J, et al. Serologic evidence of *Lyssavirus* infection in bats, Cambodia. Emerg Infect Dis. 2004 Dec; 10(12):2231-4. PMID: 15663870

39. Arguin PM, Murray-Lillibridge K, Miranda ME, Smith JS, Calaor AB, Rupprecht CE. Serologic evidence of *Lyssavirus* infections among bats, the Philippines. Emerg Infect Dis. 2002 Mar; 8(3):258-62. PMID: 11927022

40. Tortosa P, Pascalis H, Guernier V, Cardinale E, Le Corre M, Goodman SM, et al. Deciphering arboviral emergence within insular ecosystems. Infect Genet Evol. 2012 Aug; 12(6):1333-9. Epub 2012 Apr 4. doi: 10.1016/j.meegid.2012.03.024 PMID: 22504353

41. Gould AR, Hyatt AD, Lunt R, Kattenbelt JA, Hengstberger S, Blacksell SD. Characterisation of a novel lyssavirus isolated from Pteropid bats in Australia. Virus Res. 1998 Apr; 54(2):165-87. PMID: 9696125

42. Hanlon CA, Kuzmin IV, Blanton JD, Weldon WC, Manangan JS, Rupprecht CE. Efficacy of rabies bio-logics against new lyssaviruses from Eurasia. Virus Res. 2005 Jul; 111(1):44-54. PMID: 15896401

43. Brookes SM, Parsons G, Johnson N, McElhinney LM, Fooks AR. Rabies human diploid cell vaccine elicits cross-neutralising and cross-protecting immune responses against European and Australian bat lyssaviruses. Vaccine. 2005 Jul; 23(32):4101-9. PMID: 15964478

44. Brookes SM, Healy DM, Fooks AR. Ability of rabies vaccine strains to elicit cross-neutralising antibodies. Dev Biol (Basel). 2006 Jan; 125:185-93.

45. Malerczyk C, Freuling C, Gniel D, Giesen A, Selhorst T, Müller T. Cross-neutralization of antibodies induced by vaccination with purified Chick Embryo Cell Vaccine (PCECV) against different Lyssavirus species. Hum. Vaccin. Immunother. 2014 Oct; 10(10):2799-804. doi: 10.4161/21645515.2014.972741 PMID: 25483634

46. Benmansour A, Leblois H, Coulon P, Tuffereau C, Gaudin Y, Flamand A, et al. Antigenicity of rabies virus glycoprotein. J Virol. 1991 Aug; 65(8):4198-203. PMID: 1712859

47. van Thiel PP, de Bie RM, Eftimov F, Tepaske R, Zaaijer HL, van Doornum GJ, et al. Fatal human rabies due to *Duvenhage virus* from a bat in Kenya: Failure of treatment with coma-induction, ketamine, and antiviral drugs. PLoS Negl Trop Dis. 2009 Jul; 3(7):e428. doi: 10.1371/journal.pntd.0000428 PMID: 19636367

48. Guidelines for the Public Health Management of Rabies and Other Lyssavirus (including Australian Bat Lyssavirus) Exposures and Infections, 2014.