Wide Diversity of Coronaviruses in Frugivorous and Insectivorous Bat Species: A Pilot Study in Guinea, West Africa

Zoonoses can constitute a threat for public health that can have a global importance, as seen with the current COVID-19 pandemic of severe acute respiratory syndrome coronavirus (SARS-CoV2). Bats have been recognized as an important reservoir of zoonotic coronaviruses (CoVs). In West Africa, where there is a high diversity of bat species, little is known on the circulation of CoVs in these hosts, especially at the interface with human populations. In this study, in Guinea, we tested a total of 319 bats belonging to 14 genera and six families of insectivorous and frugivorous bats across the country, for the presence of coronaviruses. We found CoVs in 35 (11%) of the tested bats—in three insectivorous bat species and five fruit bat species that were mostly captured close to human habitat. Positivity rates varied from 5.7% to 100%, depending on bat species. A wide diversity of alpha and beta coronaviruses was found across the country, including three sequences belonging to SarbeCoVs and MerbeCoVs subgenera known to harbor highly pathogenic human coronaviruses. Our findings suggest that CoVs are widely spread in West Africa and their circulation should be assessed to evaluate the risk of exposure of potential zoonotic CoVs to humans.


Introduction
Emerging infectious diseases (EIDs) represent a significant challenge to global health. The frequency of EID outbreaks has dramatically increased over the last decades. Most EIDs have a zoonotic origin; more than 70% of them are caused by pathogens originating from wildlife [1]. The continuing encroachment of human populations into areas initially only occupied by wildlife, such as pristine forest-together with environmental and climatic changes-has led to increasing contacts between humans and animals and as a consequence, exposure to the pathogens they harbor [2][3][4]. The ongoing outbreak with the new severe acute respiratory syndrome coronavirus (SARS-CoV-2), responsible for the COVID-19 pandemic, is a perfect illustration of the potential impact and the global consequences resulting from a single cross-species transmission from an animal virus. The epidemic was first recognized in December 2019 in China and has spread across the entire globe in a few months. By the end of May 2020, more than five million people have been infected; more than 350,000 have died and a lockdown of four billion people across the globe has been enforced to reduce the spread of the virus [5]. It is thus important to be prepared for new EID outbreaks. A major step in understanding the risk for zoonotic infections, is to characterize the diversity of pathogens at the interface between humans and animals. Knowing the animal sources, prevalence, geographic distribution and the genetic diversity or evolutionary history of pathogens in wildlife is critical to evaluate the risk for potential emergence or reemergence of disease.
It is now well established that CoVs are widespread in bats with relative high prevalence and important genetic diversity providing thus multiple opportunities for the potential emergence of novel HCoVs in areas where contact between humans and bats are frequent [26]. Nevertheless, despite recent efforts, our knowledge on CoVs in bats is far from complete, because most studies are conducted on a limited number of species and in restricted geographic or ecological environments. Bats represent 20% of all mammals with more than 1400 different species [27]. More studies are thus still needed, especially in areas where the risk for spillover to animals or humans is high, like southeast Asia or West and Central Africa, where numerous bat species are hunted for food or medicine [28][29][30]. Moreover, bushmeat hunting is recognized as a major threat for conservation of bats in these areas [31,32]. In addition, humans can also be indirectly exposed to bat viruses through contacts with bat guano or fruit contaminated by their saliva, urine or feces [32,33].
Several studies have shown a large diversity of CoVs in bats across Africa, but limited data are currently available on bat CoVs in Sub-Saharan Africa and especially from West Africa, where we have currently only information on bats from Ghana and Nigeria [15,[34][35][36][37][38][39][40][41][42]. Nevertheless, new alpha (α-CoVs) and betacoronaviruses (β-CoVs) have been described in Africa, including some closely related to human CoVs such as HCoV-229E, or SARS-CoVs [15,34,35]. Recently, the closest relatives in bats to the highly pathogenic MERS-CoV were documented in Neoremicia and Pipistrellus bats from South-Africa and Uganda, respectively, illustrating that African bats could also play a role in outbreaks of CoVs in humans [23,43]. Here we conducted a pilot study on the presence of CoVs in bats in Guinea, West Africa, a country with a high likelihood for EID outbreaks and where humans are frequently directly or indirectly exposed to bats (A.K. Keita, personal observation).

Study Sites, Sample Collection and Ethical Permits
Samples were collected from free-ranging frugivorous and insectivorous bats in Guinea between February 2016 and January 2017 as previously described [44]. We obtained permission to conduct research and to collect samples from the National Ethics Committee for Health Research (CNERS) from Guinea (approval reference 074/CNERS /15, 26 November 2015). Briefly, bats were captured using mist nets or harp traps in roosting and foraging sites. Bats were then kept in cotton bags and released immediately after sampling. Whole blood was taken by venipuncture of the propatagial or brachial vein and blood drops were directly transferred onto Whatman 903 filter paper (GE-Healthcare, Feasterville-Trevose, PA, USA). Samples were air-dried and preserved individually in plastic bags containing silica desiccant and stored in a hermetic box as dried blood spots (DBS). DBS were stored at −40 • C in the laboratory (2 to 3 weeks after collection) and kept frozen until analysis. Rectal and oral swabs were also collected in 500 µL RNA-later to preserve RNA from degradation and stored in the field at ambient temperature for maximum two weeks and subsequently at −40 • C in the laboratory. For a subset of animals, fresh fecal samples were also collected and stored in RNA-later, i.e., feces dropped when bats were in cotton bags. For each sampled bat, information on capture sites (GPS coordinates, ecological environment), capture method, morphology (body measurements, weight, color), sex, age class (adult, juvenile) and visual species identification were recorded.

Nucleic Acid Extraction
Total DNA and RNA were extracted from swabs and feces using the NucliSENS EasyMAG platform (BioMérieux, Marcy-l'Etoile, France). Briefly, 250 µL of sample was incubated with 2 mL of lysis buffer for 15 min and extraction was performed using manufacturer's instructions. Total nucleic acids were resuspended in 60 µL elution buffer.

Molecular Confirmation of Bat Species
Species identification recorded in the field was molecularly confirmed on a subset of bats using DBS as in our previous study on Ebolaviruses in bats [44]. In addition, the species was also molecularly confirmed for each bat in which a coronavirus was identified. A fragment of approximately 800 base pairs (bp) of the mitochondrial cytochrome b (CytB) region was amplified using previously described primers to identify mammal species including bats, Cytb-L14724 (forward) and Cytb-H15506 (reverse), [45]. However, the forward primer was replaced by a newly designed primer Cytb-L1 5 -ATG ACC AAC ATC CGA AAA TCN CAC-3 or Cytb-L2 5 -ATY TCY TCM TGA TGA AAY TTY GGM T-3 to increase PCR specificity for certain bat species. PCR products were diluted and directly sequenced with BigDye Terminator version 3.1 sequencing kit (Life Technologies, Courtaboeuf, France). Electrophoresis and data collection were processed on an Applied Biosystems 3500 Genetic Analyzer (Thermo Fisher Scientific, Foster City, CA, USA). Sequences from both strands were reconstituted using SeqMan Pro tool from the package DNAStar v17.0.2 (Lasergene, Madison, WI, USA). Sequences were pasted in the NCBI BLAST web interface (https://blast.ncbi.nlm.nih.gov/Blast.cgi) to identify the most similar bat species. For samples with no or low similarity (<97%) hits with species in Genbank, a phylogenetic tree was constructed using maximum likelihood methods implemented in PhyML with reference sequences in order to obtain genus identification. The new sequences were deposited in the GenBank under accession numbers MT586773-MT586805.

RT-PCR Screening for Detection of Coronavirus RNA
First, cDNA was synthesized from denatured RNA (70 • C for 10 min) using a Reverse Transcription System kit with random primers (Promega, Madison, WI, USA), following manufacturer's instructions. PCR screening was done with a broadly reactive nested PCR approach in the highly conserved RdRp region. Our assay targeted a 440 bp fragment using degenerate primers from a modified protocol described by Chu et al. (2011). In the first round, a 602 bp fragment was targeted following previously described PCR methods for which amplification conditions were slightly adapted [46]. Briefly, cDNA was amplified using the GoTaq Hot Start Master Mix PCR kit (Promega, Madison, WI, USA) as follows: 40 cycles of 92 • C for 30 s, 48 • C for 30 s and 72 • C for 50 s for the first round of PCR. In the second round, a 440 bp product was amplified and cycling conditions were modified for the second round of PCR, using a touch-down technique to reduce nonspecific amplifications (10 cycles of 92 • C for 30 s, 53 • C for 30 s with −0.5 • C/cycle and 72 • C for 30 s, 35 cycles of 92 • C for 30 s, 53 • C for 30 s and 72 • C for 30 s). PCR products were first analyzed by agarose gel electrophoresis and positive amplicons were sequenced as indicated above and analyzed on the Applied Biosystems 3500 genetic analyzer platform. Sequences from both strands were reconstituted using SeqMan Pro tool from the package DNAStar v17.0.2 (Lasergene, Madison, WI, USA). Sequences were pasted in the NCBI BLAST web interface (https://blast.ncbi.nlm.nih.gov/Blast.cgi) to identify similar coronavirus sequences.
The PCR assay was first tested on six nasopharyngeal swabs from patients positive for SARS-CoV-2, which are anonymous and untraceable diagnostic leftover samples, kindly provided by the Virology department of Montpellier University Hospital in March 2020. Amplification products and sequences were obtained, demonstrating that the assay is also able to detect the novel SARS-CoV-2 strain (data not shown).

Phylogenetic Analyses
The new sequences were aligned with representatives of alpha and beta-coronavirus sequences that have been reported mainly in Africa and elsewhere (Table S1). A multiple sequence alignment (MSA) was obtained by using MAFFT v7 (https://mafft.cbrc.jp/alignment/server/). The alignment was manually checked and end-trimmed to match with the newly obtained RdRp sequences and to remove the PCR primers sequences. The final MSA was used for maximum likelihood (ML) phylogenetic analysis with GTR + F + I + 4Γ as the best-fit model of nucleotide substitution according to BIC and 1000 bootstrap resampling by using IQ-Tree server (http://iqtree.cibiv.univie.aC.at) [47,48]. Consensus trees were edited with FigTree v1.4.4. The new sequences were deposited in GenBank under accession numbers MT586830-MT586867.

Statistical Analyse
To explore the impact of age, sex and environment of the collected samples on the detection of coronavirus RNA, we performed Kruskal-Wallis test on bat field data. Significant results were considered for a p-value of <0.05.

Bat Species and Sampling
Samples from a total of 319 wild bats, captured at 14 sites in seven different prefectures in Guinea, were studied ( Figure 1, Table 1). Bats were captured in different ecological settings including the following environments: villages (38.6%), forest sites (30.1%), urban sites (23.2%), plantations (4.1%) and diverse other settings for the remainder (4.1%). Species identification was confirmed by CytB sequence analysis on a total of 171 (53.6%) bats to confirm at least one sample per sampling date, per capture method and per morphologic description at each site. Species identification was then extrapolated for the remaining samples by combining molecular and field data. For certain insectivorous bats (Molossidae, Rhinolophidae, Hipposideridae and Nycteridae) identification was Viruses 2020, 12, 855 5 of 17 only possible at the genus level, mostly due to the lack of reference sequences in Genbank. Because discrimination between Epomophorus gambiensis and Micropteropus pusillus was not possible based on CytB sequences only, we used morphologic details on forearm and weight measurements to discriminate the species for 93 bats, as previously documented by others [49]. Details on bat families, genera and species are shown in Table 2. We collected samples from 274 (85.9%) frugivorous bats (Pteropodidae family) representing at least eight species and 45 (14.1%) insectivorous bats from five families and at least six genera. Overall, 195 (61.1%) bats were female, 121 (37.9%) were male and for three (1.0%) sex classification was not available. The vast majority were adults (300; 94.1%), 11 (3.4%) were juveniles and for eight bats (2.5%), age could not be identified, or information was not recorded. Total 319 Sites are indicated with a two-letter code as in Figure

PCR Screening and Presence of α and β Coronaviruses
We screened a total of 634 samples from 319 bats by RT-PCR for the presence of coronavirus RNA: for 286 (89.7%) bats oral and rectal swabs were tested, 19 (6.6%) bats were tested using oral swabs and fecal samples, for five (1.6%) oral and rectal swabs and feces, for two and seven bats only rectal or oral swabs were tested, respectively. A total of 38 samples from 35 of the 319 tested bats (11%), were positive for the presence of CoVs; five in oral swabs only, 20 in rectal swabs only, seven in feces only, one in both feces and oral swab, two in rectal and oral swabs ( Table 3). The rate of positivity varied according to the sample type; the proportion of positive samples was 7.5% (22/293) in rectal swabs, 2.5% (8/317) in oral swabs and 33.3% (8/24) in fresh fecal samples. Among the three bats that were positive in two different samples, two were infected with the same virus in the rectal and oral swab or in feces and oral swab, but for one bat, an α-CoV and a β-CoV were detected in oral and rectal swab, respectively.

Genetic Diversity of Bat Coronaviruses
The phylogenetic analysis based on 374 bp nucleotides of the new sequences and reference sequences is shown in Figure 2. High genetic diversity was observed among the newly identified α and β-CoVs. α-CoVs were detected in 7/45 (15.6%) of insectivorous bats; five Hipposideros sp. and two Rhinolophus spp. (Figure 2A). All α-CoV strains from Hipposideros bats in Guinea clustered with previously reported CoV sequences from bats from this genus in Zimbabwe, Gabon or Ghana and comprising also the HCoV-229E virus. The new α-CoV strains from Rhinolophus bats fall in a cluster that groups sequences from other Rhinolophus bats captured in Mozambique and Kenya. Among frugivorous bats, 8/274 (2.9%) were infected with α-CoV; seven R. aegyptiacus and one L. angolensis. α-CoVs from R. aegyptiacus formed a separate lineage for which no close relatives (>90% homology) were detected in Genbank. The new α-CoV sequence from a frugivorous L. angolensis clustered with CoV sequences mainly obtained from insectivorous Molossidae bats from Cameroon, Tanzania, Southern Africa and the Western Indian ocean. Micropteropus pusillus or Rousettus sp. from central and Eastern Africa were also observed in this clade. Interestingly, the new β-CoV from the Nanonycteris bat in Guinea fell in the same cluster. The new strains from Lissonycteris angolensis were closely related to a recently reported β-CoV from the same species in Rwanda. Among the new β-CoVs from Rousettus aegyptiacus, two strains formed a separate lineage with no close relatives (>90% homology) in Genbank and the four remaining strains clustered with β-CoVs identified in Rousettus sp. from Cameroon or Kenya. The β-CoVs observed in insectivorous bats belonged to two subgenera. The two β-CoV strains from Hipposideros bats belonged to the SarbeCoVs subgenus, that also harbor the human SARS-CoV strains as well as other closely related CoVs from insectivorous bats, civets and pangolins. Our sequences were highly diverse and fell in a clade grouping a wide diversity of CoVs, all from Hipposideros bats in diverse African countries. Finally, the new β-CoV strain from the Nycteris bat, trapped in a cave, belongs to the MerbeCoV genus which includes the CoVs responsible for the human Middle East respiratory syndrome and closely related CoVs detected in camels and Vespertilionidae bats (Pipistrellus sp. from China, Uganda and Neoromica sp. from South Africa). Our sequence clustered specifically in a distinct clade with other strains from Nycteris bats from Mozambique and Ghana that was basal to the MERS-CoV related viruses.
Overall, Rousettus and Hipposideros bats from our study are infected with at least three different CoVs (one α-CoV and two β-CoVs) and L. angolensis bats with α-CoV and β-CoVs.  β-CoVs were observed in 18/274 (6.6%) frugivorous bats and 3/45 (4.4%) insectivorous bats. None of the new viruses was closely related to the human β-CoVs ( Figure 2B). Among β-CoVs, our sequences fell into three subgenera, recently named Sarbecovirus, Merbecovirus and Nobecovirus [50]. All the 18 β-CoV sequences detected in fruit bats clustered in the NobeCoV subgenus comprising only CoVs from Pteropodid bats from various places in Asia and Africa (i.e., Kenya, Cameroon, Rwanda, Madagascar, Congo, Tanzania). New β-CoVs strains from E. helvum clustered with previously reported CoV strains form this species observed in bats from Cameroon, Kenya and even Saudi Arabia. However, a β-CoV from a Rousettus bat from Madagascar fell also in this clade. Similarly, the five new β-CoV strains from Epomophorus bats clustered with previously documented strains from this genus from other regions in Africa, but strains from other bat species like Micropteropus pusillus or Rousettus sp. from central and Eastern Africa were also observed in this clade. Interestingly, the new β-CoV from the Nanonycteris bat in Guinea fell in the same cluster. The new strains from Lissonycteris angolensis were closely related to a recently reported β-CoV from the same species in Rwanda. Among the new β-CoVs from Rousettus aegyptiacus, two strains formed a separate lineage with no close relatives (>90% homology) in Genbank and the four remaining strains clustered with β-CoVs identified in Rousettus sp. from Cameroon or Kenya.
The β-CoVs observed in insectivorous bats belonged to two subgenera. The two β-CoV strains from Hipposideros bats belonged to the SarbeCoVs subgenus, that also harbor the human SARS-CoV strains as well as other closely related CoVs from insectivorous bats, civets and pangolins. Our sequences were highly diverse and fell in a clade grouping a wide diversity of CoVs, all from Hipposideros bats in diverse African countries. Finally, the new β-CoV strain from the Nycteris bat, trapped in a cave, belongs to the MerbeCoV genus which includes the CoVs responsible for the human Middle East respiratory syndrome and closely related CoVs detected in camels and Vespertilionidae bats (Pipistrellus sp. from China, Uganda and Neoromica sp. from South Africa). Our sequence clustered specifically in a distinct clade with other strains from Nycteris bats from Mozambique and Ghana that was basal to the MERS-CoV related viruses.
Overall, Rousettus and Hipposideros bats from our study are infected with at least three different CoVs (one α-CoV and two β-CoVs) and L. angolensis bats with α-CoV and β-CoVs.

Discussion
With the recent emergence of the COVID-19 pandemic caused by SARS-CoV-2 and the previous highly pathogenic SARS-CoV-1 and MERS-CoV outbreaks, it is important to better document the diversity and evolution of coronaviruses, particularly in regions at high-risk for emergence of infectious diseases of zoonotic origin. All CoVs in humans have a zoonotic origin with a direct or indirect link with bat CoVs. Here we documented for the first time the presence and genetic diversity of alpha and beta coronaviruses in bats in Guinea, West Africa, where contacts between humans and bats are frequent. In this pilot study, we analyzed more than 300 bats from at least eight frugivorous and six insectivorous species, collected at 14 different areas in seven prefectures across the country. We detected and characterized coronaviruses in 11% of the bats, at almost all sites and observed a wide diversity of α and β-CoVs with potential new viral lineages. Some strains belonged to the subgenera of MerbeCoV and SarbeCoVs known to harbor highly pathogenic HCoVs like MERS-CoV and SARS-CoV, respectively. We also report for the first time, the presence of a coronavirus in a bat from the genus Nanonycteris.
The overall rate of CoV detection (11%) in our study is consistent with the rates found in previous studies in Kenya and Mozambique [40,51], but infection rates in studies in central and southern Africa and the western Indian Ocean were generally lower [15,35,[38][39][40]42,52]. Positivity rates of coronavirus detection can vary by species and season, therefore comparison between studies is limited. In Guinea, we also observed differences in positivity rates according to species. Although we tested only a few samples per species, high rates were observed in three insectivorous genera from three families (7/21 Hipposideros sp., 1/3 Nycteris sp. and 2/5 Rhinolophus sp.) and two frugivorous species (4/9 Eidolon helvum and 1/1 Nanonycteris sp.). However, it cannot be excluded that this could be related to the season when samples were collected, for example several studies have shown temporal variations of CoV infections associated with parturition [37,41,[53][54][55]. Longer follow-up should be done on these bat species in Guinea to assess temporal dynamics of CoV infection. Most samples have been obtained in cities or villages close to human habitats, it is thus important to know whether exposure to bat CoVs vary during certain periods of the year. Differences in positivity rates can also be explained by the nature of samples that have been tested, for example in our study we observed higher rates in feces and rectal swabs than in saliva. This was also observed in other reports and confirms that coronaviruses are highly shed in the intestinal tract [56]. Therefore, non-invasive collection of fresh bat feces, under roosting or feeding sites, may be more efficient to study CoVs, especially for species that are difficult to access.
Phylogenetic analyses revealed a wide diversity of α and β-CoVs with potential examples of co-evolution between CoVs and their bat hosts and examples of host diversity within CoV lineages. The genetic diversity observed in our study is most likely associated with the high diversity of bat species tested, which is in accordance with previous findings [40]. It also confirms the observations from the large study on CoVs in more than 12,000 bats across Africa and Asia from Anthony et al. (2017) who demonstrated that increasing sampling efforts leads to increased identification of new lineages. Although we studied only very few samples from Hipposideros bats, we did not only find a high infection rate (7/21, 33%), but observed also a high genetic diversity among their CoVs in accordance with most previous studies on Hipposideros species in Africa [15,34,35,[39][40][41][42]. The human 229E α-CoV falls in a clade of bat CoVs that is widely present in Hipposideros bats from Kenya, Zimbabwe, Mozambique, Gabon, Ghana and also in our study in Guinea. Hipposiderids roost in hollow trees, small caves and abandoned buildings, and indirect and direct exposure to infected bats and to the wide diversity of CoVs is thus very likely [57,58]. Moreover, the β-CoVs observed in Hipposideros bats belong to the SarbeCoV subgenus where they form a sister-clade to the clade of human and animal CoVs related to SARS-CoV-1 and -2. Within this viral clade, the genetic diversity is also very high, and viruses from this clade are present in Hipposideros species across Africa [15,35,39,41,59]. CoVs have also been described in Hipposideros bats in China and Thailand, suggesting a potential new clade [55,60,61]. The β-CoV observed in the Nycteris bat from Guinea falls with other CoVs from Nycteris bats in the MerbeCoV subgenus, forming a sister clade to the MERS-CoV related viruses from humans and animals. This clade is also widespread in Nycteris species across Africa [37,40]. All β-CoVs from frugivorous bats fall within the NobeCoV subgenus, in which for certain species potential new viral lineages are seen. However, this could be related to the low number of CoVs available among GenBank references. For α-CoVs we also observed mainly association between CoVs and their bat hosts, except for α-CoV virus identified in the frugivorous L. angolensis bat (CCGU550) clustering with α-CoVs from insectivorous bats.
As in other studies, we also showed that the same bat species can be infected with different α and/or β-CoV lineages and species; this was seen in L. angolensis, R. aegyptiacus and Hipposideros sp. bats [15,34,39,51,56,62]. We observed a co-infection in a single Rousettus aegyptiacus bat with an alpha and betacoronavirus. Co-infections have been previously reported among insectivorous bats in Asia, but not yet in Africa [61][62][63].
We confirm that different viruses co-circulate in the same species and that certain viral strains may be maintained in more than one host and that cross-species transmission may occur, as previously reported [34,64]. These results are important because viral sharing and cross-species transmission may lead to recombination events and the emergence of novel bat coronaviruses or novel coronaviruses in other mammal species naturally infected with CoVs [64,65]. Recent studies demonstrated that frequent cross-species transmission seem to explain the evolutionary history of CoVs rather than co-divergence [56,66]. Moreover, viruses with a higher taxonomically and ecologically diverse host range seem to be more likely to adapt to humans and spread on a broad geographic scale and thus with a higher pandemic potential [65].
With this study in Guinea, we brought additional knowledge on coronaviruses in bats in Africa and especially in West Africa where only information was available for bats from Ghana and for small numbers of bats and species in Senegal and Nigeria [15,35,37,59]. The fact that we observed a high genetic diversity with potential new viral variants illustrates clearly that knowledge on bat coronaviruses is still incomplete and that more bat CoVs will be discovered as sampling efforts will increase. Anthony et al. (2017) even estimate that more than 3000 bat CoVs circulate in the more than 1400 known bat species [56]. Given the impact of zoonotic transmissions from CoVs to humans, it is important to document the CoVs that circulate in bats, but also in other wildlife species because for certain HCoVs, intermediate hosts were involved. For example, the closest relative to SARS-CoV-2 is the SARSr-Ra-BatCoV-RaTG13 from Rhinolophus bats from China, sharing more than 96% identity at the whole-genome level, but the receptor binding domain (RBD) exhibits stronger similarity to CoVs identified in pangolins, suggesting recombination [24,25]. However, it cannot be excluded that a not yet identified CoV in bats or another animal species is at the origin of the current SARS-CoV-2 epidemic. A prerequisite for cross-species transmission is compatibility between virus and human receptors, therefore future studies should not only focus on partial RdRp sequences but should thus also include characterization of the RBD of CoVs. In addition, as reminder, our findings are based on a small fragment of the RdRp gene, which is a highly conserved region. Therefore, complete genome sequencing and analysis of some key genes such as the Spike and ORF8 genes that have been shown to be hypervariable, should allow to assess recombination events [61,64]. This will bring additional data to retrieve phylogenetic and temporal relationships between these bat CoVs and other pathogenic strains and thus give us more insight into their zoonotic potential.

Conclusions
In this pilot study, we have demonstrated high rates and important genetic diversity of CoVs in frugivorous and insectivorous bats in Guinea. More studies are needed on CoVs in bats and wildlife in Guinea, but also in other geographic areas of Africa. This will certainly lead to the identification of more novel CoV species and maybe also new genera. Although we now have data for several thousands of bats worldwide, this still remains relatively small given the high numbers of bats that constitute colonies and the high number of bat species worldwide [67]. Future studies should focus on areas and species where contacts between humans are frequent or could become frequent due to environmental and climate changes. The high genetic diversity and complex evolutionary history of co-evolution, host switching, and recombination illustrates the complexity to track the animal reservoir of CoVs and to predict the areas at highest risk and species infected with CoVs with higher adaptive probability to infect humans. Efforts must continue not only on bats, but also on other animals to elucidate where CoVs circulate in wildlife and to clarify the role of intermediate species.