# Genetic diversity and ecology of coronaviruses hosted by cave-dwelling bats in Gabon

Gael Darren Maganga, Anaïs Pinto, Illich Manfred Mombo, Mankomra Madjitobaye, Antoine Mitte Mbeang Beyeme, Larson Boundenga, Meriadeg Ar Gouilh, Nadine N'dilimabaka, Jan Felix Drexler, Christian Drosten, et al.

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

## To cite this version:

Gael Darren Maganga, Anaïs Pinto, Illich Manfred Mombo, Mankomra Madjitobaye, Antoine Mitte Mbeang Beyeme, et al.. Genetic diversity and ecology of coronaviruses hosted by cave-dwelling bats in Gabon. Scientific Reports, 2020, 10 (1), 7314 [13p.]. 10.1038/s41598-020-64159-1. hal-04991954

HAL Id: hal-04991954

https://hal.science/hal-04991954v1

Submitted on 14 Mar 2025

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Genetic diversity and ecology of coronaviruses hosted by cave-dwelling bats in Gabon

Gael Darren Maganga<fnref n="1" />,²✉, Anaïs Pinto<fnref n="1" />, Illich Manfred Mombo<fnref n="1" />, Mankomra Madjitobaye<fnref n="1" />, Antoine Mitte Mbeang Beyeme², Larson Boundenga<fnref n="1" />, Meriadeg Ar GouilhD³, Nadine N’Dilimabaka<fnref n="1" />, Jan Felix Drexler⁴,⁵, Christian Drosten⁴,⁵ & Eric Maurice Leroy<fnref n="1" />,⁶

Little research on coronaviruses has been conducted on wild animals in Africa. Here, we screened a wide range of wild animals collected in six provinces and five caves of Gabon between 2009 and 2015. We collected a total of 1867 animal samples (cave-dwelling bats, rodents, non-human primates and other wild animals). We explored the diversity of CoVs and determined the factors driving the infection of CoVs in wild animals. Based on a nested reverse transcription-polymerase chain reaction, only bats, belonging to the *Hipposideros gigas* (4/156), *Hipposideros cf. ruber* (13/262) and *Miniopterus inflatus* (1/249) species, were found infected with CoVs. We identified alphacoronaviruses in *H. gigas* and *H. cf. ruber* and betacoronaviruses in *H. gigas*. All *Alphacoronavirus* sequences grouped with *Human coronavirus 229E* (HCoV-229E). Ecological analyses revealed that CoV infection was significantly found in July and October in *H. gigas* and in October and November in *H. cf. ruber*. The prevalence in the Faucon cave was significantly higher. Our findings suggest that insectivorous bats harbor potentially zoonotic CoVs; highlight a probable seasonality of the infection in cave-dwelling bats from the North-East of Gabon and pointed to an association between the disturbance of the bats' habitat by human activities and CoV infection.

Coronaviruses (CoVs) belonging to the *Coronaviridae* family are viruses known to infect a wide range of animals and humans. In humans, CoVs are responsible for mild to severe respiratory illnesses such as the Severe Acute Respiratory Syndrome (SARS) and the Middle-East Respiratory Syndrome (MERS). The epidemic of SARS which started in 2002 in China, and spread on various other continents such as North America and Europe, reached a mortality rate of 9%<fnref n="1" />. The other severe epidemic illness (MERS) appeared more recently in the Middle East, and like the SARS, it also spread in other countries in Africa, America and Europe with case fatality rates of 35% (reviewed by de Wit *et al.*²). Currently the world is facing an ongoing pandemic caused by a new coronavirus (SARS-CoV-2) that emerged in China in December 2019 and caused disease (COVID-19) resulting in 118 326 confirmed cases and 4292 deaths as of March 11, 2020³. Since the emergence of these respiratory diseases, CoVs have been considered a real public health problem because of their ability to become epidemic and pandemic⁴. Even though human to human transmission of SARS-CoV and MERS-CoV occurs mainly through nosocomial transmission, some human cases recorded were occasionally the results of a zoonotic transmission⁵. This raised an interest in the identification of animal reservoirs. Indeed, it appears that almost each human CoV have zoonotic origins or otherwise have a close relative that circulate in wild animals (bats)⁶⁻⁸ and domestic animals (camels and cattle)⁹,¹⁰. Besides SARS-CoVs and MERS-CoVs, numerous other CoVs have been detected in bats in Africa, Asia, Europe and America and are classified into the genera *Alphacoronavirus* and *Betacoronavirus*¹¹,¹². For SARS-CoV-2, a recent study revealed that the whole genome sequence of SARS-CoV-2 was most closely related with a bat coronavirus detected in bats from Yunnan Province in China¹³. Consequently, bats were defined as reservoirs of ancestral coronaviruses which were transmitted to humans¹⁴,¹⁵. Other wildlife species appear to

<footnote n="1">¹Centre International de Recherches Médicales de Franceville (CIRMF), BP, 769, Franceville, Gabon. ²Université des Sciences et Techniques de Masuku (USTM), Institut National Supérieur d'Agronomie et de Biotechnologies (INSAB), BP, 913, Franceville, Gabon. ³Normandie Université, EA2656, Groupe de Recherche sur l'Adaptation Microbienne, 14000, Caen, France. ⁴Charité-Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt-Universität zu, Berlin, Germany. ⁵German Centre for Infection Research (DZIF), Heidelberg, Germany. ⁶UMR (IRD 224/CNRS 5290/UM1-UM2), Institut de Recherche pour le Développement, Montpellier, France. ✉e-mail: gael maganga@yahoo.fr</footnote>
play an important role in the chain of transmission and the emergence of these viruses in humans, namely palm civets⁵ and wild rodents, in which infections by an alphacoronavirus and a betacoronavirus have been recently identified¹⁶.

Gabon, a country located in Central Africa, displays a large diversity of wildlife species. Aside from bats in which CoVs have been detected in the cave-dwelling species *Hipposideros cf. ruber*¹⁷, to our knowledge, there is no information on the carriage of other mammalian species. CoVs are prone to host switching¹⁸ and could be a current and future threat to public health.

The aim of this study was to explore the genetic diversity and the ecology of CoVs circulating among several wild mammals in Gabon in order to determine the potential reservoir species of these viruses and the risk for zoonotic emergence, and identify the factors driving CoVs infection in bats.

## Results

**Molecular identification of coronaviruses.** A total of 1867 animal samples, including 1066 from bats of 5 different species, 494 from rodents of 16 different species, 33 from 4 species of primates, and 274 from 14 species of other wild animals, were screened for CoVs (Tables 1 and 2).

Using nested RT-PCR, coronaviruses were detected in 18 intestine samples, giving an overall detection rate of 0.96% (18/1867). Positive samples originate from bats only, and the infection rate was 1.69% (18/1066). No CoVs were found in rodents, primates, or other terrestrial wild mammals. The 18 bat samples identified as CoV positive included 13 among the 262 *Hipposideros cf. ruber* (4.96%), 4 among the 156 *Hipposideros gigas* (2.56%), and 1 among the 249 *Miniopterus inflatus* (0.4%). All the positive bats came from the Ogooue-Ivindo province, in northeastern Gabon (Fig. 1), including 15 from the Faucon cave and 3 from the Batouala cave (Table 3). Besides CoVs detected in the Ogooué-Ivindo province and in the Faucon and Batouala caves, no other CoVs were detected in another province or in another cave. Bats positive for CoVs were captured in November and December 2009, October 2010 and July 2013. No individuals captured in 2011 and 2014 were found infected with CoVs. Both females and males were infected regardless of age.

**Molecular characterization of identified CoVs.** Fifteen bat nucleotide sequences were obtained and compared to those available in the public database using the algorithm “Blastn” of NCBI BLAST. The results indicated that all the nucleotide sequences from this study matched the deposited CoVs sequences. In order to determine the genetic relationships between the 18 bat CoVs from this study and previously described CoVs, phylogenetic analysis was performed based on 327-bp RdRp truncated sequences (Fig. 2a). The results showed that 12 out of 15 sequences belonged to the *Alphacoronavirus* (α-CoV) genus (Fig. 2b), whereas the remaining three belonged to the *Betacoronavirus* (β-CoV) genus (Fig. 2c). Among the 12 sequences of α-CoVs, 7 were obtained from individuals of *H. cf. ruber* caught in 2009 in the Faucon cave, 2 from other individuals of the same species but caught the same year in the Batouala cave; the other 3 sequences were obtained from 2 *H. cf. ruber* and 1 *H. gigas*, all caught in 2010 in the Faucon cave. The 3 sequences of β-CoVs were detected only in *H. gigas*.

The sequences of *Alphacoronavirus* were 91.9–100% identical with each other. The 7 sequences from 2009 from *H. cf. ruber* individuals of the Faucon cave showed 100% identity with each other. Regarding the 2 sequences of *Alphacoronavirus* from *H. cf. ruber* of the Batouala cave, the 09GB0761 (MG963198) sequence displayed either 96.5% or 96.8% with the 7 other sequences, while the 09GB0809 (MG963189) sequence showed 100% identity with the sequences from *H. cf. ruber* of the Faucon cave. Phylogenetic analysis supported this finding (Fig. 2b) and also showed that these 7 sequences clustered with 2 bat α-CoVs sequences (JX174639 and JX174640) from 2 *H. cf. ruber* from the Faucon cave caught in 2009. These two bat sequences showed 100% identity at nucleotide level with the 7 others. The 3 sequences of α-CoV (10GB0354, 10GB0309 and 10GB0318) from 2010 were 93.2–99.7% identical with each other. The 10GB0309 (MG963199) and 10GB0318 (MG963201) sequences, identical at 99.7%, from individuals of the same species (*H. cf. ruber*) caught in the Faucon cave, diverged respectively by 6.4% and 6.8% from the 10GB0354 (MG963200) sequence from an *H. gigas* captured in the same cave.

Furthermore, all α-CoV sequences grouped with the human coronavirus 229E (HCoV-229E) strain (AF304460). The nucleotide identity with HCoV-229E was 90.5–93.6%. Finally, the 10GB0354 sequence, from *H. cf. ruber* caught in the Faucon cave, showed 93.8% nucleotide identity with a CoV strain named Alpaca CoV (JQ410000), detected in *Vicugna pacos* commonly known as alpaca (family Camelidae).

The 3 sequences of β-CoV, 13GB0214 (MG963186), 13GB0215 (MG963187) and 13GB0273 (MG963188) were obtained from 3 *H. gigas* caught in 2013 in the Faucon cave. The 13GB0214 and 13GB0215 sequences displayed 99.7% identity at nucleotide level. These two sequences shared a nucleotide identity of 96.5% and 96.9% respectively with the 13GB0273 strain. In addition, these 3 sequences formed a distinct cluster with a bat β-CoV named the Zaria bat coronavirus (ZCoV) (accession number HQ166910), detected in Nigeria in 2008 in a species of bat genetically close to *H. gigas*, *Hipposideros commersoni*. The nucleotide identities between ZCoV and 13GB0214, 13GB0215 and 13GB0273 were 96.8%, 96.5% and 97.8%, respectively.

Bats of the species *H. cf. ruber* and *H. gigas* positive for alphacoronaviruses shared the same caves with the species *Coleura afra* and *Rousettus aegyptiacus*, which were negative for alphacoronaviruses. Moreover, the three *H. gigas* bats from the Faucon cave caught in 2013 (Table 3), shared the same cave with bats of the species *C. afra* and *M. inflatus* in which no *Betacoronavirus* was detected. Only the *H. gigas* species was infected with a *Betacoronavirus*.

High throughput sequencing yielded about 9.9 million reads (Fig. 3). Quality filtering discarded 86 955 (0.87%) sequences. No reads related to coronavirus could be detected with EDGE's taxonomy classification module, but BLAST alignments (98–100% and 90–95%, for respectively query coverage and identity percentage) revealed 2 small coronavirus contigs (210 bp and 271 bp) from one sample. Those contigs have been determined as being parts of replicase polyprotein 1ab.
<table><thead><tr><th rowspan="2">Group</th><th rowspan="2">Province</th><th rowspan="2">City</th><th rowspan="2">Species</th><th colspan="3">Year of sampling</th></tr><tr><th>2012</th><th>2013</th><th>2015</th></tr></thead><tbody><tr><td colspan="7">Rodents</td></tr><tr><td rowspan="10">Estuaire</td><td rowspan="2">Libreville</td><td>Mus musculus</td><td>—</td><td>23</td><td>—</td></tr><tr><td>Mus nannomys</td><td>—</td><td>3</td><td>—</td></tr><tr><td rowspan="3">Owendo</td><td>Mus musculus</td><td>—</td><td>54</td><td>—</td></tr><tr><td>Rattus rattus</td><td>—</td><td>15</td><td>—</td></tr><tr><td>undetermined</td><td>—</td><td>1</td><td>—</td></tr><tr><td rowspan="8">Haut-Ogooue</td><td rowspan="8">Franceville</td><td>Doemys ferrugiens</td><td>6</td><td>—</td><td>—</td></tr><tr><td>Lemniscomys striatus</td><td>9</td><td>—</td><td>—</td></tr><tr><td>Lophuromys ssp</td><td>3</td><td>—</td><td>—</td></tr><tr><td>Mus nannomys</td><td>4</td><td>—</td><td>—</td></tr><tr><td>Praomys sp</td><td>12</td><td>—</td><td>—</td></tr><tr><td>Rattus rattus</td><td>68</td><td>—</td><td>—</td></tr><tr><td>Stochomys longicaudatus</td><td>3</td><td>—</td><td>—</td></tr><tr><td rowspan="8">Leconi</td><td>Cricetomys emini</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Lemniscomys striatus</td><td>—</td><td>3</td><td>—</td></tr><tr><td>Lophuromys ssp</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Mus musculus</td><td>—</td><td>8</td><td>—</td></tr><tr><td>Mus nannomys</td><td>—</td><td>16</td><td>—</td></tr><tr><td>Praomys sp</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Rattus rattus</td><td>—</td><td>38</td><td>—</td></tr><tr><td rowspan="10">Ogooue-Ivindo</td><td rowspan="10">Makokou</td><td>Cricetomys emini</td><td>—</td><td>7</td><td>—</td></tr><tr><td>Grammmys poensis</td><td>1</td><td>—</td><td>—</td></tr><tr><td>Hybomys univittatus</td><td>4</td><td>—</td><td>—</td></tr><tr><td>Hylomyscus sp</td><td>20</td><td>4</td><td>—</td></tr><tr><td>Lophuromys nudicaudus</td><td>1</td><td>2</td><td>—</td></tr><tr><td>Malacomys longipes</td><td>3</td><td>1</td><td>—</td></tr><tr><td>Mus musculus</td><td>23</td><td>15</td><td>—</td></tr><tr><td>Mus nannomys setulosus</td><td>7</td><td>—</td><td>—</td></tr><tr><td>Mus monticola</td><td>1</td><td>—</td><td>—</td></tr><tr><td>Oenomys hypoxantus</td><td>—</td><td>2</td><td>—</td></tr><tr><td>Praomys sp.</td><td>40</td><td>15</td><td>—</td></tr><tr><td>Rattus rattus</td><td>19</td><td>48</td><td>—</td></tr><tr><td rowspan="2">Ogooue-Ivindo</td><td rowspan="2">Makokou</td><td>Stochomys longipes</td><td>—</td><td>1</td><td>—</td></tr><tr><td>undetermined</td><td>7</td><td>—</td><td>—</td></tr><tr><td colspan="7">Non-human primates</td></tr><tr><td rowspan="3">Haut-Ogooue</td><td rowspan="3">Ogooue-Ivindo</td><td>Cercopithecus pogonias</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Cercopithecus cephus</td><td>—</td><td>4</td><td>—</td></tr><tr><td>Cercocebus torquatus</td><td>—</td><td>2</td><td>—</td></tr><tr><td rowspan="3">Ogooue-Ivindo</td><td rowspan="3">Ogooue-Ivindo</td><td>Cercopithecus cephus</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Cercopithecus nictitans</td><td>—</td><td>—</td><td>15</td></tr><tr><td>Colobus satanus</td><td>—</td><td>—</td><td>24</td></tr><tr><td rowspan="3">Ngounie</td><td rowspan="3">Ogooue-Ivindo</td><td>Cercopithecus cephus</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Cercopithecus nictitans</td><td>—</td><td>2</td><td>—</td></tr><tr><td>Cercocebus torquatus</td><td>—</td><td>9</td><td>—</td></tr><tr><td colspan="7">Other wild animals</td></tr><tr><td rowspan="8">Haut-Ogooue</td><td rowspan="8">Ogooue-Ivindo</td><td>Eudorcas rififons</td><td>—</td><td>5</td><td>—</td></tr><tr><td>Hystrix cristata</td><td>—</td><td>4</td><td>—</td></tr><tr><td>Loxodont africana</td><td>—</td><td>—</td><td>15</td></tr><tr><td>Manis javanica</td><td>—</td><td>3</td><td>—</td></tr><tr><td>Strix aluco</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Syncerius caffer nanus</td><td>—</td><td>—</td><td>49</td></tr><tr><td>Sylvicapra grimmia</td><td>—</td><td>—</td><td>1</td></tr><tr><td>Cephalophus monticola</td><td>—</td><td>—</td><td>8</td></tr><tr><td rowspan="3">Ogooue-Ivindo</td><td rowspan="3">Ogooue-Ivindo</td><td>Cephalophus silvicultor</td><td>—</td><td>—</td><td>3</td></tr><tr><td>Cephalophus dorsalis</td><td>—</td><td>—</td><td>20</td></tr><tr><td>Civettictis civetta</td><td>—</td><td>—</td><td>7</td></tr><tr><td colspan="7">Continued</td></tr></tbody></table>
<table><caption>**Table 1.** Wild animal species analyzed by year and by locality.</caption><thead><tr><th rowspan="2">Group</th><th rowspan="2">Province</th><th rowspan="2">City</th><th rowspan="2">Species</th><th colspan="3">Year of sampling</th></tr><tr><th>2012</th><th>2013</th><th>2015</th></tr></thead><tbody><tr><td></td><td></td><td></td><td>Eudorcas rififrons</td><td>—</td><td>6</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Genetta genetta</td><td>—</td><td>—</td><td>1</td></tr><tr><td></td><td></td><td></td><td>Hypotragus equinus</td><td>—</td><td>4</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Hystrix cristata</td><td>—</td><td>1</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Loxodonta africana</td><td>—</td><td>—</td><td>20</td></tr><tr><td></td><td></td><td></td><td>Manis javanica</td><td>—</td><td>4</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Potamochoerus porcus</td><td>—</td><td>1</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Syncerius caffer nanus</td><td>—</td><td>—</td><td>15</td></tr><tr><td></td><td></td><td></td><td>Tragulus javanicus</td><td>—</td><td>—</td><td>1</td></tr><tr><td rowspan="2">Ogooue-Lolo</td><td rowspan="2"></td><td rowspan="2"></td><td>Eudorcas rififons</td><td>—</td><td>2</td><td>—</td></tr><tr><td>Hypotragus equinus</td><td>—</td><td>2</td><td>—</td></tr><tr><td rowspan="2">Ngounie</td><td rowspan="2"></td><td rowspan="2"></td><td>Hystrix cristata</td><td>—</td><td>2</td><td>—</td></tr><tr><td>Manis javanica</td><td>—</td><td>1</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Eudorcas rififrons</td><td>—</td><td>6</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Osteolaemus tetraspis</td><td>—</td><td>1</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Hypotragus equinus</td><td>—</td><td>2</td><td>—</td></tr><tr><td rowspan="2">Woleu-Ntem</td><td rowspan="2"></td><td rowspan="2"></td><td>Manis javanica</td><td>—</td><td>1</td><td>—</td></tr><tr><td>Eudorcas rififrons</td><td>—</td><td>2</td><td>—</td></tr><tr><td></td><td></td><td></td><td>Hypotragus equinus</td><td>—</td><td>1</td><td>—</td></tr></tbody></table>

<table><caption>**Table 2.** Species of bats analyzed by year and by locality.</caption><thead><tr><th rowspan="2">Province</th><th rowspan="2">Cave</th><th rowspan="2">Species</th><th colspan="5">Number of samples collected per year</th><th rowspan="2">Total</th><th rowspan="2">Number of bats infected with CoVs</th></tr><tr><th>2009</th><th>2010</th><th>2011</th><th>2013</th><th>2014</th></tr></thead><tbody><tr><td rowspan="3">Haut-Ogooue</td><td rowspan="3">Djibilong</td><td>Hipposideros cf. ruber</td><td>—</td><td>—</td><td>—</td><td>5</td><td>—</td><td>5</td><td>0</td></tr><tr><td>Miniopterus inflatus</td><td>—</td><td>—</td><td>—</td><td>69</td><td>—</td><td>69</td><td>0</td></tr><tr><td>Rousettus aegyptiacus</td><td>—</td><td>—</td><td>—</td><td>8</td><td>—</td><td>8</td><td>0</td></tr><tr><td rowspan="10">Ogooue-Ivindo</td><td rowspan="10">Batouala</td><td>Coleura afra</td><td>10</td><td>—</td><td>—</td><td>—</td><td>—</td><td>10</td><td>0</td></tr><tr><td>Hipposideros cf. ruber</td><td>58</td><td>—</td><td>—</td><td>—</td><td>—</td><td>58</td><td>2</td></tr><tr><td>Hipposideros gigas</td><td>17</td><td>—</td><td>—</td><td>—</td><td>—</td><td>17</td><td>0</td></tr><tr><td>Miniopterus inflatus</td><td>102</td><td>—</td><td>—</td><td>6</td><td>—</td><td>108</td><td>1</td></tr><tr><td>Rousettus aegyptiacus</td><td>4</td><td>—</td><td>—</td><td>6</td><td>—</td><td>10</td><td>0</td></tr><tr><td rowspan="3">Zadie</td><td>Hipposideros cf. ruber</td><td>46</td><td>16</td><td>25</td><td>1</td><td>—</td><td>88</td><td>0</td></tr><tr><td>Hipposideros gigas</td><td>7</td><td>5</td><td>25</td><td>3</td><td>—</td><td>40</td><td>0</td></tr><tr><td>Rousettus aegyptiacus</td><td>37</td><td>28</td><td>50</td><td>35</td><td>—</td><td>150</td><td>0</td></tr><tr><td rowspan="5">Faucon</td><td>Coleura afra</td><td>—</td><td>34</td><td>50</td><td>18</td><td>—</td><td>102</td><td>0</td></tr><tr><td>Hipposideros cf. ruber</td><td>57</td><td>13</td><td>25</td><td>16</td><td>—</td><td>111</td><td>11</td></tr><tr><td>Hipposideros gigas</td><td>24</td><td>30</td><td>25</td><td>16</td><td>—</td><td>95</td><td>4</td></tr><tr><td>Miniopterus inflatus</td><td>—</td><td>7</td><td>50</td><td>11</td><td>—</td><td>68</td><td>0</td></tr><tr><td>Rousettus aegyptiacus</td><td>—</td><td>—</td><td>—</td><td>—</td><td>4</td><td>4</td><td>0</td></tr><tr><td rowspan="3">Ogooue-Lolo</td><td rowspan="3">Kessipoughou</td><td>Hipposideros gigas</td><td>—</td><td>—</td><td>—</td><td>—</td><td>4</td><td>4</td><td>0</td></tr><tr><td>Miniopterus inflatus</td><td>—</td><td>—</td><td>—</td><td>—</td><td>4</td><td>4</td><td>0</td></tr><tr><td>Rousettus aegyptiacus</td><td>—</td><td>—</td><td>—</td><td>—</td><td>119</td><td>119</td><td>0</td></tr><tr><td></td><td></td><td></td><td></td><td>362</td><td>133</td><td>250</td><td>194</td><td>127</td><td><strong>1066</strong></td><td><strong>18</strong></td></tr></tbody></table>

**Factors driving CoV positivity in bats.** The comparison of the proportions of infected bats by category for each variable showed that the differences observed between bat genera (Fig. 4a) were statistically significant (p-value = 0.00002804), as well as the variations within species (p-value = 0.00003354) (Fig. 4b), ages (p-value = 0.0000355) (Fig. 4c), months of capture (p-value = 0.00003423) (Fig. 4e) and sites (p-value = 0.0004276) (Fig. 3f). There was no significant difference between the 2 sexes (Fig. 4d). Furthermore,
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Location of animal sampling sites in Gabon.</figcaption></figure>

<table><caption>**Table 3.** Characteristics of bats infected with Coronaviruses.</caption><thead><tr><td>Bat ID</td><td>Genus</td><td>Species</td><td>Sex</td><td>Age</td><td>Location</td><td>Date of capture</td><td>Detected CoVs</td><td>Genbank accession number</td></tr></thead><tbody><tr><td>09GB0274</td><td>Hipposideros</td><td>cf. ruber</td><td>F</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Alphacoronavirus</td><td>MG963194</td></tr><tr><td>09GB0296</td><td>Hipposideros</td><td>cf. ruber</td><td>F</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Alphacoronavirus</td><td>MG963195</td></tr><tr><td>09GB0323</td><td>Hipposideros</td><td>cf. ruber</td><td>M</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Alphacoronavirus</td><td>MG963192</td></tr><tr><td>09GB0328</td><td>Hipposideros</td><td>cf. ruber</td><td>F</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Undetermined</td><td>NA</td></tr><tr><td>09GB0329</td><td>Hipposideros</td><td>cf. ruber</td><td>M</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Alphacoronavirus</td><td>MG963193</td></tr><tr><td>09GB0376</td><td>Hipposideros</td><td>cf. ruber</td><td>M</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Alphacoronavirus</td><td>MG963196</td></tr><tr><td>09GB0379</td><td>Hipposideros</td><td>cf. ruber</td><td>M</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Alphacoronavirus</td><td>MG963197</td></tr><tr><td>09GB0380</td><td>Hipposideros</td><td>cf. ruber</td><td>F</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Undetermined</td><td>NA</td></tr><tr><td>09GB0383</td><td>Hipposideros</td><td>cf. ruber</td><td>F</td><td>SA</td><td>Faucon cave</td><td>November 2009</td><td>Alphacoronavirus</td><td>MG963191</td></tr><tr><td>09GB0761</td><td>Hipposideros</td><td>cf. ruber</td><td>M</td><td>A</td><td>Batouala cave</td><td>December 2009</td><td>Alphacoronavirus</td><td>MG963198</td></tr><tr><td>09GB0809</td><td>Hipposideros</td><td>cf. ruber</td><td>M</td><td>A</td><td>Batouala cave</td><td>December 2009</td><td>Alphacoronavirus</td><td>MG963189</td></tr><tr><td>09GB0812</td><td>Miniopterus</td><td>inflatus</td><td>M</td><td>A</td><td>Batouala cave</td><td>December 2009</td><td>Undetermined</td><td>NA</td></tr><tr><td>10GB0309</td><td>Hipposideros</td><td>cf. ruber</td><td>M</td><td>Juv</td><td>Faucon cave</td><td>October 2010</td><td>Alphacoronavirus</td><td>MG963199</td></tr><tr><td>10GB0318</td><td>Hipposideros</td><td>cf. ruber</td><td>F</td><td>Juv</td><td>Faucon cave</td><td>October 2010</td><td>Alphacoronavirus</td><td>MG963201</td></tr><tr><td>10GB0354</td><td>Hipposideros</td><td>gigas</td><td>M</td><td>A</td><td>Faucon cave</td><td>October 2010</td><td>Alphacoronavirus</td><td>MG963200</td></tr><tr><td>13GB0214</td><td>Hipposideros</td><td>gigas</td><td>M</td><td>A</td><td>Faucon cave</td><td>Jully 2013</td><td>Betacoronavirus</td><td>MG963186</td></tr><tr><td>13GB0215</td><td>Hipposideros</td><td>gigas</td><td>M</td><td>SA</td><td>Faucon cave</td><td>Jully 2013</td><td>Betacoronavirus</td><td>MG963187</td></tr><tr><td>13GB0273</td><td>Hipposideros</td><td>gigas</td><td>M</td><td>A</td><td>Faucon cave</td><td>Jully 2013</td><td>Betacoronavirus</td><td>MG963188</td></tr></tbody></table>

the MCA on all the bats revealed that the most important variables were the month of capture and the site followed by the species. In infected bats, the most important variables were also the month of capture, followed by the site and species (Fig. 5a). Figure 5b and Table 4 show the association between the categories of the variables and the characteristics of the different groups, respectively. According to these results, it appears that CoV infection in *H. gigas* is correlated with the months of July and October. In addition, infection with CoVs in the *H. cf. ruber* species can be associated with the months of November but also October. It appears that CoV infection is particularly associated with the Faucon cave.

## Discussion

We studied 1867 animals from 39 species including rodents, bats, primates and other wildlife species, from 6 provinces (Estuaire, Haut-Ogooué, Ogooué-Ivindo, Ogooué-Lolo, Ngounié and Woleu-Ntem) in Gabon. This study identified coronaviruses (CoVs) in 18 bats belonging to the species *H. cf. ruber*, *H. gigas* and *M. inflatus* from two caves located in the Ogooue-Ivindo province located northeast of Gabon. The nucleotides sequences obtained showed that these CoVs belonged to the *Alphacoronavirus* and *Betacoronavirus* genera. In a previous
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Phylogenetic relationship of coronaviruses based on a fragment of 495-bp of RdRP gene (a). Phylogenetic relationship between Alphacoronavirus (b), Betacoronavirus (c) and bat coronaviruses from this study are highlighted.</figcaption></figure>

study, the identification of CoVs was already reported in *H. cf. ruber* in Gabon¹⁷. Taken together, these results indicated that genetically distinct coronaviruses co-circulate among bats in the Ogooue-Ivindo province. These findings supported that bats carry a great genetic diversity of CoVs and Woo *et al.*¹⁹ suggested that bats were ideal hosts for both alphacoronaviruses and betacoronaviruses and could play an important role in the ecology and evolution of coronaviruses. Although no coronavirus was identified in other animal species, alpha- and betacoronaviruses were already found naturally in various wild mammals such as civets⁵, buffalos²⁰ and rodents²¹⁻²³. In their studies in rodents from China, Lau *et al.*²¹ and Wang *et al.*²² analyzed 725 and 1465 rodents, respectively, whereas we analyzed only 494 rodents. The lack of detection of CoVs in our study could be linked to the small sample size, unlike the aforementioned studies. However, in their study, Ge *et al.*²³ found 23 infected rodents out of the 177 studied. Unfortunately the rodents in our study come mainly from urban localities, rodent captures should also be considered in rural forest areas around caves inhabited by bats. Moreover, experimental infections have been described in non-human primates²⁴,²⁵ but no natural infection of CoV in NHP have yet been reported. For other wild animals, the small size of samples investigated by species could explain the lack of CoVs found in this study.

Alphacoronaviruses and betacoronaviruses were detected only in *H. cf. ruber* and *H. gigas*, 2 among the 5 species tested. These results from the Ogooue-Ivindo province suggest that these 2 species could be both natural reservoirs of and specific for the alpha- and the betacoronavirus, respectively. A close co-evolutionary association could be suspected between CoVs, particularly betacoronaviruses, and bats of this family²⁶. Interestingly, in our
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Identification of viral sequences from raw data generated by high-throughput sequencing. From the raw data, BLAST was performed against viral sequences from RefSeq to identify potential viruses. Then to discard false positive calls, reads having aligned with viral references were again aligned with NCBI Nucleotide database (nr).</figcaption></figure>

study, we detected a specific *Hipposideridae Betacoronavirus* clustering with a *Betacoronavirus* previously detected in Thailand in 2007²⁶. It was suggested that phylogenetic relationships between betacoronaviruses are mainly driven by the host phylogeny²⁶. However, in a study conducted in Africa, Quan *et al.*²⁷ suggested that there was no strict association between bat species and betacoronaviruses. Besides this study, previous studies on CoV host species in Asia have suggested that there was no strict association between bat species and alphacoronaviruses²⁸,²⁹. Indeed, many studies have described coronaviruses in a wide range of bat species³⁰,³¹. Nevertheless, despite the report of alphacoronaviruses in numerous species of Miniopteridae and Vespertilionidae, Ar Gouilh *et al.*³² suggested a strict association between *Alphacoronavirus* EPI6 and Rhinolophidae, a family usually associated with betacoronaviruses.

No CoVs were detected in chiropteran species sharing the same caves with the positive species (*H. cf. ruber* and *H. gigas*), suggesting a specificity of CoVs for these species or an absence of contact or interaction between some species. It was observed, for example, in the Faucon cave, that colonies of *H. cf. ruber* and *H. gigas* were found separately, in specific areas of the cave, the same with the other species (*M. inflatus* and *C. afra*), except for *R. aegyptiacus* which was not present in this cave. This spatial segregation could limit direct contact between species and represent a behavioral barrier to interspecific transmission, as previously described for lentiviruses and felids³³,³⁴ and more recently for bats and coronaviruses³².

A strain of *Alphacoronavirus* was identified in *7 H. cf. ruber* nesting in the same cave (Faucon cave) and caught the same year (2009). Likewise, 2 other *H. cf. ruber* from the Batouala cave but caught in 2010 were infected by the same alphacoronavirus infecting the *7 H. cf. ruber* from the Faucon cave. Corman *et al.*¹⁴ had already identified a closely related strain of *Alphacoronavirus* in *2 H. cf. ruber* from the Faucon cave in 2009. Altogether, these findings first suggest the transmission of CoVs between individuals of the same species through close contacts probably increased by social and roosting behavior, large colony size and reduced space in these caves²⁶; this strain of *Alphacoronavirus* could be enzootic in *H. cf. ruber* in this area of Gabon and roosts harboring a large number of bats nesting together promotes viral diffusion in the colony²⁹,³⁵, and second, local movements of bats individuals in the region imputable to seasonal/reproductive or opportunistic roosting behaviors, as the two caves being only 67 km apart. The migratory behavior of some bats provides an opportunity for pathogens to spread over long distances, as reported for flying foxes of the species *Eidolon helvum* able to migrate more than 1000 km³⁶.

No individuals caught in 2011 and 2014 were found infected with CoVs. In 2011, only the guano of *C. afra*, *H. cf. ruber*, *H. gigas*, *M. inflatus* and *R. aegyptiacus* species was collected, unlike other years where organs were taken from bats. In this cave, the *R. aegyptiacus* colony lives in areas separated from other species including *H. cf. ruber* and *H. gigas*. The lack of detection in 2011 may be due to the poor quality of the samples having undergone several thawing. In 2014, only the *R. aegyptiacus* species was sampled in the Kessipoughou cave. In this cave, the *R. aegyptiacus* populations live in areas separated from other species including *H. cf. ruber* and *H. gigas*. However, Tong *et al.*³⁷ and Lau *et al.*³⁰ detected alphacoronaviruses in *Rouettus* bats in Kenya and China.

*Alphacoronavirus* strains detected in *H. cf. ruber* from the Faucon cave displayed a variable nucleotide identity of 91 to 93% with human coronavirus 229E. This virus causes a respiratory disease and is the cause of seasonal epidemics in humans³⁸. In Ghana, Pfefferle *et al.*⁶ identified a strain of *Alphacoronavirus* in *H. cf. ruber* closely related to the human coronavirus 229E. Our finding argues that close relatives of the human coronavirus 229E exist in African bats¹⁴. The phylogenetic proximity of the CoV strains identified in this study with the human CoV-229E strain suggests a risk of emergence of this virus in humans in Africa. Furthermore, the 10GB0354 alphacoronavirus strain (MG963188), from *H. gigas* bats caught in 2010 in the Faucon cave, was phylogenetically close (93.8% nucleotide identity) to the Alpaca CoV strain detected in the *Vicugna pacos* species, commonly known as alpaca (JQ410000). Maganga *et al.*¹⁷ and Corman *et al.*¹⁴ already described two Gabonese CoV strains from hipposiderid bats caught in 2009 in the Faucon cave closely related to Alpaca CoV. Indeed, this coronavirus was described in 2007 during the epidemic of severe respiratory disorders associated with abortions in alpacas. It could be the direct progeny of its common ancestor with HCoV-229E, sharing with the latter a nucleotide
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Comparison of the proportions of infected bats by genus (a), species (b), age (c), sex (d), month (e) and site (f). The error bars represent the confidence interval.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Multiple correspondence analysis on infected bats. Projection of the variables (a); projection of individuals (b).</figcaption></figure>

<table><caption>**Table 4.** Characteristics of individuals in groups.</caption><thead><tr><th rowspan="2">Groups</th><th colspan="6">Categories</th></tr><tr><th>Genus</th><th>Species</th><th>Sex</th><th>Age</th><th>Site</th><th>Month</th></tr></thead><tbody><tr><td>Group 1</td><td>Hipposideros</td><td>caffer</td><td>F</td><td>SA</td><td>Faucon cave</td><td>November</td></tr><tr><td>Group 2</td><td>Hipposideros</td><td>caffer</td><td>F/M</td><td>J</td><td>Faucon cave</td><td>October</td></tr><tr><td>Group 3</td><td>Hipposideros</td><td>gigas</td><td>M</td><td>A/SA</td><td>Faucon cave</td><td>July/October</td></tr><tr><td>Group 4</td><td>Hipposideros/Miniopterus</td><td>caffe/inflatus</td><td>M</td><td>A</td><td>Batouala cave</td><td>December</td></tr></tbody></table>

identity of 92.2%³⁹,⁴⁰, on 440 bases of the nsp12 (polymerase). The phylogenetic link between these 2 CoVs could originate from an inter-species transmission, as a result of the opportunistic evolution of their common ancestor. As a result of genetic mutations and recombinations, the Alpaca CoV would have passed from bats to alpaca and *vice versa*. Phylogenetic analysis also revealed that two betacoronaviruses from 2 *H. gigas*, caught in 2013 in the Faucon cave, were phylogenetically close to a bat betacoronavirus named the Zaria bat coronavirus (ZCoV), detected in Nigeria in 2008 in a species of chiroptera genetically close to *H. gigas*, *H. commersoni*²⁷, suggesting that phylogenetic proximity between host species may promote inter-species transmission⁴¹.

BLAST succeeded in finding 2 contigs linked to coronavirus while no reads from this species could be detected by all EDGE's taxonomy tools. To identify coronavirus reads having being assembled into those 2 contigs, further analyses using BLAST against NCBI RefSeq virus database were carried, allowing for the identification of 7
coronavirus reads. Since BLAST proved to be more sensitive, the same was done for the 2 other samples in order to verify that no coronavirus reads could have been missed. No more reads related to coronavirus were found. Moreover, this analysis overall showed that viral reads were extremely under-represented. One could argue that not finding the targeted virus might be caused by low DNA concentration in sequencing libraries. However, even in the sample where the highest concentration has been recorded, no coronavirus reads were found. Furthermore despite having achieved a satisfying quality sequencing, the biased amplification (majority species will be over-represented) and the low viral load explain why coronavirus was identified with such weak coverage.

Knowledge on the diversity of coronaviruses in bats and the factors promoting this diversity is important to prevent the risks of emergence in humans. Based on the analysis of data on life history traits of bats collected in the field, our findings suggest that the diversity of bat coronaviruses was associated with the seasonality and the site of capture, as well as the species of bats. It appeared that CoV infection in *H. gigas* correlated with the months of July and October, while CoV infection in *H. cf. ruber* was associated with the months of October and November.

In tropical ecosystems, some bat species seem to breed year-round, whereas others have annual or biannual birth peaks⁴². In Gabon, the parturitions in *H. cf. ruber* occur synchronously but some colonies give birth in March and others in October. In *H. gigas* from the Republic of Congo, a neighboring country, the young are born at the beginning of the rainy season. In Gabon, the months of October and November correspond to the beginning of the rainy season. A study demonstrated a high prevalence of CoVs in Chinese bats during the beginning of the rainy season⁴³. The same finding was reported for the shedding of astroviruses in insectivorous bats in Asia (Borneo)⁴². Moreover, these months also correspond to the parturition season for *H. gigas* and *H. cf. ruber*. Indeed, lactating females of *H. cf. ruber* and juveniles of the same species were observed in October in northeastern Gabon. Regarding *H. gigas*, these observations are more difficult, Brosset and Saint Girons⁴⁴ reported that the seasonal events of mating and parturition in *H. gigas* were confined in large caves where these large bats are entirely dependent to raise their young. According to Plowright *et al*.⁴⁵, reproductive stress was an important driver of Hendra virus antibody prevalence. The prevalence of antibodies against the Hendra virus was higher in pregnant and lactating female of *Pteropus scapulatus* contrary to non-reproductive females. Different studies showed that parturition and lactation were important risk factors for CoVs shedding in insectivorous bats⁴⁶,⁴⁷. These two physiological states are energetically costly and thus lead to a depression of the immune system.

Moreover, our findings suggested that CoV infection was particularly associated with the Faucon cave. Of all the caves studied, the Faucon cave was the largest, with many cavities sheltering a diversity of bat species living in sympatry. This cave is regularly visited by villagers to hunt bats for human consumption. This hunting pressure could cause habitat disturbance and stress to these animals. Although Seltmann *et al*.⁴² found no association between habitat disturbance and the detection rate of CoVs in insectivorous bats, numerous studies in bats reported that pathogen prevalence was associated with habitat disturbance⁴⁸. Many authors argued that habitat disturbance in some bats species may cause chronic stress and consequently immunosuppression⁴²,⁴⁹⁻⁵¹, increasing the susceptibility of animals to acquire and shed viruses.

In summary, the results of our study highlight a genetic diversity of CoVs in insectivorous bats in northeastern Gabon and show that the beginning of the rainy season and mainly the parturition season constitute a period of risk for CoV infection of insectivorous bats *H. cf. ruber* and *H. gigas*. These results also suggest an association between the disturbance of the bats' habitat by human activities and the infection of these animals with these viruses. This study highlights a probable seasonality of the infection of insectivorous cave bats of the north-east of Gabon by CoVs. To conclude, a longitudinal monitoring should be carried out at different months of the year in order to determine the reproduction periods of the different bat species and therefore the periods favorable to viral shedding. The effect of the habitat disturbance on viral shedding should also be studied further. Bats harbor a large number of potentially zoonotic CoVs and pose a threat to public health for populations that are in contact with reservoir species of these viruses. The knowledge of CoVs ecology in bats is essential to prevent emergence in human populations.

## Methods

**Ethical approval.** All experiments were performed in accordance with relevant guidelines and regulations. Experimental protocol was approved by the National Center of Scientific and Technologic Research (CENAREST) (permission number AR0031/09) for the use of non-human primate fecal samples. Fecal sample collection from non-human primates (NHPs) in protected areas was approved by the National Parcs Agency (ANPN). Bat and rodent captures and sampling were conducted with the permission of the Wildlife and Hunting Department of the Gabonese Ministry of Water and Forestry (N°003/MEFE-PA/SG/DGEF/DCF, N°0021/MEFE-PA/SG/DGEF/DCF and N°0031/MEFDD/SG/DGEF/DFC). All the capture events, animal handling, euthanasia and transfer of samples across the country were performed in accordance with the guidelines of the American Society of Mammalogists (http://www.mammalsociety.org/committees/animal-care-and-use). Safe and humane euthanasia was achieved through the use of halothane prior to autopsy. We performed the autopsies under a microbiological safety station equipped with personal protective equipment, mask and visor.

**Study area and sample collection.** We screened several samples from bushmeat, bats, rodents and NHPs collected from 2009 to 2011 and from 2013 to 2015 in six provinces of Gabon (Fig. 1).

*Bushmeat sampling.* Bushmeat was collected from officers of the provincial direction of the Gabon Ministry of Water Affairs and Forestry after seizures of illegal hunting, as well as from salesmen during the hunting season. The sampling was performed in 2013 in four provinces of Gabon (Haut-Ogooué, Ngounié, Ogooué-Ivindo and Woleu-Ntem). Species identification was performed as previously described⁵². Samples, including whole blood,
liver, spleen and gut, were then kept in liquid nitrogen until their delivery at the CIRMF (*Centre International de Recherches Médicales de Franceville*) and stored at -80°C for further laboratory investigations.

**Bat and rodent sampling.** Bat capture took place in three provinces of Gabon, namely Ogooué-Ivindo, Ogooué-Lolo and Haut-Ogooué in 2009, 2010, 2011, 2013 and 2014. In Ogooué-Ivindo, capture was performed in three caves, the Faucon (1°07 N 13°20 E), Zadie (0°98 N 13°19 E) and Batouala (0°82 N 13°45 E) caves, located in the heart of the Belinga mountains. The distances between the three different caves located in Belinga were as follows: 33.7 km between the Batouala and Zadie caves, 11.7 km between the Zadie and Faucon caves and 37.7 km between the Batouala and Faucon caves⁵³.

In the Ogooué-Lolo province, bats were caught in the Kessipoughou cave (0°86 S 12°77 E); and in the Haut-Ogooué province, bats were captured in the Djibilong cave (1°36 S 13°46 E), located north of Franceville. For bats, sex, stage of maturity or age (adult, subadult, juvenile) and capture season were recorded. Rodents were sampled in 2012 and 2013, using Tomahawk and Sherman traps as previously described⁵⁴. The traps were set in houses and their vicinities in three provinces of Gabon: Estuaire (Libreville and Owendo), Ogooué-Ivindo (Makokou) and Haut-Ogooué (Franceville and Leconi). After autopsy, organs were collected in bats and rodents, frozen and then stored at -80°C until analysis. Bat and rodent species were identified by trained field biologists.

**Feces collection.** Fecal samples were collected from NHPs, bats and wild animals. Fresh bat droppings were collected on plastic film arranged below roost sites in the caves⁵⁵. Additionally, fecal pellets were directly collected in bags in which the bats were put individually after capture. Fecal samples from NHPs were collected in the field in the Ogooué-Ivindo province (Lope National Park) in 2015. Feces from other wild animals were sampled in 2015 in the Haut-Ogooué (Lekedi Park in Bakoumba) and Ogooué-Ivindo (Lope National Park and Mopia) provinces. Fresh fecal samples were preserved in approximately 10 ml of RNAlater (Ambion) or frozen and transferred to the CIRMF laboratory, where they were stored at -80°C until analysis.

**Specimen preparation and RNA extraction.** Total RNAs were extracted from the intestines of bats and rodents coming from bushmeat, and from the feces of bats, PNHs and other wild animals. Around 100 mg of intestine was pooled according to the species and crushed in 600 µl of cold phosphate buffer saline (PBS) (Biological Diagnostic Supplies Ltd) in a ball-mill tissue grinder (Geno/Grinder 2000, Spex Centripep). Total RNA was extracted from 100 µl of homogenate supernatant using a Biorobot EZ1 and the EZ1 RNA tissue mini kit (Qiagen) according to the manufacturer's guidelines in a BSL3. For feces preserved in RNAlater solution, RNAlater was taken out before feces was suspended in PBS and centrifuged at 1500 rpm for 5 min. RNAs were then extracted from 400 µl of supernatant, using the EZ1 RNA Viral Mini kit (Qiagen) according to the manufacturer's recommendations.

The RNA was then tested for the detection of CoVs using a nested reverse transcription-polymerase chain reaction (nRT-PCR)⁵⁶ using generic primers targeting a 494 bp fragment of a conserved region of the RNA-dependant RNA polymerase (RdRp) gene. The first-round PCR was carried out using the Qiagen One-step RT PCR kit (Qiagen). Amplification involved 30 min at 50 °C; 15 sec at 95 °C, followed by a touchdown element of 10 cycles of 20 sec at 94 °C, 30 sec starting at 60 °C with a decrease of 0.5 °C per cycle, and 40 sec at 72 °C; and 40 cycles of 20 sec at 95 °C, 30 sec at 54 °C, 40 sec at 72 °C, with a final elongation step of 5 min at 72 °C. The second-round PCR was performed using the *Platinum Taq* DNA polymerase kit (Life Technologies). Amplification involved 3 min at 94 °C followed by 40 cycles of 20 sec at 94', 30 sec at 60 °C, 40 sec at 72 °C, with a final elongation step of 1 min at 72 °C. Individual samples from pools which gave positive results were tested separately by the same nRT-PCR.

**Hight-throughput sequencing.** The hight-throughput sequencing was performed using a MiSeq machine. The DNA library was prepared from 3 samples (13GB0214, 13GB0215 and 09GB0323) using the NEBNext Ultra RNA Library Prep kit for Illumina (New England BioLabs, UK) following the manufacturer's recommandations. Briefly, after RNA quantification and fragmentation, cDNAs were generated with SuperScript III reverse transcriptase. Agencourt AMpure XP beads (Beckmann Coulter, The Netherlands) was used for size selection purification and resulted in a total library size 300–500 bp. Library quality was evaluated on an Agilent BioAnalyser 2100 (Agilent Technologies, USA). Sequencing was performed on Illumina's MiSeq (V2 reagents) to obtain 2 ×150 bp paired-end reads. Bioinformatic analyses were conducted with EDGE v2.3.1 pipeline⁵⁷. EDGE's quality module (FastQC v2.08 tool) filtered sequencing reads based on a 25 PHRED score threshold, no ambiguous bases and a 50 bp minimum size. Then taxonomy classification was conducted on both trimmed reads and contigs after assembly (SPAdes v3.13.0). The EDGE's taxonomy module allows for a broad exploration using minimap2 v2.16, Kraken2 v2.0.7, Metaphlan v2.7.7, BWA v0.7.12 and GOTTCHA2 v2.1.5. Unclassified contigs were aligned with BLAST⁵⁸ and NCBI RefSeq virus database (September 16, 2019 release).

**Phylogenetic analyses.** The partial RdRp sequences obtained were cleaned and the contigs assembled using ChromasPro 1.5. The nucleotide sequences were compared to those available in GenBank using the "Blastn" algorithm of NCBI BLAST. Multiple sequence alignments were performed using the ClustalW algorithm implemented in the MEGA 7 software⁵⁹. Phylogenetic analyses were performed by maximum likelihood algorithm using 1,000 bootstraps replications with PhyML⁶⁰ and also under a Bayesian statistical framework implemented in BEAST (version 1.8.3)⁶¹ using the general time reversible model of substitution, with a gamma distribution and a proportion of invariant sites (GTR + I + G).

**Nucleotide sequence accession numbers.** All the sequences obtained in this study have been deposited in Genbank under the accession numbers MG963186-MG963189 and MG963191-MG963201.
<table><caption>**Table 5.** Description of risk factors for infection of cave-dwelling bats with CoVs.</caption><thead><tr><th>Variables</th><th>Description</th><th>No. of animals</th></tr></thead><tbody><tr><td><strong>Species</strong></td><td></td><td></td></tr><tr><td></td><td>Coleura afra</td><td>112</td></tr><tr><td></td><td>Hipposideros cf. ruber</td><td>262</td></tr><tr><td></td><td>Hipposideros gigas</td><td>156</td></tr><tr><td></td><td>Miniopterus inflatus</td><td>249</td></tr><tr><td></td><td>Rousettus aegyptiacus</td><td>287</td></tr><tr><td><strong>Sex</strong></td><td></td><td></td></tr><tr><td></td><td>Male</td><td>620</td></tr><tr><td></td><td>Female</td><td>443</td></tr><tr><td></td><td>NA</td><td>3</td></tr><tr><td><strong>Sites/Caves</strong></td><td></td><td></td></tr><tr><td></td><td>Batouala</td><td>203</td></tr><tr><td></td><td>Djibilong</td><td>82</td></tr><tr><td></td><td>Faucon</td><td>376</td></tr><tr><td></td><td>Kessipoughou</td><td>127</td></tr><tr><td></td><td>Zadie</td><td>278</td></tr><tr><td><strong>Months</strong></td><td></td><td></td></tr><tr><td></td><td>January</td><td>65</td></tr><tr><td></td><td>February</td><td>127</td></tr><tr><td></td><td>March</td><td>93</td></tr><tr><td></td><td>April</td><td>129</td></tr><tr><td></td><td>June</td><td>101</td></tr><tr><td></td><td>July</td><td>112</td></tr><tr><td></td><td>October</td><td>77</td></tr><tr><td></td><td>November</td><td>74</td></tr><tr><td></td><td>December</td><td>288</td></tr><tr><td><strong>Age</strong></td><td></td><td></td></tr><tr><td></td><td>Adult</td><td>823</td></tr><tr><td></td><td>Subadult</td><td>205</td></tr><tr><td></td><td>Juvenile</td><td>35</td></tr><tr><td></td><td>NA</td><td>3</td></tr></tbody></table>

**Statistical analysis.** Bat data were analyzed for factors associated with CoV positivity. The methodology is based on the combined use of multiple correspondence analysis and hierarchical clustering. First, multiple correspondence analysis (MCA) was done on all the individuals, and then only on the positive bats to better identify these factors. For each variable included in the analyses (age, sex, bat genus and species, month of capture, site) (Table 5), the proportions of the infected animals were compared using the Chi square test. A *p* value of <0.05 was considered for significant results. A characterization of the infected bats was made through a hierarchical clustering by the k-means method. The data were analyzed using R version 3.4.3.

Received: 17 September 2018; Accepted: 6 April 2020;

Published online: 30 April 2020

**References**

1. Peiris, J. S., Yuen, K. Y., Osterhaus, A. D. & Stöhr, K. The severe acute respiratory syndrome. *New England Journal of Medicine* **349**, 2431–2441 (2003).

2. de Wit, E., van Doremalen, N., Falzarano, D. & Munster, V. J. SARS and MERS: recent insights into emerging coronaviruses. *Nature Reviews Microbiology* **14**, 523 (2016).

3. WHO. Coronavirus disease 2019 (COVID-19). Situation Report-51. Available at, https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200311-sitrep-51-covid-19.pdf. Accessed March 13. 9 (2020).

4. Assiri, A. Hospital Outbreak of Middle East Respiratory Syndrome Coronavirus (vol 369, pg 407, 2013). *New England Journal of Medicine* **369**, 886–886 (2013).

5. Wang, M. et al. SARS-CoV infection in a restaurant from palm civet. *Emerging infectious diseases* **11**, 1860 (2005).

6. Pfefferle, S. et al. Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats, Ghana. *Emerging infectious diseases* **15**, 1377 (2009).

7. Donaldson, E. F. et al. Metagenomic analysis of the viromes of three North American bat species: viral diversity among different bat species that share a common habitat. *Journal of virology* **84**, 13004–13018 (2010).

8. Huynh, J. et al. Evidence supporting a zoonotic origin of human coronavirus strain NL63. *Journal of virology* **86**, 12816–12825 (2012).

9. Sabir, J. S. et al. Co-circulation of three camel coronavirus species and recombination of MERS-CoVs in Saudi Arabia. *Science* **351**, 81–84 (2016).
10. Vijgen, L. et al. Complete genomic sequence of human coronavirus OC43: molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event. *Journal of virology* 79, 1595–1604 (2005).

11. Drexler, J. F., Corman, V. M. & Drosten, C. Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS. *Antiviral research* 101, 45–56 (2014).

12. King, A. M., Lefkowitz, E., Adams, M. J. & Carstens, E. B. Virus taxonomy: ninth report of the International Committee on Taxonomy of Viruses. Vol. 9 (Elsevier, 2011).

13. Zhou, P. et al. Discovery of a novel coronavirus associated with the recent pneumonia outbreak in humans and its potential bat origin. *Nature Reviews Microbiology* 579 (2020).

14. Corman, V. M. et al. Evidence for an Ancestral Association of Human Coronavirus 229E with Bats. *Journal of Virology* 89, 11858–11870, https://doi.org/10.1128/jvi.01755-15 (2015).

15. Hu, B., Ge, X., Wang, L.-F. & Shi, Z. Bat origin of human coronaviruses. *Virology journal* 12, 221 (2015).

16. Monchatre-Leroy, E. et al. Identification of alpha and beta coronavirus in wildlife species in France: bats, rodents, rabbits, and hedgehogs. *Viruses* 9, 364 (2017).

17. Maganga, G. D. et al. Bat distribution size or shape as determinant of viral richness in african bats. *PLoS One* 9 (2014).

18. Woo, P. C., Lau, S. K., Huang, Y. & Yuen, K.-Y. Coronavirus diversity, phylogeny and interspecies jumping. *Experimental Biology and Medicine* 234, 1117–1127 (2009).

19. Woo, P. C. et al. Discovery of seven novel Mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus. *Journal of virology* 86, 3995–4008 (2012).

20. Lau, S. et al. First genome sequences of buffalo coronavirus from water buffaloes in Bangladesh. *New microbes and new infections* 11, 54–56 (2016).

21. Lau, S. K. et al. Discovery of a novel coronavirus, China Rattus coronavirus HKU24, from Norway rats supports the murine origin of Betacoronavirus 1 and has implications for the ancestor of Betacoronavirus lineage A. *Journal of virology* 89, 3076–3092 (2015).

22. Wang, W. et al. Discovery, diversity and evolution of novel coronaviruses sampled from rodents in China. *Virology* 474, 19–27 (2015).

23. Ge, X.-Y. et al. Detection of alpha-and betacoronaviruses in rodents from Yunnan, China. *Virology journal* 14, 98 (2017).

24. Munster, V. J., de Wit, E. & Feldmann, H. Pneumonia from human coronavirus in a macaque model. The New England journal of medicine 368 (2013).

25. Johnson, R. F. et al. 3B11-N, a monoclonal antibody against MERS-CoV, reduces lung pathology in rhesus monkeys following intratracheal inoculation of MERS-CoV Jordan-n3/2012. *Virology* 490, 49–58 (2016).

26. Gouilh, M. A. et al. SARS-Coronavirus ancestor's foot-prints in South-East Asian bat colonies and the refuge theory. *Infection, Genetics and Evolution* 11, 1690–1702 (2011).

27. Quan, P.-L. et al. Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria. *MBio* 1, e00208–00210 (2010).

28. Chu, D. et al. Coronaviruses in bent-winged bats (Miniopterus spp.). *Journal of General Virology* 87, 2461–2466 (2006).

29. Tang, X. et al. Prevalence and genetic diversity of coronaviruses in bats from China. *Journal of virology* 80, 7481–7490 (2006).

30. Lau, S. K. et al. Recent transmission of a novel alphacoronavirus, bat coronavirus HKU10, from Leschenault's rousettes to pomona leaf-nosed bats: first evidence of interspecies transmission of coronavirus between bats of different suborders. *Journal of virology* 86, 11906–11918 (2012).

31. Goffard, A. et al. Alphacoronaviruses detected in french bats are phylogeographically linked to coronaviruses of european bats. *Viruses* 7, 6279–6290 (2015).

32. Gouilh, M. A. et al. SARS-CoV related Betacoronavirus and diverse Alphacoronavirus members found in western old-world. *Virology* 517, 88–97 (2018).

33. Troyer, J. L. et al. FIV cross-species transmission: an evolutionary prospective. *Veterinary immunology and immunopathology* 123, 159–166 (2008).

34. VandeWoude, S., Troyer, J. & Poss, M. Restrictions to cross-species transmission of lentiviral infection gleaned from studies of FIV. *Veterinary immunology and immunopathology* 134, 25–32 (2010).

35. Woo, P. C. et al. Molecular diversity of coronaviruses in bats. *Virology* 351, 180–187 (2006).

36. Epstein, J. H. et al. Pteropus vampyrus, a hunted migratory species with a multinational home-range and a need for regional management. *Journal of Applied Ecology* 46, 991–1002 (2009).

37. Tong, S. et al. Detection of novel SARS-like and other coronaviruses in bats from Kenya. *Emerging infectious diseases* 15, 482 (2009).

38. Su, S. et al. Epidemiology, genetic recombination, and pathogenesis of coronaviruses. *Trends in microbiology* 24, 490–502 (2016).

39. Crossley, B. M. et al. Identification of a novel coronavirus possibly associated with acute respiratory syndrome in alpacas (Vicugna pacos) in California, 2007. *Journal of veterinary diagnostic investigation* 22, 94–97 (2010).

40. Crossley, B. M., Mock, R. E., Callison, S. A. & Hietala, S. K. Identification and characterization of a novel alpaca respiratory coronavirus most closely related to the human coronavirus 229E. *Viruses* 4, 3689–3700 (2012).

41. Baric, R. S., Sullivan, E., Hensley, L., Yount, B. & Chen, W. Persistent infection promotes cross-species transmissibility of mouse hepatitis virus. *Journal of virology* 73, 638–649 (1999).

42. Seltmann, A. et al. Seasonal fluctuations of astrovirus, but not coronavirus shedding in bats inhabiting human-modified tropical forests. *EcoHealth* 14, 272–284 (2017).

43. Lau, S. K. et al. Ecoepidemiology and complete genome comparison of different strains of severe acute respiratory syndrome-related Rhinolophus bat coronavirus in China reveal bats as a reservoir for acute, self-limiting infection that allows recombination events. *Journal of virology* 84, 2808–2819 (2010).

44. Brosset, A. & Saint Girons, H. Cycles de reproduction des microchiroptères troglophiles du nord-est du Gabon. *Mammalia* 44, 225–232 (1980).

45. Plowright, R. K. et al. Reproduction and nutritional stress are risk factors for Hendra virus infection in little red flying foxes (Pteropus scapulatus). *Proceedings of the Royal Society B: Biological Sciences* 275, 861–869 (2008).

46. Gloza-Rausch, F. et al. Detection and prevalence patterns of group I coronaviruses in bats, northern Germany. *Emerging infectious diseases* 14, 626 (2008).

47. Drexler, J. F. et al. Amplification of emerging viruses in a bat colony. *Emerging infectious diseases* 17, 449 (2011).

48. Turmelle, A. S. et al. Ecology of rabies virus exposure in colonies of Brazilian free-tailed bats (Tadarida brasiliensis) at natural and man-made roosts in Texas. *Vector-Borne and Zoonotic Diseases* 10, 165–175 (2010).

49. Sapolsky, R. M., Romero, L. M. & Munck, A. U. How do glucocorticoids influence stress responses? Integrating permissive, suppressive, stimulatory, and preparative actions. *Endocrine reviews* 21, 55–89 (2000).

50. Suorsa, P. et al. Effects of forest patch size on physiological stress and immunocompetence in an area-sensitive passerine, the Eurasian treecreeper (Certhia familiaris): an experiment. *Proceedings of the Royal Society of London. Series B: Biological Sciences* 271, 435–440 (2004).

51. Wikelski, M. & Cooke, S. J. Physiological assessment of environmental stressors. *Trends in Ecology & Evolution* 1, 38–46 (2006).

52. Boundenga, L. et al. Haemosporidian parasites of antelopes and other vertebrates from Gabon, Central Africa. *Plos One* 11 (2016).

53. Rougeron, V. et al. Characterization and phylogenetic analysis of new bat astroviruses detected in Gabon, Central Africa. *Acta Virol* 60, 386–392 (2016).
54. Duplantier, J.-M. Les rongeurs myomorphes forestiers du nord-est du Gabon: structure du peuplement, démographie, domaines vitaux. *Revue d'écologie* (1989).

55. Drexler, J. F. et al. Henipavirus RNA in African bats. *PloS one* 4 (2009).

56. de Souza Luna, L. K. et al. Generic detection of coronaviruses and differentiation at the prototype strain level by reverse transcription-PCR and nonfluorescent low-density microarray. *Journal of clinical microbiology* 45, 1049–1052 (2007).

57. Li, P.-E. et al. Enabling the democratization of the genomics revolution with a fully integrated web-based bioinformatics platform. *Nucleic acids research* 45, 67–80 (2017).

58. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. *Journal of molecular biology* 215, 403–410 (1990).

59. Kumar, S., Stecher, G. & Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. *Molecular biology and evolution* 33, 1870–1874 (2016).

60. Dereeper. A. et al. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. *Nucleic Acids Res* 36, W465–469, https://doi.org/10.1093/nar/gkn180 (2008).

61. Drummond, A. J. & Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7, 214, https://doi.org/10.1186/1471-2148-7-214 (2007).

## Acknowledgements

We acknowledge Heidi Lançon for the English revision, and Andy Nkili Meyong and Ingrid Labouba for high-throughput sequencing and bioinformatics analysis. We thank all the persons involved in sample collection. The CIRMF is supported by the Government of Gabon, Total-Fina-Elf Gabon, and the Ministère de la Coopération Française.

## Author contributions

G.D.M. and E.M.L. designed and conceived the work. G.D.M., A.P. and M.M. carried out the experiments, G.D.M., I.M.M., A.M.M.B., L.B., M.A.G. and N.N. analyzed the data, J.F.D., C.D. and E.M.L provided essential material. G.D.M. and I.M.M. wrote the main manuscript text. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to G.D.M.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

Open Access This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020