# Investigating the Circulation of Ebola Viruses in Bats during the Ebola Virus Disease Outbreaks in the Equateur and North Kivu Provinces of the Democratic Republic of Congo from 2018

Audrey Lacroix, Placide Mbala Kingebeni, Simon Pierre Ndimbo Kumugo, Guy Lempu, Christelle Butel, Laetitia Serrano, Nicole Vidal, Guillaume Thaurignac, Amandine Esteban, Daniel Mukadi Bamuleka, et al.

► To cite this version:

Audrey Lacroix, Placide Mbala Kingebeni, Simon Pierre Ndimbo Kumugo, Guy Lempu, Christelle Butel, et al. Investigating the Circulation of Ebola Viruses in Bats during the Ebola Virus Disease Outbreaks in the Equateur and North Kivu Provinces of the Democratic Republic of Congo from 2018. Pathogens, 2021, 10 (5), pp.557. 10.3390/pathogens10050557. [https://doi.org/10.36525/p pathogens.10050557](https://doi.org/10.36525/p pathogens.10050557)

HAL Id: hal-03652563

https://hal.umontpellier.fr/hal-03652563

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

Article

# Investigating the Circulation of Ebola Viruses in Bats during the Ebola Virus Disease Outbreaks in the Equateur and North Kivu Provinces of the Democratic Republic of Congo from 2018

Audrey Lacroix <fnref n="1" />, Placide Mbala Kingebeni <fnref n="2" />,<fnref n="3" />, Simon Pierre Ndimbo Kumugo <fnref n="2" />, Guy Lempu <fnref n="2" />, Christelle Butel <fnref n="1" />, Laetitia Serrano <fnref n="1" />, Nicole Vidal <fnref n="1" />, Guillaume Thaurignac <fnref n="1" />, Amandine Esteban <fnref n="1" />, Daniel Mukadi Bamuleka <fnref n="3" />,<fnref n="4" />, Jacques Likofata <fnref n="5" />, Eric Delaporte <fnref n="1" />, Jean-Jacques Muyembe Tamfum <fnref n="2" />,<fnref n="3" />, Ahidjo Ayouba <fnref n="1" />, Martine Peeters <fnref n="1" />,* and Steve Ahuka Mundeke <fnref n="2" />,<fnref n="3" />,*

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

check for
updates

**Citation:** Lacroix, A.; Mbala Kingebeni, P.; Ndimbo Kumugo, S.P.; Lempu, G.; Butel, C.; Serrano, L.; Vidal, N.; Thaurignac, G.; Esteban, A.; Mukadi Bamuleka, D.; et al. Investigating the Circulation of Ebola Viruses in Bats during the Ebola Virus Disease Outbreaks in the Equateur and North Kivu Provinces of the Democratic Republic of Congo from 2018. *Pathogens* **2021**, *10*, 557.
https://doi.org/10.3390/pathogens10050557

Academic Editor:
Tatjana Vilibic-Cavlek

Received: 31 March 2021

Accepted: 30 April 2021

Published: 4 May 2021

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

<figure><img src="image_4.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:** With 12 of the 31 outbreaks, the Democratic Republic of Congo (DRC) is highly affected by Ebolavirus disease (EVD). To better understand the role of bats in the ecology of Ebola viruses, we conducted surveys in bats during two recent EVD outbreaks and in two areas with previous outbreaks. Dried blood spots were tested for antibodies to ebolaviruses and oral and rectal swabs were screened for the presence of filovirus using a broadly reactive semi-nested RT-PCR. Between 2018 and 2020, 892 (88.6%) frugivorous and 115 (11.4%) insectivorous bats were collected. Overall, 11/925 (1.2%) to 100/925 (10.8%) bats showed antibodies to at least one Ebolavirus antigen depending on the positivity criteria. Antibodies were detected in fruit bats from the four sites and from species previously documented to harbor Ebola antibodies or RNA. We tested for the first time a large number of bats during ongoing EVD outbreaks in DRC, but no viral RNA was detected in the 676 sampled bats. Our study illustrates the difficulty to document the role of bats as a source of Ebolaviruses as they might clear quickly the virus. Given the increasing frequency of EVD outbreaks, more studies on the animal reservoir are urgently needed.

**Keywords:** Ebola virus; bats; the Democratic Republic of the Congo

## 1. Introduction

Since the first recognized Ebolavirus disease (EVD) outbreak in humans in 1976 in the Democratic Republic of Congo (DRC), a total of 31 outbreaks have now been recognized in Africa [1–3]. The majority remained limited in the number of infected individuals and in the geographic spread, but the outbreaks in West Africa (December 2013 to March 2016) and in Eastern DRC (August 2018 to June 2020) clearly showed that they can also become epidemic and infect thousands of individuals over large geographic areas [4]. In addition, the frequency of EVD outbreaks also increased over the last decades. With twelve outbreaks, DRC is one of the most affected countries. Since May 2018, the country has been continuously facing the circulation of the virus during four consecutive outbreaks [3,5–7].

<footnote n="1">¹ TransVIHMI (Recherches Translationnelles sur VIH et Maladies Infectieuses), Université de Montpellier, Institut de Recherche pour le Développement, INSERM, 34394 Montpellier, France; audrey.lacroix@ird.fr (A.L.); christel.butel@ird.fr (C.B.); laetitia.serrano@ird.fr (L.S.); nicole.vidal@ird.fr (N.V.); guillaume.thaurignac@ird.fr (G.T.); amandine.esteban@ird.fr (A.E.); eric.delaporte@ird.fr (E.D.); ahidjo.ayouba@ird.fr (A.A.)</footnote>

<footnote n="2">² Institut National de Recherche Biomédicale (INRB), 1197 Kinshasa, Democratic Republic of the Congo; placidembala@inrb.net (P.M.K.); simonp_ndimbok@yahoo.fr (S.P.N.K.); kajosmpiana@gmail.com (G.L.); jjmuyembet@gmail.com (J.-J.M.T.)</footnote>

<footnote n="3">³ Service de Microbiologie, Cliniques Universitaires de Kinshasa, 1197 Kinshasa, Democratic Republic of the Congo; drmukadi@gmail.com</footnote>

<footnote n="4">⁴ Institut National de Recherche Biomédicale (INRB), Goma, Democratic Republic of the Congo</footnote>

<footnote n="5">5 Laboratoire Provincial de Mbandaka, Equateur, Democratic Republic of the Congo; jacqueslikofata@gmail.com</footnote>

<footnote>* Correspondence: martine.peeters@ird.fr (M.P.); amstev04@yahoo.fr (S.A.M.)</footnote>
The largest EVD outbreak occurred in the eastern provinces of North Kivu and Ituri over almost two years (July 2018–June 2020). It infected 3470 people and 2299 (66.2%) died [6,7]. The eleventh outbreak occurred in the Equateur province between June and November 2020, close to areas of the 2018 outbreak in the same province [2,7]. On the 7th February 2021, a new resurgence has been declared in North Kivu, involving 11 cases to date [3,8].

Except for a few recent outbreaks linked to a resurgence in humans [9,10], EVD outbreaks are most likely the result of independent zoonotic transmission event(s) and bats are considered as reservoir species [11–13]. Although there is no direct evidence of exposure to infected bats, two EVD outbreaks have been suspected to be linked to bats; i.e., in Luebo (DRC) in 2007 and the major outbreak in West Africa in 2013 [12,14]. Other findings increase the likelihood of the role of bats in the filovirus ecology. For example, the detection of another filovirus, Marburg virus (MARV), in bats across Africa (Uganda, RDC, Kenya, South Africa, Gabon, Zambia, Sierra Leone) [15–25]. The detection of other filoviruses in bats, such as Lloviu virus in Europe [26,27], filoviruses in bats from China [28,29] and the Philippines [30,31] or Bombali virus (BOMV) in insectivorous bats in Africa provides additional evidence that filoviruses, including Ebola virus, have a Chiropteran origin [32–34]. So far, MARV is the only pathogenic filovirus isolated from bats [18], and Egyptian fruit bats (*Rousettus aegyptiacus*), identified as the reservoir, seem to support virus replication with no apparent disease [35,36].

Zaire Ebola virus (EBOV) RNA and antibodies were detected in three frugivorous bat species (*Epomops franqueti*, *Hypsignathus monstrosus* and *Myonycteris torquata*) during EVD outbreaks in 2003 in Gabon and the Republic of Congo [13]. Surveys of bats in West, Central and East Africa have revealed the presence of antibodies to Ebolaviruses in at least eight frugivorous and one insectivorous species (*Mops condylurus*) [13,16,37–40]. Only a few studies have been performed on bats and ebolaviruses in DRC [40–43]. Our previous serological survey in the Western province of Bas Congo where no previous EVD outbreaks have been reported, investigated a total of 830 bats, including 428 frugivorous and 402 insectivorous bats, [40]. We identified antibodies in *E. helvum* but no viral RNA was detected.

Bats may represent a source of pathogen spillover into human populations in many African countries through their hunting and butchering for bushmeat or through indirect exposure to fruits contaminated by infected saliva, urine or faeces [44–46]. Nevertheless, intermediate amplifying hosts may also play a major role in the emergence of EVD outbreaks, as illustrated in Gabon and Ivory Coast where apes have been confirmed as the source of infection in several outbreaks [11,47]. Similar to humans, bats can also transmit their viruses to other mammals through contact with fruits contaminated with secretions or by hunting, for example, non-human primates (monkeys and bonobos) are reported to hunt bats [48,49].

Today, the ecology of Ebolaviruses is still poorly understood and the role of bats in outbreaks needs to be further clarified. In order to increase the chances to detect viral RNA and Ebola antibodies, we investigated the presence of RNA from filoviruses and Ebola antibodies in bats from DRC collected during the EVD outbreaks in the Equateur and North Kivu provinces and in other areas that have already experienced outbreaks.

**2. Results**

2.1. *Samples and Species per Collection Site*

Between February 2018 and October 2019, a total of 1007 bats were sampled in four provinces with previous and ongoing Ebola outbreaks (Figure 1, Table S1). Specifically, 288 bats were sampled during the ninth EVD outbreak in May–July 2018 in the Equateur province in three sites were Ebola cases were reported, i.e., Bikoro (BK), Ingende (IG) and Iboko (IB). A total of 453 bats were sampled between December 2018 and June 2019 during the tenth EVD outbreak in the Eastern provinces of North Kivu and Ituri, i.e., Komanda (KM), Mangina (MG), Beni (BN), Butembo (BT). Additionally, 98 and 168 bats
were sampled around Kikwit (KK) in February 2018 and in Boende (BD) during September-October 2019, where the third and the seventh EVD outbreaks occurred in 1995 and 2014, respectively. Bats were captured in different ecological environments including 839 (83.3%) in anthropized sites, i.e., villages (14.4%) or cities (68.9%), and 168 in forest sites (16.6%).

The identification of the species was confirmed on 903 (90%) of the 1007 collected bats, 807 cytb and 96 12S sequences were obtained. The species identification of the remaining samples was extrapolated as described in the methods. For some insectivorous bats (Molossidae, Hipposideridae, Vespertilionidae, Nycteridae), identification was only possible at the genus level due to the lack of reference sequences in Genbank. Details on bat families, genera and species collected at each site are shown in Table 1. A total of 892 (88.6%) samples were from frugivorous bats belonging to 11 genera and 115 (11.4%) were from insectivorous bats, representing seven genera from five families. Overall, 48.4% of the samples were from female bats and 51.6% were from males, 97.4% were from adult bats and 2.6% were from juveniles.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Collection sites of bat samples. Sites where bat samples were collected are highlighted with yellow circles. Grey stars represent filovirus outbreaks in DRC. Sites are abbreviated as follows: BK: Bikoro; IB: Iboko; IG: Ingende; BN: Beni; BT: Butembo; KM: Komanda; MG: Mangina; KK: Kikwit; BD: Boende.</figcaption></figure>
<table><caption>**Table 1.** Details and numbers of species that were collected at each site. Two letter code of sites corresponds to codes on Figure 1, -means that no samples were collected. Sites are abbreviated as follows: BK: Bikoro; IB: Iboko; IG: Ingende; BN: Beni; BT: Butembo; KM: Komanda; MG: Mangina; KK: Kikwit; BD: Boende.</caption><thead><tr><th rowspan="2"></th><th colspan="8">Equateur</th><th>North-Kivu/Ituri</th><th>Kikwit</th><th>Boende</th><th rowspan="2">Total</th></tr><tr><th>BK</th><th>IB</th><th>IG</th><th>subtotal</th><th>BN</th><th>BT</th><th>KM</th><th>MG</th><th>subtotal</th><th>KK</th><th>subtotal</th><th>BD</th><th>subtotal</th></tr></thead><tbody><tr><td colspan="15"><strong>FRUGIVOROUS BATS</strong></td></tr><tr><td colspan="15"><strong>Family PTEROPODIDAE</strong></td></tr><tr><td>Casinycteris argynnis</td><td>1</td><td>11</td><td>-</td><td>12</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>4</td><td>4</td><td>16</td></tr><tr><td>Eidolon helvoum</td><td>2</td><td>-</td><td>3</td><td>5</td><td>2</td><td>-</td><td>2</td><td>10</td><td>14</td><td>-</td><td>130</td><td>130</td><td>149</td></tr><tr><td>Epomophorus sp.</td><td>2</td><td>-</td><td>-</td><td>2</td><td>185</td><td>34</td><td>4</td><td>43</td><td>266</td><td>-</td><td>-</td><td>-</td><td>268</td></tr><tr><td>Epomops franqueti</td><td>20</td><td>57</td><td>27</td><td>104</td><td>-</td><td>-</td><td>21</td><td>35</td><td>56</td><td>2</td><td>2</td><td>17</td><td>17</td></tr><tr><td>Hypsignathus monstrosus</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>5</td><td>6</td><td>-</td><td>5</td><td>5</td><td>11</td></tr><tr><td>Lissonycteris angolensis</td><td>-</td><td>-</td><td>-</td><td>-</td><td>23</td><td>1</td><td>2</td><td>-</td><td>26</td><td>-</td><td>-</td><td>-</td><td>26</td></tr><tr><td>Megaloglossus woermanni</td><td>-</td><td>4</td><td>-</td><td>4</td><td>-</td><td>-</td><td>1</td><td>-</td><td>1</td><td>5</td><td>5</td><td>1</td><td>11</td></tr><tr><td>Micropteropus pusillus</td><td>118</td><td>-</td><td>-</td><td>118</td><td>3</td><td>-</td><td>3</td><td>-</td><td>6</td><td>33</td><td>33</td><td>3</td><td>160</td></tr><tr><td>Myonycteris torquata</td><td>11</td><td>16</td><td>2</td><td>29</td><td>-</td><td>-</td><td>1</td><td>1</td><td>2</td><td>30</td><td>30</td><td>5</td><td>66</td></tr><tr><td>Rousettus aegyptiacus</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>4</td><td>-</td><td>5</td><td>-</td><td>-</td><td>-</td><td>5</td></tr><tr><td>Scotonycteris bergmansi</td><td>-</td><td>1</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td><strong>Subtotal frugivorous bats</strong></td><td>154</td><td>89</td><td>32</td><td>275</td><td>213</td><td>36</td><td>39</td><td>94</td><td>382</td><td>70</td><td>70</td><td>165</td><td>892</td></tr><tr><td colspan="15"><strong>INSECTIVOROUS BATS</strong></td></tr><tr><td colspan="15"><strong>Family EMBALLONURIDAE</strong></td></tr><tr><td>Taphozous mauritianus</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>7</td><td>-</td><td>7</td><td>-</td><td>-</td><td>-</td><td>7</td></tr><tr><td colspan="15"><strong>Family HIPPOSIDERIDAE</strong></td></tr><tr><td>Hipposideros caffer</td><td>-</td><td>-</td><td>-</td><td>-</td><td>61</td><td>-</td><td>-</td><td>-</td><td>61</td><td>-</td><td>-</td><td>-</td><td>61</td></tr><tr><td>Hipposideros sp.</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td><td>-</td><td>1</td></tr><tr><td colspan="15"><strong>Family MOLOSSIDAE</strong></td></tr><tr><td>Chaerephon sp.</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>6</td><td>6</td><td>1</td><td>7</td></tr><tr><td>indeterminate</td><td>2</td><td>-</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>Mops sp.</td><td>11</td><td>-</td><td>-</td><td>11</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>11</td></tr><tr><td colspan="15"><strong>Family NYCTERIDAE</strong></td></tr><tr><td>Nycteris arge</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td></tr><tr><td>Nycteris sp.</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td></tr><tr><td colspan="15"><strong>Family VESPERTILIONIDAE</strong></td></tr><tr><td>Neoromicia nanus</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td><td>2</td><td>-</td><td>2</td></tr><tr><td>Neoromicia sp.</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>4</td><td>4</td><td>-</td><td>4</td></tr><tr><td>Scotophilus dinganii</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>3</td><td>-</td><td>3</td><td>14</td><td>14</td><td>-</td><td>17</td></tr><tr><td>Scotophilus nux</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td><td>-</td><td>1</td></tr><tr><td><strong>Subtotal insectivorous bats</strong></td><td>13</td><td>-</td><td>-</td><td>13</td><td>61</td><td>-</td><td>10</td><td>-</td><td>71</td><td>28</td><td>28</td><td>3</td><td>115</td></tr><tr><td><strong>TOTAL</strong></td><td>167</td><td>89</td><td>32</td><td>288</td><td>274</td><td>36</td><td>49</td><td>94</td><td>453</td><td>98</td><td>98</td><td>168</td><td>1007</td></tr></tbody></table>
## 2.2. Molecular Screening for Filoviruses

A total of 1121 swab samples from 676 bats, all from areas of ongoing outbreaks, were screened for the presence of viral RNA. Depending on the availability, oral and rectal swabs were both tested for 425 (63%) bats, and only oral and rectal swabs were tested for 56 (8%) and 193 (29%) bats, respectively (Table 2). Almost all the samples collected during the two EVD outbreaks were tested; i.e., 288/288 (100%) and 388/453 (86%) of the samples from the outbreak in the Equateur and North-Kivu/Ituri provinces, respectively. For 65 bats from the North Kivu outbreak, including all samples from Komanda (KM), swab samples were not available. None of the 1121 swabs tested positive for filovirus viral RNA. For 11 samples, PCR products at the expected size were observed, but MinIon sequencing did not reveal the presence of Ebolavirus sequences. The 11 samples were also tested using the VP35 ebolavirus specific RT-PCR, and were negative.

<table><caption>**Table 2.** Numbers of samples from each species tested for presence of viral RNA by RT-PCR as described in methods, details are shown on number of samples that were tested for oral or rectal swabs only or for both. -; not tested because not available.</caption><thead><tr><th></th><th>N Tested</th><th>Oral Swab</th><th>Oral+Rectal Swab</th><th>Rectal Swab</th></tr></thead><tbody><tr><td colspan="5"><strong>FRUGIVOROUS BATS</strong></td></tr><tr><td colspan="5"><strong>Family PTEROPODIDAE</strong></td></tr><tr><td>Casinycteris argynnis</td><td>12</td><td>-</td><td>-</td><td>12</td></tr><tr><td>Eidolon helvum</td><td>17</td><td>-</td><td>15</td><td>2</td></tr><tr><td>Epomophorus sp.</td><td>263</td><td>-</td><td>261</td><td>2</td></tr><tr><td>Epomops franqueti</td><td>126</td><td>2</td><td>63</td><td>61</td></tr><tr><td>Hypsignathus monstrosus</td><td>5</td><td>-</td><td>5</td><td>-</td></tr><tr><td>Lissonycteris angolensis</td><td>24</td><td>-</td><td>22</td><td>2</td></tr><tr><td>Megaloglossus woermanni</td><td>4</td><td>-</td><td>-</td><td>4</td></tr><tr><td>Micropteropus pusillus</td><td>121</td><td>-</td><td>48</td><td>73</td></tr><tr><td>Myonycteris torquata</td><td>30</td><td>-</td><td>5</td><td>25</td></tr><tr><td>Rousettus aegyptiacus</td><td>1</td><td>-</td><td>1</td><td>-</td></tr><tr><td>Scotonycteris bergmansi</td><td>1</td><td>-</td><td>-</td><td>1</td></tr><tr><td><strong>Subtotal frugivorous bats</strong></td><td><strong>604</strong></td><td><strong>2</strong></td><td><strong>420</strong></td><td><strong>182</strong></td></tr><tr><td colspan="5"><strong>INSECTIVOROUS BATS</strong></td></tr><tr><td colspan="5"><strong>Family EMBALLONURIDAE</strong></td></tr><tr><td>Taphozous mauritianus</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="5"><strong>Family HIPPOSIDERIDAE</strong></td></tr><tr><td>Hipposideros caffer</td><td>59</td><td>56</td><td>3</td><td>-</td></tr><tr><td>Hipposideros sp.</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="5"><strong>Family MOLOSSIDAE</strong></td></tr><tr><td>Chaerephon sp.</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td>indeterminate</td><td>2</td><td>-</td><td>2</td><td>-</td></tr><tr><td>Mops sp.</td><td>11</td><td>-</td><td>-</td><td>11</td></tr><tr><td colspan="5"><strong>Family NYCTERIDAE</strong></td></tr><tr><td>Nycteris arge</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Nycteris sp.</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="5"><strong>Family VESPERTILIONIDAE</strong></td></tr><tr><td>Neoromicia nanus</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Neoromicia sp.</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Scotophilus dinganii</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Scotophilus nux</td><td>0</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Subtotal insectivorous bats</strong></td><td><strong>72</strong></td><td><strong>56</strong></td><td><strong>5</strong></td><td><strong>11</strong></td></tr><tr><td><strong>TOTAL</strong></td><td><strong>676</strong></td><td><strong>58</strong></td><td><strong>425</strong></td><td><strong>193</strong></td></tr></tbody></table>

## 2.3. Serological Screening for Presence of Antibodies against Ebolaviruses

We tested dried blood spots (DBS) from 925/1007 (91.9%) bats for Ebolavirus antibodies. All samples from Boende and Kikwit and 99.6% and 82.1% of those collected during the outbreaks in Equateur Province and eastern provinces of DRC, respectively,
were tested. As in our previous study, we used a range to express the number and percentages of reactive samples according to different cutoff calculations as described in the methods and shown in Supplementary Table S2 and Supplementary material 4. According to our previously defined positivity criteria, i.e., simultaneous reactivity for NP and GP for the same Ebolavirus species, only one *Epomophorus* bat collected in Mangina (MG) in the North-Kivu province, tested positive for Sudan virus (SUDV) using the mean + 4SD method (Figure 1, Table 3). As shown in Table 3, reactivity to glycoprotein antigens of the different Ebolavirus species was highest, ranging between 0.2% and 0.7% to 4.0% and 5.7% for Ebola (EBOV), Sudan (SUDV) or Bundibugyo (BDBV). None of the samples reacted with GP from Reston virus (RESTV). Reactivities were lower with the VP40 antigens and weakest with nucleoprotein (Table 3). In addition, simultaneous reactivity to the same antigen from different virus lineages was observed for the GPs for 0.4% to 5.2% (4–48) of samples (data not shown). Another simultaneous reactivity to at least 2 antigens from the same virus lineage was observed only for the SUDV lineage using the less stringent mean+4SD method; i.e., two *Rousettus aegyptiacus* and two *Eidolon helvum* were reactive to both GP and VP40. All insectivorous bat samples were negative using both cutoffs for all antigens.

<table><caption>**Table 3.** Numbers and percentages of samples with antibodies for each Ebolavirus antigen in the Luminex assay are shown at each site and in total. Numbers and percentages are expressed as a range, corresponding to values obtained by stringent or less stringent cut-off calculations as described in methods. The assay used recombinant proteins of Nucleoprotein (NP), Viral Protein-40 (VP40) or Glycoprotein (GP) for different Ebola virus lineages: Zaire (EBOV), Sudan (SUDV), Bundibugyo (BDBV) and Reston (RESTV). GP proteins from the Mayinga (GP-M) and the Kissidougou (GP-K) strain were used for EBOV.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th>Kikwit</th><th>Boende</th><th>Equateur</th><th>North Kivu/Ituri</th><th>Total</th></tr><tr><th>n = 98<br>N (%) pos</th><th>n = 168<br>N (%) pos</th><th>n = 287<br>N (%) pos</th><th>n = 372<br>N (%) pos</th><th>n = 925<br>N (%) pos</th></tr></thead><tbody><tr><td rowspan="5">EBOV</td><td>NP</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>1-1 (0.3-0.3)</td><td>0-0 (0-0)</td><td>1-1 (0.1-0.1)</td></tr><tr><td>GP-K</td><td>1-3 (1-3)</td><td>1-21 (0.6-30.9)</td><td>0-4 (0-1.4)</td><td>3-19 (0.8-5.1)</td><td>5-47 (0.5-5.1)</td></tr><tr><td>GP-M</td><td>0-4 (0-4)</td><td>1-8 (0.6-4.8)</td><td>0-2 (0.0.7)</td><td>2-7 (0.5-1.9)</td><td>3-21 (0.3-2.3)</td></tr><tr><td>VP40</td><td>0-0 (0-0)</td><td>0-4 (0.0-2.4)</td><td>0-2 (0.0.7)</td><td>0-0 (0-0)</td><td>0-6 (0-0.6)</td></tr><tr><td>NP+GP</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td></tr><tr><td rowspan="4">SUDV</td><td>NP</td><td>0-1 (0-1)</td><td>0-0 (0-0)</td><td>0-1 (0-0.3)</td><td>1-2 (0.3-0.5)</td><td>1-4 (0.1-0.4)</td></tr><tr><td>GP</td><td>0-3 (0-3)</td><td>4-42 (2.4-25)</td><td>0-3 (0-1.1)</td><td>3-5 (0.4-1.3)</td><td>7-53 (0.7-5.7)</td></tr><tr><td>VP40</td><td>0-2 (0-2)</td><td>0-1 (0-0.6)</td><td>0-1 (0-0.3)</td><td>0-4 (0-1.1)</td><td>2-8 (0.2-0.8)</td></tr><tr><td>NP+GP</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-1 (0-0.3)</td><td>0-1 (0-0.1)</td></tr><tr><td rowspan="2">BDBV</td><td>GP</td><td>0-2 (0-2)</td><td>1-18 (0.6-10.7)</td><td>0-3 (0-1.1)</td><td>1-14 (0.3-3.7)</td><td>2-37 (0.2-4.0)</td></tr><tr><td>VP40</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-2 (0-0.2)</td></tr><tr><td>RESTV</td><td>GP</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td><td>0-0 (0-0)</td></tr></tbody></table>

Because the simultaneous presence of antibodies against NP and GP recombinant proteins might be too stringent as criteria, we also examined the number of samples that were reactive with at least one antigen from one of the Ebolavirus species. A total of eleven (1.4%) to 100 (12.3%) samples from frugivorous bats belonging to 4 and 6 species tested positive for at least one antigen (Table 4, Supplementary Table S3). Reactive samples were observed in the four sites studied; 2 (2%) to 8 (8.2%) in Kikwit, 1 (0.3%) to 8 (2.8%) at the Equateur Province, 3 (0.8%) to 39 (10.5%) in the Eastern region and 5 (3%) to 45 (26.8%) in Boende. The reactive samples were from fruit bat species that were previously documented to harbor Ebola antibodies, i.e., *Eidolon helvum*, *Epomophorus* sp., *Epomops franqueti*, *Micropteropus pusillus*, *Myonycteris torquata*, *Rousettus aegyptiacus* (Table 4). Due to the heterogeneity of species per site and differences in numbers per site, it is difficult to draw meaningful comparisons and conclusions about whether antibody levels may differ between sites and sample collection period.
<table><caption>**Table 4.** Numbers and percentages of samples with antibodies to at least one Ebolavirus antigen in the Luminex assay are shown for each species at each site and in total. Numbers and percentages are expressed as a range, corresponding to values obtained by stringent or less stringent cut-off calculations as described in methods. na, means not applicable because no samples of the species were collected at the given site. Details on reactivity to each Ebola antigen are provided in the Supplementary Table S3.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Equateur</th><th colspan="2">North Kivu/Ituri</th><th colspan="2">Kikwit</th><th colspan="2">Boende</th><th colspan="2">Total</th></tr><tr><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th></tr></thead><tbody><tr><td colspan="11"><strong>FRUGIVOROUS BATS</strong></td></tr><tr><td colspan="11"><strong>Family PTEROPODIDAE</strong></td></tr><tr><td>Casinycteris argynnis</td><td>11</td><td>0 (0)</td><td>na</td><td></td><td>na</td><td></td><td>4</td><td>0 (0)</td><td>15</td><td>0 (0)</td></tr><tr><td>Eidolon helvum</td><td>5</td><td>0 (0)</td><td>14</td><td>2-6<br>(14.3-42.9)</td><td>na</td><td></td><td>130</td><td>5-44<br>(3.8-33.8)</td><td>149</td><td>7-50 (4.7-33.6)</td></tr><tr><td>Epomophorus sp.</td><td>2</td><td>0 (0)</td><td>187</td><td>0-22 (0-11.8)</td><td>na</td><td></td><td>na</td><td></td><td>189</td><td>0-22 (0-11.6)</td></tr><tr><td>Epomops franqueti</td><td>104</td><td>0-3 (0-2.9)</td><td>56</td><td>0-5 (0-8.9)</td><td>2</td><td>0 (0)</td><td>17</td><td>0 (0)</td><td>179</td><td>0-8 (0-4.5)</td></tr><tr><td>Hypsignathus monstrosus</td><td>na</td><td></td><td>6</td><td>0 (0)</td><td>na</td><td></td><td>5</td><td>0 (0)</td><td>11</td><td>0 (0)</td></tr><tr><td>Lisonycteris angolensis</td><td>na</td><td></td><td>26</td><td>0 (0)</td><td>na</td><td></td><td>na</td><td></td><td>26</td><td>0 (0)</td></tr><tr><td>Megaloglossus woermanni</td><td>4</td><td>0 (0)</td><td>1</td><td>0 (0)</td><td>5</td><td>0 (0)</td><td>1</td><td>0 (0)</td><td>11</td><td>0 (0)</td></tr><tr><td>Micropteropus pusillus</td><td>118</td><td>1-4 (0.8-3.4)</td><td>4</td><td>0-2 (0-50)</td><td>33</td><td>1-4 (3-12.1)</td><td>3</td><td>0 (0)</td><td>158</td><td>2-10 (1.2-6.3)</td></tr><tr><td>Myonycteris torquata</td><td>29</td><td>0-1 (0-3.4)</td><td>2</td><td>0 (0)</td><td>30</td><td>1-4 (3.3-13.3)</td><td>5</td><td>0-1 (0-20)</td><td>66</td><td>1-6 (1.5-9.1)</td></tr><tr><td>Rouettus aegyptiacus</td><td>na</td><td></td><td>5</td><td>1-4 (20-80)</td><td>na</td><td></td><td>na</td><td></td><td>5</td><td>1-4 (20-80)</td></tr><tr><td>Scotonycteris bergmansi</td><td>1</td><td>0 (0)</td><td>na</td><td></td><td>na</td><td></td><td>na</td><td></td><td>1</td><td>0 (0)</td></tr><tr><td><strong>Subtotal frugivorous bats</strong></td><td><strong>274</strong></td><td><strong>1-8 (0.4-2.9)</strong></td><td><strong>301</strong></td><td><strong>3-39 (1-13)</strong></td><td><strong>70</strong></td><td><strong>2-8 (2.9-11.4)</strong></td><td><strong>165</strong></td><td><strong>5-45 (3-27.3)</strong></td><td><strong>810</strong></td><td><strong>11-100<br>(1.4-12.3)</strong></td></tr><tr><td colspan="11"><strong>INSECTIVOROUS BATS</strong></td></tr><tr><td colspan="11"><strong>Family EMBALLONURIDAE</strong></td></tr><tr><td>Taphozous mauritianus</td><td>na</td><td></td><td>7</td><td>0 (0)</td><td>na</td><td></td><td>na</td><td></td><td>7</td><td>0 (0)</td></tr><tr><td colspan="11"><strong>Family HIPPOSIDERIDAE</strong></td></tr><tr><td>Hipposideros caffer</td><td>na</td><td></td><td>61</td><td>0 (0)</td><td>na</td><td></td><td>na</td><td></td><td>61</td><td>0 (0)</td></tr><tr><td>Hipposideros sp.</td><td>0</td><td></td><td>na</td><td></td><td>1</td><td>0 (0)</td><td>na</td><td></td><td>1</td><td>0 (0)</td></tr><tr><td colspan="11"><strong>Family MOLOSSIDAE</strong></td></tr><tr><td>Chaerephon sp.</td><td>na</td><td></td><td>na</td><td></td><td>6</td><td>0 (0)</td><td>1</td><td>0 (0)</td><td>7</td><td>0 (0)</td></tr><tr><td>Indeterminate</td><td>2</td><td>0 (0)</td><td>na</td><td></td><td>na</td><td></td><td>na</td><td></td><td>2</td><td>0 (0)</td></tr><tr><td>Mops sp.</td><td>11</td><td>0 (0)</td><td>na</td><td></td><td>na</td><td></td><td>na</td><td></td><td>2</td><td>0 (0)</td></tr></tbody></table>
<table><caption>Table 4. *Cont.*</caption><thead><tr><th rowspan="2"></th><th colspan="2">Equateur</th><th colspan="2">North Kivu/Ituri</th><th colspan="2">Kikwit</th><th colspan="2">Boende</th><th colspan="2">Total</th></tr><tr><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th><th>N Tested</th><th>N (%) pos</th></tr></thead><tbody><tr><td><strong>Family NYCTERIDAE</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Nycteris arge</td><td>na</td><td></td><td>na</td><td></td><td>na</td><td></td><td>1</td><td>0 (0)</td><td>1</td><td>0 (0)</td></tr><tr><td>Nycteris sp.</td><td>na</td><td></td><td>na</td><td></td><td>na</td><td></td><td>1</td><td>0 (0)</td><td>1</td><td>0 (0)</td></tr><tr><td><strong>Family VESPERTILIONIDAE</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Neoromicia nanus</td><td>na</td><td></td><td>na</td><td></td><td>2</td><td>0 (0)</td><td>na</td><td></td><td>2</td><td>0 (0)</td></tr><tr><td>Neoromicia sp.</td><td>na</td><td></td><td>na</td><td></td><td>4</td><td>0 (0)</td><td>na</td><td></td><td>4</td><td>0 (0)</td></tr><tr><td>Scotophilus dinganii</td><td>na</td><td></td><td>3</td><td>0 (0)</td><td>14</td><td>0 (0)</td><td>na</td><td></td><td>17</td><td>0 (0)</td></tr><tr><td>Scotophilus nux</td><td>na</td><td></td><td>na</td><td></td><td>1</td><td>0 (0)</td><td>na</td><td></td><td>1</td><td>0 (0)</td></tr><tr><td><strong>Subtotal insectivorous bats</strong></td><td>13</td><td>0 (0)</td><td>71</td><td>0</td><td>28</td><td>0 (0)</td><td>3</td><td>0 (0)</td><td>115</td><td>0 (0)</td></tr><tr><td><strong>Total</strong></td><td>287</td><td>1–8 (0.3–2.8)</td><td>372</td><td>3–39<br>(0.8–10.5)</td><td>98</td><td>2–8 (2–8.2)</td><td>168</td><td>5–45 (3–26.8)</td><td>925</td><td>11–100<br>(1.2–10.8)</td></tr></tbody></table>
### 3. Discussion

Although there is evidence for the involvement of bats in the ecology of Ebola viruses, their exact role is still not elucidated. Several studies reported presence of antibodies but only one study in Gabon identified viral RNA in addition to antibodies in a handful of animals from three species of fruit bats [13]. To clarify the meaning of Ebola virus antibodies, it is essential to document the extent to which viral RNA and shedding can be detected in species with antibodies.

In order to increase our chances to identify viral RNA in bats, we focused our efforts on studying bats as early as possible during outbreaks and tested samples from 741 bats collected during two recent EVD outbreaks and 266 from two regions that have already experienced EVD outbreaks in DRC. Overall, samples were thus collected from 1007 bats, for 903 the species identification was confirmed by sequence analyses, swabs of 676 bats were tested by RT-PCR for the presence of viral RNA, and 925 bats were tested for the presence of antibodies to Ebolavirus. We also focused on frugivorous bats because previous studies on more than 8000 bats sampled in Africa, showed higher rates of Ebolavirus antibodies in frugivorous than in insectivorous bats. Today, antibodies to Zaire Ebolavirus have been detected in eight species of frugivorous bats (*E. franqueti*, *H. monstrosus*, *M. torquata*, *E. helvum*, *E. gambianus*, *R. aegyptiacus*, *M. pusillus*, *L. angolensis*) and only in one genus (*Mops* sp) of insectivorous bats [13, 16, 37–40, 50]. Unfortunately, we were unable to document the presence of viral RNA in oral and/or rectal swabs from the 604 frugivorous and 72 insectivorous bats sampled during two EVD outbreaks.

Our study includes 197 samples from the three species previously shown to harbor viral RNA, i.e.; *Epomops franqueti* (n = 160), *Myonycteris torquata* (n = 31) and *Hypsignathus monstrosus* (n = 6). Compared to the study conducted during EVD outbreaks in 2003 in Gabon and the Republic of Congo, the total number of samples that we tested here is comparable for *E. franqueti* (160 versus 117 in Gabon), but lower for *H. monstrosus* (6 versus 21) and *M. torquata* (31 versus 141) [13]. Unlike studies in Gabon, where viral RNA was detected in 4.6% (13/279) of the samples from these three species, all of our 197 samples were negative. It is important to note that oral and rectal swabs were tested from bats that were released after sampling, while in the study from Gabon, organs (i.e., liver, spleen, kidney.) of euthanized bats were tested [13]. It is most likely that viral loads are lower in swabs than in organs and influence the result of the PCR tests. So far, infectious ebolavirus has never been isolated from bats.

Ebolaviruses are most likely cleared from their hosts and can therefore only be detected for a limited period of time. Despite the fact that we collected samples during outbreaks, sampling started several weeks after the zoonotic transmission events to the index case, either directly from bats to humans or via an intermediate mammal host. It is therefore possible that the bats sampled in the study have already cleared their viruses. For example, experimental inoculation of *R. aegyptiacus* bats with Zaire Ebolavirus showed antibody development but the detection of viral RNA or shedding was infrequent [36]. In contrast, inoculation of the same species with the Marburg virus showed viremia in organs and oral and rectal shedding. The study showed also that MARV can be horizontally transmitted between bats through direct or indirect contact with infectious body fluids [35,36,51]. This also suggests that the virus may be transmitted to other animals, including humans, by the same mechanisms.

It is also important to note that in the study conducted by Leroy and colleagues in Gabon, no single bat was simultaneously positive for antibodies and RT-PCR [13]. This strengthens the case that active viral replication leading to PCR positivity is transient in infected bats and that virus is quickly cleared, with minimal, if any, overlap with seropositivity. Although, in the study of Forbes et al. on *Mops* bats in Kenya, antibodies to EBOV VP40 antigen were detected in bats that tested positive by PCR for the presence of Bombali virus [52]. Whether the rates of viral replication and infectivity differ between the different bat species in which antibodies to Ebola viruses have been identified needs to be further investigated.
The presence of Ebola antibodies confirms that bats are able to handle infection without serious signs of disease, which has already been experimentally documented, as well as for MARV [36,53]. Despite their remarkable immunity, in the case of Ebolavirus, it is still questionable whether bats are the natural host or if they are maybe an intermediate host. It would thus be worth studying the bat biology more in detail and consider for example signs of the Ebolaviruses in their food or interaction with other susceptible mammals, as putative sources of the presence of ebolavirus in bats.

Our study focused on sampling sessions in areas during outbreaks or with previous outbreaks. We did not take into account the potential variation of virus infection and shedding over the seasons and reproduction period. For example, a study on MARV demonstrated pulses of virus circulation in juvenile Egyptian fruit bats correlating to breeding cycles. Juvenile bats are thus the predominant animals within a population that are susceptible to new infections, and significant increases in MARV infection have been demonstrated within this naïve group [19]. In contrast, the majority of bats in our survey were adults (97.4%) potentially reducing the likelihood to detect virus.

Studying antibodies provides information on a previous acute infection. Using our previously defined positivity criteria, i.e., the simultaneous presence of antibodies against the GP and NP proteins, we only observe one bat (0.2%) out of the 575 frugivorous bats tested for antibodies in outbreak areas that were positive for Sudan Ebolavirus, while the Zaire Ebolavirus was responsible for the current outbreaks. The sample was collected in the North Kivu province. In our previous studies using the same serological assay, we showed that the extent of antibodies to Ebolaviruses varied among species. The highest rates were observed in *E. helvum*, *H. monstrosus*, *M. pusillus* and *R. aegyptiacus*. Except for *M. pusillus* (n = 124), we only had few samples of these species in outbreak areas, 19 *E. helvum*, 6 *H. monstrosus*, and 5 *R. aegyptiacus*, thus reducing the likelihood to identify positive samples. Nevertheless, we identified positive antibody samples in *E. helvum* in our previous study in the western part of DRC and in other areas in Cameroon and Guinea where no EVD outbreak has been documented, illustrating that Ebolaviruses can circulate over a large geographic area across Africa between EVD outbreaks [40]. Previous studies have also shown antibodies against the Sudan virus in several species of frugivorous bats [39,40].

Although 12 of the 31 documented EVD outbreaks occurred in DRC, there are still a relatively limited number of studies reporting Ebola viruses in bats in the country. With this new study, more than 3000 samples from bats have now been studied in DRC at eight different areas with seven that are in locations with known outbreaks [40–43]. In all studies, insectivorous bats were negative for antibodies or viral RNA when tested. Unfortunately, the early studies conducted after original outbreaks around Tandala and Yambuku in 1979 and 1980 and in Kikwit in 1995, only collected a small number of frugivorous bats; 26/426 and 123/539 and all were negative for antibodies and/or virus isolation [41,42]. We report here the first study which tested, by molecular and serological assays, a large number of bats collected during ongoing EVD outbreaks in DRC. Today, only a few other studies have been conducted during outbreaks, i.e., one in Gabon and one in southeast Guinea, where the west African outbreak started [13,54].

These studies in DRC and in other areas of Africa illustrate the difficulty of documenting the role of bats as a reservoir of Ebolaviruses and their role in the ecology of these viruses. The challenges relate to the difficulties of sampling bats in their natural environment, identifying viral RNA in a bat since it is very likely that bats can clear the virus and to the different assays used for antibody testing and to the different antibody detection tests and positivity criteria used. By analogy with observations in EVD survivors, we used the simultaneous presence of antibodies against NP and GP as a criterion of positivity in order to increase specificity as described previously [40]. However, it cannot be excluded that this criterion is too strict, because the dynamics of antibodies directed against different antigens is not known in bats.

Future studies should probably focus on regular surveys of bat colonies suspected of being infected in order to increase the likelihood of detecting viral RNA. The longitudinal
following may be even more relevant regarding the previously described seasonality of MARV, another filovirus [19]. Given the increasing frequency of EVD outbreaks and their potential impact, studies on the ecology and animal reservoir of Ebolaviruses are now urgently needed.

# 4. Materials and Methods

## 4.1. *Collection Sites*

Samples were collected in four provinces between February 2018 and October 2019 from free-ranging bats in various locations during ongoing EVD outbreaks and in areas where outbreaks have already been documented in DRC. Samples were collected as described previously [40]. Briefly, bats were captured using mist nets or harp traps in roosting and foraging sites. Bats were released immediately after sampling for ethical and conservation purposes. Dried blood spots (DBS), rectal and oral swabs were collected and preserved as described previously [40,55]. Samples were stored in the field at ambient temperature then at -20 °C in the laboratory. Data on capture sites (GPS coordinates, ecological environment), capture method, morphology (body measurements, weight, color), sex, age class (adult, juvenile) and visual identification of species were recorded for each bat sampled.

Permission to conduct research and to collect samples was obtained from the Ministry of Health and Ministry of Environment and the national ethics committee from the DRC (ESP/CE/009/2016).

## 4.2. *Screening for Ebola Virus Antibodies*

Dried blood spots (DBS) were tested with a Luminex-based serological assay adapted for bats as previously described, with minor modifications [40]. The assay included a total of 10 recombinant Ebola virus proteins, i.e., glycoprotein (GP), nucleoprotein (NP), or viral protein 40 (VP40) for four Ebola virus species: Zaire (EBOV), Sudan (SUDV), Bundibugyo (BDBV) and Reston (RESTV). The Taï Forest Ebola species was not included in the assay due to a lack of available recombinant protein. Briefly, whole blood from DBS was reconstituted as previously described [40,56] and 100 µL of sample adjusted at a final plasma dilution of 1:2000, taking into account the hematocrit, were incubated with 50 µL of magnetic beads coated with recombinant protein (2 µg protein/1.25 × 106 beads) in 96-well flat-bottom chimney plates (Greiner bio one, Frickenhausen, Germany) on a plate shaker at 300 rpm for 16 h at 4 °C in the dark. After washing, 0.1 µg/mL of goat anti-bat biotin-labeled IgG (Euromedex, Souffelweyersheim, France) was added to each well and incubated for 30 min at 300 rpm at room temperature. After washing, we added 50 µL of 4 µg/mL streptavidin-R-phycoerythrin (Fisher Scientific/Life Technologies, Illkirch, France) per well and incubated for 10 min at 300 rpm at room temperature. Reactions were read with BioPlex-200 (BioRad, Marnes-la-Coquette, France) or MagPix (Luminex, Austin, TX, USA). At least 100 events were read for each bead set, and results were expressed as median fluorescence intensity (MFI) per 100 beads. Samples that showed positive signals were repeated in order to validate results.

## 4.3. *Cut-Off Determination*

For each of the 10 antigens, we first used the general formula which used the MFI values of negative control samples as described previously, i.e., mean of 145 negative samples plus 4 × SD (standard deviation) [40]. Negative control samples were from insectivorous and frugivorous bat species born in captivity in Europe [40]. Second, we determined consensus cutoffs as the mean of three cutoffs obtained using three previously described statistical methods used in the absence of well-documented positive controls [16,40,56,57]. The cutoffs were calculated using the same methods as those described by De Nys et al. [40] and are described in details in Supplementary material 4; i.e., change point analysis method and fitted univariate binomial and exponential distributions to define the cutoff at a 0.001 risk for error. For these calculations, we used a larger panel of 8741 bats from Guinea, DRC and Cameroon, which included the samples from our previous study [40] and samples
from subsequent surveys in the same countries. The large number of bat samples used for
the cutoff calculation allowed the robustness of the statistical methods (Supplementary
Table S2, Supplementary material 4). Analyses were performed with R version 4.0.2 soft-
ware (https://www.r-project.org/, accessed on 30 March 2021). We considered a sample
reactive with an antigen if the MFI values were above the cutoff value. We defined the
positivity of Ebola virus antibodies as reactivity to glycoprotein (GP) and nucleoprotein
(NP) of the same lineage, as in our previous studies [40,56,58].

4.4. *Molecular Confirmation of Bat Species*

The species identification recorded in the field was confirmed at the molecular level on a subset of bats using DNA extracted from DBS as described in our previous studies [40,55]. Briefly, a fragment of approximately 800 base pairs (bp) of the mitochondrial cytochrome b (cytb) was amplified using primers adapted from Irwin et al. [59]. In addition, for certain species of bats, primers targeting the 12 S region of ribosomal RNA were used [60]. PCR products were directly sequenced on an ABI 3500 sequencer (Applied Biosystems, Courtaboeuf, France). Sequences were submitted to NCBI for BLAST analysis to identify the most similar bat species. For sequences with no or low similarity (<97%) hits with species in Genbank, phylogenetic tree analysis with reference sequences was performed using maximum likelihood methods with RAxMLv8 [61] implemented in MegAlignPro version 17.2 (DNASTAR. Madison, WI, USA) in order to determine the genus. At least one sample per species and per sampling event was confirmed at the molecular level. Species identification was then extrapolated to the remaining samples by combining the molecular and morphological data records in the field. We also used morphological details on the forearm and weight measurements to distinguish *Epomophorus* sp. and *Micropteropus pusillus*, as previously documented [55,62].

4.5. Nucleic Acid Extraction, RT-PCR Screening and Nanopore Sequencing for Detection of Filoviruses

Total DNA and RNA were extracted from oral and rectal swabs using the NucliSens MiniMAG® system (Biomerieux, Marcy-l'Etoile, France). Briefly, 250 µL of RNAlater swab sample was incubated with 1 mL of lysis buffer for 10 min and extraction was performed using manufacturer's instructions. cDNA was first synthesized from denatured RNA using the Reverse Transcription System kit with random primers (Promega, Madison, WI, USA), following the manufacturer's instructions. PCR screening was performed with a broadly reactive semi-nested PCR targeting a 630 bp fragment of the L gene using degenerated primers that detect *Filoviridae* family level, as previously described [15,32]. Briefly, cDNA was amplified using the GoTaq Hot Start Master Mix PCR kit (Promega, Madison, WI, USA) as follows for first and second PCR rounds: 10 cycles of 92 °C for 20 s, 50 °C for 30 s with -0.5 °C/cycle and 72 °C for 1 min, 35 cycles of 92 °C for 20 s, 50 °C for 30 s and 72 °C for 1 min. After agarose gel electrophoresis, PCR products of the expected size were sequenced by Nanopore sequencing technology. Briefly, PCR products were barcoded and pooled using the Native Barcoding Kit (Oxford Nanopore Technologies, Oxford, UK). Sequencing libraries were generated from the barcoded products using the Genomic DNA Sequencing Kit EXP-NBD104-114/SQK-LSK109 (Oxford Nanopore Technologies) and were loaded onto an R9.4 flow cell. Base-calling, adapter removal and demultiplexing were achieved using MinKnow software, and the obtained Fastq sequences were analyzed using Nanopipe pipeline (http://www.bioinformatics.uni-muenster.de, accessed on 30 March 2021), and compared to other Filovirus genus references. Additionally, samples with bands of expected size with universal primers were also tested using a previously described Ebola Zaire specific PCR assay targeting a 184 bp fragment on the VP35 gene [5].

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/pathogens10050557/s1, Table S1: Collection sites, date of collection, and number of bats sampled in the Democratic Republic of the Congo, and number of serological and molecular tests performed; Table S2: MFI cut-off values calculated with different methods as described in Methods
section for the different Ebolavirus strains [Zaire (EBOV), Sudan (SUDV), Bundibugyo (BDBV) and Reston (RESTV)] and the different antigens used per virus. Table S3: Numbers of samples with antibodies to each Ebolavirus antigen, to each Ebolavirus species and simultaneous NP+GP EBOV and SUDV antigens in the Luminex assay for each tested bat species. Numbers are obtained by stringent (Mean of the cutoffs; CO 1) or less stringent cutoff (Mean + 4SD; CO 2) calculations as described in methods. Supplementary material 4: Determination of Cutoffs.

**Author Contributions:** Conceptualization, S.A.M., A.A. and M.P.; methodology, A.L., A.A., P.M.K., M.P.; formal analysis, A.L., A.A., M.P.; investigation, P.M.K., S.P.N.K., G.L., A.L., G.T., A.E., C.B., L.S., N.V., D.M.B., J.L., A.A.; writing—original draft preparation, A.L., A.A., M.P.; writing—review and editing, P.M.K., A.A., S.A.M., M.P.; supervision, S.A.M., P.M.K., J.-J.M.T., A.A. and M.P.; funding acquisition, A.A., M.P., J.-J.M.T. and E.D. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported in part by grants from Institut National de la Santé et de la Recherche Médicale (INSERM)/the Ebola Task Force/REACTing, EBO-SURSY project funded by the European Union (FOOD/2016/379-660) and Christophe Mérieux Prize 2015 (to J.-J.M.T.). P.M.K. was supported by fellowships from IRD (ARTS, Allocation de Recherche pour une thèse au Sud). BIODIV-AFREID project (BiodivERsA ERA-Net COFUND programme, 2018-2019 joint call for proposals (grant number ANR-19-EBI3-0004).

**Institutional Review Board Statement:** The study was conducted in compliance with the Ministry of Health and Ministry of Environment and the national ethics committee from the DRC (ESP/CE/009/2016).

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data presented in this study are available upon request, and data on bat sampling are available on the EBO-SURSY website (https://rr-africa.oie.int/fr/projets/ebo-sursy-fr/, accessed on 30 March 2021).

**Acknowledgments:** We thank the Ministries of Health, Environment and the National Ethics committee from DRC, for permission to perform this study. We thank all the field staff from DRC for the collection of bat samples. We also thank the provincial authorities of Equateur, North Kivu and Ituri for the field missions.

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

## References

1. Languon, S.; Quaye, O. Filovirus Disease Outbreaks: A Chronological Overview. *Virol. Res. Treat.* **2019**, *10*. [CrossRef] [PubMed]

2. Mahase, E. New Ebola Outbreak in DRC Has Killed 20 People so Far. *BMJ* **2020**, *370*, m2867. [CrossRef]

3. Adepoju, P. Ebola Returns to Guinea and DR Congo. *Lancet* **2021**, *397*, 781. [CrossRef]

4. LaBrunda, M.; Amin, N. The Emerging Threat of Ebola. *Glob. Health Secur.* **2019**, *10*, 103–139. [CrossRef]

5. Mbala Kingebeni, P.; Villabona-Arenas, C.-J.; Vidal, N.; Likofata, J.; Nsio-Mbeta, J.; Makiala-Mandanda, S.; Mukadi, D.; Mukadi, P.; Kumakamba, C.; Djokolo, B.; et al. Rapid Confirmation of the Zaire Ebola Virus in the Outbreak of the Equateur Province in the Democratic Republic of Congo: Implications for Public Health Interventions. *Clin. Infect. Dis.* **2018**, *68*, 330–333. [CrossRef]

6. Wise, J. Ebola: DRC's 10th Outbreak Is Declared Over. *BMJ* **2020**, *369*, m2601. [CrossRef] [PubMed]

7. End of the 11th Ebola Outbreak-Democratic Republic of the Congo. Available online: https://www.afro.who.int/health-topics/ebola-virus-disease/end-ebola-11-drc (accessed on 6 December 2020).

8. Epidémie de La Maladie à Virus Ebola Province Du Nord Kivu, République Démocratique Du Congo Sitrep N°50/2021 (Date de Publication: 29 Mars 2021)-Democratic Republic of the Congo. Available online: https://reliefweb.int/report/democratic-republic-congo/epid-mie-de-la-maladie-virus-ebola-province-du-nord-kivu-r-30 (accessed on 30 March 2021).

9. Keita, A.K.; Düx, A.; Diallo, H.; Calvignac-Spencer, S.; Sow, M.S.; Keita, M.B.; Sidibe, Y.; Ayouba, A.; Lemarcis, F.; Marí Saéz, A.; et al. Resurgence of Ebola Virus in Guinea after 5 Years Calls for Careful Attention to Survivors without Creating Further Stigmatization. Available online: https://virological.org/t/guinea-2021-ebov-genomes/651 (accessed on 30 March 2021).

10. Mbala-Kingebeni, P.; Pratt, C.; Mutafali-Ruffin, M.; Pauthner, M.G.; Bile, F.; Nkuba-Ndaye, A.; Black, A.; Kinganda-Lusamaki, E.; Faye, M.; Aziza, A.; et al. Ebola Virus Transmission Initiated by Relapse of Systemic Ebola Virus Disease. *N. Engl. J. Med.* **2021**, *384*, 1240–1247. [CrossRef]

11. Leroy, E.M.; Rouquet, P.; Formenty, P.; Souquière, S.; Kilbourne, A.; Froment, J.-M.; Bermejo, M.; Smit, S.; Karesh, W.; Swanepoel, R.; et al. Multiple Ebola Virus Transmission Events and Rapid Decline of Central African Wildlife. *Science* **2004**, *303*, 387–390. [CrossRef] [PubMed]
12. Leroy, E.M.; Epelboin, A.; Mondonge, V.; Pourrut, X.; Gonzalez, J.-P.; Muyembe-Tamfum, J.-J.; Formenty, P. Human Ebola Outbreak Resulting from Direct Exposure to Fruit Bats in Luebo, Democratic Republic of Congo, 2007. *Vector Borne Zoonotic Dis.* **2009**, *9*, 723–728. [CrossRef] [PubMed]

13. Leroy, E.M.; Kumulungui, B.; Pourrut, X.; Rouquet, P.; Hassanin, A.; Yaba, P.; Délicat, A.; Paweska, J.T.; Gonzalez, J.-P.; Swanepoel, R. Fruit Bats as Reservoirs of Ebola Virus. *Nature* **2005**, *438*, 575–576. [CrossRef]

14. Marí Saéz, A.; Weiss, S.; Nowak, K.; Lapeyre, V.; Zimmermann, F.; Düx, A.; Kühl, H.S.; Kaba, M.; Regnaut, S.; Merkel, K.; et al. Investigating the Zoonotic Origin of the West African Ebola Epidemic. *EMBO Mol. Med.* **2015**, *7*, 17–23. [CrossRef] [PubMed]

15. Amman, B.R.; Bird, B.H.; Bakarr, I.A.; Bangura, J.; Schuh, A.J.; Johnny, J.; Sealy, T.K.; Conteh, I.; Koroma, A.H.; Foday, I.; et al. Isolation of Angola-like Marburg Virus from Egyptian Rousette Bats from West Africa. *Nat. Commun.* **2020**, *11*, 1–9. [CrossRef] [PubMed]

16. Pourrut, X.; Souris, M.; Towner, J.S.; Rollin, P.E.; Nichol, S.T.; Gonzalez, J.-P.; Leroy, E. Large Serological Survey Showing Cocirculation of Ebola and Marburg Viruses in Gabonese Bat Populations, and a High Seroprevalence of Both Viruses in Rousettus Aegyptiacus. *BMC Infect. Dis.* **2009**, *9*, 159. [CrossRef] [PubMed]

17. Towner, J.S.; Pourrut, X.; Albariño, C.G.; Nkogue, C.N.; Bird, B.H.; Grard, G.; Ksiazek, T.G.; Gonzalez, J.-P.; Nichol, S.T.; Leroy, E.M. Marburg Virus Infection Detected in a Common African Bat. *PLoS ONE* **2007**, *2*, e764. [CrossRef] [PubMed]

18. Towner, J.S.; Amman, B.R.; Sealy, T.K.; Carroll, S.A.R.; Comer, J.A.; Kemp, A.; Swanepoel, R.; Paddock, C.D.; Balinandi, S.; Khristova, M.L.; et al. Isolation of Genetically Diverse Marburg Viruses from Egyptian Fruit Bats. *PLoS Pathog.* **2009**, *5*, e1000536. [CrossRef]

19. Amman, B.R.; Carroll, S.A.; Reed, Z.D.; Sealy, T.K.; Balinandi, S.; Swanepoel, R.; Kemp, A.; Erickson, B.R.; Comer, J.A.; Campbell, S.; et al. Seasonal Pulses of Marburg Virus Circulation in Juvenile Rousettus Aegyptiacus Bats Coincide with Periods of Increased Risk of Human Infection. *PLoS Pathog.* **2012**, *8*, e1002877. [CrossRef] [PubMed]

20. Amman, B.R.; Nyakarahuka, L.; McElroy, A.K.; Dodd, K.A.; Sealy, T.K.; Schuh, A.J.; Shoemaker, T.R.; Balinandi, S.; Atimnedi, P.; Kaboyo, W.; et al. Marburgvirus Resurgence in Kitaka Mine Bat Population after Extermination Attempts, Uganda. *Emerg. Infect. Dis.* **2014**, *20*, 1761–1764. [CrossRef]

21. Swanepoel, R.; Smit, S.B.; Rollin, P.E.; Formenty, P.; Leman, P.A.; Kemp, A.; Burt, F.J.; Grobbelaar, A.A.; Croft, J.; Bausch, D.G.; et al. Studies of Reservoir Hosts for Marburg Virus. *Emerg. Infect. Dis.* **2007**, *13*, 1847–1851. [CrossRef]

22. Kuzmin, I.V.; Niezgoda, M.; Franka, R.; Agwanda, B.; Markotter, W.; Breiman, R.F.; Shieh, W.J.; Zaki, S.R.; Rupprecht, C.E. Marburg Virus in Fruit Bat, Kenya. *Emerg. Infect. Dis.* **2010**, *16*, 352–354. [CrossRef] [PubMed]

23. Storm, N.; Jansen Van Vuren, P.; Markotter, W.; Paweska, J.T. Antibody Responses to Marburg Virus in Egyptian Rousette Bats and Their Role in Protection against Infection. *Viruses* **2018**, *10*. [CrossRef]

24. Kajihara, M.; Hang'ombe, B.M.; Changula, K.; Harima, H.; Isono, M.; Okuya, K.; Yoshida, R.; Mori-Kajihara, A.; Eto, Y.; Orba, Y.; et al. Marburgvirus in Egyptian Fruit Bats, Zambia. *Emerg. Infect. Dis.* **2019**, *25*, 1577–1580. [CrossRef] [PubMed]

25. Paweńska, J.T.; Jansen van Vuren, P.; Kemp, A.; Storm, N.; Grobbelaar, A.A.; Wiley, M.R.; Palacios, G.; Markotter, W. Marburg Virus Infection in Egyptian Rousette Bats, South Africa, 2013–20141. *Emerg. Infect. Dis.* **2018**, *24*, 1134–1137. [CrossRef] [PubMed]

26. Negredo, A.; Palacios, G.; Vázquez-Morón, S.; González, F.; Dopazo, H.; Molero, F.; Juste, J.; Quetglas, J.; Savji, N.; de la Cruz Martínez, M.; et al. Discovery of an Ebolavirus-Like Filovirus in Europe. *PLoS Pathog.* **2011**, *7*, e1002304. [CrossRef]

27. Kemenesi, G.; Kurucz, K.; Dallos, B.; Zana, B.; Földes, F.; Boldogh, S.; Görföl, T.; Carroll, M.W.; Jakab, F. Re-Emergence of Lloviu Virus in Miniopterus Schreibersii Bats, Hungary, 2016. *Emerg. Microbes. Infect.* **2018**, *7*, 66. [CrossRef]

28. Yuan, J.; Zhang, Y.; Li, J.; Zhang, Y.; Wang, L.-F.; Shi, Z. Serological Evidence of Ebolavirus Infection in Bats, China. *Virol. J.* **2012**, *9*, 236. [CrossRef]

29. Yang, X.-L.; Zhang, Y.-Z.; Jiang, R.-D.; Guo, H.; Zhang, W.; Li, B.; Wang, N.; Wang, L.; Waruhiu, C.; Zhou, J.-H.; et al. Genetically Diverse Filoviruses in Rousettus and Eonycteris Spp. Bats, China, 2009 and 2015. *Emerg. Infect. Dis.* **2017**, *23*, 482–486. [CrossRef]

30. Taniguchi, S.; Watanabe, S.; Masangkay, J.S.; Omatsu, T.; Ikegami, T.; Alviola, P.; Ueda, N.; Iha, K.; Fujii, H.; Ishii, Y.; et al. Reston Ebolavirus Antibodies in Bats, the Philippines. *Emerg. Infect. Dis.* **2011**, *17*, 1559. [CrossRef]

31. Jayme, S.I.; Field, H.E.; de Jong, C.; Olival, K.J.; Marsh, G.; Tagtag, A.M.; Hughes, T.; Bucad, A.C.; Barr, J.; Azul, R.R.; et al. Molecular Evidence of Ebola Reston Virus Infection in Philippine Bats. *Virol. J.* **2015**, *12*, 107. [CrossRef]

32. Goldstein, T.; Anthony, S.J.; Gbakima, A.; Bird, B.H.; Bangura, J.; Tremeau-Bravard, A.; Belaganahalli, M.N.; Wells, H.L.; Dhanota, J.K.; Liang, E.; et al. The Discovery of Bombali Virus Adds Further Support for Bats as Hosts of Ebolaviruses. *Nat. Microbiol.* **2018**, *3*, 1084–1089. [CrossRef] [PubMed]

33. Forbes, K.M. Bombali Ebola Virus in Mops Condylurus Bat, Kenya. *Emerg. Infect. Dis.* **J. CDC** **2019**, *25*. [CrossRef]

34. Karan, L.S.; Makenov, M.T.; Korneev, M.G.; Sacko, N.; Boumbaly, S.; Yakovlev, S.A.; Kourouma, K.; Bayandin, R.B.; Gladysheva, A.V.; Shipovalov, A.V.; et al. Bombali Virus in Mops Condylurus Bats, Guinea. *Emerg. Infect. Dis.* **2019**, *25*. [CrossRef]

35. Amman, B.R.; Jones, M.E.B.; Sealy, T.K.; Uebelhoer, L.S.; Schuh, A.J.; Bird, B.H.; Coleman-McCray, J.D.; Martin, B.E.; Nichol, S.T.; Towner, J.S. Oral Shedding of Marburg Virus in Experimentally Infected Egyptian Fruit Bats (Rousettus Aegyptiacus). *J. Wildl. Dis.* **2015**, *51*, 113–124. [CrossRef]

36. Jones, M.E.B.; Schuh, A.J.; Amman, B.R.; Sealy, T.K.; Zaki, S.R.; Nichol, S.T.; Towner, J.S. Experimental Inoculation of Egyptian Rousette Bats (Rousettus Aegyptiacus) with Viruses of the Ebolavirus and Marburgvirus Genera. *Viruses* **2015**, *7*, 3420–3442. [CrossRef] [PubMed]
37. Pourrut, X.; Délicat, A.; Rollin, P.E.; Ksiazek, T.G.; Gonzalez, J.-P.; Leroy, E.M. Spatial and Temporal Patterns of Zaire Ebolavirus Antibody Prevalence in the Possible Reservoir Bat Species. *J. Infect. Dis.* **2007**, *196*, S176–S183. [CrossRef] [PubMed]

38. Hayman, D.T.S.; Yu, M.; Crameri, G.; Wang, L.-F.; Suu-Ire, R.; Wood, J.L.N.; Cunningham, A.A. Ebola Virus Antibodies in Fruit Bats, Ghana, West Africa. *Emerg. Infect. Dis.* **2012**, *18*, 1207–1209. [CrossRef]

39. Ogawa, H.; Miyamoto, H.; Nakayama, E.; Yoshida, R.; Nakamura, I.; Sawa, H.; Ishii, A.; Thomas, Y.; Nakagawa, E.; Matsuno, K.; et al. Seroepidemiological Prevalence of Multiple Species of Filoviruses in Fruit Bats (Eidolon Helvum) Migrating in Africa. *J. Infect. Dis.* **2015**, *212*, S101–S108. [CrossRef]

40. De Nys, H.M.; Kingebeni, P.M.; Keita, A.K.; Butel, C.; Thaurignac, G.; Villabona-Arenas, C.-J.; Lemarcis, T.; Geraerts, M.; Vidal, N.; Esteban, A.; et al. Survey of Ebola Viruses in Frugivorous and Insectivorous Bats in Guinea, Cameroon, and the Democratic Republic of the Congo, 2015–2017. *Emerg. Infect. Dis.* **2018**, *24*. [CrossRef] [PubMed]

41. Breman, J.G.; Teams, E.V.S.; Johnson, K.M.; Teams, E.V.S.; van der Groen, G.; Teams, E.V.S.; Robbins, C.B.; Teams, E.V.S.; Szczeniowski, M.V.; Teams, E.V.S.; et al. A Search for Ebola Virus in Animals in the Democratic Republic of the Congo and Cameroon: Ecologic, Virologic, and Serologic Surveys, 1979–1980. *J. Infect. Dis.* **1999**, *179*, S139–S147. [CrossRef]

42. Leirs, H.; Mills, J.N.; Krebs, J.W.; Childs, J.E.; Akaibe, D.; Woollen, N.; Ludwig, G.; Peters, C.J.; Ksiazek, T.G.; Members, O.S.G. Search for the Ebola Virus Reservoir in Kikwit, Democratic Republic of the Congo: Reflections on a Vertebrate Collection. *J. Infect. Dis.* **1999**, *179*, S155–S163. [CrossRef]

43. Gryseels, S.; Mbala-Kingebeni, P.; Akonda, I.; Angoyo, R.; Ayouba, A.; Baelo, P.; Mukadi, D.B.; Bugenthö, E.; Bushmaker, T.; Butel, C.; et al. Role of Wildlife in Emergence of Ebola Virus in Kaigbono (Likati), Democratic Republic of the Congo, 2017. *Emerg. Infect. Dis.* **2020**, *26*, 2205–2209. [CrossRef] [PubMed]

44. Kamins, A.O.; Rowcliffe, J.M.; Ntiamoa-Baidu, Y.; Cunningham, A.A.; Wood, J.L.N.; Restif, O. Characteristics and Risk Perceptions of Ghanaians Potentially Exposed to Bat-Borne Zoonoses through Bushmeat. *EcoHealth* **2015**, *12*, 104–120. [CrossRef] [PubMed]

45. Han, H.-J.; Wen, H.; Zhou, C.-M.; Chen, F.-F.; Luo, L.-M.; Liu, J.; Yu, X.-J. Bats as Reservoirs of Severe Emerging Infectious Diseases. *Virus Res.* **2015**, *205*, 1–6. [CrossRef] [PubMed]

46. Baudel, H.; de Nys, H.; Mpoudi Ngole, E.; Peeters, M.; Desclaux, A. Understanding Ebola Virus and Other Zoonotic Transmission Risks through Human-Bat Contacts: Exploratory Study on Knowledge, Attitudes and Practices in Southern Cameroon. *Zoonoses Public Health* **2019**, *66*, 288–295. [CrossRef]

47. Formenty, P.; Boesch, C.; Wyers, M.; Steiner, C.; Donati, F.; Dind, F.; Walker, F.; Le Guenno, B. Ebola Virus Outbreak among Wild Chimpanzees Living in a Rain Forest of Côte d'Ivoire. *J. Infect. Dis.* **1999**, *179*, S120–S126. [CrossRef]

48. Tapanes, E.; Detwiler, K.M.; Cords, M. Bat Predation by Cercopithecus Monkeys: Implications for Zoonotic Disease Transmission. *EcoHealth* **2016**, *13*, 405–409. [CrossRef]

49. Bermejo, M.; Rodríguez-Teijeiro, J.D.; Illera, G.; Barroso, A.; Vilà, C.; Walsh, P.D. Ebola Outbreak Killed 5000 Gorillas. *Science* **2006**, *314*, 1564. [CrossRef]

50. Hayman, D.T.S.; Emmerich, P.; Yu, M.; Wang, L.-F.; Suu-Ire, R.; Fooks, A.R.; Cunningham, A.A.; Wood, J.L.N. Long-Term Survival of an Urban Fruit Bat Seropositive for Ebola and Lagos Bat Viruses. *PLoS ONE* **2010**, *5*, e11978. [CrossRef]

51. Schuh, A.J.; Amman, B.R.; Jones, M.E.B.; Sealy, T.K.; Uebelhoer, L.S.; Spengler, J.R.; Martin, B.E.; Coleman-McCray, J.A.D.; Nichol, S.T.; Towner, J.S. Modelling Filovirus Maintenance in Nature by Experimental Transmission of Marburg Virus between Egyptian Rousette Bats. *Nat. Commun.* **2017**, *8*, 14446. [CrossRef]

52. Forbes, K.M.; Webala, P.W.; Jääskeläinen, A.J.; Abdurahman, S.; Ogola, J.; Masika, M.M.; Kivistö, I.; Alburkat, H.; Plyusnin, I.; Levanov, L.; et al. Bombali Virus in Mops Condylurus Bat, Kenya. *Emerg. Infect. Dis.* **2019**, *25*, 955–957. [CrossRef]

53. Swanepoel, R.; Leman, P.A.; Burt, F.J.; Zachariades, N.A.; Braack, L.E.; Ksiazek, T.G.; Rollin, P.E.; Zaki, S.R.; Peters, C.J. Experimental Inoculation of Plants and Animals with Ebola Virus. *Emerg. Infect. Dis.* **1996**, *2*, 321–325. [CrossRef] [PubMed]

54. Leendertz, S.A.J.; Gogarten, J.F.; Düx, A.; Calvignac-Spencer, S.; Leendertz, F.H. Assessing the Evidence Supporting Fruit Bats as the Primary Reservoirs for Ebola Viruses. *EcoHealth* **2016**, *13*, 18–25. [CrossRef]

55. Lacroix, A.; Vidal, N.; Keita, A.K.; Thaurignac, G.; Esteban, A.; de Nys, H.; Diallo, R.; Toure, A.; Goumou, S.; Soumah, A.K.; et al. Wide Diversity of Coronaviruses in Frugivorous and Insectivorous Bat Species: A Pilot Study in Guinea, West Africa. *Viruses* **2020**, *12*. [CrossRef] [PubMed]

56. Ayouba, A.; Touré, A.; Butel, C.; Keita, A.K.; Binetruy, F.; Sow, M.S.; Foulongne, V.; Delaporte, E.; Peeters, M. Development of a Sensitive and Specific Serological Assay Based on Luminex Technology for Detection of Antibodies to Zaire Ebola Virus. *J. Clin. Microbiol.* **2017**, *55*, 165–176. [CrossRef]

57. Peel, A.J.; McKinley, T.J.; Baker, K.S.; Barr, J.A.; Crameri, G.; Hayman, D.T.S.; Feng, Y.-R.; Broder, C.C.; Wang, L.-F.; Cunningham, A.A.; et al. Use of Cross-Reactive Serological Assays for Detecting Novel Pathogens in Wildlife: Assessing an Appropriate Cutoff for Henipavirus Assays in African Bats. *J. Virol. Methods* **2013**, *193*, 295–303. [CrossRef]

58. Ayouba, A.; Ahuka-Mundeke, S.; Butel, C.; Mbala Kingebeni, P.; Loul, S.; Tagg, N.; Arenas, C.-J.V.; Lacroix, A.; Ndimbo-Kumugo, S.-P.; Keita, A.K.; et al. Extensive Serological Survey of Multiple African Non-Human Primate Species Reveals Low Prevalence of IgG Antibodies to Four Ebola Virus Species. *J. Infect. Dis.* **2019**. [CrossRef]

59. Irwin, D.M.; Kocher, T.D.; Wilson, A.C. Evolution of the Cytochromeb Gene of Mammals. *J. Mol. Evol.* **1991**, *32*, 128–144. [CrossRef] [PubMed]
60. Kocher, T.D.; Thomas, W.K.; Meyer, A.; Edwards, S.V.; Pääbo, S.; Villablanca, F.X.; Wilson, A.C. Dynamics of Mitochondrial DNA Evolution in Animals: Amplification and Sequencing with Conserved Primers. *Proc. Natl. Acad. Sci. USA* **1989**, *86*, 6196–6200.
[CrossRef]

61. Stamatakis, A. RAxML Version 8: A Tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. *Bioinformatics* **2014**, *30*, 1312–1313. [CrossRef] [PubMed]

62. Nesi, N.; Nakouné, E.; Cruaud, C.; Hassanin, A. DNA Barcoding of African Fruit Bats (*Mammalia, Pteropodidae*). The Mitochondrial Genome Does Not Provide a Reliable Discrimination between Epomophorus Gambianus and Micropteropus Pusillus. *Comptes Rendus Biol.* **2011**, *334*, 544–554. [CrossRef] [PubMed]